xprime.m

这是一个求杜芬方程解的小程序。开发环境为matlab。通过仿真曲线可以很容易的看到解的曲线。

function xdot=xprime(t,x)
xdot=zeros(2,1);
xdot(1)=x(2);
% a=0.3;b=33;x1(0)=0.3;x2(0)=0.5;
xdot(2)=-0.3*x(2)-x(1)^3+31.5*cos(t);

% xdot(1)=(1-x(2)^2)*x(1)-x(2);
% xdot(2)=x(1);