Difference between revisions of "Упрaжнение 2. Matlab"

From Ilianko
Line 241: Line 241:
 
</pre></code>
 
</pre></code>
  
 
+
Намиране корени на полином по метода на Newton
>> x = 0:0.01:1.5;
+
<code><pre>
>> plot(x,newton(x));
+
>> x = 0:0.01:1.5;;
>>  grid on
 
>> plot(x,newton(x));
 
>> plot(x,newton(x));
 
 
>> plot(x,newton(x)), grid on;
 
>> plot(x,newton(x)), grid on;
 
>>  
 
>>  
>> a = 0.1, b = 0
+
>> a = 0.1, b = 0;
 
 
a =
 
 
 
    0.1000
 
 
 
 
 
b =
 
 
 
    0
 
 
 
 
>> eps = 1.e-10
 
>> eps = 1.e-10
  
Line 270: Line 257:
 
a = a-newton(a)/fp(a)
 
a = a-newton(a)/fp(a)
 
end
 
end
 
+
</pre></code>
 
 
>> 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)
 
  
 
[[Category:Matlab]]
 
[[Category:Matlab]]
end
 

Revision as of 10:19, 15 April 2011

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

Намиране корени на полином по метода на Newton

>> x = 0:0.01:1.5;;
>> plot(x,newton(x)), grid on;
>> 
>> a = 0.1, b = 0;
>> eps = 1.e-10

eps =

   1.0000e-10

>> while  abs (b-a) >=eps
b=a;
a = a-newton(a)/fp(a)
end