杜博凱,曾向陽,洪 汐
(西北工業(yè)大學(xué)航海學(xué)院,陜西 西安 710072)
聲場重放(Sound Field Reproduction)旨在重現(xiàn)一個(gè)完整的聲學(xué)環(huán)境或者合成一個(gè)期望的聲學(xué)場景??臻g聲場重放技術(shù)是一種通過耳機(jī)或者揚(yáng)聲器陣列為聽者提供空間立體聲體驗(yàn)的方法。使用揚(yáng)聲器陣列的重放,可以構(gòu)建較大區(qū)域的聲學(xué)環(huán)境以及構(gòu)建相對(duì)較廣的聽音區(qū)域,近年來受到了廣泛的關(guān)注。
目前主流的揚(yáng)聲器聲場重放方法為波場合成[1-3](Wave Field Synthesis, WFS)與高階高保真立體聲重放[4-6](Higher Order Ambisonics, HOA)。WFS 基于的是惠更斯原理(Huygens principle)以及 Kirchhoff-Helmholtz積分方程,即如果已知封閉曲面上聲壓和法向質(zhì)點(diǎn)振速分布,那么可以在其表面布置一系列單極子和偶極子點(diǎn)聲源或平面波模型推導(dǎo)聲源的驅(qū)動(dòng)函數(shù)。但是為了避免空間混疊,WFS通常需要幾十到幾千個(gè)揚(yáng)聲器來實(shí)現(xiàn)聲場的重放,實(shí)用性與可操作性低,目前沒有商業(yè)性大規(guī)模推廣。HOA是另外一種常見的聲場重放方法,即將聲場在球坐標(biāo)系下分解成一組球諧函數(shù)(基函數(shù))的疊加,利用初級(jí)聲場和重放聲場之間的球諧函數(shù)模態(tài)進(jìn)行匹配求解出次級(jí)聲源的驅(qū)動(dòng)函數(shù)。但是 HOA技術(shù)也存在一些缺陷,為了提高重放聲場的精度,需要使用大量傳聲器采集聲場,且重放揚(yáng)聲器需要特殊的排列形狀如球形或圓形等,實(shí)用性受到了一定的限制。
由于WFS與HOA系統(tǒng)在信號(hào)錄制與重放過程中的限制,一些研究也關(guān)注基于聲壓匹配法[7-8]的聲場重放方法。不同于WFS與HOA是基于潛在的空間連續(xù)物理問題的解析方法,聲壓匹配法(Pressure Matching, PM)是基于空間離散策略的方法。通常選取位于重放目標(biāo)區(qū)域周圍或在該區(qū)域內(nèi)的控制點(diǎn),設(shè)計(jì)揚(yáng)聲器驅(qū)動(dòng)函數(shù)使得在這些點(diǎn)上盡可能準(zhǔn)確地合成目標(biāo)聲源的聲壓場。但是為了提高重放精度,需要大量的測量,此方法也存在一定限制。
等效源法最早是在聲場再現(xiàn)(Sound field reconstruction)的近場聲全息(Nearfield Acoustic Holography, NAH)中應(yīng)用的[9],其基本思想是利用一定形式的傳聲器陣列在距離聲源很近的表面采集聲場信息,然后通過空間聲場變換計(jì)算聲源和整個(gè)聲場的聲壓、質(zhì)點(diǎn)振速及聲強(qiáng)等分布。1989年Koopmann等[10]提出了波疊加法,即等效源法,該方法將單極子或偶極子等簡單源作為等效源配置在振動(dòng)體內(nèi),通過近場測量求解等效源的源強(qiáng),再現(xiàn)整個(gè)聲場分布。該方法適用于任意形狀的聲場,且不存在邊界元法(Boundary Element Method, BEM)的缺陷。其團(tuán)隊(duì)隨后[11-12]研究了等效源的位置、數(shù)目等對(duì)聲場再現(xiàn)精度的影響。在之前的工作中,我們利用多層等效源再現(xiàn)封閉空間聲場,在600 Hz以上比等效源法獲得更低的再現(xiàn)誤差[13]。隨后,文獻(xiàn)[14]利用等效源方法進(jìn)行了多區(qū)域聲場重放問題的研究,相對(duì)于傳統(tǒng)的聲壓匹配法保持了重放精度的同時(shí)實(shí)現(xiàn)了更安靜的暗區(qū)。然而,作者之前的工作對(duì)于等效源的數(shù)目選取,位置確定以及混響條件對(duì)方法性能影響等關(guān)鍵內(nèi)容缺乏明確的討論。同時(shí),相對(duì)于傳統(tǒng)的聲壓匹配法,也未提升在重放區(qū)域內(nèi)的重放精度。
本文基于之前的工作,提出一種基于等效源法插值的聲場重放方法。此方法可以有效降低傳聲器采樣工作量,同時(shí)在目標(biāo)區(qū)域內(nèi)實(shí)現(xiàn)較低誤差的聲場重放。首先給出等效源法插值聲場重放的理論推導(dǎo)。作為進(jìn)一步討論,本文對(duì)等效源數(shù)目、等效源距離以及不同混響條件下的重放誤差進(jìn)行了詳細(xì)的比較與分析。最后分別對(duì)比了低采樣與高采樣聲壓匹配法與本文提出方法的重放誤差。
假設(shè)目標(biāo)聲場由一個(gè)傳聲器陣列采樣,傳聲器測得的聲壓為:
其中,M為傳聲器的數(shù)量。揚(yáng)聲器到重放區(qū)域中的傳聲器位置處的聲傳遞函數(shù)表示為
傳聲器點(diǎn)處的聲壓一般為實(shí)際聲場景中的測量的聲壓或期望合成聲場的聲壓,為已知量。當(dāng)?shù)刃г吹奈恢门c傳聲器的位置固定時(shí),聲傳遞函數(shù)矩陣G確定。通過求解式(7)可以得到等效源的權(quán)重q。通常情況下,設(shè)置的等效聲源的數(shù)量大于傳聲器的測量數(shù)量,則求解過程具有不適定性,需要借助正則化方法來消除矩陣中較小的奇異值引起的不穩(wěn)定擾動(dòng)。通常使用的正則化方法為截?cái)嗥娈愔捣纸?Truncated Singular Value Decompositon, TSVD)與Tikhonov正則化等。對(duì)于上述問題,等效源權(quán)重q的求解可以等效為2范數(shù)優(yōu)化問題:
其中,λ為正則化過程中的超參數(shù)。通常選取參數(shù)的方法包括L曲線法、廣義交叉驗(yàn)證法(Generalized Cross Validation, GCV)法等[12]。
為清晰描述本方法,圖1為聲壓匹配法與等效源法插值聲場重放方法的流程。首先,揚(yáng)聲器播放的聲信號(hào)通過傳聲器陣列采樣。對(duì)于聲壓匹配法,可以通過直接計(jì)算得到揚(yáng)聲器的驅(qū)動(dòng)函數(shù)。對(duì)于等效源法插值聲場重放,需要先計(jì)算揚(yáng)聲器至重放區(qū)域內(nèi)虛擬點(diǎn)的聲傳遞函數(shù)插值,得到插值后的聲場,從而計(jì)算揚(yáng)聲器驅(qū)動(dòng)函數(shù)。
圖1 聲壓匹配法與等效源法插值聲場重放方法流程圖Fig.1 The flowchart of the pressure matching method and the equivalent source reproduction method
為驗(yàn)證所提出的等效源(Equivalent Source Method, ESM)插值重放方法,以3D矩形封閉混響房間內(nèi)的平面聲場重放為例進(jìn)行研究。如圖 2所示,房間的三維尺寸為4 m×4 m×3 m。坐標(biāo)原點(diǎn)為揚(yáng)聲器陣列、傳聲器陣列以及等效聲源的中心點(diǎn)。目標(biāo)聲源的位置為rs(0,2,0)。揚(yáng)聲器陣列為一個(gè)均勻分布、L=32個(gè)揚(yáng)聲器組成的圓陣,其中揚(yáng)聲器都為單極子聲源。重放區(qū)域在揚(yáng)聲器陣列內(nèi)部,為0.5 m×0.5 m 的正方形區(qū)域,其中傳聲器采樣點(diǎn)共M=16個(gè),且均勻分布在重放區(qū)域的邊界處。傳聲器陣列測量的聲信號(hào)為Pmic。重放區(qū)域內(nèi)部設(shè)置100個(gè)均勻分布的虛擬點(diǎn)。揚(yáng)聲器外圍設(shè)置 N=48個(gè)等效源,呈圓形均勻分布。鑒于揚(yáng)聲器陣列的空間奈奎斯特(Nyquist)頻率為876 Hz,本文所關(guān)注的聲場重放頻率為100~1 000 Hz。自由場中聲傳遞函數(shù)由自由場格林(Green)函數(shù)給出,混響環(huán)境下的聲傳遞函數(shù)則由鏡像源方法(Image Source Method,ISM)[15]給出。另外,傳聲器的測量聲壓添加信噪比為30 dB的高斯白噪聲。
圖2 揚(yáng)聲器、傳聲器、等效聲源和目標(biāo)聲源位置圖Fig.2 The positions of loudspeakers, microphones, equivalent sources and the target sound source
其中:S為聲場重放區(qū)域的面積;pre為重放聲壓;pdes為目標(biāo)重放聲壓;E( f)的單位為dB。
在計(jì)算等效源法插值聲場重放的過程中,等效源的數(shù)目、等效源的距離以及聲場內(nèi)部虛擬點(diǎn)的數(shù)目與距離,均對(duì)重放結(jié)果產(chǎn)生影響。在傳聲器采樣位置、數(shù)量,目標(biāo)聲場內(nèi)部虛擬點(diǎn)以及等效源的距離均相等時(shí),增加等效源的數(shù)量,對(duì)于不同頻率下的重放結(jié)果都有所提升,如圖3所示。
圖3 不同頻率下等效源數(shù)目對(duì)重放誤差的影響Fig.3 Effect of the number of equivalent sources on reproduction errors at different frequencies
對(duì)于300 Hz,當(dāng)?shù)刃г磾?shù)目增加到14時(shí),重放誤差降低,而后,隨著等效源的數(shù)目增加,重放誤差趨于穩(wěn)定。而對(duì)于800 Hz,當(dāng)?shù)刃г磾?shù)目低于26時(shí),重放誤差存在波動(dòng),當(dāng)?shù)刃г磾?shù)目多于 26時(shí),重放誤差趨于穩(wěn)定。對(duì)于不同的頻率,在保證重放效果穩(wěn)定的前提下,所需的最小等效源數(shù)目是不一定的。在本文中,對(duì)于所關(guān)注的頻率區(qū)間,選取的等效源數(shù)目48個(gè)。
圖4所示為不同等效源距離下的重放誤差。隨著等效源距離的增加,聲場的重放誤差降低,尤其是在500 Hz以下的低頻段部分。但是,當(dāng)?shù)刃г淳嚯x為揚(yáng)聲器陣列距離的兩倍即2 m時(shí),對(duì)于頻率在500 Hz以下的部分,重放效果最好,500 Hz以上部分,重放誤差與等效源距離2.3 m的重放誤差接近,而等效源距離為2.3 m時(shí),低于500 Hz頻率的重放誤差比等效源距離2.0 m稍高。本文中,對(duì)于之后的仿真驗(yàn)證,均選擇等效源距離為2.0 m。
圖4 不同等效源距離下的重放誤差對(duì)比Fig.4 Comparison of reproduction errors fort different equivalent source distances
在封閉空間中的聲場重放過程,會(huì)受到房間壁面反射聲的影響,從而降低重放效果。圖5對(duì)比了自由場條件下與封閉空間混響條件下的重放效果。圖5中封閉空間的吸聲系數(shù)為0.6。在200~300 Hz這一頻率段,封閉空間混響條件下的重放誤差較大,在部分頻點(diǎn)處超過了-15 dB。在600 Hz以下頻段處,自由場中的重放效果穩(wěn)定,波動(dòng)小,而封閉空間混響條件下的重放效果波動(dòng)大。結(jié)果表明,對(duì)于低頻部分,房間壁面的反射更易影響聲場重放效果。在600 Hz以上頻率部分,混響條件下的重放誤差與自由場下的誤差走勢基本相同,重放誤差大。封閉空間中,不同的吸聲系數(shù),對(duì)于聲場重放的效果也有不同的影響。圖6對(duì)比了不同吸聲系數(shù)下的聲場重放效果。在 4 m×4 m×3 m的矩形房間中,0.2、0.4、0.6和0.8的吸聲系數(shù)對(duì)應(yīng)的混響時(shí)間分別為0.48、0.24、0.16和0.12 s。由圖6中的結(jié)果可知,吸聲系數(shù)越大,混響時(shí)間越短,重放效果越好。但是,混響存在的環(huán)境中,對(duì)于不同的吸聲系數(shù),除了個(gè)別頻點(diǎn),重放誤差相互之間的差距在10 dB內(nèi)。即使對(duì)于吸聲系數(shù)為0.8的情況,重放效果整體仍比自由場條件下的聲場重放效果差。但是對(duì)于整個(gè)頻率區(qū)間,不同混響環(huán)境下的重放誤差均低于-10 dB在可以接受的重放誤差區(qū)間內(nèi)。
圖5 不同聲學(xué)環(huán)境中重放誤差對(duì)比Fig.5 Comparison of reproduction errors under different acoustic conditions
圖6 不同混響條件下重放誤差對(duì)比Fig.6 Comparison of reproduction errors under different reverberant conditions
為對(duì)比驗(yàn)證所提出的等效源法插值的聲場重放方法,分別計(jì)算了聲壓匹配法在低采樣與高采樣情況下的重放誤差,同時(shí)對(duì)比了本文提出的等效源法插值方法,結(jié)果如圖7所示。圖7中低采樣為目標(biāo)聲場由 16個(gè)傳聲器測量的聲場信號(hào),高采樣為在低采樣測量聲場的基礎(chǔ)上,目標(biāo)聲場增加100個(gè)測量點(diǎn)獲得的聲場信號(hào)。由結(jié)果可以得出,在低采樣的情況下,重放誤差較大,部分頻率處的重放誤差超過-10 dB。在500 Hz以下,高采樣的聲壓匹配法與本文提出的等效源法插值方法得到的重放誤差幾乎相等。這意味著使用本文提出的方法可以大大降低聲傳遞函數(shù)測量的工作量同時(shí)獲得良好的重放效果。在500 Hz以上的頻率,等效源法插值方法得到的重放誤差大于高采樣的聲壓匹配法的重放誤差。這是因?yàn)殡S著頻率的升高,在虛擬點(diǎn)的插值聲壓存在誤差,然而相對(duì)于僅采用低采樣的聲壓匹配法的結(jié)果,本文提出的方法仍然獲得了更好的重放效果。
圖7 PM法低采樣、高采樣和ESM插值聲場重放誤差對(duì)比Fig.7 Comparison of reproduction errors for the PM method of low and high sampling and the interpolation method
本文提出了一種等效源法插值聲場重放方法。相對(duì)于傳統(tǒng)的聲壓匹配法,在傳聲器低采樣的條件下,利用等效源法對(duì)揚(yáng)聲器陣列聲傳遞函數(shù)進(jìn)行插值,提升了重放效果。在傳聲器陣列、揚(yáng)聲器陣列不變的條件下,通過仿真實(shí)驗(yàn),選擇合適的等效源數(shù)目與等效源距離,對(duì)比了本文提出的方法與傳統(tǒng)聲壓匹配法的重放效果。結(jié)果表明,本文提出的方法比聲壓匹配法的重放誤差小,提升了重放效果。并且對(duì)比了自由場與封閉空間混響條件下的重放誤差。由于受到封閉空間內(nèi)壁面的聲反射影響,重放效果降低,尤其是對(duì)600 Hz以下低頻段,在進(jìn)一步的研究中,會(huì)針對(duì)封閉空間的房間補(bǔ)償進(jìn)行了展開,降低混響條件下的重放誤差。