clear?all;
clc;
close?all;
N=512;%序列長度
Fs=50000000;%采樣率為50Mhz
n=0:N-1;
d=20;
F=800000;%信號頻率0.8Mhz
a1=5;
x1=a1*cos（2*pi*F*n/Fs）;
x2=zeros（1N）;
for?i=（d+1）:N??
x2（i）=x1（i-d）;
end
%加噪聲
%randn（‘state‘sum（100*clock））;
x1n=awgn（x15）;
x2n=awgn（x25）;
subplot（411）;
plot（nx1）;
title（‘信號1‘）;
grid?on;
subplot（412）;
plot（nx1n）;
title（‘加噪聲的信號1‘）;
grid?on
subplot（413）;
plot（nx2）;
title（‘信號2‘）;
grid?on;
subplot（414）;
plot（nx2n）;
title（‘加噪聲的信號2‘）;
grid?on


X1=fft（x1n2*N-1）;???%x1=fft（x）
X2=fft（x2n2*N-1）;???%x2=fft（y）

%求基本互相關函數
G12=conj（X1）.*X2;???%信號x1與x2的互功率譜密度函數
R=fftshift（ifft（G12））;%對互功率譜反FFT變換得到互相關函數

%求出最大值max及最大值所在的位置
[maxklocationk]=max（R）;%求最大值，及最大值對應的位置
locationk
%d=location-N/2-1
dk=locationk-N%算出延遲了幾個點
Delayk=dk/Fs%求得延遲時間

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2799??2013-09-21?17:38??GCCbijiao5db.m

-----------?---------??----------?-----??----

?????????????????2799????????????????????1