%用希爾伯特算法對二次相關函數進行加權
clear;
m=100;
t=0:m-1;
s1=10*sin（0.5*t）;
D=4;???????????
%產生延時信號
s2=zeros（1m）;
for?i=（D+1）:m
s2（i）=s1（i-D）;
end
x=s1;
y=s2;
x1=fft（x2*m-1）;???%x1=fft（x）
x2=fft（y2*m-1）;???%x2=fft（y）
G12=conj（x1）.*x2;??
G11=conj（x1）.*x1;?
G=conj（G11）.*G12;
r1=real（fftshift（ifft（G12）））;?
r=real（fftshift（ifft（G）））;?
R=hilbert（r）;???%H12進行希爾伯特變換
Rn=imag（R）;?????%hilbert（r）計算出的結果其虛部是希爾伯特變化值
c=r-abs（Rn）;
[c1p]=mapminmax（c）;
[c2p]=mapminmax（r）;
[dl]=max（c1）
subplot（211）
plot（c1）
xlabel（‘Hilbert變換后的相關函數‘）;
subplot（212）
plot（c2）
xlabel（‘二次相關函數‘）;

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

?????文件????????657??2013-11-25?09:28??3_二次相關法程序\ercixiangguan.asv

?????文件????????648??2013-11-27?14:26??3_二次相關法程序\ercixiangguan.m

?????文件?????256082??2013-04-26?17:40??3_二次相關法程序\基于二次相關的時延估計方法研究.caj

?????目錄??????????0??2013-11-25?09:28??3_二次相關法程序

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

???????????????257387????????????????????4