简介

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

版权声明：本文为eat-too-much原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://www.cnblogs.com/eat-too-much/p/13444998.html