请各位学长看看我这个程序错在哪里了？

美好

请各位学长看看我这个程序错在哪里了？画达芬方程分岔图的程序：

M文件：
function xdot=duffing2(t,x,omega)
global omega;
xdot=[x(2);x(1)-x(1)^3-0.3*x(2)+0.25*cos(omega*t)];

计算程序：
clear;
global omega;
omega=1:0.01:2;
for h=1:length(omega)
T=2*pi/omega(h);
[t,x]=ode45('duffing2',[0:T/100:120*T],[0 0], omega(h));
plot(omega(h),x(10000:100:end,2),'k.');hold on
end

谢谢！

咕噜噜

function xdot=duffing2(t,x,omega)?
改成function xdot=duffing2(t,x)
[t,x]=ode45('duffing2',[0:T/100:120*T],[0 0], omega(h));也少了点东西，搜索论坛，相关程序很多

无水1324

应该改为：
function xdot=duffing2(t,x,flag,omega)

美好

还是不行呀学长，还是有错误，可是哪里出了问题呢？

无水1324

把错误的信息发出来吧
还有把两个global omega去掉

美好

这是提示的错误信息！
??? Error using ==> funfun\private\odearguments
Correct syntax is ode45(DUFFING2,tspan,y0,options).

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...
请学长帮忙看看！谢谢！

无水1324

你按照上面的修改了吗？这好像不是第一次给你修改了吧

octopussheng

M文件：
function xdot=duffing2(t,x,omega)
xdot=[x(2);x(1)-x(1)^3-0.3*x(2)+0.25*cos(omega*t)];

计算程序：
clear;
omega=1:0.01:2;
for h=1:length(omega)
T=2*pi/omega(h);
[t,x]=ode45('duffing2',[0:T/100:120*T],[0 0], omega(h));
plot(omega(h),x(1000:100:end,2),'k.');hold on
end

试试吧


[ 本帖最后由 octopussheng 于 2007-10-25 22:41 编辑 ]

yuxinting

clear;
omega=1:0.01:2;
for h=1:length(omega)
T=2*pi/omega(h);
[t,x]=ode45('duffing2',[0:T/100:120*T],[0 0],[], omega(h));%加个'[]',就可以运行了.
plot(omega(h),x(1000:100:end,2),'k.');hold on
end

美好

非常感谢各位学长的热心帮助！只是按照各位的改法依旧运行不通，着实很郁闷，不知哪个小环节出了问题。恳请各位有时间的时候再帮我看看！还是非常感谢大家伸出援助之手！谢谢各位学长！

yuxinting

按照8楼学长给的程序，只需改动[],就可以运行了。