chua.m

使用四阶龙泽库塔法进行蔡氏电路仿真matlab程序

function [xq,bi]=chua(x,aaa)
h=0.01;a=h/2;aa=h/6;
R=2122.98;
C2=100*10^(-9);
symbol_no=250;
xx=[];
for j=1:symbol_no
   k0=chua_map(x,1,aaa);
   x1=x+k0*a;
   k1=chua_map(x1,1,aaa);
   x1=x+k1*a;
   k2=chua_map(x1,1,aaa);
   x1=x+k2*h;
   k3=chua_map(x1,1,aaa);
   x=x+aa*(k0+2*(k1+k2)+k3);%+10^(-11)*sin(2*pi*j*0.01*R*C2/(100*0.01)+O);
   xx=[xx x];
end
for j=1:250-1
    if (xx(2,j)>0)&&(xx(2,j+1)<0)
       x2q=0;
       x1q=(xx(1,j)*xx(2,j+1)-xx(2,j)*xx(1,j+1))/(xx(2,j+1)-xx(2,j));
       x3q=(xx(3,j)*xx(2,j+1)-xx(2,j)*xx(3,j+1))/(xx(2,j+1)-xx(2,j));
       xq=[x1q;x2q;x3q];
       bi=xx(2,j)/(-xx(2,j+1)+xx(2,j));
      % Q=2*pi*(j+bi)*0.01*R*C2/(100*0.01)+O;
    end
end