duffingpinlv.m

来自「利用duffing振子检测微弱信号频率的matlab程序。可以检测与周期策动力频」· M 代码 · 共 16 行

%Duffing检测频率系统方程
%方程为 x”(t)+kx'(t)-x(t)+x(t)^3=r*cos(t)+a*cos((1+w)*t+p)
%其中系统周期策动力的频率设为1，故待测信号频率为1+w，w表示频率差，w不大于0.03时应该能产生间歇性混沌
%r表示周期策动力的幅度，在检测时应设为临界值，当步长为0.04时，r=0.0815
%a表示待测信号的幅度，q是待测信号与周期策动力的相位差，因为很小，所以对系统影响不大，这里设为0

function dx=f(t,x)
k=0.5;
r=0.7375;
w=0.03;
a=0.02;
q=0;
dx(1)=x(2);
dx(2)=-k*x(2)+x(1)-x(1)^3+r*cos(t)+a*cos(t+w*t+q)+wgn(1,1,-13.9794);
dx=[dx(1);dx(2)];