摘要：傳統(tǒng)有源噪聲控制系統(tǒng)中的音頻編解碼器只能使用高速率的單一采樣率，因此系統(tǒng)中 的建模濾波器和自適應(yīng)濾波器需要的長度較長，而且預(yù)留給每次迭代計算的時間很短，這導(dǎo)致 了計算量大、難以實現(xiàn)實時在線計算的問題。針對以上困難，首先提出一種多采樣率的方法， 即采用抽取和插值的方法為有源噪聲控制系統(tǒng)提供更低的采樣率，從而在減少計算量的同時增 加可用于迭代計算的時間，并基于最小均方算法對次級聲通道進行建模，推演多采樣率雙通道 反饋有源噪聲控制系統(tǒng)模型；其次搭建基于座椅頭靠和數(shù)字信號處理器的實驗平臺，對多采樣 率的雙通道有源噪聲反饋控制系統(tǒng)的有效性進行實驗驗證。驗證結(jié)果表明，針對低頻噪聲，在 相應(yīng)的頻段上本文降噪系統(tǒng)可以達到約 15 dB 的降噪量，降噪效果明顯。

關(guān)鍵詞：有源噪聲控制；多采樣率；雙通道；反饋控制；低頻噪聲

中圖分類號：TB535

文獻標志碼：A

隨著工業(yè)化和城市化的發(fā)展，噪聲污染的問題 日益突出，如何降低生活環(huán)境中的噪聲成為提高生 活質(zhì)量的重大問題之一。傳統(tǒng)的降低噪聲的方法主 要是通過安裝吸音材料等手段阻隔噪聲的傳播，該 手段在很多場合都有廣泛的應(yīng)用，但是實踐證明，這 種方法雖然對高頻噪聲有很好的阻隔效果，但是對 于中低頻噪聲的降噪效果不好[1]。有源噪聲控制 （Active Noise Control， ANC） 方法的出現(xiàn)使得低頻噪 聲得到了較好的控制[2]。

ANC 方法主要是利用電聲器件 （如揚聲器） 發(fā) 出和原始噪聲幅度相同且相位相反的聲波，利用聲 波的相互干涉，從而達到降噪的效果。對于管道這 種一維的聲場來說，利用單通道 ANC 系統(tǒng)便可以獲 得良好的降噪效果。但是對于三維空間這種更復(fù)雜 的聲場環(huán)境，要想獲得良好的降噪效果，就必須拓展 通道的數(shù)量[3] ，從而解決單通道降噪在三維空間中存 在的問題。

上世紀 80 年代以來，國內(nèi)外學(xué)者針對三維空間 噪聲的有源控制進行了大量研究。文獻 [4] 研究了 三維空間聲場中有源靜區(qū)的范圍同噪聲波長的關(guān) 系，結(jié)果表明，有源靜區(qū)的半徑大約為聲波波長的 1/10；Nelson 等[5] 采用多通道前饋控制方法對飛機艙 內(nèi)的窄帶噪聲進行控制，取得了一定的降噪效果；李 卓林等[6] 采用前饋控制系統(tǒng)設(shè)計了一種雙通道室內(nèi) 有源降噪器，并通過實驗驗證了該系統(tǒng)的有效性；張 振超等[7] 提出一種綜合最優(yōu)控制器的有源降噪耳機 的設(shè)計方案，優(yōu)化改善噪聲從不同角度入射情況下 的降噪效果； Rafaely 等[8] 在忽略次級通道間相互耦 合的情況下，設(shè)計了基于 H2 /H∞的有源降噪頭靠系 統(tǒng)，達到了一定的降噪效果；雷成友[9] 采用增加次級 聲源的方法對多通道系統(tǒng)進行解耦合，并通過仿真 實驗驗證了算法的效果。

綜合以上研究可以看出，目前學(xué)者對于空間噪 聲有源控制的探究主要集中在大規(guī)?？臻g范圍的有 源降噪、多通道前饋控制、有源降噪算法等方面，對 于多通道反饋控制、減少計算量等的研究相對較少， 其他的一些算法研究大多受限于算法的復(fù)雜度或者 硬件的計算能力而僅僅局限于仿真研究。另外，傳統(tǒng)的多通道前饋控制系統(tǒng)需要安裝參考傳感器以獲 得參考信號，然而在某些場景里參考信號是不易獲 得甚至是無法獲得的，這就導(dǎo)致采用前饋控制的降 噪系統(tǒng)無法工作。

本文以座椅頭靠為研究對象，設(shè)計并實現(xiàn)一種 只有兩個誤差麥克風(fēng)和兩個次級揚聲器的雙通道反 饋控制的主動降噪系統(tǒng)，基于最小均方 （Least Mean Square， LMS） 算法對次級聲通道進行建模，推演雙通 道反饋 ANC 模型；采用抽取和插值的方法為 ANC 系統(tǒng)提供低的采樣率，并通過插值濾波將控制信號 通過次級揚聲器播放，以減少算法的總計算量；最后 搭建實驗平臺，開展基于數(shù)字信號處理器 （Digital Signal Processor， DSP） 硬件系統(tǒng)進行窄帶噪聲的降噪 實驗，驗證多采樣率控制系統(tǒng)的有效性與實用性。

1""" 雙通道反饋控制原理

1.1 雙通道反饋 ANC 基本結(jié)構(gòu)

雙通道反饋 ANC 系統(tǒng)由兩個誤差麥克風(fēng) M1 和 M2 以及兩個次級揚聲器 L1 和 L2 組成，ANC 系統(tǒng) 經(jīng)過自適應(yīng)濾波器計算產(chǎn)生兩個 n 時刻的“反噪聲” 信號 y1 （n） 和 y2 （n），并由這兩個信號驅(qū)動次級揚聲器 發(fā) 出 反 相 聲 波 與 麥 克 風(fēng) 處 的 原 始 噪 聲 d1 （n） 和 d2 （n） 疊加，麥克風(fēng) M1 和 M2 得到誤差信號 e1 （n） 和 e2 （n），將誤差信號傳送到 DSP 中計算并更新自適應(yīng) 濾波器的參數(shù)。如此迭代循環(huán)，當(dāng)殘余噪聲達到最 小時降噪系統(tǒng)進入穩(wěn)定狀態(tài)。雙通道反饋 ANC 基本結(jié)構(gòu)如圖 1 所示，其中，S11、S12、S21、S22 為 4 個 次級通道的模型。

1.2 雙通道反饋 ANC 算法

雙通道 ANC 算法是在傳統(tǒng)單通道 ANC 算法的 基礎(chǔ)上拓展而來。但是需要注意的是，雙通道并不 是在單通道的基礎(chǔ)上簡單地增加一個通道，而且還 增加了兩個耦合通道 S12 和 S21，因此在算法中也要加 入這兩個耦合通道的影響[10]。

雙通道反饋 ANC 算法的結(jié)構(gòu)如圖 2 所示。圖 中 W1 和 W2 為主控制器，x1 和 x2 為兩個麥克風(fēng)處的 參考信號，z 表示系統(tǒng)中的 z 變換。其中 W1 和 W2 采用 自適應(yīng)濾波器實現(xiàn)，為了達到自適應(yīng)濾波的目的，基于雙 通道 FxLMS（Filtered-x LMS） 算法更新其參數(shù)。FxLMS 算法是基于最小均方差算法并考慮了次級通道影響 的自適應(yīng)濾波算法，通過調(diào)節(jié)濾波器的權(quán)系數(shù)使得 誤差信號趨向于 0，從而實現(xiàn)降噪的目的[11-13]。

3""" 實驗平臺的搭建

3.1 硬件系統(tǒng)

雙通道有源噪聲反饋控制系統(tǒng)的硬件結(jié)構(gòu)如 圖 5 所示，主要由兩個麥克風(fēng)、話筒放大器、揚聲 器、功率放大器以及 DSP 和音頻芯片組成 ，其中 DSP 采用德州儀器公司的 TMS320C6748 高性能浮 點 型 DSP，音頻芯片采 用 AIC3106， DSP 通 過 I2C （Inter-Integrated Circuit）協(xié)議配置音頻芯片的各寄存 器，從而使音頻芯片的模數(shù)轉(zhuǎn)換器 （ADC）、數(shù)模轉(zhuǎn)換 器 （DAC）、左右聲道等功能。麥克風(fēng)采集到的音頻 模擬信號通過音頻芯片的 ADC 采樣為數(shù)字信號，并 通過 I2S（Inter-IC Sound） 協(xié)議傳輸?shù)?DSP 的 MCASP （Multichannel Audio Serial Port） 外設(shè)，從而使信號在 DSP 中進行 ANC 算法的計算。同樣，計算好的控制 信號經(jīng)過 MCASP 傳輸?shù)揭纛l芯片，并通過 DAC 轉(zhuǎn) 化為模擬信號驅(qū)動揚聲器發(fā)聲。

雙通道有源噪聲反饋控制系統(tǒng)實驗裝置如圖 6 所示，兩個誤差麥克風(fēng)和次級揚聲器分列于座椅頭 靠的兩側(cè)，對應(yīng)于人的雙耳。上位機通過仿真器將 代碼燒寫到 DSP 中，揚聲器和麥克風(fēng)分別通過 Line out 和 Line in 孔連接到 DSP 開發(fā)板。

3.2 軟件系統(tǒng)

軟件設(shè)計大致分為兩個部分：次級通道的辨識 和 ANC 算法的實現(xiàn)。設(shè)置 DSP 主頻為 400 MHz 和 音頻編解碼器的采樣率為 8000 Hz，并通過抽取和插 值的方式將 ANC 系統(tǒng)的采樣率設(shè)置為 2000 Hz。

程序采用直接存儲器訪問 （DMA） 觸發(fā)中斷的方式， 在中斷服務(wù)程序中進行自適應(yīng) FIR 濾波器沖激響應(yīng) 的迭代更新。圖 7為軟件程序的流程框圖，當(dāng)系統(tǒng)接 通電源后，先對 C6748DSP 和 AIC3106 音頻編解碼器 的有關(guān)功能外設(shè)以及相關(guān)變量進行初始化；然后對 4 個次級通道進行辨識，獲得次級通道的模型參數(shù) 后，啟動雙通道 ANC 程序調(diào)用次級通道的模型即可 實現(xiàn)降噪。

4""" 實驗驗證與分析

4.1 次級通道離線建模結(jié)果

使用 128 階的 FIR 濾波器對次級通道進行建 模。DSP 生成白噪聲，通過揚聲器進行播放，同時左 右聲道的兩個麥克風(fēng)采集參考信號，算法收斂后即 可完成對一條主通道和一條耦合通道的建模。圖 8 和圖 9 分別示出了 S11 和 S12 兩個次級通道模型的沖 激響應(yīng)，S21 和 S22 與此類似。由圖 8 和圖 9 可知，每 個次級通道模型的沖激響應(yīng)均約在時刻為 70 時衰減到 0 左右，因此次級通道模型的階數(shù)完全可以保證建 模的準確度。直接次級通道 S11 模型的沖激響應(yīng)幅 值大約是耦合次級通道 S12 模型沖激響應(yīng)幅值的 2 倍多，與實驗平臺的實際情況吻合。

得到次級通道模型的參數(shù)后，利用 MATLAB 的 Filter Designer 工具分析次級通道 S11 和 S12 的頻率響 應(yīng)，結(jié)果分別如圖 10 和圖 11 所示。由圖可知，次級 通道在 0～100 Hz和 800～1000 Hz 時幅值響應(yīng)很低，這是由于實驗的電聲器件無法對 100 Hz以下的信號 做出響應(yīng)，以及抗混疊濾波器的截止頻率略低于 1000 Hz。但是次級通道可以對 100～800 Hz 頻段的 聲信號做出很好的響應(yīng)。

4.2 ANC 實驗結(jié)果分析

通過對次級通道頻率響應(yīng)的分析可知，本系統(tǒng) 可以對 100～800 Hz 的噪聲做出較好的響應(yīng)，因此初 級噪聲的頻率范圍不應(yīng)該超出該范圍。分別對 300 Hz 和 400 Hz 的窄帶噪聲進行降噪，以驗證該系統(tǒng)的有 效性。實驗方法為：通過 MATLAB 分別編寫 300 Hz和 400 Hz的窄帶噪聲作為初級噪聲的發(fā)聲程序，并 通過電腦連接音響進行播放。設(shè)置雙通道 FxLMS 算法步長 μ=0.01，啟動 ANC 系統(tǒng)進行在線計算和降 噪，通過誤差麥克風(fēng)采集降噪開啟前后的誤差信號。 圖 12 和圖 13 分別示出了開啟降噪后誤差麥克 風(fēng)采集到的信號。由圖可知，當(dāng)開啟降噪后誤差麥 克風(fēng)接收到的聲信號幅值迅速減小，算法很快就收 斂，殘余噪聲的幅值進入平穩(wěn)狀態(tài)。300 Hz 和 400 Hz 下平穩(wěn)后的殘余噪聲的幅值均大約是原噪聲幅值 的 50%。

功率譜密度 （Power Spectral Density， PSD） 可以 很好地反映信號在一定頻率范圍內(nèi)的平均功率分 布。圖 14 和圖 15 分別示出了 300 Hz 和 400 Hz 噪 聲環(huán)境下原始噪聲和計算收斂后殘余噪聲信號的 PSD 對比圖，為了顯示更加直觀，截取了噪聲分布的 主要頻段，可以看到相應(yīng)頻段的噪聲被明顯削弱， 300 Hz 頻段上實現(xiàn)了約 18 dB 的降噪量，400 Hz 頻 段上實現(xiàn)了約 15 dB 的降噪量。

實驗數(shù)據(jù)在時域和頻域的分析均表明，采用多 采樣率的雙通道反饋式有源降噪系統(tǒng)可以有效消減 其頻響范圍內(nèi)的窄帶噪聲。

4.3 DSP 運算負荷分析

2.3 節(jié)分別對多采樣率和單采樣的雙通道有源 噪聲反饋控制系統(tǒng)在每次迭代時的計算量 （主要是 浮點型加法和浮點型乘法） 進行了對比分析和統(tǒng) 計。對本系統(tǒng)而言，fs=8000 Hz，LSD=LWD=128，所需 的乘法次數(shù)為 1860，加法次數(shù)為 1846，如若采用單 采樣率系統(tǒng)，則需 7172 次乘法運算和 7162 次加法 運算。不同架構(gòu)的 DSP 處理相同計算量時的計算負 荷會有所不同。對于本文所使用的 TI 公司 C674x 系 列 DSP 來說，其每個時鐘周期可進行 4 次單精度浮 點加法和 2 次單精度浮點乘法。對于 C674x 系列 DSP 來說，乘法所需的時間是加法所需時間的 2 倍， 故影響計算時間最主要的因素是乘法。表 2 對比總 結(jié)了多采樣率和單采樣率的雙通道有源噪聲反饋控制系統(tǒng)在 C6748 DSP 上進行實現(xiàn)時每次迭代計算所 需的時鐘周期數(shù)，可以看出，采用 4 倍抽取的多采樣 系統(tǒng)的周期數(shù)大約是單采樣系統(tǒng)周期數(shù)的 1/4。

5""" 結(jié)束語

本文針對單采樣率雙通道反饋式 ANC 系統(tǒng)數(shù) 據(jù)量多、計算復(fù)雜度高的問題，提出一種利用抽取與 插值實現(xiàn)多采樣率的反饋式雙通道有源降噪系統(tǒng)， 并推導(dǎo)了相應(yīng)的算法。以空間窄帶噪聲為實驗對 象，搭建了以座椅和 DSP 為平臺的實驗系統(tǒng)，進行降 噪實驗，對多采樣率雙通道有源噪聲反饋控制系統(tǒng) 的有效性進行驗證，得到了以下結(jié)論：

（1） 本系統(tǒng)通過 4 倍的抽取和插值處理使得主動 降噪系統(tǒng)的采樣率降低為音頻編解碼器采樣率的 1/4，為 DSP 減少了約 3/4 的計算量，極大地減輕了系 統(tǒng)的計算負擔(dān)。

（2） 采用多采樣率的雙通道反饋式 ANC 系統(tǒng)在 300 Hz 和 400 Hz 的窄帶噪聲環(huán)境下分別取得了 18 dB 和 15 dB 的降噪效果。

為了增加本系統(tǒng)的實用性，未來可以考慮對本 系統(tǒng)進行改進使其適用于更加復(fù)雜噪聲 （如混頻、寬 帶噪聲） 的降噪；也可以將音頻娛樂系統(tǒng)嵌入到本系 統(tǒng)中，以實現(xiàn)在消除環(huán)境噪聲的同時能夠播放高品 質(zhì)的音頻以供學(xué)習(xí)和娛樂。

參考文獻：

ELLIOTT S J， NELSON P A. Active noise control[J]. IEEE Signal Processing Magazine， 1993， 10（4）： 12-35.

KUO" S" M，" MORGAN" D." Active" Noise" Control" Systems： Algorithms" and" DSP" Implementations" [M]." Manhattan： John Wiley amp; Sons Inc， 1996.

陳克安. 有源噪聲控制 [M]. 北京： 國防工業(yè)出版社， 2014.

GARCIA-BONITO" J，" ELLIOTT" S" J，" BONILHA" M. Active" cancellation" of" pressure" at" a" point" in" a" pure" tone diffracted" diffuse" sound" field[J]. Journal" of" Sound" and Vibration， 1997， 201（1）： 43-65.

ELLIOTT" S" J，" NELSON" P" A." The" active" control" of sound[J]." Electronics" amp;" Communications" Engineering Journal， 1990， 46（4）： 127-136.

李卓林，馬令坤，馮智睿. 雙通道室內(nèi)主動降噪器設(shè)計[J]. 科學(xué)技術(shù)與工程， 2019， 19（31）： 108-115.

張振超， 安峰巖， 吳鳴， 等. 一種數(shù)字有源降噪耳機控制器 設(shè)計方法[J]. 應(yīng)用聲學(xué)， 2017， 36（4）： 317-323.

RAFAELY B， ELLIOTT S J. H-2/H-infinity active control of sound in a headrest： Design and implementation[J]. IEEE Transactions" on" Control" System" Technology，" 1999，" 7（1）： 79-84.

雷成友. 局部有源噪聲控制的理論和實驗研究 [D]. 北京： 中國科學(xué)院大學(xué)， 2015.

TOKHI O， VERES S. Active Sound and Vibration Control Theory" and" Applications" [M]." London：" The" Institution" of Electrical Engineers，2002.

BURGESS J C. Active adaptive sound control in a duct： A computer simulation[J]. The Journal of the Acoustical Soci[1]ety of America， 1998， 70（2）： 715-726.

MORGAN" D" R." History，" applications，" and" subsequent development" of" the" FXLMS" algorithm[J]. IEEE" Signal Processing Magazine， 2013， 30（3）： 172-176.

Morgan D R. An analysis of multiple correlation cancella[1]tion loops with a filter in the auxiliary path[J]. IEEE Trans[1]action on Acoustics ， Speech and Signal Processing， 1980， 28（4）： 454-467.

ELLIOTT S J， SUTTON T J. Performance of feedforward and feedback systems for active control[J]. IEEE Transac[1]tions" on" Speech" and" Audio" Processing，" 1996，" 4（3）：" 214- 223.

ERIKSSON L J， ALLIE M C. Use of random noise for on[1]line" transducer" modeling" in" an" adaptive" active" attenuation system[J]. The" Journal" of" the" Acoustical" Society" of" Ame[1]rica， 1989， 85（2）： 797-802.

KUO" S" M，" BOB" H" L，" WENSHUN" T." Real-Time" Digital Signal Processing-Implementations， Applications [M]. UK： Wiley： West Sussex， 2013.

LEE K G W， KUO S M. Subband Adaptive Filtering： The[1]ory" and" Implementation" [M]." UK：" Wiley：" West" Sussex， 2009.

WANG L， GAN W， KHONG A， et al. Convergence analy[1]sis" of" narrowband" feedback" active" noise" control" system with imperfect secondary path estimation[J]. IEEE Transac[1]tions" on" Audio，" Speech，" and" Language" Processing，" 2013， 21： 2403-2411.