变形蔡氏分岔图，图形不太对劲，求帮助。。谢谢。。。
function dx=bianxing(t,x)global alfa;
dx=zeros(3,1);
dx(1)=alfa*(x(1)/3-x(1)*abs(x(1))/3+x(2));
dx(2)=x(1)-x(2)+x(3);
dx(3)=-16*x(2);
clear;
global alfa
%取参数alfa变化范围及增长幅度
range=;
period=2*pi;
k=0;
YY1=[];
step=2*pi/100;
for alfa=range
y0=;
alfa
k=k+1;
%除去前面60个周期的数据，并将最后的结果作为下一次积分的初值。
tspan=;
=ode45('bianxing',tspan ,y0);
y0=Y(end,:);
j=1;
for i=60:200
tspan=;
=ode45('bianxing',tspan ,y0);
YY1(k,j)=Y(end,1);
j=j+1;
%取出每一个周期内的第一个解的最后一个值。
y0=Y(end,:);
end
end
plot(range,YY1,'k','markersize',1);
xlabel('alfa');
ylabel('Y');
title('变形蔡氏电路分岔图');
页:
[1]