I m not getting the exact bifurcation plot for matlab code. my bifurcation plot for the particular system is not showing the divergence to nonlinearity from linear region. what should be the error in my codes?

1 次查看（过去 30 天）

Rizwana Junaid 2011-10-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/19140-i-m-not-getting-the-exact-bifurcation-plot-for-matlab-code-my-bifurcation-plot-for-the-particular-s

% mlcmainfile

 clf
 clc
 global f 
 for f=0:0.01:50
   initialx(1)=0.8;initialx(2)=0.04;
 [t,x]=ode45('mlceq',[0 450],[initialx(1),initialx(2)]);
 plot(f,x(440:450,1),'r');
 hold on
 end
 xlabel('f');ylabel('x')
 axis ([0 50 -1.2 1.2])

--------------------------------------------------------------------------------------------------------------------------------------

 % filename:mlceq(t,x)
 function s=mlceq(t,x)
 w=0.7;
 global f
 xdot=1.65* mlcg(x(1))-0.55*x(1)+x(2);
 ydot=-x(2)-x(1)+f*sin(w*t);
 s=[xdot;ydot];
 % % xdot =1.65*G(x)-0.55*x+y;
 % % ydot =-y-x+fsin(w*t);
-------------------------------------------------------------------
 % Filename: mlcg(x)
 function s = mlcg(x) 
 c=0.52;
 if (x(1))>c
          s=c;
 else
     if (x(1))<-c
           s=-c;
    else
           s=x(1);
     end
 end
 end