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

xt=x+n;



i=0;
for t=0.01:0.01:20
    i=i+1;
    y(i)=cos(t);%信號為幅度為2的正弦信號
    %%信號延時為幅度為2的余弦信號
end

%xttao=xtao+n;
%[Ry,c]=xcorr(y,y)%做y(t)的自相關(guān)
%[Rn,c]=xcorr(n,n);
[Rxy,c]=xcorr(xt,y);%做x(t)的自相關(guān)
RR=Rxy;

plot(c/100,Rxy);
xlabel('時間差值');
ylabel('相關(guān)幅度')
title('信號的互相關(guān)函數(shù)');