此问题为非刚性问题,故需用ode45求解器,不能用ode15s微分方程函数为: function dx = vdp0.2(t,x)%u=0.2dx = zeros(2,1); dx(1) = -x(2);dx(2) = -0.2*(1 - x(1)^2)*x(2) + x(1);指定时间间隔为[0,20],设初值为x(1)=2,x(2)=0;利用ode45求解[t,x] = ode45(@vdp0.2,[0 20],[2 0]);plot(t,x(:,1),'-',t,x(:,2),'--')希望你能对你有所帮助! |