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');