姜曉云，李文靖

（山西大同大學(xué)物理與電子科學(xué)學(xué)院，山西大同037009）

數(shù)字混響器及其改進(jìn)

姜曉云，李文靖

（山西大同大學(xué)物理與電子科學(xué)學(xué)院，山西大同037009）

混響時間是衡量音質(zhì)的重要參量。為了獲得最佳混響時間，采用數(shù)字混響器對聲源信號進(jìn)行處理是目前經(jīng)常使用的方法之一。本文對兩種混響器模型進(jìn)行了詳細(xì)分析，并對其進(jìn)行了改進(jìn)。仿真驗(yàn)證了理論的正確性。

混響時間；數(shù)字混響器；濾波器

混響是室內(nèi)常見的一種聲音現(xiàn)象。它是室內(nèi)聲源發(fā)出的聲音在傳播過程中被室內(nèi)各種界面不斷反射和積累的結(jié)果。室內(nèi)的聲源連續(xù)發(fā)聲，當(dāng)達(dá)到平衡時停止發(fā)聲，此時室內(nèi)聲源發(fā)射出的聲能與被吸收的聲能相等，這時室內(nèi)還會有“殘留”的聲音，這種現(xiàn)象就是混響。

混響時間是指室內(nèi)聲源停止發(fā)聲后，聲能密度衰減60 dB所經(jīng)歷的時間。混響時間是衡量房間音質(zhì)的重要參數(shù)之一，為了達(dá)到預(yù)期的聽音效果就要確保室內(nèi)具有合適的混響時間，即最佳混響時間。

房間要達(dá)到一個最佳的混響時間并不容易，影響混響時間的因素很多，還有許多不確定因素是我們無法控制的。以前人們只能通過改變室內(nèi)筑材或者通過室內(nèi)裝修來改變室內(nèi)混響，有它的局限性?，F(xiàn)在，隨著數(shù)字信號處理技術(shù)的快速發(fā)展，處理器的處理和運(yùn)算能力不斷提高，可以對聲源發(fā)出的聲音信號進(jìn)行數(shù)字化處理從而改善其混響時間。首先對室內(nèi)聲源產(chǎn)生的沖激響應(yīng)進(jìn)行測量來確定該房間的混響特性，然后通過某種數(shù)字濾波器對原始聲音進(jìn)行處理來達(dá)到人們所預(yù)期的混響效果，完成這種功能的數(shù)字濾波器稱為數(shù)字混響器。房間的脈沖響應(yīng)如圖1所示。

數(shù)字混響器的沖激響應(yīng)具有很長的長度，原始的聲音信號經(jīng)過數(shù)字混響器處理后就會產(chǎn)生一定混響效果，理論上這個數(shù)字混響的能產(chǎn)生如圖1所示的房間脈沖響應(yīng)形式。

圖1 房間的脈沖響應(yīng)

從20世紀(jì)70年代開始，數(shù)字混響器才被人們所熟知，到現(xiàn)在已衍生出了很多種模型，主要包括Moorer混響模型、Schroeder混響模型、梳狀濾波器混響模型、全通濾波器混響模型等。其中Schroeder的數(shù)字混響模型影響比較大，人們在它的基礎(chǔ)上對數(shù)字混響模型進(jìn)行了更進(jìn)一步的發(fā)展。下面我們著重討論兩種數(shù)字混響模型對混響的處理。

1 梳狀濾波器混響模型

梳狀濾波器的內(nèi)部有一些按一定頻率間隔相同排列的通帶和阻帶，只有在一定頻率范圍的信號才能通過。由于梳狀濾波器形成的特性曲線呈梳狀形，于是稱其為梳狀濾波器。

根據(jù)理論上對混響生成的分析，對于后期的混響聲，我們可以得出簡單的數(shù)字混響器模型：

式（1）中，x（n），y（n）為分別為輸入輸出信號，D為延遲時間，a為衰減系數(shù)（a﹤l）。那么該混響器的系統(tǒng)函數(shù)為：

用等比數(shù)列求和公式可將（2）式變?yōu)椋?/p>

由式（3）可以得到該混響器的系統(tǒng)框圖如圖2所示。

圖2 梳妝濾波器系統(tǒng)框圖

在Matlab平臺下對其進(jìn)行仿真，可得到它的頻率特性和脈沖響應(yīng)（取m=8，g=0.5），如圖3、圖4所示。

圖3 梳狀濾波器頻率特性

從頻率特性圖3來看，其幅頻特性呈明顯梳狀，這樣系統(tǒng)接收到不同頻率的信號時就會產(chǎn)生波動情況，出現(xiàn)聲染色現(xiàn)象，使聲音信號失真，經(jīng)過處理得到的聲音聽起來就會很不自然。由圖4可見，該脈沖響應(yīng)幅度隨時間呈間斷的逐漸遞減的情形，實(shí)際的混響情形與上面圖1所示的走向基本相似，但是可以明顯看到與實(shí)際極不相符的情況，就是該混響系統(tǒng)產(chǎn)生的聲音譜線是平均分布的，且非常稀疏。這意味著它產(chǎn)生的回聲密度幾乎是不變的且非常低，而且它不會隨著時間的延長而加強(qiáng)，這樣不能達(dá)到我們所期望的混響效果。

圖4 梳狀濾波器的脈沖響應(yīng)

2 全通濾波器混響模型

全通濾波器具有平坦的頻率響應(yīng)，它不能衰減任何頻率的信號，也就是說它沒有濾波作用。雖然它不能改變輸入信號的頻率特性，但它會改變輸入信號的相位。

根據(jù)全通濾波器的上述特性人們設(shè)計出一種全通濾波器的混響模型來改善梳狀濾波器混響模型出現(xiàn)的問題，其系統(tǒng)框圖如圖5所示。

圖5 全通濾波器系統(tǒng)框圖

由圖5可以看出其改進(jìn)的方法就是在梳狀濾波器混響模型的基礎(chǔ)上加了一個反饋支路。圖中Z-m表示延時，g（g﹤l）表示反饋增益。圖5所示混響模型的系統(tǒng)函數(shù)為：

為了更清楚的了解全通濾波器的混響特性，我們在Matlab平臺下對該混響模型的頻率特性和脈沖響應(yīng)進(jìn)行仿真（取m=8，g=0.5）結(jié)果如圖6、圖7所示。

圖6 全通濾波器的頻率特性

圖7 全通濾波器的脈沖響應(yīng)

由圖6所示的頻率特性可見，幅度譜呈一水平直線，即使輸入不同頻率的聲音信號也不會有波動的現(xiàn)象，聲音信號也不會失真，也就不會出現(xiàn)聲染色現(xiàn)象；它的相位譜是呈線性波浪狀的，隨時間逐漸遞減，這樣輸入的聲音信號就會有不同的相位，使聲音信號疊加或者抵消，在頻譜中產(chǎn)生許多波峰和波谷，從而產(chǎn)生不同的聲音，達(dá)到預(yù)期混響效果。

從圖7所示的脈沖響應(yīng)來看，與梳狀濾波器的脈沖響應(yīng)基本相似，問題仍是系統(tǒng)產(chǎn)生的回聲密度不夠高，且不會隨時間的延長而加強(qiáng)。這樣還是不能達(dá)到預(yù)期的混響效果。

現(xiàn)在我們要解決的問題就明顯凸現(xiàn)出來了，由于全通濾波器比梳狀濾波器的優(yōu)勢在于能解決聲染色現(xiàn)象，每一個濾波器頻譜都是全通的，現(xiàn)在如果把幾個全通濾波器串聯(lián)在一起形成一個整體，理論上這個整體的頻譜也是全通的，這時各個濾波器產(chǎn)生的回聲就能疊加在一起從而使回聲密度加大。這樣既能保證沒有聲染色現(xiàn)象發(fā)生，又能解決回聲密度低的問題。在MATLAB平臺下對此改進(jìn)的濾波器模型進(jìn)行仿真，結(jié)果如圖8所示。由于它的頻率特性與全通濾波器模型一致，現(xiàn)在只考慮它脈沖響應(yīng)的變化。

圖8 模型改進(jìn)后的脈沖響應(yīng)

由圖8可以很清楚的看到這種方法得到脈沖響應(yīng)圖與上面的幾乎一樣，也就是說這樣系統(tǒng)模型產(chǎn)生的回聲密度還是不夠大，與實(shí)際混響效果還是相差甚遠(yuǎn)，為了解決其回聲密度低的問題還需要繼續(xù)分析研究。

在實(shí)際房間內(nèi)人們聽到的是混響音效的多次反射波，這些聲波隨著多次反射而逐漸衰減，從仿真結(jié)果可以看到每隔相同的時間脈沖和相位就會遞減，也就是說該混響模型產(chǎn)生的聲波在房間內(nèi)每次反射時間間隔相同。如果把幾個全通濾波器嵌套在一起，即將圖5中的Z-m部分換成一個全通濾波器。為了與實(shí)際情況更加貼近，先把三個全通濾波器嵌套在一起，每個濾波器都取不同參數(shù)，然后在MATLAB平臺下進(jìn)行仿真，得到它的頻譜特性、脈沖響應(yīng)如圖9、圖10所示。3個全通濾波器的參數(shù)分別取m=8，g= 0.5；m=6，g=0.6；m=8，g=0.8。

圖9 三個全通濾波器嵌套后系統(tǒng)的頻率特性

圖10 三個全通濾波器嵌套后系統(tǒng)的脈沖響應(yīng)

由全通濾波器的特性可知，嵌套在一起后形成的新的濾波器模型應(yīng)該也是全通的，這與圖9正好吻合，可見該混響模型系統(tǒng)也不會產(chǎn)生聲染色現(xiàn)象。由圖10可以明顯看到隨著時間推移，線條越來越密集，也就是回聲密度逐漸變大，這正好符合實(shí)際情況，從而解決了回聲密度小的問題。這是因?yàn)楫?dāng)替換原系統(tǒng)內(nèi)部的全通濾波器生成回聲時，回聲就可以通過外部濾波器的反饋環(huán)路重新作為輸入信號輸入，生成的回聲數(shù)目遠(yuǎn)比原全通濾波器模型的多，回聲密度就會變大，使得信號產(chǎn)生的混響時間變長，聽音效果就會顯得比較自然了。

3 總結(jié)

本文詳細(xì)分析了兩種數(shù)字混響模型，它們的模型結(jié)構(gòu)都比較簡單，并且容易實(shí)現(xiàn)，但都存在著不足，因此對其進(jìn)行了改進(jìn)。改進(jìn)后的混響模型結(jié)構(gòu)一定程度上解決了前者的一些缺陷，但這種結(jié)構(gòu)的缺點(diǎn)是太復(fù)雜，不容易實(shí)現(xiàn)。所以對數(shù)字混響模型必須進(jìn)一步的進(jìn)行研究，使其在保證室內(nèi)聽音效果最佳的同時，又能簡單方便的運(yùn)用到實(shí)踐生活中。

[1]李正周.MATLAB數(shù)字信號處理與應(yīng)用[M].北京：清華大學(xué)出版社，2008：89-127.

[2]李承智.人工混響算法[J].電聲技術(shù)，2002，26（10）：11-13.

[3]周婧潔，陳建平，吉曉東.幾種數(shù)字混響模型的比較與分析[J].南通大學(xué)學(xué)報：自然科學(xué)版，2010，9（2）：23-28.

[4]王興國，蔣偉峰.數(shù)字混響器設(shè)計[J].電子學(xué)報，2000，28（7）：84-87.

[5]彭妙顏，張承云.人工混響的設(shè)計方法[J].電子技術(shù)，2006，30（1）：11-13.

[6]GardnerWG.EfficientConvolutionwithout Input-OutputDelay[J].Journalof theAudio Engineering Society，1995，43：127-136.

〔責(zé)任編輯 李海〕

Analysis and I mprovement of the D igital S ound R reverberatol

JIANG Xiao-yun，LIWen-jing
（School of Physics and Electronics Science，Shanxi Datong University，Datong Shanxi，037009）

Reverberation time is the important parameters tomeasure quality.In order to get the optimal reverberation time，adopting digital reverb units to sound signal processing is currently one of the regular use of themethod.In this paper，the two kinds of reverb units on detailed analysismodel，and analyses the improvement.Simulation results prove the validity of the theoretical analysis.

r everberation time；digital sound reverberator；filter

TN713

A

1674-0874（2012）04-0028-04

2012-05-02

山西省自然科學(xué)基金資助項目[59637050]

姜曉云（1969-），男，山西大同人，碩士，講師，研究方向：聲信號處理。