[發(fā)明專(zhuān)利]一種基于MDCT系數(shù)進(jìn)行正弦頻率估計(jì)的方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201410325387.8
• 申請(qǐng)日: 2014-07-09
• 公開(kāi)（公告）號(hào): CN104111374B
• 公開(kāi)（公告）日: 2017-01-04
• 發(fā)明（設(shè)計(jì)）人: 頓玉潔;劉貴忠
• 申請(qǐng)（專(zhuān)利權(quán)）人: 西安交通大學(xué)
• 主分類(lèi)號(hào): G01R23/02
• 分類(lèi)號(hào): G01R23/02
• 代理公司: 西安通大專(zhuān)利代理有限責(zé)任公司61200
• 代理人: 陸萬(wàn)壽

1.一種基于MDCT系數(shù)進(jìn)行正弦頻率估計(jì)的方法，其特征在于，包括以下步驟：

1)構(gòu)建兩個(gè)關(guān)于頻率偏移與系數(shù)比值關(guān)系的查找表R_α(ξ)和R_β(ξ)；

2)獲取正弦信號(hào)一個(gè)分析幀的MDCT系數(shù)，記做X(k)，其中，k＝0,1,…,N，N為MDCT分析的子帶數(shù)；

3)查找所述一個(gè)分析幀的MDCT系數(shù)X(k)絕對(duì)值的極大值及其所在的子帶位置，分別記做X_max,i和p_i，其中，i＝1,...,l,l為檢測(cè)到的極大值的個(gè)數(shù)；

4)根據(jù)所述每個(gè)極大值X_max,i所對(duì)應(yīng)MDCT系數(shù)X(p_i)，及其相鄰的MDCT系數(shù)值X(p_i-2)、X(p_i-1)、X(p_i+1)以及X(p_i+2)，使用決策算法得到對(duì)該極大值所對(duì)應(yīng)分量進(jìn)行頻率估計(jì)的分支決策參數(shù)，其中，分支決策參數(shù)包括小系數(shù)跳轉(zhuǎn)因子a_i、正向比值I_0,i、反向比值I_1,i和小系數(shù)判定因子λ_i；

5)使用所述得到的分支決策參數(shù)a_i、I_0,i、I_1,i、λ_i以及MDCT變換所使用的窗函數(shù)類(lèi)型，通過(guò)查表和內(nèi)插的方法估算每一個(gè)頻率分量的正弦頻率值。

2.根據(jù)權(quán)利要求1所述的基于MDCT系數(shù)進(jìn)行正弦頻率估計(jì)的方法，其特征在于，步驟1)中，構(gòu)建兩個(gè)關(guān)于頻率偏移與系數(shù)比值關(guān)系的查找表包括：

(1)設(shè)置頻率偏移的步長(zhǎng)Δξ；

(2)根據(jù)所述頻率偏移的步長(zhǎng)Δξ，計(jì)算系數(shù)比值R_α(ξ)|_ξ＝nΔξ和R_β(ξ)|_ξ＝nΔξ，其計(jì)算公式分別如下：

R_α(ξ)|_ξ＝nΔξ＝A_H(nΔξ-1)/A_H(nΔξ+1)????????????(1)

R_β(ξ)|_ξ＝nΔξ＝A_H(nΔξ-2)/A_H(nΔξ+2)??????????(2)

其中：A_H(ξ)是MDCT變換所使用的窗函數(shù)h(n)的中心離散傅里葉變換系數(shù)的實(shí)部，ξ的取值范圍為[-1.5,-0.5]，n為頻率偏移數(shù)，按所需精度設(shè)定步長(zhǎng)Δξ，通過(guò)控制頻率偏移數(shù)n，獲得在ξ的取值范圍內(nèi)的系數(shù)比值R_α(ξ)|_ξ＝nΔξ和R_β(ξ)|_ξ＝nΔξ。

3.根據(jù)權(quán)利要求1所述的基于MDCT系數(shù)進(jìn)行正弦頻率估計(jì)的方法，其特征在于，步驟2)中，具體步驟如下：

(1)當(dāng)給定信號(hào)為時(shí)域信號(hào)樣本x(n)時(shí)，取連續(xù)2N個(gè)點(diǎn)作為一幀信號(hào)，對(duì)所述一幀信號(hào)加窗h(n)后計(jì)算其MDCT系數(shù)X(k)；

(2)當(dāng)給定信號(hào)為MDCT域壓縮數(shù)據(jù)，結(jié)合壓縮域數(shù)據(jù)解碼的過(guò)程，獲得分析子帶從低到高依次排布的MDCT系數(shù)X(k)。