Difference between revisions of "Упрaжнение 2. Matlab"
| Line 226: | Line 226: | ||
== Въвеждане на стойности на функции в диалогов режим == | == Въвеждане на стойности на функции в диалогов режим == | ||
| + | <code><pre> | ||
% безкраен цикъл | % безкраен цикъл | ||
| − | |||
while 1 | while 1 | ||
x = input('x = '); | x = input('x = '); | ||
| Line 238: | Line 238: | ||
break; | break; | ||
end | end | ||
| + | end | ||
| + | </pre></code> | ||
| − | |||
| − | |||
>> x = 0:0.01:1.5; | >> x = 0:0.01:1.5; | ||
>> plot(x,newton(x)); | >> plot(x,newton(x)); | ||
| Line 299: | Line 299: | ||
b=a; | b=a; | ||
a = a-newton(a)/fp(a) | a = a-newton(a)/fp(a) | ||
| + | |||
| + | [[Category:Matlab]] | ||
end | end | ||
Revision as of 10:15, 15 April 2011
Contents
- 1 Help
- 2 Индексиране на вектори
- 3 Хистограма
- 4 Двумерна графика
- 5 Въвеждане на специални символи
- 6 Скриптови файлове
- 7 Създаване на директория
- 8 Достъп до променливи
- 9 Пример скрипт файл
- 10 Извикване на коментар на функция от команден ред
- 11 Смяна на показваната точност
- 12 Предаване на данни
- 13 Условни оператори
- 14 Въвеждане на стойности на функции в диалогов режим
Help
>>help
>>help symbolic toolbox
>>demo
Индексиране на вектори
v = [1 2 3 Inf 0 2 NaN]
Inf <=> bezkrajnost NaN <=> nesyshtestwuwa
функция 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])
Създаване на директория
Конфигуриране на собствена директория и даване на пътя
file -> set path ...
тестване: >>path
Достъп до променливи
Локални Глобални променливи
Workspace - обхваща цялото пространство
Всяка функция си има собствено пространство
f(x,y) с аргументите на функцията се предават променливи от основния воркспаце към функцията
>> global x,y; прави глобални променливите, трябва да се извика и във функцията;
Пример скрипт файл
Създайте файл със следните редове
A = input( 'A = [x1 y1 z1] , wekto s tri koordinati:')
B = input( 'B = [x2 y2 z2] , wekto s tri koordinati:')
d = sqrt( (A(1)-B(1))^2 + (A(2)-B(2))^2 + (A(3)-B(3))^2 );
disp(['sumata e ', num2str(d)])
Изпълнете в работаната среда
>> ab
A = [x1 y1 z1] , wekto s tri koordinati:[0 0 0]
A =
0 0 0
B = [x2 y2 z2] , wekto s tri koordinati:[1 1 1]
B =
1 1 1
sumata e 1.7321
Извикване на коментар на функция от команден ред
function d = ab_f(A,B)
% Prograма за намиране разстянието между две триизмерни точки
% Входни аргументи
% Вектор А = [a b c] с координатите на т.А
% Вектор B = [a b c] с координатите на т.B
d = sqrt( (A(1)-B(1))^2 + (A(2)-B(2))^2 + (A(3)-B(3))^2 );
%
end
Иzwevdaне не komentaрите до първия изпълним ред
>>help ab_f
Prograма за намиране разстянието между две триизмерни точки
Входни аргументи
Вектор А = [a b c] с координатите на т.А
Вектор B = [a b c] с координатите на т.B
ab_f([0 0 0], [1 1 1])
ans =
1.7321
Смяна на показваната точност
>>format long
Предаване на данни
чрез обявяването им като глобални
function d = ab_g(A,B)
% Prograма за намиране разстянието между две триизмерни точки
% Входни аргументи
% Вектор А = [a b c] с координатите на т.А
% Вектор B = [a b c] с координатите на т.B
global p;
d = sqrt( (A(1)-B(1))^2 + (A(2)-B(2))^2 + (A(3)-B(3))^2 );
a = d/p
%
end
[s , x] = ab_g([0 0 0], [1 1 1])
Условни оператори
for Втори вариант
A = magic(7)
A =
30 39 48 1 10 19 28
38 47 7 9 18 27 29
46 6 8 17 26 35 37
5 14 16 25 34 36 45
13 15 24 33 42 44 4
21 23 32 41 43 3 12
22 31 40 49 2 11 20
>> for i = A i', sum(i) end
while
while <expression>
statements
end
Въвеждане на стойности на функции в диалогов режим
% безкраен цикъл
while 1
x = input('x = ');
y = sin(x);
disp(['resultat',num2str(y)]);
disp('press any button');
pause;
an = input('prodylzhi (y/n)', 's'); % s - ще въвеждаме стринг
if ( (an == 'y' || an == 'Y'))
break;
end
end
>> x = 0:0.01:1.5;
>> plot(x,newton(x));
>> grid on
>> plot(x,newton(x));
>> plot(x,newton(x));
>> plot(x,newton(x)), grid on;
>>
>> a = 0.1, b = 0
a =
0.1000
b =
0
>> eps = 1.e-10
eps =
1.0000e-10
>> while abs (b-a) >=eps b=a; a = a-newton(a)/fp(a) end
>> x = 0:0.01:1.5;
>> plot(x,newton(x));
>> grid on
>> plot(x,newton(x));
>> plot(x,newton(x));
>> plot(x,newton(x)), grid on;
>>
>> a = 0.1, b = 0
a =
0.1000
b =
0
>> eps = 1.e-10
eps =
1.0000e-10
>> while abs (b-a) >=eps b=a; a = a-newton(a)/fp(a) end
