main.m

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

% 把Henon映射产生的混沌现象控制到周期-1轨道上一个不动点上
% 不动点：xF=yF=(B-1+((1-B).^2+4*A).^0.5)/0.5
%  -(1-1/ju)*hu*g*d<xu<(1-1/ju)*hu*g*d  % 控制启动条件
% 画出x随着n的变化图示

R1=3.3*10^3;
R2=22*10^3;
R3=22*10^3;
R4=2.2*10^3;
R5=220;
R6=220;

Ga=-1/R1-1/R4;
Gb=1/R3-1/R4;

L=23.3*10^(-3);
C1=10.25*10^(-9);
C2=100*10^(-9);
R=2122.98;
a=C2/C1;
b=(R^2)*C2/L;
m0=Ga*R;
m1=Gb*R;
A=[-a*(m1+1) a 0;1.0 -1.0 1.0;0 -b 0];
[V,D]=eig(A);
for j=1:3
    if real(D(j,j))<-1
       ju=D(j,j);
       eu=V(:,j);
       one=[zeros(1,j-1) 1 zeros(1,3-j)];
    end
end
hu=real(one/V);
g=det(eye(3,3)-A)^(-1)*[1;1;1];
x1f=5.7795;
x2f=0;
x3f=-5.9766;
x0=0.05*randn(3,1);
x=chua_initial(x0,-b);
p=zeros(3,300);
for n=1:300
      if n==1
          Q=0;xq=[x1f;x2f;x3f];
      end
      xpre=xq;
      diag=[x(1,1)-x1f;x(2,1)-x2f;x(3,1)-x3f];              
      p(n)=(ju*hu*diag)/((ju-1)*(hu*g)); 
      [xq,Q,bi]=chua(x,-b,Q);
      xq=xq+[p(n);0;0];
      [x,Q]=chua_chua(xq,-b,Q,bi);
end