鞍结分岔图的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; 请问分岔图用什么做、 分岔图

gghhjj

简单来说分岔图就是不同分岔参数下庞加莱点的集合，论坛有很多类似的程序可以参考