下面的分岔图是否正确

咕噜噜

这一行function ff=www(t,x,flag,w,b)
把w去掉，把w=2.0;写到主程序里面
ff=[dx(1);dx(2);dx(3);dx(4)]我觉得还是写成ff＝zeros(4,1)的好，写在dx前面
就没问题了

中原

为什么我发新帖时出现，“对不起，您填写的内容(如签名、帖子、短消息等)包含不良内容而无法提交，请返回修改。
”，我晕，帖子的名称为“什么是分岔”，怎么不良了，请版主说明一下:@L

chuandong418

按你所说的做了修改，但是还是不行,居然matlab没反应了!

无水1324

没有定义b就是在你画# 25 楼图的那个点‘.’ 里面加一个k就全部是黑色的了

chuandong418

555,还是不知道问题在那里,还是出现同样的错误，望指点,谢谢!

咕噜噜

:@o 修改之后我这里能运行啊
b=2.0:0.1:15.0;
options = odeset('RelTol',1e-6,'AbsTol',[1e-6 1e-6 1e-6 1e-6]);
for  n=1:length(b);
w=2.0;
T=2*pi/w;
ts=[0:T/100:100*T];
x0=[0.001 0 0.001 0];
[t,X]=ode45('www',ts,x0,options,b(n));
figure(1)
plot(b(n),X(5000:100:10000,1),'.');
xlabel('\fontsize{18}\omega');
ylabel('\fontsize{18}x');grid
hold on
figure(2)
plot(b(n)*30/pi,X(5000:100:10000,2),'.');
xlabel('\fontsize{18}\omega');
ylabel('\fontsize{18}y');grid
hold on
end


function ff=www(t,x,flag,b)
ff=zeros(4,1);
g=9800;
E=0.12;
f=0.12;
a=0.5;
%b=10.8;
q=0.16;
f0=25;
u=0.136/q;
w=2.0;
e=sqrt(x(1)^2+x(2)^2);
G=g/((2*pi*f0)^2*q);
ff=zeros(4,1);
dx(1)=x(3);
dx(2)=x(4);
dx(3)=-2*E*x(3)-x(1)-a*(x(1)^2+x(2)^2)*x(1)-b*(1-1/e)*(x(1)-f*x(2))+u*w^2*cos(w*t);
dx(4)=-2*E*x(4)-x(2)-a*(x(1)^2+x(2)^2)*x(2)-b*(1-1/e)*(f*x(1)+x(2))+u*w^2*sin(w*t)-G;
%ff=[dx(1);dx(2);dx(3);dx(4)]

octopussheng

我试了一下，算是能算，不过结果有问题的哦！呵呵

chuandong418

为什么不需要ff=[dx(1);dx(2);dx(3);dx(4)]？

无水1324

应该是书写错误，必须要ff=[dx(1);dx(2);dx(3);dx(4)]的

octopussheng

这条语句确实是可要可不要的，因为已经在这里定义好了：
dx(1)=x(3);
dx(2)=x(4);
dx(3)=-2*E*x(3)-x(1)-a*(x(1)^2+x(2)^2)*x(1)-b*(1-1/e)*(x(1)-f*x(2))+u*w^2*cos(w*t);
dx(4)=-2*E*x(4)-x(2)-a*(x(1)^2+x(2)^2)*x(2)-b*(1-1/e)*(f*x(1)+x(2))+u*w^2*sin(w*t)-G;

chuandong418

刚刚又计算了下以参数Ｕ为变量的分叉图，得出的结果也是一条直线，现在我都在怀疑这种方法是不是只是使用与w为变量的分叉，其他的不行，是不是要换一种方法来计算呢？还望大家帮忙，好苦恼！

咕噜噜

看看你的初始条件改一下

chuandong418

修改初值也不行哪，急死我了，大家帮忙哦！５５５５５

octopussheng

有没有可能系统不会产生分岔啊？如果是这样的话，只有改系统参数了！

咕噜噜

即便是不分叉也不应该是一直线，奇怪