把以下內容copy到edit中,運行
%=====================
function hhh
clear
ff=0.222; %force的值
options=odeset('RelTol',1e-7); %定義誤差精度的,系統默認1e-3,如果改為1e-3,X將等於0
tt=2*pi/1.2 %定義步長的
[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],ff)
[ttt,xxx]=ode45(@dafen,[0:tt/100:80*tt],[0,0],[],[],ff)%按照系統默認值1e-3
plot(x(2000:end,1),x(2000:end,2),'-')
X=x-xxx %檢驗options的
%pojialai
hold on
i=2000:100:3000
plot(x(i,1),x(i,2),'*')
function df=dafen(t,x,flag,force)
df=[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];
%===========================