matlab蔡氏电路相图

clear all;
hold off;

x（1）=0.02;y（1）=0.01;z（1）=0.03;
m0=-1/7;m1=2/7;a=10;b=15;
T=1e-2;
N=5e4;
for i=1:N
    i
    t（i）=i;
    f（i）=m1*x（i）+0.5*（m0-m1）.*（abs（x（i）+1）-abs（x（i）-1））;
    x（i+1）=x（i）+T*a*y（i）-T*a*f（i）;
    y（i+1）=T*x（i）+（1-T）*y（i）+T*z（i）;
    z（i+1）=-T*b*y（i）+z（i）;
end

n1=N/2;n2=N;
figure（1）
plot（t（n1:n2）x（n1:n2）‘b-‘）;
xla