大家帮忙看看这分岔图对吗

ymsncx

根据liliangbo 的duffing方程改的蔡氏电路的分岔图，大家看看给点意见吧。
function dx=yang(t,X)
global a;
x=x(1);
y=x(2);
z=x(3);
dx=zeros(3,1);
dx(1)=6.3*y-1.8*x-a*(abs(x-1)-abs(x+1));
dx(2)=0.7*x-0.7*y+z;
dx(3)=-7.0*y;


clear;
global a
range=[0.7: 0.05: 1.5];
period=2*pi;
k=0;
YY1=[];
step=2*pi/100;
for a=range
  y0=[0 0.01 0];

a

k=k+1;
%出去前面60个周期的数据，并将最后的结果作为下一次积分的初值。

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

[t,Y]=ode45('yang1',tspan ,y0);

y0=Y(end,:);

j=1;

for i=60:200

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

[t,Y]=ode45('yang1',tspan ,y0);

YY1(k,j)=Y(end,1);

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

y0=Y(end,:);

end
end
bifdata=YY1(:,end-51:end);
plot(range,bifdata,'k','markersize',1);
根据上面的程序做的分岔图











大家看看对吗 ？？
















[ 本帖最后由 ymsncx 于 2008-6-3 18:20 编辑 ]

liliangbiao

Why not take these points plotted?

xueer

对           应该画点图