自相关法信号检测.docxVIP

自相关法信号检测根据要求我们设定：1、余弦信号的频率为2HZ，抽样频率为100HZ；2、设定N(0,0.25)高斯白噪声，及噪声功率为0.25W；求Y（t）的自相关函数我们可得到：RY()=RS()+RN()+2RSN()因为S(t)与平稳高斯噪声N(t)相互独立。所以RSN()=0并且RS()=a*a/2*cos(ω+θ) RN()=N0/2信噪比s=a*a/(2*0.25)，其中a为原信号的振幅理论上我们可以看出Y（t）的自相关函数的图像在零点有一个冲激，其他都与原信号的自相关函数相一致，所以通过观察Y（t）的自相关函数的图像我们可以检测出原余弦信号的存在。在matlab仿真中得到了六个图：y=a*cos(2*pi*fc*x+p)N(0，0.25)的高斯白噪声叠加了高斯白噪声的cosx原余弦的信号自相关函数高斯白噪声的自相关函数叠加了高斯白噪声后的自相关函数通过改变振幅a的大小来改变信噪比的大小当信噪比s=0.005时：当信噪比s=0.02时：当信噪比s=0.5时：由以上三个matlab仿真图形，通过改变信噪比的大小最后检测出余弦信号。代码：fs=100;fc=2;a=0.1;s=a*a/(2*0.25); %信噪比为sp=unifrnd(0,2*pi);x=(0:1/fs:2);n=201;y1=a*cos(2*pi*fc*x+p); %原余弦信号，频率为2a=0;b=0.5; %均值为a，方差为b^2subplot(3,2,1);plot(x,y1,r);title(y=a*cos(2*pi*fc*x+p));ylabel(y);xlabel(x/20pi);grid;y2=a+b*rand(1,n); %高斯白噪声subplot(3,2,2);plot(x,y2,b);title(N(0，0.25)的高斯白噪声);ylabel(y);xlabel(x/20pi);grid;y=y1+y2; %加入噪声之后的信号subplot(3,2,3);plot(x,y,c);title(叠加了高斯白噪声的cosx);ylabel(y);xlabel(x/20pi);grid; [c,d]=xcorr(y1,unbiased);subplot(3,2,4);plot(d*0.01,c)title(原信号自相关函数); [e,f]=xcorr(y2,unbiased);subplot(3,2,5);plot(f*0.01,e)title(高斯白噪声的自相关函数); [j,m]=xcorr(y,unbiased);subplot(3,2,6);plot(m*0.01,j)title(叠加了高斯白噪声后的自相关函数);