IMAG0152 (11)

V, cizyA"-^:	?ib*dy/
icho ^1,11
fclC glob*^1 bb-inp^ul ^	^1 podaj
to - 0/ *	( i na l
yO * li*'	5 0]; !“ 31
1 ol ■
t raC© *	L;
tfinal ■	t f i na 1

opisane w pliku rciezar.m

% Współczynnik tłumienia % Zadany przedział % Warunki początkowe % Dokładność obliczeń % Krok

I 11,y) • od*23('roherm ',tO,tfinal,yO);

% procedura rozwiązanie r. rozn. II rzędu metoda Runge-Kutty |t,y] * ode23('reiezaz',tO,tfinal,yO,tol, tracę) ;

% sporządzenie wykresu fazowego plot (y(i /1) / y(i/2)), titłe('Równanie ciężarka na sprężynie trajektoria fazowa') iclabel (’ y¹), y\ abel ('dy/dt*) kle

rciezar. m

% T. Sałacinakł & G. Sycz, Warszawa 1997

% rciezar.m funkcja opisuje równanie ciężarka zawieszonego na eprezynie

■ d2y/dt2+b*dy/dt:+y»0

/% tciezar.m jest wołane przez funkcje ODE23 w pliku ciężar.ra l;unction yp • rciezar(t/y)

o)obfll b    % Deklaracja globalna współczynnika tłumienia

yp zero# (2,1)    % Zerov/y wektor dwuelementowy

yp(1) - y(2)    % y(l)*y y(2)«dy/dt

(2) - -l*y (1) - b*y (2);

Rfekt wykonania powyźizago tri pUku znajduje się na rys.5.