現(xiàn)實(shí)的語(yǔ)音通信可能發(fā)生在嘈雜的噪聲環(huán)境中，例如工廠中的手機(jī)通訊會(huì)受到機(jī)器轟鳴聲的影響；火車(chē)駕駛室中的語(yǔ)音通信會(huì)受到電機(jī)運(yùn)行和鐵軌碰撞聲的干擾。噪聲的統(tǒng)計(jì)特性跟現(xiàn)場(chǎng)密切相關(guān)，即使同一場(chǎng)合的噪聲統(tǒng)計(jì)特性也可能會(huì)隨時(shí)間發(fā)生某種變化，這都要求消噪設(shè)備必須具有噪聲跟蹤的自適應(yīng)能力。自適應(yīng)信號(hào)處理的理論經(jīng)過(guò)40多年的發(fā)展和完善，已經(jīng)在許多領(lǐng)域中得到應(yīng)用[1,2]。本文將自適應(yīng)信號(hào)處理技術(shù)應(yīng)用于語(yǔ)音通信的噪聲對(duì)消，通過(guò)仿真研究和樣機(jī)實(shí)驗(yàn)對(duì)軟硬件參數(shù)進(jìn)行了優(yōu)化設(shè)計(jì)，研制了一種適用于手機(jī)麥克風(fēng)、火車(chē)駕駛室中使用的自適應(yīng)對(duì)消系統(tǒng)。

1噪聲對(duì)消原理

自適應(yīng)噪聲對(duì)消系統(tǒng)的原理如圖1所示。他有兩個(gè)輸入：原始輸入和參考輸入，參考輸入為噪聲源ν1(n)，原始輸入為受噪聲污染的信號(hào)x(n)=s(n)+ν0(n)。當(dāng)噪聲成分ν0(n)與信號(hào)s(n)不相關(guān)、與噪聲源ν1(n)相關(guān)時(shí)，自適應(yīng)濾波器AF可以根據(jù)誤差信號(hào)ej來(lái)調(diào)整自身濾波器的系數(shù)，使其輸出yj趨于原始輸入中的ν0(n)，從而使誤差信號(hào)0趨于信號(hào)s(n)。

最小均方(LMS)算法以其算法簡(jiǎn)單、運(yùn)算量小、實(shí)現(xiàn)容易等優(yōu)點(diǎn)在眾多自適應(yīng)信號(hào)處理算法中占有舉足輕重的地位[3,4],本文的系統(tǒng)設(shè)計(jì)采用遞推方法實(shí)現(xiàn)LMS算法。

自適應(yīng)濾波器AF在j時(shí)刻的輸出表示為：

這里的Wj為j時(shí)刻的濾波器系數(shù)，Xj為j時(shí)刻的濾波器輸入。下一時(shí)刻(j+1)的權(quán)系數(shù)要根據(jù)當(dāng)前時(shí)刻j的誤差信號(hào)ej進(jìn)行調(diào)整，調(diào)整的遞推算法為：



其中μ為步長(zhǎng)因子，他影響系統(tǒng)的收斂性，系統(tǒng)收斂的充分條件為：



μ值要針對(duì)語(yǔ)音通信的實(shí)際折衷考慮：μ太小時(shí)，權(quán)系數(shù)收斂慢，不能適應(yīng)噪聲不平穩(wěn)的場(chǎng)合；μ太大時(shí)，對(duì)消效果差，甚至引起系統(tǒng)發(fā)散。

2 系統(tǒng)仿真

2．1聲音樣本的獲取

Matlab 7．0的信號(hào)處理工具箱提供了從計(jì)算機(jī)標(biāo)準(zhǔn)音頻設(shè)備實(shí)時(shí)讀取音頻數(shù)據(jù)的Simulink模塊--From Wave Device，如圖2所示。

圖2中的Signal To workspace將采集的音頻數(shù)據(jù)保存到Matlab的工作空間。作為示例，本文采集了約0．3 s的人聲作為原始信號(hào)，采樣率為16 kHz，其時(shí)域波形如圖3所示。

2．2 步長(zhǎng)因子優(yōu)化

用250 Hz的正弦波模擬現(xiàn)場(chǎng)噪聲，考慮到現(xiàn)場(chǎng)噪聲情況，噪聲幅度的取值應(yīng)與之相當(dāng)，這里取0．01。將聲音樣本和噪聲疊加作為系統(tǒng)的原始輸入，如圖4所示。

取μ=0．5，采用128階的自適應(yīng)濾波器進(jìn)行語(yǔ)音消噪，由圖5可以看出經(jīng)過(guò)0．03 s對(duì)消系統(tǒng)的輸出與原始信號(hào)相似。

自適應(yīng)濾波器的Wj，需要迭代多次才能達(dá)到理想值即E[ej2]E趨勢(shì)于[Eej2)min需要一個(gè)過(guò)程，鑒于步長(zhǎng)因而會(huì)顯著影響系統(tǒng)的對(duì)消效果，下面重點(diǎn)考察了μ取不同時(shí)的學(xué)習(xí)曲線，如圖所示。

由圖6可以看出：
(1)隨著μ值變大，系統(tǒng)的學(xué)習(xí)速度明顯變快；

(2)信號(hào)非平穩(wěn)時(shí)，過(guò)大的μ值容易引起系統(tǒng)失調(diào)，0．1處的均方誤差明顯變大。因而μ的取值應(yīng)當(dāng)折衷考率:一方面當(dāng)噪聲非平穩(wěn)時(shí)，系統(tǒng)的學(xué)習(xí)時(shí)間應(yīng)小于一個(gè)音的平穩(wěn)時(shí)間(約0．1 s)；另一方面，盡量減小系統(tǒng)失調(diào)。細(xì)對(duì)比可以發(fā)現(xiàn)：μ=0．25時(shí)的學(xué)習(xí)時(shí)間已經(jīng)小于一個(gè)節(jié)的平穩(wěn)時(shí)間，這比μ=0．1時(shí)已經(jīng)大大改善，所以將長(zhǎng)因子取為0．25就可以滿足系統(tǒng)需要。

3 系統(tǒng)實(shí)現(xiàn)

3．1 電路設(shè)計(jì)

自適應(yīng)對(duì)消系統(tǒng)的信號(hào)處理器采用TI公司的定點(diǎn)DSP-TMS320VC5509，其內(nèi)核時(shí)鐘可以穩(wěn)定工作在200 MHz，完成兩次乘累加(MAC)運(yùn)算僅需要一個(gè)指令周期(5 ns)。

系統(tǒng)有兩個(gè)麥克風(fēng)通道，其中一個(gè)用于采集原輸入，另一個(gè)用于采集參考輸入。兩通道的電路形式完全一致，信號(hào)調(diào)理與模數(shù)轉(zhuǎn)換電路如圖7所示。其中16 b Codec AD73311通過(guò)同步串行接口與DSP相連。AD73311的參考輸出經(jīng)過(guò)AD8058的Buffer 作為麥克風(fēng)的直流偏置；駐極體麥克風(fēng)信號(hào)采用交流耦合輸入，放大電路的增益為5