回复 #90 咕噜噜 的帖子
前辈好！我的算了一晚上，到现在机子出了点问题，报废了:@L ．也许是因为别人给的ＭＡＴＬＡＢ７．０是盗版的？盗版的回影响以后应用吗？
还有，plot(omega(h),x(25000:100:end,2),'k.');hold on中，是说从２５０００开始每隔１００取值，那从接近３００００取不是更好?
请指教！ 你总共才30000个点
回复 #92 无水1324 的帖子
请教前辈：那摸这个问题去掉瞬态响应还有好办法吗?我前面画出的图太乱了！现在该成２５０００看起来也很慢．还在算．顺便问一下：盗版的７．０会因影响以后的计算吗？:@) 盗版的７．０一般不会影响计算，不行你就改成28000 :@L :@L 汉！！取太慢了，我就取了100×T，结果是这样的，这个分岔图看不出效果
回复 #95 octopussheng 的帖子
前辈好！你的这个算了多长时间？我前面给的图大约算了１．５个小时．:lol 6.5以上版本的计算速度相对而言要慢些
改成6.5程序不要变。 function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=;
clear;clc;close all;
omega=0.1:0.01:4;
for h=1:length(omega)
T=2*pi/omega(h);
=ode45('dddd',,,[],omega(h));
plot(omega(h),x(21000:100:end,2),'k.');hold on
end
回复 #98 无水1324 的帖子
谢谢前辈！不过程序里有个小变化才行．如红色标记：
function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=-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);
=ode45('dddd',,,[],omega(h));
plot(omega(h),x(21000:100:end,2),'k.');hold on
end 差不多吧，1小时左右
回复 #99 mechanic05 的帖子
负号的话是他自己的方程的问题了。我没有看到方程，就这能对程序作出判断 这个确实应该是方程的问题了，按照她这个方程，分岔图看来是效果不好的！我拿我那个方法再给她试试看！
回复 #102 octopussheng 的帖子
前辈好！我的方程是：
dx1/dt=x2
dx2/dt=－a*x1－b*x1^3－eta*x2+gamma*cos(omega*t)．本来就有负号的．谢谢！ 呵呵，看到你的方程了，我用我的方法把分岔图做了一下，第一个图是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)
第二个分岔图看的比第一个要好很多！
也请无水、咕噜噜这些高手看看，让他们帮你分析分析吧
回复 #104 octopussheng 的帖子
多谢前辈！我的程序是：function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=;
clear;clc;close all;
omega=0.01:0.01:4;
for h=1:length(omega)
T=2*pi/omega(h);
=ode45('dddd',,,[],omega(h));
plot(omega(h),x(10000:100:end,2),'k.');hold on
end