chaios.m

来自「应用于混沌电路研究」· M 代码 · 共 7 行

lnlnit=simget('chaios','initialstate')
[t,x,y]=sim('chaios',100);
alfa=9;beta=14+2/7;k1=0.25;k2=0.1*3;        % 给出4个电路参数alfa,beta,k1,k2
A=[0 alfa 0;1 -1 1;0 -beta 0];              % 初始化矩阵A;
opts=simset('chiaos',a);
[tt,xx,yy]=sim('chiaos',100,opts);
plot(t,x(:,1),':b',tt,xx(:,1),'r')