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