chua_map.m

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

function [x]=chua_map(xx,time_variable,aaa)
m0=-1.6083;
m1=-0.8685;
if xx(1)>=1
   hx=m1*xx(1)+m0-m1;
elseif abs(xx(1))<=1
   hx=m0*xx(1);
else
   hx=m1*xx(1)-m0+m1;
end
A=[-9.7561 9.7561 0;1.0 -1.0 1.0;0 aaa 0];
x=A*xx;
x=x+[-9.7561*hx 0 0]';