Упрaжнение 2. Matlab
Contents
Help
>>help
>>symbolic toolbox
>>demo
- Индексиране на вектори
v = [1 2 3 Inf 0 2 NaN]
функция isfinite() - връща само крайните стойности v(isfinite(v) ) = []
v(v>4 | v<1) = []
вектор с равномерно разпределение
a + (b-a)*rand(1,7)
Хистограма
Хистограма на нормално разпределение
>> y = round(100*randn(1,10000)); >> hist(y,20)
Двумерна графика
Цвят на линията, ..., включване на грид >> hist(y,20), grid on
Много прозорци plot в един
subplot(3,2,1)
m = 3 -> редове n = 2 -> koloni a = 1 -> ацтивен
plot
Въвеждане на специални символи
Синтаксис от Latex
Скриптови файлове
поредица от команди записани във файл
Разделяме прозореца 4 части
>> subplot(2,2,1)
затихваща синусоида с различен коеф. на затихваме
y = e^(-k*x.sin(x))
>>k = [0.1, 0.2, 0.3, 0.4] >>x = 0:pi/100:6*pi
Създаване на потребителски функции. Два входящи елемента x,k
Създаваме нов m файл!
- Името на файла трябва да отговаря на името на функцията
- оператор точка за по елементно умножение
for - при незадаване на стъпка, тя е единица
Управляващи оператори на plot(x,y, 'r:')
function y = f( x,k )
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
y = 2.^( (-k*x).*sin(x) );
end
>>for i =1:4 subplot(2,2,i); plot(x, f(x,k(i)), 'r' ); end
>>ezsurf('sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)', [-6*pi 6*pi])