實(shí)驗(yàn)二嘯叫檢測抑制系統(tǒng)

1、實(shí)驗(yàn)二嘯叫檢測與控制系統(tǒng)一、實(shí)驗(yàn)?zāi)康?、認(rèn)識擴(kuò)音系統(tǒng)中的嘯叫現(xiàn)象；2、認(rèn)識嘯叫產(chǎn)生的條件；3、解析嘯叫信號的頻譜；4、掌握嘯叫信號檢測的方法；5、掌握嘯叫控制的原理；6、實(shí)現(xiàn)嘯叫控制。二、實(shí)驗(yàn)設(shè)備、音頻功率放大系統(tǒng)；、錄音機(jī)；、計(jì)算機(jī)；、matlab軟件三、實(shí)驗(yàn)內(nèi)容1、認(rèn)識嘯叫及其危害聲反響現(xiàn)象在平常生活中非常常有。在多動能報(bào)告廳，KTV等同時出現(xiàn)揚(yáng)聲器和麥克風(fēng)的場合，由于揚(yáng)聲器和麥克風(fēng)之間存在電聲耦合，必然會以致聲反響現(xiàn)象的產(chǎn)生。聲反響會在反響回路中產(chǎn)生再生混響，使講話、唱歌的聲音嚴(yán)重失真，音質(zhì)碰到破壞，清楚度大大降低。嚴(yán)重時甚至?xí)a(chǎn)生自激嘯叫，限制擴(kuò)聲系統(tǒng)傳聲增益的提升，使整個系統(tǒng)的正常工

2、作碰到影響。還可能會燒毀系統(tǒng)放大器、揚(yáng)聲器中的高音單元，甚至?xí)θ说穆犃υ斐蓚?。因此，在?yīng)用到擴(kuò)聲系統(tǒng)的場合，嘯叫的檢測和控制特別重要。2、嘯叫信號產(chǎn)生的原理聲音信號第一從麥克風(fēng)拾入，經(jīng)過擴(kuò)聲系統(tǒng)的功率放大器放大后由揚(yáng)聲器送出，經(jīng)過各種阻擋物的多次反射后，又被麥克風(fēng)拾入，從而形成一個閉合環(huán)路。若是傳聲器對某些頻點(diǎn)的拾音矯捷度過高，以致聲音在這些頻點(diǎn)的增益是正當(dāng)，就形成了一個正反響過程，聲音信號經(jīng)過多次屢次循環(huán)放大后，在某些頻點(diǎn)的聲音強(qiáng)度高出必然的增益上限，就會發(fā)生自激振蕩，從而產(chǎn)生嘯叫。圖1聲反響原理圖嘯叫的產(chǎn)生必定同時具備以下三個條件：1）揚(yáng)聲器和拾音設(shè)備（麥克風(fēng)）要處于同一聲場中，從而保

3、證揚(yáng)聲器輸出的信號能被拾音設(shè)備再度采集，形成正反響；2）拾音設(shè)備的拾音矯捷度高，系統(tǒng)的傳聲增益大；3）聲場存在弊端共振，即擴(kuò)聲系統(tǒng)的頻譜特色不平坦，在某些頻點(diǎn)上簡單出現(xiàn)共振。3、嘯叫信號的特色為了檢測和控制嘯叫信號，需要對嘯叫信號的基本特色做必然的解析，對嘯叫信號在時域和頻域上進(jìn)行解析，最后得出以下結(jié)論：嘯叫信號的時域波形是一個頻率恒定的正弦波，其幅值隨著時間的增加迅速增大，直到高出了功放放大區(qū)，進(jìn)入飽和區(qū)和截止區(qū)時，產(chǎn)生削波現(xiàn)象，如圖2所示。嘯叫信號的頻譜圖中存在單一的，而且固定不變的嘯叫頻率點(diǎn)，而且嘯叫頻率點(diǎn)對應(yīng)的幅值遠(yuǎn)遠(yuǎn)大于語音信號中其他頻率點(diǎn)的幅值，如圖2所示。圖2嘯叫時時域波形圖3嘯

4、叫時頻域波形該實(shí)驗(yàn)要做的主要工作就是在聲反響建立的初期，在盡可能少的時間內(nèi)檢測出嘯叫對應(yīng)的的頻率點(diǎn)，并設(shè)計(jì)相應(yīng)的陷波器將其控制掉，保證系統(tǒng)工作在正常的線性放大區(qū)，預(yù)防諧波信號的產(chǎn)生。4、嘯叫控制的方法及原理對嘯叫的控制方法大體上分為被動型預(yù)防和主動型控制兩大類。被動型預(yù)防是依照聲場特色，從擴(kuò)聲系統(tǒng)的系統(tǒng)設(shè)計(jì)、設(shè)備選型、聲場布局以及聲場調(diào)整等環(huán)節(jié)著手做好預(yù)防嘯叫的工作；主動型控制力求在控制嘯叫的同時提升擴(kuò)聲系統(tǒng)的傳聲增益。主要有以下幾種方法：基于建筑聲學(xué)特色在音響界，最初解決回聲嘯叫的方法是經(jīng)過降低擴(kuò)聲系統(tǒng)的增益來保證系統(tǒng)的牢固工作，但是，這樣就無法完好滿足擴(kuò)聲音量的需要。因此，工程技術(shù)人員開始

5、在室內(nèi)建筑聲學(xué)、結(jié)構(gòu)聲學(xué)以及室內(nèi)裝飾裝修方面下功夫，來解決聲反響問題和嘯叫現(xiàn)象。但是，這種方法可是被動地控制嘯叫，造價及施工難度太大，不擁有寬泛實(shí)行意義。（2）均衡法圖4均衡器均衡法的原理就是經(jīng)過衰減或切除聲音信號中某些過強(qiáng)的頻率來控制住聲反響。均衡器第一將聲音頻率進(jìn)行分段，爾后經(jīng)過調(diào)治各頻率段增益來填充聲場合揚(yáng)聲器的弊端，達(dá)到控制某一頻點(diǎn)嘯叫的目的。但是，這個方法需專業(yè)音響師手動調(diào)治，要求音響師正確判斷嘯叫頻點(diǎn)并迅速做出反響；其次，在實(shí)質(zhì)操作中，由于均衡器對嘯叫頻點(diǎn)的定位不夠精確，會嚴(yán)重傷害音質(zhì)。移頻法移頻就是搬動頻率，移頻器的工作原理是對聲音信號的頻點(diǎn)向上或向下移幾個或幾十個頻點(diǎn)，破壞構(gòu)成

6、聲反響的條件，最后達(dá)到防范嘯叫的目的。但是，移頻法以致了整個聲音頻率范圍內(nèi)的頻率失真，對音質(zhì)傷害較大；而且，移頻器的工作狀態(tài)與現(xiàn)場有無嘯叫點(diǎn)沒有關(guān)系，這樣就會造成新的失真。因此，應(yīng)該盡量少用這種方法。壓限器它的工作原理是依照輸入聲音信號的大小自動改變輸出信號的增益。當(dāng)聲音信號強(qiáng)度高出設(shè)定的電平值時，壓縮器就將壓縮比調(diào)到：1，語音信號強(qiáng)度就不會再連續(xù)增大，嘯叫也就自但是然地被除掉去了。但是，這種控制方法會對聲音帶來動向損失，不適合要求較高的場合。5）窄帶陷波器利用頻譜解析和判斷條件找出聲音信號的嘯叫頻率點(diǎn)，設(shè)計(jì)參數(shù)可調(diào)的數(shù)字陷波器對聲反響頻點(diǎn)進(jìn)行控制。由于窄帶陷波器在某個很窄的頻帶內(nèi)有很大的衰減

7、，因此對嘯叫頻點(diǎn)定位精度高，陷波器帶寬窄，可以較為正確地控制嘯叫頻點(diǎn)，而且對音質(zhì)的傷害比較低，適合于對音質(zhì)要求較高的場合。嘯叫控制系統(tǒng)基本性能指標(biāo)設(shè)計(jì)的嘯叫控制系統(tǒng)要具備以下條件：第一，陷波器的阻帶帶寬要越窄越好，這樣可以減少聲音的動向失真。第二，對嘯叫頻點(diǎn)進(jìn)行控制時，衰減量的陷波深度越深越好，保證對嘯叫的完好控制；第三，嘯叫頻率出現(xiàn)在任何頻點(diǎn)，系統(tǒng)都可以自動檢測到該頻點(diǎn)，并對其進(jìn)行控制；該實(shí)驗(yàn)的性能指標(biāo)以下：辦理字長：16bit二階IIR數(shù)字陷波器陷波濾波器頻率范圍：100Hz-5KHz陷波濾波器帶寬：1/10倍頻程頻率定位精度：1Hz響應(yīng)時間：1s6、系統(tǒng)算法及其實(shí)現(xiàn)嘯叫的檢測和判斷算法簡

8、單的檢判算法傳統(tǒng)的嘯叫檢測是利用FFT算法實(shí)現(xiàn)的，用FFT算法對每幀信號xn,n0,1,.,N1進(jìn)行頻譜分析，求得對應(yīng)的功率譜Xk,k0,1,.,N/21，爾后計(jì)算出判決值R,計(jì)算公式以下：R10log10(ppeak/pav)（1）在1式中，ppeak是信號功率譜中的最大幅值，所對應(yīng)得頻率點(diǎn)是Fpeak，pav為每幀信號中所有功率幅值的平均值。若是在連續(xù)的一段時間，判決值都比預(yù)設(shè)門限值T高,且最大幅值對應(yīng)的頻點(diǎn)Fpeak都為同一個值，則判斷Fpeak為嘯叫頻點(diǎn)。細(xì)化頻譜解析算法用上面的算法對聲音信號的頻譜進(jìn)行解析時，分辨率為：DfftFs/N48KHz/102446.875Hz（2）不問可知

9、，該算法的頻率分辨率過低，我們可以采用CZT算法來細(xì)化頻譜解析。下面簡單的介紹一下CZT算法:CZT（ChirpZ-transform）算法的實(shí)質(zhì)就是對z變換采用螺線抽樣。z變換的定義式為：N1XzxnZn（3）n0k其中，xn0,1,.,N1是有限長序列。若令抽樣點(diǎn)zAW,k0,1,.,M1，沿z平面nk上的一段螺線作均分角抽樣，就可以獲取抽樣點(diǎn)處的z變換：N1XzkxnAnWnk,k0,1,.M1（4）其中，AA0ej20,WW0ej2n00,A0為初步抽樣點(diǎn)Z0的矢量半徑長度，0為初步抽樣點(diǎn)的相角，W0為螺線的伸展率（當(dāng)W01，為螺線內(nèi)縮；當(dāng)W01，為螺線外伸），0為相鄰抽樣點(diǎn)的角度差（

10、當(dāng)00，為逆時針；當(dāng)00，為順時針），M為要解析的復(fù)頻譜點(diǎn)數(shù)。令0Fl/Fs,0(FhFl)/Fs，就可以獲取Fl:Fh范圍的頻譜，CZT算法的頻率分辨率為(FhFl)/Fs?？梢?，CZT算法可以任意采用FlFh的頻率范圍進(jìn)行解析，可以依照所需的頻率分辨率來確定頻率取樣點(diǎn)數(shù)M。CZT算法的原理表示圖如圖5所示：圖5CZT算法原理圖CZT算法細(xì)化頻譜收效如圖6所示，第一幅圖是直接利用FFT變換后的頻譜圖，第二幅圖是利用CZT變換后的細(xì)化頻譜圖。圖6CZT細(xì)化頻譜收效圖在對嘯叫信號進(jìn)行頻域解析時，若是使用CZT對整個頻域進(jìn)行頻譜解析，那么需要N*M次復(fù)乘，（N-1）M次復(fù)加，運(yùn)算量特別大。依照CZ

11、T算法頻率分辨率的公式Dczt(FhFl)/M，其中Fh:Fl為整個頻帶范圍，要想提升分辨率，只能經(jīng)過增加頻率采樣點(diǎn)數(shù)M來實(shí)現(xiàn)，而嘯叫檢測系統(tǒng)要求的頻率分辨率又很高，需要的M值就會很大，使運(yùn)算量大大增加。頻譜解析方案設(shè)計(jì)綜上所述，F(xiàn)FT算法頻率分辨率太低，達(dá)不到系統(tǒng)設(shè)計(jì)的要求，同時CZT算法運(yùn)算量又太大。綜合考慮FFT算法和CZT算法的優(yōu)弊端，我們采用一個折中的方法：先利用FFT算法對整個頻域進(jìn)行粗解析，一旦檢測到疑似聲反響頻率點(diǎn)Fc存在，就利用Fc設(shè)定CZT算法的Fl和Fh，其中FlFc512，F(xiàn)hFc512，爾后對FlFh頻帶范圍進(jìn)行頻譜細(xì)化解析，如圖7所示。這種算法在保證較高頻率分辨率的

12、同時減小了CZT算法的運(yùn)算量,綜合性能較好。圖7FFT+CZT算法表示圖7、聲反響控制算法實(shí)現(xiàn)方案簡述本實(shí)驗(yàn)采用參數(shù)可調(diào)的IIR陷波器來除掉聲反響頻點(diǎn)的自激嘯叫。陷波器實(shí)質(zhì)是一種帶阻濾波器，它的阻帶很窄，近似于一個特定的頻點(diǎn)，因此又稱為點(diǎn)阻濾波器。當(dāng)聲反響檢判算法檢測到一個聲反響頻點(diǎn)，就設(shè)計(jì)一其中心頻率為該頻點(diǎn)的陷波器對其進(jìn)行控制。圖8陷波器頻率特色圖陷波器設(shè)計(jì)數(shù)字陷波濾波器的設(shè)計(jì)要從濾波器的種類、模擬濾波器的設(shè)計(jì)準(zhǔn)則以及模擬濾波器到數(shù)字濾波器的照射方法三個方面來考慮。陷波器的設(shè)計(jì)可以采用有限沖擊響應(yīng)FIR（FiniteImpulseResponse）濾波器也許無量沖擊響應(yīng)IIR（Infini

13、teImpulseResponse）濾波器。先比較一下兩者的優(yōu)弊端：FIR濾波器有嚴(yán)格線性相位，牢固，適合定點(diǎn)DSP實(shí)現(xiàn)等優(yōu)點(diǎn)，但同時它運(yùn)算量太大且頻率選擇性不高；對照之下，IIR濾波器有優(yōu)異的頻率選擇性，在相同的技術(shù)指標(biāo)下，可用比FTR濾波器少得多的階數(shù)實(shí)現(xiàn)要求，計(jì)算量小，有利于陷波器設(shè)計(jì)的實(shí)時性。綜合考慮兩類濾波器的優(yōu)弊端，我們最后采用采用IIR濾波器設(shè)計(jì)陷波器。IIR濾波器設(shè)計(jì)的基本步驟：第一，將數(shù)字濾波器技術(shù)指標(biāo)，比方通帶截止頻率，阻帶截止頻率等，變換成對應(yīng)的模擬濾波器技術(shù)指標(biāo)；第二，設(shè)計(jì)模擬濾波器Ha(s)使其滿足上述技術(shù)指標(biāo)；第三，將設(shè)計(jì)出的模擬濾波器系統(tǒng)函數(shù)Ha(s)按某種方法變

14、換成數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。在整個設(shè)計(jì)過程中存在著兩種不相同的變換：一是頻帶變換，立刻低通濾波器變換成高通、帶通、帶阻濾波器：二是變換域變換，立刻模擬濾波器變換成數(shù)字濾波器。依照這兩種變換的先后序次，引出以模擬濾波器為基礎(chǔ)的兩種設(shè)計(jì)方法7，如圖9和10所示。圖9IIR數(shù)字濾波器設(shè)計(jì)方法1圖10IIR數(shù)字濾波器設(shè)計(jì)方法2模擬低通濾波器幾種常有的種類有貝塞爾濾波器、巴特沃思濾波器、切比雪夫I型、II型和橢圓濾波器。它們的性能比較如表1所示。表1五類模擬低通濾波器性能比較圖從模擬濾波器到數(shù)字濾波器，有兩種常用的照射方法：脈沖響應(yīng)不變法和雙線性變換法。脈沖響應(yīng)不變法是從濾波器的脈沖響應(yīng)出發(fā)，使數(shù)

15、字濾波器的單位脈沖響應(yīng)h(n)模擬模擬濾波器的沖擊響應(yīng)ha(t)。只有當(dāng)模擬濾波器的頻率響應(yīng)是帶限的，且?guī)抻谡郫B頻率以內(nèi)，才能使數(shù)字濾波器的頻率響應(yīng)在折疊頻率以內(nèi)重現(xiàn)模擬濾波器的頻率響應(yīng)，保證不產(chǎn)生頻率失真。脈沖響應(yīng)不變法最大的弊端是有頻率響應(yīng)的混疊現(xiàn)象，因此只適用于限帶的模擬濾波器，比方衰減特色很好的帶通和帶阻濾波器，而不適用于高通和帶阻濾波器，因此本文采用雙線性變換法。雙線性變換法的基根源理：它是使數(shù)字濾波器的頻率響應(yīng)模擬模擬濾波器頻率響應(yīng)的一種濾波器設(shè)計(jì)方法。雙線性變換法把整個s平面照射到整個z平面，s平面的左半平面照射到z平面的單位圓內(nèi)，s平面的右半平面照射到z平面的單位圓外，s平面

16、的虛軸照射到z平面的單位圓上。這樣就使s平面與z平面是一一對應(yīng)的關(guān)系，它戰(zhàn)勝了多值照射頻譜混疊的弊端，同時也滿足了因果牢固的模擬濾波器可以變換成因果牢固的數(shù)字濾波器的要求。用雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器的設(shè)計(jì)步驟是：1）確定參數(shù)T。雙線性變換法中的參數(shù)T的選擇和最后設(shè)計(jì)出的數(shù)字濾波器沒關(guān)，因此可以取實(shí)際關(guān)系中的值，也可以為設(shè)計(jì)簡單起見，取T=2。由于參數(shù)T的取值會影響計(jì)算精度，因此一般地說，使cT1左右比較適合，不宜取太大或太小的數(shù)值。（2）將數(shù)字濾波器的界線頻率p、s變換為模擬濾波器的界線頻率p、s，變換公式為2tan(p2tan(sp)s)T2，T2。（3）依照模擬濾波器技術(shù)指標(biāo)p,p

17、,s,s，設(shè)計(jì)模擬濾波器H(s)。a（4）用雙線性變換法將模擬濾波器Ha(s)變換成數(shù)字濾波器H(z)，即H(z)Ha(s)|s21z11T1z本課題中采用雙線性變換設(shè)計(jì)二階IIR數(shù)字陷波器的步驟，假設(shè)二階模擬陷波器的傳輸函數(shù)為22s0（5）H(s)s2Bs02其中B為3dB帶寬。用雙線性變換，可以獲取二階數(shù)字濾波器的傳輸函數(shù)為：令：112020H(z)H(s)|1Z1s1Z1B120B120(102)2(102)Z1(102)Z2(1B)2(12)Z(1(6)020102B)Z2則上式可改為：0.5(1)2(1)Z1(1)Z2H(z)(1)Z1Z2（7）1tan()依照雙線性變換法中數(shù)字頻率

18、與模擬頻率的對應(yīng)關(guān)系2有：11則求得IIR陷波器差分方程系數(shù)以下：202cos0（8）000.5*(1)1(1)20.5*(1)（9）b1(1)b2（10）二階IIR數(shù)字陷波器經(jīng)過錯散時間系統(tǒng)來實(shí)現(xiàn)，用差分方程描述以下：y(n)2x(ni)2i)aiby(n（11）i四、實(shí)驗(yàn)過程及結(jié)果解析i0i1該實(shí)驗(yàn)在matlab下完成，對IIR濾波器設(shè)計(jì)全過程中的各個步驟，Matlab都供應(yīng)了相應(yīng)的工具箱函數(shù)，使IIR數(shù)字濾波器設(shè)計(jì)變得特別簡單。IIR數(shù)字濾波器的設(shè)計(jì)步驟和所用的工具箱函數(shù)如圖11所示。圖11IIR數(shù)字濾波器設(shè)計(jì)步驟流程圖1、運(yùn)行算法程序，給出以下結(jié)果：嘯叫時語音信號的時域波形圖嘯叫時語音

19、信號的頻譜圖嘯叫控制后的時域波形圖嘯叫控制后的頻譜圖2、對上述結(jié)果進(jìn)行解析，說明濾波后的結(jié)果。附錄：fs=10000;%設(shè)置采樣頻率N=fs*58;%設(shè)置信息長度%讀入音頻x,fs=wavread(,N);%將音頻輸入到程序中%顯示語音信號時域波形figure(1);subplot(2,1,1);plot(x);xlabel(n);ylabel(幅度/V);title(帶嘯叫的語音信號時域波形);%顯示語音信號x的頻譜；n=0:1:N-1;f=n*fs/N;%頻譜圖中的頻率軸temp=fft(x,N);%對語音信號x作fft變換mag=abs(temp);%求功率譜figure(2);subp

20、lot(2,1,1);plot(f(1:N/2),mag(1:N/2);%畫語音信號x的頻率譜title(帶嘯叫的語音信號頻譜);xlabel(頻率/Hz);ylabel(幅度);axis(0,6000,0,10000);gridon;%對嘯叫點(diǎn)進(jìn)行檢測maxsindex=max(mag(1:N/2);%找出所有頻點(diǎn)中幅值最大的點(diǎn)A=f(index);%利用CZT做細(xì)化頻譜解析f1=A-512;f2=A+512;Fs=1000;m=1024;w=exp(-j*2*pi*(f2-f1)/(m*Fs);a=exp(j*2*pi*f1/Fs);z=czt(x,m,w,a);%細(xì)化頻率范圍fz=(0:length(z)-1)*(f2-f1)/length(z)+f1;figure(4);plot(fz,abs(z);axis(f1f2010000);xlabel(Hz);ylabel(Magnitude);title(CZT細(xì)化頻譜圖);maxsindex=max(abs(z);B=fz(index)%從聲卡播放x數(shù)組，聽到帶有嘯叫的語音sound(x,fs);%設(shè)計(jì)陷波器%模擬指標(biāo)阻帶下截止頻率fs1=A-200,阻帶上截止頻率fs1=A+200,阻帶最小衰減as=50db%通帶下截止頻率fp1=A-300,通帶上截止頻率fp2=A+300,阻帶最小衰減ap=rp=;%通帶衰減dBrs