duffing_success.asv

利用混沌振子检测微弱正弦信号

clear all
clc
k=0.3;%阻尼比
r=0.6;%周期摄动力幅度,0.535和0.536的时候效果比较明显，
         %而提高灵敏度为0.5352和05356的时候，混沌不是很复杂，有一段时间的周期
%D=0.1;
w=1;%周期摄动力频率
w2b=0;w2c=0;%积分器初值
xd=-1;x=-1;%初始条件
tfinal=500;%仿真时间
fs=1000;%取样频率
delt=1/fs;%取样间隔
npts=fs*tfinal;%仿真次数
xdv=zeros(1,npts);%dx/dt取样矢量
xv=zeros(1,npts);%x的取样矢量
tv=zeros(1,npts);
%load bpsk_tool
%for n=1:2
 %   phase=[0,pi];
for i=1:npts
    t=(i-1)*delt;%仿真时刻
    A=randn;
    xdd=-k*w*xd+w^2*(x-x^3+r*cos(w*t)+);
    w1b=xdd+w2b;
    w2b=xdd+w1b;%第一次梯形积分
    xd=w1b/(2*fs);%第一次积分的输出
    w1c=xd+w2c;
    w2c=xd+w1c;%第二次积分
    x=w1c/(2*fs);%第二次积分的输出
    xdv(1,i)=xd;%将第一次积分的结果存储
    xv(1,i)=x;%将第一次积分的结果存储
    tv(1,i)=t;
end
plot(tv,xv)
% xlabel('x(v)')
% ylabel('xd(v/s)')
% % legend('k=0.535,a=0,sigma=1',4)