8416073118

Przykład:

W poniższym przykładzie szukamy rozwiązania numerycznego równania różniczkowego drugiego stopnia van der Pola

Przyjmując podstawienie

yl⁼ y 2

otrzymujemy układ równań pierwszego stopnia w postaci:

y[=y₂

function dydt = f(t,y,mu) % funkcja w której zapisaliśmy równanie różniczkowe dydt = [y(2); mu*(1-y(1)'‘2)*y(2)-y(1) ];

» Mu=100; % przyjmujemy parametr p

»tspan = [0; max(20,3*MU)]; % przedział, dla którego szukamy rozwiązań » yO = [2; 0]; %warunki początkowe » [ty] = ode15s(@f,tspan,y0); % metoda ode15s » plot(t,y(:,1)); % zobrazowanie wyników »title(['Solution of van der Pol Equation, \mu = ' num2str(MU)]);

» xlabel('time t');

» ylabel('solution y_1 ’);

» axis([tspan(1) tspan(end) -2.5 2.5]);

17