matlab画虫口模型分岔图,Matlab画图片碰磨转子系统中的分岔图，poincare映射图？？？...

图片中是转子碰磨故障的动力学模型，应用龙格库塔法求解，在画分岔图和poincare截面图，画出来不对，不知道为什么，下面是我自己编的程序，请大侠们帮忙看看我错在哪里了？

global omega;              %以转速为变化画分岔图

range=[1.1:0.001:1.35];

k=0;

YY1=[];

step=2*pi/100;  %步长。

for omega=range

period=2*pi/omega;

x0=[0;0;0.001;0];

omega

k=k+1;

% discard the first 60 periodic data;

%除去前面60个周期的数据，并将最后的结果作为下一次积分的初值

tspan=[0:step:60*period];

[t,x]=ode45(@fangcheng,tspan,x0);

y0=x(end,:);

j=1;

for i=60:200

tspan=[i*period:step:(i+1)*period];

[t,x]=ode45(@fangcheng,tspan,x0);

YY1(k,j)=x(end,1);   % get the omega data from every period end

j=j+1;               %取出每一个周期内的第一个解的最后一个值。

y0=x(end,:);

end

<