liueq.m

来自「研究了利用混沌相变进行弱信号检测的理论及仿真试验( 对基于)振子初值敏感性检测弱」· M 代码 · 共 29 行

function OUT=liueq(t,X)

a=10; b=40; c=2.5; k=1; h=4;

x(1)=X(1);x(2)=X(2);x(3)=X(3);

Q=[ X(4),  X(7),  X(10);
    X(5),  X(8),  X(11);
    X(6),  X(9),  X(12)];
    
    
    

dx(1)=a*(x(2)-x(1));
dx(2)=b*x(1)-k*x(1)*x(3);
dx(3)=-c*x(3)+h*x(1)^2;




DX1=[dx(1);dx(2);dx(3)];	

 J=[-a,        a,           0;
    b-k*x(3),  0,         -k*x(1);
    2*h*x(1),  0,      -c;];
  
F=J*Q;

OUT=[DX1; F(:)];