Instrukcje- opracowanie

Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
Tryb interpretera
2 * 3
 
ans =
6
2 * 3;
pi
 
 
ans =
3.1416
sin(pi / 2)
 
ans =
1
[1 4; 6 8]
 
ans =
1
6
4
8
6: 2: 20
 
ans =
6
8
10
12
14
16
18
20
Wprowadzenie do Matlaba, folia 1
Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
b
 
? Undefined function or variable ’b’
A = [1 4; 6 8]
 
A =
1
6
4
8
A - 1
 
ans =
0
5
3
7
A * A
 
ans =
25
54
36
88
sin(A)
 
ans =
0.8415
-0.2794
-0.7568
0.9894
Wprowadzenie do Matlaba, folia 2
Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
Programowanie
Zadanie 1. Napisa´ funkcj˛ znajdujaca
c
e
˛ ˛
pierwiastek równania liniowego.
function x = rown1(a, b)
if a == 0
if b ~= 0
x = [];
else
x = NaN;
end;
else
x = -b / a;
end;
Uwaga: T˛ funkcj˛ zapisuje si˛ w pliku rown1.m !
e
e
e
Efekt działania:
rown1(1, 2)
 
ans =
-2
Wprowadzenie do Matlaba, folia 3
Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
Zadanie 2. Napisa´ funkcj˛ rozwiazujaca
c
e
˛ ˛ ˛
równania liniowe i kwadratowe.
function x = rown21(a, b, c)
if nargin == 2
x = rown1(a, b);
elseif a == 0
x = rown1(b, c);
else
delta = b * b - 4 * a * c;
if delta = 0
if b 0
x = (-b-sqrt(delta))/(2*a);
else
x = (-b+sqrt(delta))/(2*a);
end;
if x == 0
x = [x 0];
else
x = [x (c/a)/x];
end;
end;
end;
Wprowadzenie do Matlaba, folia 4
Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
Wskazana jest równie˙ diagnostyka błedów:
z
if nargin

(…)

…
z´
Newtona.
metoda
˛
©
¤ ¤
. (Łatwo
2. Podzieli´ trójmian przez
c
pokaza´ , ze odpowiedni iloraz ma posta´
c ˙
c
, gdzie:
,
.)
£
¨©¤  ¡
§
¨
¤
¨ § ¡
¤
¨ § ¨ ¦ ¤
¤©
 
¨
3. Rozwiaza´ otrzymane równanie kwadratowe
˛ c
otrzymujac pozostałe dwa pierwiastki.
˛
Wprowadzenie do Matlaba, folia 19
Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
function x1 = newton1(a, b, c, d)
epsi = 1.0e…
…
z
 
¡
 
¡
,
£
,
©
£  %# ¡¢¢¡¢# ©   ¡¢ 
¡
$ £
¥
¤
    ¡  
#
¤   #  %¤ ¡  
¥
¡
¡
 
¤
¥
gdzie:
Metoda Eulera
¦
¡
¢
¨ !  %# ! ¤
¥
¡
¡
¨ !¤ ¡ ©$ !¤
¨ !  ¡ © $ ! 
gdzie:
Wprowadzenie do Matlaba, folia 44
Wprowadzenie do Matlaba, folia 45
function yh = euler(f, x, y, h)
yh = y + h * feval(f, x, y);
Zadanie. Zaimplementowa´ wymienione metody.
c…
… jest palindromem.
function r = palindrom(s)
r = strcmp(s, s(end: -1: 1));
Wprowadzenie do Matlaba, folia 28
Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
Zadanie. Napisa´ funkcj˛ zamieniajaca pierwsze
c
e
˛ ˛
litery słów danego ciagu znaków na du˙ e litery.
˛
z
function s1 = duze_lit(s)
s = [’ ’ s];
ind = findstr(s, ’ ’);
s = [s ’ ’];
s(ind+1) = upper(s(ind+1));
s1 = s(2: end-1…
…
0
−0.2
−0.4
Euler
rk2
−0.6
rk4
rozw. analityczne
−0.8
−1
ode45
0
1
2
3
4
5
6
7
8
9
10
Wprowadzenie do Matlaba, folia 50
Instytut Sterowania i Systemów Informatycznych, Politechnika Zielonogórska
¡
¢ 
£
Zadanie. W zbiorniku umieszczono moli
i
1 mol benzenu. Oznaczmy odpowiednio przez , ,
, , i ilo´ci
s
, benzenu, chlorobenzenu,
dwuchlorobenzenu i trójchlorobenzenu w chwili .
Dla obj˛ to´ci produktu…
... zobacz całą notatkę