0 引言

    噪聲污染廣泛存在于工業(yè)生產(chǎn)和日常生活中，而其中的周期性噪聲，如發(fā)動機、壓縮機、風扇及螺旋槳等旋轉或往復設備發(fā)出的噪聲也十分常見^[1]。由于其噪聲能量集中分布于特定的頻帶，所以被稱為窄帶噪聲。隨著電子技術的進步，主動噪聲控制技術逐漸發(fā)展，相比于傳統(tǒng)的被動噪聲控制，主動噪聲控制具有成本低、靈活性高等優(yōu)點^[2]。而且窄帶噪聲不同于寬帶噪聲，其頻率數(shù)值往往可以通過非聲學傳感器采集得到，進而避免了聲反饋問題，更適合于采用主動控制的方法^[3]。

    目前主動噪聲控制系統(tǒng)工程實現(xiàn)時主要采用的是DSP處理器，其優(yōu)點在于相關驅(qū)動比較成熟并且芯片成本較低^[4]。但是隨著窄帶噪聲頻率的提高，需要的音頻采樣率也要提高，因此計算量會相應增加，DSP的計算速度就會成為瓶頸，所以目前的主動噪聲控制系統(tǒng)所應用的場景都是低頻噪聲^[5]。為了能夠?qū)Ω哳l噪聲進行主動控制，同時也為了應對多頻率、多通道情況下計算量成倍增加的情況，本系統(tǒng)采用并行計算能力較強的FPGA作為核心處理器，并對并聯(lián)結構的窄帶前饋FxLMS算法進行了工程實現(xiàn)與實際測試。此外考慮到并聯(lián)結構算法對于乘法器資源的消耗大，本文采用乘法器資源復用技術^[6]，極大降低了系統(tǒng)成本，便于系統(tǒng)的實際工程應用。

1 FxLMS算法分析與改進

1.1 窄帶前饋FxLMS算法

    前饋ANC系統(tǒng)有兩種實現(xiàn)方法，其主要區(qū)別在于參考信號的來源^[2]。雖然用聲學傳感器采集初級噪聲作為參考信號的算法可以針對寬帶和窄帶噪聲進行主動控制，但是聲反饋問題會嚴重影響算法的穩(wěn)定性。所以本文設計的系統(tǒng)只針對窄帶噪聲信號，從而可利用加速度計、轉速計等非聲學傳感器采集并計算出噪聲信號主要頻率，進而在系統(tǒng)內(nèi)部生成參考信號。窄帶前饋FxLMS算法的結構如圖1所示。算法主要包括參考信號的產(chǎn)生、次級聲通道的離線辨識以及LMS自適應濾波3個部分^[7]。

    圖1中W(z)代表FIR濾波器，P(z)代表初級聲通道，S(z)代表從次級噪聲數(shù)模轉換到傳聲器采集到音頻數(shù)據(jù)之間的次級聲通道，是次級聲通道S(z)的M階FIR估計。次級聲通道會導致系統(tǒng)不穩(wěn)定^[8]，為了補償次級聲通道的影響，F(xiàn)xLMS算法需要得到次級聲通道S(z)的M階FIR估計 。假設次級聲通道是時不變系統(tǒng),則可以采用離線辨識的方法。圖1中虛線框內(nèi)即為系統(tǒng)所采用的離線辨識模塊。穩(wěn)定后得到次級聲通道的傳遞函數(shù)為：

1.2 并聯(lián)結構FxLMS算法

    在實際中，窄帶噪聲往往還包括其基頻的高次諧波或者其他能量較高的頻率分量。為了解決這個問題，可以使上一節(jié)中的正弦信號發(fā)生器包含多個對應的頻率分量，這樣做的優(yōu)點是算法結構簡單、資源消耗低。但同時會帶來一個嚴重的問題，即當噪聲中各頻率分量差距較大時，同一個步長因子μ無法使所有頻率達到最佳的收斂狀態(tài)，會產(chǎn)生振蕩或者不收斂的情況。而當各頻率分量差距較小時，就需要自適應濾波器的階數(shù)足夠長，以滿足相鄰頻率的高分辨率，這會導致收斂速度慢、量化誤差大等問題^[2]。所以，解決這個問題最好的辦法是單獨處理各個頻率分量，即采用并聯(lián)結構的FxLMS算法。算法結構如圖2所示。

    自適應算法的系數(shù)調(diào)整公式為：

    在并聯(lián)結構的FxLMS算法中，不同頻率的正弦參考信號由其對應的自適應濾波器分別處理。所有濾波器輸出信號相加后作為次級噪聲信號，用于抵消初級噪聲。通過這種方式將噪聲中不同的頻率成分分開處理，系統(tǒng)的穩(wěn)定性和收斂速度都會得到明顯提高，同時增加了系統(tǒng)應用的靈活性。相比于一般結構的FxLMS算法，并聯(lián)結構算法會成倍地增加對系統(tǒng)資源的占用，但是這個不足可以利用系統(tǒng)FPGA實現(xiàn)過程中采用的資源復用技術克服。

2 ANC系統(tǒng)的搭建與實現(xiàn)

    系統(tǒng)采用FPGA作為信號處理器的原因有兩點。首先是將算法進行硬件實現(xiàn)后，其運算速度更快，利于搭建實時處理系統(tǒng)；其次利用FPGA并行計算的優(yōu)點，通過簡單的復制，即可將現(xiàn)有ANC系統(tǒng)擴展為多通道，使系統(tǒng)可以應用于更多的實際場景。然而FPGA在實現(xiàn)數(shù)字信號處理算法，尤其是濾波器等乘加運算較多的算法時，其內(nèi)部的乘法器資源消耗會很大，而選擇乘法器數(shù)量多的FPGA型號會極大提高系統(tǒng)成本。所以為了提高ANC系統(tǒng)的實用性，需要對FPGA內(nèi)部的乘法器進行復用。

2.1 乘法器資源復用技術

    考慮到系統(tǒng)中FPGA的工作時鐘為50 MHz，而音頻模塊的采樣時鐘為32 kHz，滿足系統(tǒng)時鐘頻率遠大于數(shù)據(jù)時鐘頻率的條件。所以采用將同一時刻并行進行的乘法運算部分轉換為串行計算的方法對乘法器進行復用，即所謂的時間復用。時間復用的時序過程如圖3所示。

    假設有N個并行的階數(shù)為M的濾波器，則當數(shù)據(jù)時鐘上升沿到來時，需要計算出N×M個乘法結果，即需要N×M個乘法器。而如果采用復用技術，乘法器在數(shù)據(jù)有效后接受控制器的統(tǒng)一調(diào)度，逐一計算各濾波器在該級的乘法運算，只要滿足兩次乘法運算的時間間隔大于乘法器輸出結果的時間，則N個乘法器串行計算的輸出結果正確。這樣并行濾波器所消耗的乘法器數(shù)目就會降為M個，大大節(jié)約了乘法器資源。

    由于本系統(tǒng)是對相同類型的算法模塊進行乘法器復用，所以可以采用簡單的計數(shù)器實現(xiàn)調(diào)度控制，即在數(shù)據(jù)有效后對系統(tǒng)時鐘開始進行計數(shù)，在不同的計數(shù)值進行各算法模塊的乘法運算，同時在不同的計數(shù)值輸出各乘法運算結果。

2.2 系統(tǒng)分析

    本文所實現(xiàn)的ANC系統(tǒng)所完成的功能是對環(huán)境噪聲中3個主要頻率分量進行主動控制，達到整體降噪的目的。降噪算法采用3個并聯(lián)的FxLMS濾波器實現(xiàn)，如果不采用資源復用技術，除了次級聲通道辨識以外，算法還需要3個128階FIR濾波器以及3個64階LMS自適應濾波器，理論上共需要960個乘法器。由于采用了復用技術，所以系統(tǒng)算法包括一個128階的次級聲通道辨識模塊、一個三輸入三輸出的128階FIR濾波器和一個三輸入三輸出的64階LMS自適應濾波器，只需要320個乘法器即可，相比復用前節(jié)約了2倍的乘法器資源。系統(tǒng)采用的FPGA芯片為Altera公司Stratix III系列中的EP3SE260F1152C2，系統(tǒng)實現(xiàn)后所消耗的資源及比例如表1所示。

    由表1可以看出，F(xiàn)PGA的資源可以滿足算法要求，經(jīng)調(diào)試后系統(tǒng)整體可以正常工作。

3 結果分析

    對所搭建的ANC系統(tǒng)進行實際應用測試，測試所采用的初級噪聲為音箱發(fā)出的窄帶噪聲信號，其主要頻率為1 000 Hz、2 000 Hz和3 000 Hz。ANC系統(tǒng)工作之前，拾音器采集到的聲音為初級噪聲，ANC系統(tǒng)工作之后，拾音器采集到的聲音為初級噪聲和次級噪聲疊加以后的誤差噪聲信號。通過將ANC系統(tǒng)工作前后拾音器所采集到的信號的時域和頻域進行對比，即可以判斷預期的系統(tǒng)功能是否實現(xiàn)。

    對參考信號頻率與噪聲頻率基本相同的情況進行測試，即參考信號的頻率成分1 000 Hz、2 000 Hz和3 000 Hz。圖4顯示了ANC系統(tǒng)工作前后的噪聲信號時域和頻域的對比。

    定義收斂時間為時域波形中的誤差噪聲信號幅值降低到初級噪聲幅值的20%以下的時間，對噪聲時域波形的分析可以得出在無頻偏的情況下，ANC系統(tǒng)的收斂時間為12.5 ms。分析頻域圖可以得出，無頻偏情況下系統(tǒng)對于噪聲中1 000 Hz、2 000 Hz和3 000 Hz的頻率分量分別可以達到15.5 dB、14.3 dB和11.0 dB的衰減?；贒SP的FxLMS算法的系統(tǒng)在汽車發(fā)動機噪聲主動控制測試的結果為13.5 dB^[10]，可以看出本文所實現(xiàn)的系統(tǒng)在降噪效果和功能擴展方面具有明顯優(yōu)勢。

4 結論

    本文提出了一種基于FPGA的窄帶噪聲主動控制系統(tǒng)的設計方法。該系統(tǒng)克服了傳統(tǒng)ANC系統(tǒng)只能針對低頻窄帶噪聲控制應用的不足，使其更適于實際應用。此外，資源復用技術的應用大大減少了系統(tǒng)對于FPGA乘法器資源的消耗，降低了系統(tǒng)成本。由于FPGA的并行性特點，可以通過模塊復制將系統(tǒng)應用于多頻率多通道噪聲主動控制的環(huán)境。

    通過對系統(tǒng)的實際測試驗證了系統(tǒng)噪聲控制的有效性。通過與不同噪聲頻率分析算法配合，可以將該系統(tǒng)應用于多種窄帶噪聲的控制領域。例如對于汽車車內(nèi)發(fā)動機噪聲的主動控制，可以通過實時采集汽車發(fā)動機信號，分析出汽車發(fā)動機噪聲主要頻率，進而利用本文所提出的系統(tǒng)對其進行主動控制。但是本文系統(tǒng)自適應算法的內(nèi)部參數(shù)，包括濾波器階數(shù)以及算法步長因子等是通過仿真實驗確定的，并沒有形成一個統(tǒng)一的選擇方法，還需要后續(xù)的研究加以完善，以擴大系統(tǒng)的實際應用范圍。

參考文獻

[1] ROULET C A.Indoor environment quality in buildings and its impact on outdoor environment[J].Energy and Buildings，2001，33：183-191.

[2] ROURE A.Self-adaptive broadband active sound control system[J].Journal of Sound and Vibration，1985，101(3)：429-441.

[3] KUO S M，MORGAN D R.Active noise control systems：algorithms and DSP implementations[M].New York：Wiley，1996.

[4] 孫琳.基于DSP的窄帶主動噪聲控制系統(tǒng)實現(xiàn)[D].哈爾濱：哈爾濱工業(yè)大學，2011.

[5] 劉劍.基于FXLMS算法的窄帶主動噪聲控制系統(tǒng)性能分析研究[D].哈爾濱:哈爾濱工業(yè)大學，2011.

[6] 劉純武，黃芝平.一種多通道信號處理復用結構及其FPGA實現(xiàn)方法[J].電子測量與儀器學報，2008，22(1)：72-75.

[7] VU H S，CHEN K H，SUN S F，et al.A power-efficient circuit design of feed-forward FxLMS active noise cancellation for in-ear headphones[C].VLSI Design，Automation and Test(VLSI-DAT)，2015 International Symposium on.IEEE，2015：1-4.

[8] XIAO Y，IKUTA A，MA L，et al.Stochastic analysis of the FXLMS-based narrowband active noise control system[J].IEEE Transactions on Audio，Speech，and Language Processing，2008，16(5)：1000-1014.

[9] 劉小明，洪一.基于查找表和Taylor展開的正余弦函數(shù)的實現(xiàn)[J].現(xiàn)代電子技術，2009，32(13)：165-166.

[10] 計煒梁，凌強，朱學俊，等.基于DSP的主動降噪系統(tǒng)設計與實現(xiàn)[J].微型機與應用，2015，34(24)：32-34.

[11] 王雯婷，葛海波，王艷艷.變步長LMS算法抑制窄帶干擾技術研究[J].電子技術應用，2015，41(2)：89-92.

作者信息：

姜宇程，關  添，王佳飛

（清華大學深圳研究生院 生物醫(yī)學工程研究中心，廣東 深圳518055）