lk1.m

matlab时频分析配套源码第14章的程序

function ydot=LK1(t,n,y)
global w;
%参数值
a0=3.1815;
a1=0.2063;
a2=0.0108;
b0=3.223;
b1=0.59;
b2=0.0327;
G0=2.56*10^(-4);
n0=0.05;
L=20;
a=2*pi*n0;
b=2*pi*sqrt(G0);
temp=2*pi*sqrt(G0*v(t));
%状态方程组
ydot(1)=-a*v(t)*y(1)+temp*w(n);
ydot(2)=-a*v(t)*y(2)+temp*(a1-a2*b1/b2)*y(5)+temp*(a0-a2*b0/b2)*y(6)+temp*a2/b2*w(n);
ydot(3)=-a*v(t)*y(3)+temp*(1-L/(2*v(t)))*y(7)+temp*w(n);
ydot(4)=-a*v(t)*y(4)+temp*(a1-a2*b1/b2)*y(5)+temp*(a0-a2*b0/b2)*y(6)+...
    temp*(1-L/(2*v(t)))*y(7)+temp*a2/b2*w(n);
ydot(5)=-b1/b2*y(5)-b0/b2*y(6)+w(n)/b2;
ydot(6)=y(5);
ydot(7)=-6*v(t)/L*y(7)-12*(v(t))^2/L^2*y(8)+12*(v(t))^2/L^2*w(n);
ydot(8)=y(7);
ydot=[ydot(1) ydot(2) ydot(3) ydot(4) ydot(5) ydot(6) ydot(7) ydot(8)];
%输出
ydot=ydot';