我的分叉图为什么是一条直线呢

wxfx1988 · 发表于 2013-8-30 16:34

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

请教大家，我画的分叉图为什么会是这样的一条直线呢

外激励是变化参数

下面是我的程序：
y0=[0;0];Omega=0.8883;T1=2*pi/Omega;tt=T1/400;max=500*T1;
g=0.1:0.005:2;
  for j=1:length(g)
   [t,x]=ode45('fsm',[0:tt:max],y0,[],g(j));
   plot(g(j),x(100000:400:end,1),'k.')
   xlabel('g');ylabel('x1');

  end
function fy=fsm(t,x,flag,g)
%a是阻尼；b是固有频率，beta是非线性项系数，g是外激励幅值，Omega是外激励频率
a=0.2;b=0.888523;beta=0.394737;Omega=0.8883;
fy=[x(2);-a*x(2)-b^2*x(1)-3*beta*x(1)^2-beta*x(1)^3-g*cos(Omega*t)*x(1)-g*cos(Omega*t)];
end

lihaitao123 · 发表于 2013-9-2 16:24

y0=[0;0];Omega=0.8883;T1=2*pi/Omega;tt=T1/400;max=500*T1;
g=0.1:0.005:2;
  for j=1:length(g)
   [t,x]=ode45('fsm',[0:tt:max],y0,[],g(j));
hold on;
   plot(g(j),x(100000:400:end,1),'k.')
   xlabel('g');ylabel('x1');

wxfx1988 · 发表于 2013-9-6 17:31

lihaitao123 发表于 2013-9-2 16:24

登录/注册后可看大图

y0=[0;0];Omega=0.8883;T1=2*pi/Omega;tt=T1/400;max=500*T1;
g=0.1:0.005:2;
for j=1:length(g)
...

谢谢啦，问题解决了

账号		自动登录	找回密码
密码			我要加入

[稳定性与分岔] 我的分叉图为什么是一条直线呢

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。