考虑vanderpol方程是什么意思

问题描述

精选答案

clear allglobal mumu=1;y0=[0 1];options = odeset('AbsTol',1e-8);[T Y]=ode45('Vandelpol',[0 100],y0,options)

;plot(T,Y(:,1))%%%%%%%%%%%%%%%%%%%function dy=Vanderpol(t,y)global mudy=zeros(2,1);dy(1)=y(2)

;dy(2)=mu*(1-y(1)^2)*y(2)-y(1);