neweqq.m

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

function OUT=neweqq(t,X)

c1=1.766;c2=38.27;c3=4.914;c4=36.15;c5=7.789;c6=0.01;
a=10;b=40;c=2.5;k=1;h=4;

x(1)=X(1);x(2)=X(2);x(3)=X(3);x(4)=X(4);x(5)=X(5);x(6)=X(6);
x(7)=X(7);x(8)=X(8);x(9)=X(9);

Q=[ X(10),  X(19), X(28), X(37), X(46),X(55), X(64), X(73), X(82);
    X(11),  X(20), X(29), X(38), X(47),X(56), X(65), X(74), X(83);
    X(12),  X(21), X(30), X(39), X(48),X(57), X(66), X(75), X(84);
    X(13),  X(22), X(31), X(40), X(49),X(58), X(67), X(76), X(85);
    X(14),  X(23), X(32), X(41), X(50),X(59), X(68), X(77), X(86);
    X(15),  X(24), X(33), X(42), X(51),X(60), X(69), X(78), X(87);
    X(16),  X(25), X(34), X(43), X(52),X(61), X(70), X(79), X(88);
    X(17),  X(26), X(35), X(44), X(53),X(62), X(71), X(80), X(89);
    X(18),  X(27), X(36), X(45), X(54),X(63), X(72), X(81), X(90)];
    
    
    

dx(1)=x(2);
dx(2)=x(3);
dx(3)=-c4*x(3)-c5*x(2)-c6*x(1)+a*c1*x(6)-a*c1*x(3)+a*c2*x(5)-a*c2*x(2)+a*c3*x(4)-a*c3*x(1);

dx(4)=x(5);
dx(5)=x(6);
dx(6)=-c4*x(6)-c5*x(5)-c6*x(4)+b*c1*x(3)-k*c1*(x(3)*x(7)+2*x(2)*x(8)+x(1)*x(9))+b*c2*x(2)-k*c2*(x(2)*x(7)+x(1)*x(8))+b*c3*x(1)-k*c3*x(1)*x(7);

dx(7)=x(8);
dx(8)=x(9);
dx(9)=-c4*x(9)-c5*x(8)-c6*x(7)-c*c1*x(9)+h*c1*(2*x(1)*x(3)+2*x(2)^2)-c*c2*x(8)+2*h*c2*x(1)*x(2)-c*c3*x(7)+h*c3*x(1)^2;



DX1=[dx(1);dx(2);dx(3);dx(4);dx(5),dx(6);dx(7);dx(8);dx(9)];	

 J=[0,                 1,               0,                0,         0,      0,     0,    0,    0;
    0,                 0,               1,                0,         0,      0,     0,    0,    0;
    -c6-a*c3,          -c5-a*c2,        -c4-a*c1,         a*c3,      a*c2,   a*c1,  0,    0,    0;
    0,                 0,               0,                0,         1,      0,     0,    0,    0;
    0,                 0,               0,                0,         0,      1,     0,    0,    0;
    b*c3-k*c1*x(9)-k*c2*x(8)-k*c3*x(7),b*c2-2*k*c1*x(8)-k*c2*x(7), b*c1-k*c1*x(7), -c6,-c5,-c4,-k*c1*x(3)-k*c2*x(2)-k*c3*x(1),-2*k*c1*x(2)-k*c2*x(1),-k*c1*x(1);
    0,0,0, 0, 0, 0, 0,    1,    0;
    0,0,0, 0, 0, 0, 0,    0,    1;
    2*h*c1*x(3)+2*h*c2*x(2)+2*h*c3*x(1), 4*h*c1*x(2)+2*h*c2*x(1),2*h*c1*x(1),0,0,0,-c6-c*c3,-c5-c*c2,-c4-c*c1];

F=J*Q;

OUT=[DX1; F(:)];