关于汽车前轮转向失稳的问题（鞍结分岔图的matlab程序）

victorpn

现在附上相图及相图的程序，希望大牛们能给出平衡点分岔图的编写程序！谢谢
function bifurcation
%解方程
global deltaf
      deltaf=0.015;
%画图
options = odeset('RelTol',1e-4);
for p=-0.1:0.05:0.1;
    for q=-0.4:0.1:0.4;
    x0=[p,q];
    [t,x]=ode45(@f,[0,30],x0,options);
   plot(x(:,1),x(:,2));
       hold on
    axis([-0.15,0.15,-0.6,0.6]);
   pause(0.01);
end
end



function dx=f(t,x)
dx=zeros(2,1);
global deltaf
af=1.2;ar=1.3;Iz=3000;m=1500;u=20;Bf=11.275;Cf=1.56;Df=-2574.7;Ef=-1.9990;
Br=18.631;Cr=1.56;Dr=-1749.7;Er=-1.7908;
Jf=x(1)+atan((af*x(2)*cos(x(1)))/u)-deltaf;
Jr=x(1)-atan((ar*x(2)*cos(x(1)))/u);
Ff=Df*sin(Cf*atan(Bf*(1-Ef)*Jf+Ef*atan(Bf*Jf)));
Fr=Dr*sin(Cr*atan(Br*(1-Er)*Jr+Er*atan(Br*Jr)));
dx(1)=(Ff+Fr-m*u*x(2))/(m*u);
dx(2)=((af*Ff-ar*Fr)*cos(x(1)))/Iz;

下面的分岔图程序怎么编写

分岔图

牛小贱

坐等大仙出来！！

victorpn

大神至今没出来

牛小贱

楼主可以到【Matlab专区】问一下，看有没有人解答