chua_initial.m

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

function [x0]=chua_initial(x,aaa)
h=0.01;a=h/2;aa=h/6;
x=[-0.03;0.6;-0.01];
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);
for k=2:400
   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);
end
x0=x;