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

wxfx1988

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

外激励是变化参数

下面是我的程序：
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

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

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)     
...

谢谢啦，问题解决了