Priručnik za MATLAB - Ponude.biz

More documents

Recommendations

Info

gde je TSPAN = [TO TFINAL] interval na kome se integriše sistem diferencijalnih jednačina y’ = F(t,y), uz početni uslov Y0. ‘F’ je string koji sadrži naziv ODE datoteke. Funkcija F(T, Y) mora da vrati kolona vektor. Svaka vrsta u rešenju Y odgovara vremenu vraćenom u kolona vektoru T. Da biste dobili rešenje za određene trenutke T0, T1, ... , TFINAL, koristite TSPAN = [T0 T1 … TFINAL]. Oblik [T, Y] = ODE45 (‘F’, TSPAN, Y0, OPTIONS) radi na sličan način, samo su parametri integracije definisani preko parametra OPTIONS, koji daje funkcija ODESET. (Za detaljnije objašnjenje u komandnoj liniji treba navesti help odeset.) Obično se koriste skalarna vrednost, tolerancija relativne greške ‘RelTol’ (podrazumeva se vrednost 1e-3), i vektorska vrednost, tolerancija apsolutne greške ‘AbsTol’ (podrazumeva se vrednost 1e-6 za sve komponente). Vratimo se sada našem problemu. Prvo ćemo definisati ekvivalentni sistem diferencijalnih jednačina 1-og reda. Definišimo vektore u 1 = θ , i u 2 = θ '. Sada imamo sledeći sistem: gde je g 2 =g/L. u 1 '(t)=u 2 (t), u 1 (0) = θ 0 u 2 '(t)= - g 2 sin(u 1 (t)), u 2 (0) = 0 Definisaćemo sada m-datoteku za prosto klatno, tj. klatnou.m function uprim=klatnou(t,u) % jednacina prostog klatna % uprim(1)=u(2) % uprim(2)=-omega*sin(u(1)) omega=16;uprim=[u(2);-omega*sin(u(1))]; Definisaćemo sada m-datoteku za linearizovan model (za male uglove možemo smatrati da je sin θ = θ ), tj. klatnow.m function wprim=klatnow(t,w) % linearizovana jednacina klatna omega=16;wprim=[w(2);-omega*w(1)]; „Glavni” program je m-datoteka klatno1.m koja poziva prethodne funkcije, tj. mdatoteke. % model prostog klatna clear all,tspan=[0,2*pi];u0=[pi/4;0]; opcije=odeset('RelTol',1e-5,'AbsTol',1e-8) [t,u]=ode45('klatnou',tspan,u0,opcije); % za linearizovan model prostog klatna [s,w]=ode45('klatnow',tspan,u0,opcije); plot(t,u(:,1),'k-',s,w(:,1),'r:') xlabel('vreme')
ylabel('crveno = linearizovanu model, crno = nelinearni model') Pokrenimo MATLAB i proverimo kako ovo funkcioniše. Pripremio Dragan Marković
Page 1 and 2: 2D grafika MATLAB raspolaže bogati
Page 3 and 4: Dobijen je rezultat kao da ste sva
Page 5 and 6: for i = 1:100 xnovo = (xstart + x/x
Page 7: ilo kog funkcijskog programa za pro
Page 11 and 12: Pripremio Dragan Marković U ovom t
Page 13 and 14: Napomena: Kada radite sa velikim sk
Page 15 and 16: Ako želite da obojite male element
Page 17 and 18: Sada ćemo nacrtati eksponencijalno
Page 19 and 20: Komanda subplot Komanda subplot omo
Page 21 and 22: definiše da su podeoci na obe ose
Page 23 and 24: gde a -1 označava inverznu matricu
Page 25 and 26: x1 - 4x2 + 3x3 = -7 3x1 + x2 - 2x3
Page 27 and 28: T1_var = (std(vremtemp(:,2)))^2 daj
Page 29 and 30: y = a0 + a1x + a2x2 + a3x3 + ......
Page 31 and 32: 4x1 + 3x2 - 2x3 + 2x4 + x6 = -1 3x1
Page 33 and 34: Kao i kod svakog drugog softvera, u
Page 35 and 36: c cyan x x-mark -. dashdot r red +
Page 37 and 38: Unos matrica Matrice možete uneti
Page 39 and 40: X = A; X(3,4) = 17 X = 1 1 3 0 4 0
Page 41 and 42: sum(A(1:3,3)) izračunava sumu tre
Page 43 and 44: A = [ ... 16.0 3.0 2.0 13.0 5.0 10.
Page 45 and 46: ' Transponovanje ( ) Zadavanje redo
Page 47 and 48: Matrična izračunavanja Pošto mat
Page 49 and 50: odgovarajućim elementom druge matr
Page 51 and 52: kod koga se odgovarajuće stavke mn
Page 53: 1. a = g(:,2) 2. b = g(4,:) 3. c =

Priručnik za MATLAB - Ponude.biz

Create successful ePaper yourself

Delete template?

Save as template?