實驗二嘯叫檢測抑制系統(tǒng)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

17、,s,s，設計模擬濾波器H(s)。a（4）用雙線性變換法將模擬濾波器Ha(s)變換成數(shù)字濾波器H(z)，即H(z)Ha(s)|s21z11T1z本課題中采用雙線性變換設計二階IIR數(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、與模擬頻率的對應關系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)來實現(xiàn)，用差分方程描述以下：y(n)2x(ni)2i)aiby(n（11）i四、實驗過程及結(jié)果解析i0i1該實驗在matlab下完成，對IIR濾波器設計全過程中的各個步驟，Matlab都供應了相應的工具箱函數(shù)，使IIR數(shù)字濾波器設計變得特別簡單。IIR數(shù)字濾波器的設計步驟和所用的工具箱函數(shù)如圖11所示。圖11IIR數(shù)字濾波器設計步驟流程圖1、運行算法程序，給出以下結(jié)果：嘯叫時語音信號的時域波形圖嘯叫時語音

19、信號的頻譜圖嘯叫控制后的時域波形圖嘯叫控制后的頻譜圖2、對上述結(jié)果進行解析，說明濾波后的結(jié)果。附錄：fs=10000;%設置采樣頻率N=fs*58;%設置信息長度%讀入音頻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;%對嘯叫點進行檢測maxsindex=max(mag(1:N/2);%找出所有頻點中幅值最大的點A=f(index);%利用CZT做細化頻譜解析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);%細化頻率范圍fz=(0:length(z)-1)*(f2-f1)/length(z)+f1;figure(4);plot(fz,abs(z);axis(f1f2010000);xlabel(Hz);ylabel(Magnitude);title(CZT細化頻譜圖);maxsindex=max(abs(z);B=fz(index)%從聲卡播放x數(shù)組，聽到帶有嘯叫的語音sound(x,fs);%設計陷波器%模擬指標阻帶下截止頻率fs1=A-200,阻帶上截止頻率fs1=A+200,阻帶最小衰減as=50db%通帶下截止頻率fp1=A-300,通帶上截止頻率fp2=A+300,阻帶最小衰減ap=rp=;%通帶衰減dBrs