IMAG0149 (9)

3.

yłconywsnie wykresów fazowych

7,'rkresy fazowe można wykonywać na kilka sposobów np. metodą izoklin lub metodą Jednak ze względu na złożoność tych metod do sporządzania wykresów w ćwiczeniu „ ydzie program MA TL AB

Sposób sporządzania wykresu będzie omówiony na przykładzie klasycznego oscylatora jujp-ooicznego opisanego równaniem:

d?

+ y = 0,

dla którego przyjęto warunki początkowe: y(0)=l i y(0) = 0.

Do wykreślenia trajektorii fazowej sporządzono skrypt w postaci m-pliku o nazwie oharm.m. i si-plik funkcyjny o nazwie roharm.m.

■ t. Salacinski k G. Sycz, Warszawa 1997

p charm.si rozwiązuje numerycznie r, różniczkowe II rzędu P d2y/dt2+y=0 opisane w pliku roharm.m

Lchc off

fco - 0; tfinal = 10;	% Zadany przedział
ifl 1 [0 1];	% Warunki początkowe dla	t“0
BLol = l.e-3;	% Dokładność obliczeń
Erace - 1;	% Krok
Stflnal = tfinal*(1+eps]	; >. eps - precyzja zmiennoprzecinkowa
V lt,y) B ode23(^1roharra	',tO,t finał,yO);
h procedura rozwiązanie	r. rozn. 11 rzędu metoda	Runge-Kutty
Bfc,y] i| ode23 (' roharm',	tO,tfinał,yO,tol,tracę);
ft, Sporządzenie wykresu	fazowego

Be axis('egual')

plot(y(:,1),y(:, 2)), title(¹ Równanie oscylatora trajektoria fazowa') irlabeł (' y'), ylabel (idy/dt¹) %, grid ęlc

harmonicznego

Wszystkie teksty występujące po znaku % oznaczają komentarze. Pierwsze dwie funkcje mają itarakter porządkowy. Echo steruje wyprowadzaniem komentarzy na ekran, a dc czyści ekran