[發(fā)明專利]一種基于互譜相位擬合的時延估計方法

在審

• 申請?zhí)枺?/td>: 202010633065.5
• 申請日: 2020-07-02
• 公開（公告）號: CN111856401A
• 公開（公告）日: 2020-10-30
• 發(fā)明（設(shè)計）人: 唐嵐;李保偉;張興敢
• 申請（專利權(quán)）人: 南京大學(xué)
• 主分類號: G01S5/22
• 分類號: G01S5/22
• 代理公司: 蘇州威世朋知識產(chǎn)權(quán)代理事務(wù)所(普通合伙) 32235
• 代理人: 楊林潔

1.一種基于互譜相位擬合的時延估計方法，其特征在于，包括以下步驟：

1)利用麥克風(fēng)陣列采集聲源信號，計算每兩路接收信號相對時延的整數(shù)部分；

2)根據(jù)整數(shù)部分時延，重新抽樣，構(gòu)建延遲點(diǎn)數(shù)小于1的兩路接收信號，求出互功率譜；

3)對所述互功率譜進(jìn)行頻域?yàn)V波，去除噪聲、雜波干擾，計算濾波后的互譜相位；

4)利用互譜相位跟頻率的線性關(guān)系計算每個頻點(diǎn)的時延，擬合得到分?jǐn)?shù)部分時延估計值，綜合得到整體時延。

2.根據(jù)權(quán)利要求1所述一種基于互譜相位擬合的時延估計方法，其特征在于：所述步驟1)包括以下步驟：

11)利用麥克風(fēng)陣列采集聲源信息，在同一時刻抽樣長度為N的兩路接收信號其中x₁[n]表示第1個麥克風(fēng)接收信號，n表示信號x[n]的索引，x₂[n]表示第2個麥克風(fēng)接收信號；

12)根據(jù)傳統(tǒng)GCC時延估計算法計算相對時延整數(shù)部分m。

3.根據(jù)權(quán)利要求2所述一種基于互譜相位擬合的時延估計方法，其特征在于：所述步驟2)包括以下步驟：

21)信號不變，將信號抽樣位置順延m個采樣單元，重新抽樣長度N得到

22)對信號分別進(jìn)行FFT轉(zhuǎn)換到頻域，得到k表示X[k]的索引，再將取共軛與相乘得到互功率譜

4.根據(jù)權(quán)利要求2所述一種基于互譜相位擬合的時延估計方法，其特征在于：所述步驟3)包括以下步驟：

31)將聲源信號頻率范圍[f₁,f₂]除以頻率分辨率f₀、取整，得到相應(yīng)的正頻率部分頻點(diǎn)k值范圍[k₁,k₂]，其中f₁表示頻率下限，f₂表示頻率上限；

32)設(shè)計一個長度為N的濾波器，系數(shù)在范圍內(nèi)的值取1，其余為0，而后與互譜進(jìn)行點(diǎn)乘運(yùn)算得到新互譜，實(shí)現(xiàn)頻域?yàn)V波；

33)利用angle函數(shù)獲得濾波后的互譜相位N'為濾波后保留的互譜長度，表示第k個頻點(diǎn)的相位。

5.根據(jù)權(quán)利要求2所述一種基于互譜相位擬合的時延估計方法，其特征在于：所述步驟3)包括以下步驟：

31)將聲源信號頻率范圍[f₁,f₂]除以頻率分辨率f₀、取整，得到相應(yīng)的正頻率部分頻點(diǎn)k值范圍[k₁,k₂]，其中f₁表示頻率下限，f₂表示頻率上限；

32)區(qū)分正負(fù)頻率，直接將互譜在頻點(diǎn)范圍[k₁,k₂]、[-k₂,-k₁]以外的值置0的方式達(dá)到濾波效果；

33)利用angle函數(shù)獲得濾波后的互譜相位N'為濾波后保留的互譜長度，表示第k個頻點(diǎn)的相位。

6.根據(jù)權(quán)利要求2所述一種基于互譜相位擬合的時延估計方法，其特征在于：所述步驟4)包括以下步驟：

41)根據(jù)互譜相位與頻率的線性關(guān)系，將除以得到每個頻點(diǎn)的時延值p[k]表示第k個頻點(diǎn)的時延值；

42)對進(jìn)行最小二乘線性擬合，確定兩路接收信號分?jǐn)?shù)部分時延估計值

43)與步驟1)得到的整數(shù)部分時延m求和得到整體時延值Δn。