by zzb
- 取兩路時(shí)間序列s0和s1;
- 取s0靠中間約80%的部分(克服帶通效應(yīng)),取s1全部;
- 在FFT前,將兩路數(shù)據(jù)擴(kuò)展到一定位數(shù)(FFT內(nèi)插);

原始時(shí)間序列S0(左上)和S1(右上),以及它們對(duì)應(yīng)的FFT變換序列(下)
- 互相關(guān)函數(shù)求峰值位置p1,希爾伯特函數(shù)求取最大最小值間過零點(diǎn)的位置p2。實(shí)際上p2要比p1更準(zhǔn)確。

上述兩時(shí)間序列的互相關(guān)函數(shù)(左上)和希爾伯特變換(右上),以及它們對(duì)應(yīng)的FFT逆變換(下圖)
見程序
PLOTRPO@mbp