%信號程序
clear all
clc
j=0;
%tao=0.5;
for t=0.01:0.01:20
    j=j+1;
    x(j)=cos(t);%信號為幅度為2的余弦信號
    %%信號延時(shí)為幅度為2的余弦信號
end
ll=length(x)
n=normrnd(0,1,1,ll);%噪聲為均值為0方差為1的高斯白噪聲

xt=x+n;
%xttao=xtao+n;
[Rn,c]=xcorr(n,n);%做x(t)的自相關(guān)
[Rxn,c]=xcorr(xt,xt);%做x(t)的自相關(guān)
RR=Rxn-Rn;
subplot(2,1,1);
plot(c/100,Rxn);
xlabel('時(shí)間差值');
ylabel('相關(guān)幅度')
title('信號加噪聲的自相關(guān)函數(shù)');

subplot(2,1,2);
plot(c/100,RR)
xlabel('時(shí)間差值');
ylabel('相關(guān)幅度')
title('提取信號的自相關(guān)函數(shù)');