neweq.m

研究了利用混沌相变进行弱信号检测的理论及仿真试验( 对基于)振子初值敏感性检测弱信号的方法分 析后指出

function OUT=neweq(t,X)

a=43; b=9; c=-5.75; d=12; e=8;

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

Q=[ X(6),  X(11), X(16), X(21), X(26);
    X(7),  X(12), X(17), X(22), X(27);
    X(8),  X(13), X(18), X(23), X(28);
    X(9),  X(14), X(19), X(24), X(29);
    X(10), X(15), X(20), X(25), X(30)];
    
    
    

dx(1)=a*(x(2)-x(1))+x(2)*x(3)*x(4);
dx(2)=b*(x(1)+x(2))-x(1)*x(3)*x(4)*x(5);
dx(3)=-c*(x(3)-x(5))+x(1)*x(2);
dx(4)=-d*(x(4)-x(3))+x(1)*x(2)*x(5);
dx(5)=-e*x(4)*x(5)+x(1)*x(2)*x(3)*x(4);



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

 J=[-a,                a+x(3)*x(4),     x(2)*x(4),        x(2)*x(3),               0;
    b-x(3)*x(4)*x(5),  b,               -x(1)*x(4)*x(5),  -x(1)*x(3)*x(5),         -x(1)*x(3)*x(4);
    x(2),              x(1),            -c,               0,                       c;
    x(2)*x(5),         x(1)*x(5),       d,                -d,                      x(1)*x(2);
    x(2)*x(3)*x(4),    x(1)*x(3)*x(4),  x(1)*x(2)*x(4),   -e*x(5)+x(1)*x(2)*x(3),  -e*x(4)];
  
F=J*Q;

OUT=[DX1; F(:)];