3.分别用ode45和ode15s求解Van-del-Pol方程

2021-11-04 社会 157阅读
clearall
globalmu

mu=1;
y0=[01];

options=odeset('AbsTol',1e-8);
[TY]=ode45('Vandelpol',[0100],y0,options);

plot(T,Y(:,1))

%%%%%%%%%%%%%%%%%%%

functiondy=Vanderpol(t,y)
globalmu
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=mu*(1-y(1)^2)*y(2)-y(1);
声明:你问我答网所有作品(图文、音视频)均由用户自行上传分享,仅供网友学习交流。若您的权利被侵害,请联系fangmu6661024@163.com