語音降噪方法及語音采集設備技術

本發明專利技術提供一種語音降噪方法及語音采集設備，其中，通過一濾波組件，對一濾波形成的消噪語音信號，以及一濾波形成的參考噪聲信號進行頻域轉換并濾波，以濾波輸出、對應所述參考噪聲信號的原始噪聲信號、以及對應所述消噪語音信號的含噪語音信號，分別控制形成所述參考噪聲信號的第一抗串擾濾波單元，以及形成所述消噪語音信號的第二抗串擾濾波單元進行濾波系數更新。其技術方案的有益效果是：通過僅使用一對自適應濾波器，并通過后濾波器控制該一對自適應濾波器的濾波系數更新，從而簡化了的結構；通過對語音信號和噪聲信號的濾波系數更新權重分別進行有針對性的調整，實現優秀的抗串擾和降噪效果。

【技術實現步驟摘要】

本專利技術涉及通訊領域，尤其涉及語音降噪方法及語音采集設備。
技術介紹
隨著通訊技術的革新，人們對于語音通訊的質量要求不斷提高。其中，移動通訊終 端的抗噪性能是語音通訊質量的重要關注點。在經歷了單麥克風系統通過單通道語音增 強方案降低噪聲，提高語音質量的階段以后，越來越多的移動通訊終端開始配置雙麥克風 (Dual Mic)系統。這種雙麥克風系統通常采用主次麥克風的結構，它將一個主麥克風放置 在手機下端，靠近嘴的位置，用于接收含噪語音信號，另一個次麥克風放置在手機上端的背 部或頂部，靠近耳朵的位置，用于接收以噪聲為主的參考噪聲信號，如圖1所示。 基于這種物理結構，自適應噪聲消除（Adaptive Noise Cancellation，ANC)可以 利用含噪語音信號和參考噪聲信號這兩個信號來進行分析和計算，得到干凈的語音。但是， 采用目前的技術方案在實際應用中起到的降噪效果并不是很理想，輸出的語音還是有存在 較大的失真。 因此，如何進一步提高通話質量一直為本領域技術人員所致力研究的方向。
技術實現思路
針對現有的自適應噪聲消除技術中存在的上述問題，本專利技術旨在提供一種結構簡 單，并可有效實現有效的抗串擾以及提高語音質量的語音降噪方法及語音采集設備。 具體技術方案如下： -種語音降噪方法，其中， 通過一濾波組件，對一濾波形成的消噪語音信號，以及一濾波形成的參考噪聲信 號進行頻域轉換并濾波，以濾波輸出、對應所述參考噪聲信號的原始噪聲信號、以及對應 所述消噪語音信號的含噪語音信號，分別控制形成所述參考噪聲信號的第一抗串擾濾波單 元，以及形成所述消噪語音信號的第二抗串擾濾波單元進行濾波系數更新。 優選的，所述濾波組件采用離散傅立葉變換，或離散余弦變換，或改進余弦變換對 所述消噪語音信號進行頻域轉換得到第一幅度譜；和/或 所述濾波組件采用離散傅立葉變換，或離散余弦變換，或改進余弦變換對所述參 考噪聲信號進行頻域轉換得到第二幅度譜。 優選的，采用離散傅立葉變換得到所述第一幅度譜的公式為： 其中，Eal為所述第一幅度譜，e Jt)為所述參考噪聲信號，e為自然對數的底數，j 為虛數單位，j = (-1)α5, w_(n)為窗函數，η = 1，2,…Μ, Μ為頻域分析的時間窗長度；和 /或 采用離散傅立葉變換得到所述第二幅度譜的公式為： 其中，Ea2為第二幅度譜，e 2(t)為所述消噪語音信號，e為自然對數的底數，j為虛 數單位，j= (_l)a5，w_(n)為窗函數，n = l，2,…M，M為頻域分析的時間窗長度。 優選的，所述窗函數w_(η)為矩形窗，或者正弦窗，或者漢寧窗，或者海明窗，或者 Tukey 窗。 優選的，Μ的取值范圍為fs/100/2〈M〈0. 4fs，其中fs為采樣頻率。 優選的，所述濾波組件提供一后濾波器，所述后濾波器為具有Μ點沖激響應的頻 域維納濾波器，通過下式計算所述頻域維納濾波器的解： 其中，G為所述維納濾波器的解，Psl為所述消噪語音信號的自功率譜，Ρ "為噪聲的 自功率譜；或者 所述濾波組件提供一后濾波器，所述后濾波器為具有Μ點沖激響應的頻域維納濾 波器，通過下式計算所述頻域維納濾波器的解： SNRt= a G t JSNlUkL i+d-a )max(SNRPt-l，0) 其中，G為所述維納濾波器的解，SNR為信噪比，SNRP為先驗信噪比，Psl為所述消 噪語音信號的自功率譜，P n為噪聲的自功率譜，〇〈 α〈1。 優選的，所述消噪語音信號的自功率譜Psl通過下式計算： Pslt= |Ealt|2， 其中，Eal為所述消噪語音信號經頻域轉換后獲得的第一幅度譜。 優選的，所述噪聲的自功率譜^通過下式進行估計： 其中，η。n2, n3為平滑因子，〇〈 n p n2, n3〈i; Ps為所述消噪語音信號的自功率譜，或者P 3為所述消噪語音信號的自功率譜與所 述參考噪聲信號的自功率譜兩者之間的較小值。 優選的，Ps為所述消噪語音信號的自功率譜時，P s通過下式計算： Pst= |Ealt|2； 其中，Eal為所述消噪語音信號經頻域轉換后獲得的第一幅度譜。 優選的，Ps為所述消噪語音信號的自功率譜與所述參考噪聲信號的自功率譜兩者 之間的較小值時，ps通過下式計算： Pst=min(|Ealt| 2, |Ea2t|2); 其中， Eal為所述消噪語音信號經頻域轉換后獲得的第一幅度譜； Ea2為所述參考噪音信號經頻域轉換后獲得的第二幅度譜。 優選的，提供一語音激活檢測單元，所述語音激活檢測單元連接所述濾波組件； 所述語音激活檢測單元根據所述含噪語音信號，或者所述含噪語音信號和所述原 始噪聲器信號，判斷是否采集到語音信號，并向所述濾波組件輸出判斷結果。 優選的，Ps通過下式計算： 其中，VAD為所述語音激活檢測單元的輸出值，true為采集到語音信號，false為 未米集到語音信號； Eal為所述消噪語音信號經頻域轉換后獲得的第一幅度譜； Ea2為所述參考噪音信號經頻域轉換后獲得的第二幅度譜； β為控制因子，1〈β〈5。 優選的，所述濾波組件提供一第一自適應濾波器，及一第二自適應濾波器； 所述第一自適應濾波器根據所述頻域維納濾波器的解G、所述參考噪聲信號以及 所述含噪語音信號計算所述第一抗串擾濾波單元的濾波系數更新項并通過所述濾波 系數更新項更新所述第一抗串擾濾波單元的濾波系數glit; 所述第二自適應濾波器根據所述頻域維納濾波器的解G、所述消噪語音信號以及 所述噪聲信號，計算所述第二抗串擾濾波單元的濾波系數更新項 §1;，并采用所述濾波系數 更新項g2/更新所述第二抗串擾濾波單元的濾波系數g2it。 優選的，所述第一抗串擾濾波單元通過下式更新濾波系數glit: 所述第二抗串擾濾波單元通過下式更新濾波系數g2it: 優選的，所述第一自適應濾波器通過頻域最小均方算法計算所述第一抗串擾濾波 單元的濾波系數更新項 §1,，具體算式如下： 其中，GfG' λ為大于〇的非線性控制因子； μ為更新步長，且0 < μ < 2 ; ε為防止除零錯誤的微小正數； ?為點乘符號； 上標*為共軛操作； 下標NXΝ表示矩陣為Ν行Ν列； 下標NX (M-N)表示矩陣為N行M-N列； t為所述含噪語音信號s i (t)的Μ點傅立葉變換； E2, t為所述參考噪聲信號e 2 (t)的Μ點傅立葉變換； F-為反傅立葉變換矩陣； E為|Slit|2在t時刻的短期數學期望。 優選的，Slit通過下式計算： Sl t= FT, 其中，上標T為轉置操作；和/或 E2,t通過下式計算： E2j t= F { 1ΧΝ }τ, 其中，上標Τ為轉置操作；和/或 E通過自回歸法進行獲取，具體算式如下： E = YE + (l-y) |Slit|2,0<Y<l〇 優選的，所述第二自適應濾波器通過頻域最小均方算法計算所述第二抗串擾濾波 單元的濾波系數更新項g 2f，具體算式如下： 其中，G2= 1-G 人或G2= (1-G)入； λ為大于〇的非線性控制因子； μ為更新步長，且0 < μ < 2 ; ε為防止除零錯誤的微小正數； ?為點乘符號； 上標*為共軛操作； 下標NXN表示矩陣為N行N列； 下標N本文檔來自技高網...

【技術保護點】
一種語音降噪方法，其特征在于，通過一濾波組件，對一濾波形成的消噪語音信號，以及一濾波形成的參考噪聲信號進行頻域轉換并濾波，以濾波輸出、對應所述參考噪聲信號的原始噪聲信號、以及對應所述消噪語音信號的含噪語音信號，分別控制形成所述參考噪聲信號的第一抗串擾濾波單元，以及形成所述消噪語音信號的第二抗串擾濾波單元進行濾波系數更新。

【技術特征摘要】

【專利技術屬性】
技術研發人員：吳晟，蔣斌，林福輝，
申請(專利權)人：展訊通信上海有限公司，
類型：發明
國別省市：上海;31