请教关于分岔程序

mechanic05

前辈好！
我的算了一晚上，到现在机子出了点问题，报废了:@L ．也许是因为别人给的ＭＡＴＬＡＢ７．０是盗版的？盗版的回影响以后应用吗？
还有，plot(omega(h),x(25000:100:end,2),'k.');hold on中，是说从２５０００开始每隔１００取值，那从接近３００００取不是更好?

请指教！

无水1324

你总共才30000个点

mechanic05

请教前辈：那摸这个问题去掉瞬态响应还有好办法吗?我前面画出的图太乱了！现在该成２５０００看起来也很慢．还在算．顺便问一下：盗版的７．０会因影响以后的计算吗？:@)

咕噜噜

盗版的７．０一般不会影响计算，不行你就改成28000

octopussheng

:@L :@L 汉！！

取[0:T/100:300*T]太慢了，我就取了100×T，结果是这样的，这个分岔图看不出效果

mechanic05

前辈好！
你的这个算了多长时间？我前面给的图大约算了１．５个小时．:lol

无水1324

6.5以上版本的计算速度相对而言要慢些
改成6.5程序不要变。

无水1324

function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=[x(2);a*x(1)-b*x(1)^3-eta*x(2)+gamma*cos(omega*t)];


clear;clc;close all;
omega=0.1:0.01:4;
for h=1:length(omega)
    T=2*pi/omega(h);
    [t,x]=ode45('dddd',[0:T/100:220*T],[0 0],[],omega(h));
    plot(omega(h),x(21000:100:end,2),'k.');hold on
end

mechanic05

谢谢前辈！
不过程序里有个小变化才行．如红色标记：
function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=[x(2);-a*x(1)-b*x(1)^3-eta*x(2)+gamma*cos(omega*t)];


clear;clc;close all;
omega=0.1:0.01:4;
for h=1:length(omega)
    T=2*pi/omega(h);
    [t,x]=ode45('dddd',[0:T/100:220*T],[0 0],[],omega(h));
    plot(omega(h),x(21000:100:end,2),'k.');hold on
end

octopussheng

差不多吧，1小时左右

无水1324

负号的话是他自己的方程的问题了。我没有看到方程，就这能对程序作出判断

octopussheng

这个确实应该是方程的问题了，按照她这个方程，分岔图看来是效果不好的！

我拿我那个方法再给她试试看！

mechanic05

前辈好！
我的方程是：
dx1/dt=x2
dx2/dt=－a*x1－b*x1^3－eta*x2+gamma*cos(omega*t)．本来就有负号的．谢谢！

octopussheng

呵呵，看到你的方程了，我用我的方法把分岔图做了一下，第一个图是a前面没有负号，即
a*x(1)-b*x(1).^3-eta*x(2)+gamma*cos(omega*t)


第二个分岔图是有负号的，即
－a*x(1)-b*x(1).^3-eta*x(2)+gamma*cos(omega*t)

第二个分岔图看的比第一个要好很多！

也请无水、咕噜噜这些高手看看，让他们帮你分析分析吧

mechanic05

多谢前辈！我的程序是：
function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=[x(2);-a*x(1)-b*x(1)^3-eta*x(2)+gamma*cos(omega*t)];


clear;clc;close all;
omega=0.01:0.01:4;
for h=1:length(omega)
    T=2*pi/omega(h);
    [t,x]=ode45('dddd',[0:T/100:220*T],[0 0],[],omega(h));
    plot(omega(h),x(10000:100:end,2),'k.');hold on
end