chua_chua.m

来自「使用四阶龙泽库塔法进行蔡氏电路仿真matlab程序」· M 代码 · 共 14 行

function [x,Q]=chua_chua(x,aaa,bi)
R=2122.98;
C2=100*10^(-9);
h=0.01*(1-bi);a=h/2;aa=h/6;
   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*0.01*(1-bi)*R*C2/(100*0.01)+Q);
  % Q=2*pi*0.01*(1-bi)*R*C2/(100*0.01)+Q;