蔡氏电路的分岔图

lizhijun

各位大侠，这是我的蔡氏电路分岔图程序，为何得不到正确结果，请诸位多指教
function dx=chua(t,x)
global a；
b=15;
m0=-1/7;
m1=2/7;
f=m1*x(1)+0.5*(m0-m1)*(abs(x(1)+1)-abs(x(1)-1));
dx=[a*(x(2)-f);x(1)+x(3)-x(2);-b*x(2)];


clear;
global a
%取参数a变化范围及增长幅度
range=[8:0.01:9];
period=2*pi;
k=0;
YY1=[];
step=2*pi/100;
for a=range
  y0=[0.1 0 0];
  a
  k=k+1;
  %除去前面60个周期的数据，并将最后的结果作为下一次积分的初值。
  tspan=[0:step:100*period];
  [t,Y]=ode45('cai',tspan ,y0);
  y0=Y(end,:);
  j=1;
  for i=100:200
    tspan=[i*period:step:(i+1)*period];
    [t,Y]=ode45('cai',tspan ,y0);
    YY1(k,j)=Y(end,1);
    j=j+1;
    %取出每一个周期内的第一个解的最后一个值。
    y0=Y(end,:);
  end
end
plot(range,YY1,'.','markersize',1);
xlabel('a');
ylabel('Y');
title('蔡氏电路分岔图');