solve_vdpol.m

matlab从入门到精通第6章书后源码对初学者很有帮助

% solve_vdpol.m
% 用常微分方法求解Van der Pol微分方程组
% 积分的时间区间
tspan=[0 20];
y0=[2;0];
disp('ode45求解Van der Pol微分方程：')
[t,y]=ode45(@vdpol,tspan,y0);
% 绘制数据结果
plot(t,y(:,1),t,y(:,2),'--');
xlabel('t');
legend('x','dx/dt');
title('Van der Pol solution')