关于汽车前轮转向失稳的问题(鞍结分岔图的matlab程序)
现在附上相图及相图的程序,希望大牛们能给出平衡点分岔图的编写程序!谢谢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=;
=ode45(@f,,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;
下面的分岔图程序怎么编写
坐等大仙出来!!{:{39}:} 大神至今没出来 楼主可以到【Matlab专区】问一下,看有没有人解答{:{26}:}
页:
[1]