1、第一步，建立一个M文件，用来存贮函数，本例题以达芬方程（Duffing）为例，其中force为参数function df=dafen(t,x,flag,force)df=[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];
2、第二步，建立一个画图的M文件clearforce=0.222;options=odeset猾诮沓靥('RelTol',1e-7);%定义误差精度的，系统默认1e-3，如果观鲼视防改为1e-3，X将等于0tt=2*pi/1.2 %定义步长的[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],force);figureplot(x(2000:end,1),x(2000:end,2),'-')%X=x-xxx %检验options的%pojialaihold oni=2000:100:3000plot(x(i,1),x(i,2),'*')
