matlab 求解高階方程

簡介

van der Pol 方程

code

dy = @(t,y)[y(2); 1000 * (1-y(1)^2)*y(2)-y(1)]; % 定義匿名函數
[t,y]= ode15s(dy,[0 3000],[2;0]); % 求數值解
plot(t,y(:,1),'*');
title('solution of van der pol equ, mu=1000');
xlabel('time t');
ylabel('solution y');