基于M-估計的沖擊類噪聲有源控制方法
技術(shù)領(lǐng)域：
：本發(fā)明涉及一種基于M-估計的沖擊類噪聲有源控制方法，具體涉及一種針對沖擊類噪聲的自適應(yīng)主動控制算法，尤其是對沖擊類噪聲的前饋自適應(yīng)控制系統(tǒng)。
背景技術(shù)：
：：主動噪聲控制技術(shù)(ANC，ActiveNoiseControl)逐漸成為抑制低頻噪聲問題的主要手段，尤其是數(shù)字信號處理器(DSP)的發(fā)展極大地促進了ANC技術(shù)的商業(yè)化。ANC的基本原理是由控制系統(tǒng)放出與噪聲波形幅度相等、相位相反的次級聲波進行相互干涉抵消。經(jīng)典的ANC主動控制系統(tǒng)一般用自適應(yīng)濾波器來作為控制器，通過濾波最小均方算法(FXLMS，F(xiàn)iltered-xLeastMeanSquare)來不斷調(diào)節(jié)濾波器參數(shù)。然而，F(xiàn)XLMS算法的提出是基于噪聲信號的二階統(tǒng)計矩(均方誤差，E[e2(n)])，一般假設(shè)信號是確定性的或者統(tǒng)計上服從高斯分布。因此，F(xiàn)XLMS算法對于很多帶有沖擊性的噪聲收斂性能將會很差，控制系統(tǒng)極易失穩(wěn)。為了解決這一難題，很多學(xué)者提出了相應(yīng)的改進算法：一類是采用穩(wěn)健的誤差統(tǒng)計矩。1995年，Leahy等人提出了最小均方p范數(shù)(FXLMP，F(xiàn)iltered-xLeastMeanp-norm)算法(R.Leahy,Z.Zhou,andY.-C.Hsu,"Adaptivefilteringofstableprocessesforactiveattenuationofimpulsivenoise,"inAcoustics,Speech,andSignalProcessing,1995.ICASSP-95.,1995InternationalConferenceon,1995,pp.2983-2986vol.5.)，該算法是基于誤差信號的分數(shù)階矩E[|e(n)|p](而不是二階矩)，對于沖擊類噪聲具有很好的穩(wěn)定性。但是，不足之處是需要對噪聲信號進行預(yù)先的統(tǒng)計估計(非高斯α穩(wěn)態(tài)分布)以較精確地確定p值以及由范數(shù)計算引起的DSP大量計算量。因此，2011年南京大學(xué)聲學(xué)所的吳等人用誤差信號的對數(shù)變換的均方值作為罰函數(shù)(L.Wu,H.He,andX.Qiu,"AnActiveImpulsiveNoiseControlAlgorithmWithLogarithmicTransformation,"Audio,Speech,andLanguageProcessing,IEEETransactionson,vol.19,pp.1041-1044,2011.)，提出了基于對數(shù)變化的脈沖類噪聲自適應(yīng)控制算法，并且在2014年獲得專利授權(quán)，公告號為CN101789771B。該算法無需要預(yù)先估計脈沖信號的統(tǒng)計特征，便于實際應(yīng)用。但是，對數(shù)變換的誤差值計算也會帶來巨大的計算成本，不過吳等人提出可以預(yù)先存儲對數(shù)變化的值通過查表的方式避免在實時計算中的運算量。2007年,Thanigai等人提出了最小均方M-估計算法(FXLMM，thefiltered-xleastmeanM-estimate)，該算法使用誤差信號的穩(wěn)健M-估計罰函數(shù)E[ρ{e(n)}]來代替?zhèn)鹘y(tǒng)的最小均方(P.Thanigai,S.M.Kuo,andR.Yenduri,"Nonlinearactivenoisecontrolforinfantincubatorsinneo-natalintensivecareunits,"2007IeeeInternationalConferenceonAcoustics,Speech,andSignalProcessing,VolI,Pts1-3,Proceedings,pp.109-112,2007.)，典型的HampleM-估計通過非線性分段函數(shù)來去除沖擊信號對收斂過程的擾動。但是，當(dāng)參考信號通道的脈沖信號與誤差信號通道具有延時情況出現(xiàn)時，算法的魯棒性會大大下降。另一類改進算法是通過引入沖擊信號的閾值對經(jīng)典的FXLMS算法進行簡單的修改。2006年，Sun等人提出了改進的FXLMS算法用于沖擊類噪聲控制(X.Sun,S.M.Kuo,andG.Meng,"Adaptivealgorithmforactivecontrolofimpulsivenoise,"JournalofSoundandVibration,vol.291,pp.516-522,2006.)，主要是在誤差信號通道上引入了統(tǒng)計閾值來約束過大樣本值對算法收斂性的干擾。同樣，Akhtar等人擴展了Sun的方法(M.T.AkhtarandW.Mitsuhashi,"ImprovingperformanceofFxLMSalgorithmforactivenoisecontrolofimpulsivenoise,"JournalofSoundandVibration,vol.327,pp.647-656,2009.)，在參考信號和誤差信號通道都引入了統(tǒng)計閥值，該閥值的功能不僅是舍棄大幅值的信號量，而是用統(tǒng)計信號的上下限值來替代沖擊數(shù)據(jù)，因此可以更大程度地提高算法的穩(wěn)定性和穩(wěn)態(tài)性能。技術(shù)實現(xiàn)要素：目的：為了克服現(xiàn)有技術(shù)中存在的不足，本發(fā)明提供一種基于M-估計的沖擊類噪聲有源控制方法。技術(shù)方案：為解決上述技術(shù)問題，本發(fā)明采用的技術(shù)方案為：一種基于M-估計的沖擊類噪聲有源控制方法，該方法是以誤差信號的魯棒M-估計的均方值為代價函數(shù)，并在參考信號通道引入閾值來限制沖擊樣本對算法迭代過程的影響，具體包括以下步驟：1)定義代價函數(shù)J(n)＝E[ρ{e(n)}]，ρ{e(n)}為具有高魯棒性的M-估計函數(shù)，采用通用的Hampel三段梯度函數(shù)；2)估計M-估計函數(shù)的閥值參數(shù)；3)設(shè)定參考信號通道的閾值，確定參考信號的截斷處理結(jié)果xc(n)；4)更新濾波器系數(shù)：其中，n為采樣時間序列，u為收斂系數(shù)，q{e(n)}為誤差信號的M-估計，Xc(n)為經(jīng)過閾值處理的參考信號向量，為次通道傳遞函數(shù)的估計模型。所述M-估計函數(shù)的閥值參數(shù)通過在線估計實現(xiàn)，算法的相對計算量來自在線均值的估計上，通過平行和順序冒泡法實現(xiàn)。下面將簡單列出算法的推導(dǎo)過程：定義自適應(yīng)濾波器的代價函數(shù)為：J(n)＝E[ρ{e(n)}]，這里,ρ{e(n)}為具有高魯棒性的M-估計函數(shù)。比較通用的Hampel三段梯度M-估計函數(shù)如圖1所示，具體的表達式：這里ξ,Δ1和Δ2是閾值參數(shù)，可以通過在線或者離線的“無脈沖”干擾的信號的方差估計來確定這些閾值。接著，定義評價函數(shù)和權(quán)重函數(shù)所以，代價函數(shù)的一階梯度可表示為：接著使用最陡梯度算法可以獲得自適應(yīng)濾波系數(shù)W(n)的迭代公式：其中，u是算法的收斂系數(shù)，次通道傳遞函數(shù)的模型。式(4)可以看作是經(jīng)典FXLMS算法的通用形式。當(dāng)誤差信號的絕對值|e(n)|小于閾值ξ時，q{e(n)}等于1使得FXLMM算法與FXLMS算法等價。而當(dāng)誤差信號大于ξ時,q{e(n)}開始遞減直到當(dāng)|e(n)|≥Δ2時趨于零值。但是，傳統(tǒng)的FXLMM算法對參考信號的大樣本值不具有魯棒性，尤其是當(dāng)沖擊類噪聲的參考路徑與誤差路徑具有沖擊延時時。因此，本發(fā)明提出的改進算法為：其中，參考信號的修正值為：有益效果：本發(fā)明提供的基于M-估計的沖擊類噪聲有源控制方法，通過改進傳統(tǒng)的M-估計FXLMM算法，在參考信號通道引入沖擊信號的M-估計該去除沖擊數(shù)據(jù)，更大程度地提高算法對沖擊類噪聲的控制魯棒性。該發(fā)明提出的算法系統(tǒng)地包括了基于統(tǒng)計閥值的其他類算法，給出了更加科學(xué)合理的統(tǒng)計性描述?；隰敯艚y(tǒng)計學(xué)的M-估計函數(shù)，從根本上解決傳統(tǒng)的二階統(tǒng)計矩(基于高斯分布)對沖擊類噪聲的失穩(wěn)、發(fā)散局限，具有以下優(yōu)點：魯棒統(tǒng)計特征量M-估計具有對沖擊類噪聲的穩(wěn)定性，噪聲控制性能顯著提高；另一方面系統(tǒng)化描述閾值類算法，指導(dǎo)實際應(yīng)用。附圖說明圖1是HampelM-估計函數(shù)ρ{e(n)}以及評價函數(shù)ψ{e(n)}。圖2是基于改進型MFXLMM算法的沖擊類噪聲單通道前饋控制示意圖。圖3是窗的寬度以及濾波器階數(shù)對MFXLMM算法的計算量的影響。圖4是仿真實驗中采取的主、次通道傳遞函數(shù)的頻率響應(yīng)圖。圖5是單次數(shù)值仿真時域控制結(jié)果。圖6算法平均收斂特性的比較結(jié)果。具體實施方式下面結(jié)合附圖對本發(fā)明作更進一步的說明。如圖1至圖6所示，一種基于M-估計的沖擊類噪聲有源控制方法，該方法是以誤差信號的魯棒M-估計的均方值為代價函數(shù)，并在參考信號通道引入閾值來限制沖擊樣本對算法迭代過程的影響，具體包括以下步驟：1)定義代價函數(shù)J(n)＝E[ρ{e(n)}]，ρ{e(n)}為具有高魯棒性的M-估計函數(shù)，采用通用的Hampel三段梯度函數(shù)；2)估計M-估計函數(shù)的閥值參數(shù)；3)設(shè)定參考信號通道的閾值，確定參考信號的截斷處理結(jié)果xc(n)；4)更新濾波器系數(shù)：其中，n為采樣時間序列，u為收斂系數(shù)，q{e(n)}為誤差信號的M-估計，Xc(n)為經(jīng)過閾值處理的參考信號向量，為次通道傳遞函數(shù)的估計模型。所述M-估計函數(shù)的閥值參數(shù)通過在線估計實現(xiàn)，算法的相對計算量來自在線均值的估計上，通過平行和順序冒泡法實現(xiàn)。下面通過數(shù)值仿真實驗對本發(fā)明提出的算法進行詳細的說明。首先，分析算法的計算量，并與經(jīng)典的FXLMS算法進行比較。接著，設(shè)計單通道的ANC仿真系統(tǒng)驗證算法的有效性。算法計算量分析：對于自適應(yīng)控制算法的計算量，可以通過算法所用的加法和乘法量來描述。一般來說，F(xiàn)XLMS算法主要的計算量是由控制器FIR濾波器長度N以及次通道函數(shù)的長度M決定，總的計算量為：2(2N+M)。而本發(fā)明提出的MFXLMM算法，主要的運算量將是M-估計函數(shù)里包含的閾值的在線識別。閥值的在線識別可以通過Zou提出的方法(Y.X.Zou,S.C.Chan,andT.S.Ng,"ArobustM-estimateadaptivefilterforimpulsenoisesuppression,"Icassp'99:1999IeeeInternationalConferenceonAcoustics,Speech,andSignalProcessing,ProceedingsVolsI-Vi,pp.1765-1768,1999.)，主要是通過估計“非脈沖”信號的方差，估算公式里需要確定一段長度信號的中間值?？梢酝ㄟ^冒泡法進行估計，MFXLMM算法的總體計算量為2(2N+M)+2NwlogNw)或者2(2N+M)+2Nw)。結(jié)果如圖2所示，不同濾波器長度以及窗的寬度對算法的計算量有很大影響(歸一化的計算量，相比于FXLMS算法)。仿真系統(tǒng)構(gòu)成：圖2所示為基于MFXLMM算法的沖擊類噪聲單通道前饋控制示意圖。參考信號x(n)可以通過麥克風(fēng)傳聲器獲取，誤差麥克風(fēng)用來測量誤差信號e(n)。W(z)是自適應(yīng)濾波器控制器，一般采用有限脈沖(FIR)結(jié)構(gòu)?？刂破鞯妮敵鲂盘杫(n)是用來驅(qū)動次級聲源產(chǎn)生次級聲波來抵消主要噪聲d(n)。xc(n)是經(jīng)過閾值的參考信號，而x′(n)是經(jīng)過次通道函數(shù)濾波后得到，用來彌補實際控制路徑的相位以及幅值影響。主、次級通道傳遞P(z)和S(z)函數(shù)如圖4所示，本實驗用FIR濾波器來估計主、次通道，階數(shù)分別是256和128。并且假設(shè)次通道估計模型沖擊類參考信號通過在隨機白噪聲中加入固定數(shù)目的脈沖來生存。通過對三個算法FXLMS、FXLMM和MFXLMM算法進行比較分析，不斷地增加系統(tǒng)的收斂系數(shù)，來找到系統(tǒng)的最優(yōu)性能。仿真結(jié)果：圖5是單次仿真實驗的時域比較結(jié)果，圖5(a)是控制前沖擊噪聲信號，圖5(b)是經(jīng)典FXLMS算法的控制結(jié)果，圖5(c)是傳統(tǒng)的FXLMM算法，圖5(d)是本發(fā)明提出的MFXLMM算法。由圖結(jié)果可以看出，所提出的算法能有效地控制沖擊類噪聲，在四個突發(fā)脈沖位置，系統(tǒng)能有很好的魯棒性。然而經(jīng)典的FXLMS和FXLMM算法在第四個脈沖位置出現(xiàn)明顯的發(fā)散。較為系統(tǒng)的結(jié)果如圖6所示，進行了25次仿真實驗求出平均降噪(MeanNoiseReduction)效果。可以看出本發(fā)明提出的算法具有很好的穩(wěn)定性，在沖擊噪聲處能夠很好地保障魯棒性，并且能夠較好地控制。以上所述僅是本發(fā)明的優(yōu)選實施方式，應(yīng)當(dāng)指出：對于本
技術(shù)領(lǐng)域：
：的普通技術(shù)人員來說，在不脫離本發(fā)明原理的前提下，還可以做出若干改進和潤飾，這些改進和潤飾也應(yīng)視為本發(fā)明的保護范圍。當(dāng)前第1頁1 2 3 當(dāng)前第1頁1 2 3