趙 晨，馮丹平，楊明明，高守勇

（91388部隊(duì)95分隊(duì)，廣東 湛江 524022）

0 引 言

近場(chǎng)聲全息技術(shù)（Near-field Acoustical Holography,NAH）是噪聲源識(shí)別、定位、分析的一項(xiàng)重要技術(shù)。其基本原理是通過測(cè)量近場(chǎng)全息面的聲場(chǎng)數(shù)據(jù)（如聲壓、振速、聲強(qiáng)等），采用一定的近場(chǎng)聲全息算法，反演重建聲源面聲場(chǎng)信息或者預(yù)測(cè)遠(yuǎn)場(chǎng)輻射聲場(chǎng)信息[1]。根據(jù)測(cè)量面形狀的不同，近場(chǎng)聲全息技術(shù)分為平面近場(chǎng)聲全息、柱面近場(chǎng)聲全息和球面近場(chǎng)聲全息[2-5]。其中基于聲壓測(cè)量的平面近場(chǎng)聲全息技術(shù)，對(duì)于大型目標(biāo)的噪聲源分析更為適用，并已經(jīng)應(yīng)用于實(shí)際測(cè)量工作中。

平面近場(chǎng)聲全息算法一直是聲學(xué)領(lǐng)域研究的熱點(diǎn)。常用平面近場(chǎng)聲全息算法所采用的基本方法包括空間二維傅里葉變換法、邊界元法、等效源法、統(tǒng)計(jì)最優(yōu)法等[6-10]?；诳臻g二維傅里葉變換法的近場(chǎng)聲全息（Near-Field Acoustic Holography Based on Spatial Two-Dimensional Fourier Transform,SFT-NAH）算法是近場(chǎng)聲全息算法中的一種典型算法，通過將全息面上聲壓進(jìn)行空間二維傅里葉變換再乘以格林函數(shù)，最后進(jìn)行空間二維傅里葉逆變換得到重建面或者預(yù)測(cè)面聲場(chǎng)信息。相比較其他方法，空間二維傅里葉變換法具有算法簡(jiǎn)單、不需要占用較大內(nèi)存和計(jì)算速度快等優(yōu)點(diǎn)。

采用空間二維傅里葉變換法進(jìn)行近場(chǎng)聲全息重建時(shí)，由于環(huán)境噪聲、系統(tǒng)誤差等因素的影響，對(duì)重建過程有重要影響的高波數(shù)域的倏逝波成分很容易被各類噪聲誤差淹沒，導(dǎo)致重建結(jié)果產(chǎn)生較大誤差。因此重建過程中必須進(jìn)行波數(shù)域的低通濾波，指數(shù)濾波器是基于空間二維傅里葉變換法的近場(chǎng)聲全息算法中常用的濾波方法[11]。該方法中有兩個(gè)關(guān)鍵參數(shù)，分別是指數(shù)濾波器截止波數(shù)和指數(shù)濾波器窗函數(shù)陡度系數(shù)。這兩個(gè)關(guān)鍵參數(shù)如果設(shè)置不合理就會(huì)嚴(yán)重影響指數(shù)濾波器濾波效果，進(jìn)而影響近場(chǎng)聲全息算法聲場(chǎng)重建效果。因此對(duì)其進(jìn)行設(shè)計(jì)優(yōu)化，對(duì)于提高聲全息算法的重建精度和分辨率具有重要意義。

本文以水下聲全息技術(shù)應(yīng)用為背景，以大孔徑近場(chǎng)聲全息面為分析模型，重點(diǎn)分析對(duì)應(yīng)不同聲源頻率，指數(shù)濾波器窗函數(shù)及截止波數(shù)的取值對(duì)聲全息重建誤差的影響。相鄰兩聲源相干引起的分辨率誤差是近場(chǎng)聲全息重建誤差的一個(gè)主要來源，有效降低該誤差是提高近場(chǎng)聲全息算法精度的一項(xiàng)重要工作，因此選取兩個(gè)相鄰聲源作為檢測(cè)的對(duì)象，針對(duì)基于空間二維傅里葉變換法的近場(chǎng)聲全息算法中指數(shù)濾波器的截止波數(shù)和窗函數(shù)陡度系數(shù)這兩個(gè)關(guān)鍵參數(shù)，采用最小二乘法進(jìn)行了設(shè)計(jì)優(yōu)化。優(yōu)化結(jié)果表明，有效提高了基于空間二維傅里葉變換法的近場(chǎng)聲全息算法的重建精度和分辨率。

1 基于空間二維傅里葉變換法的近場(chǎng)聲全息算法基本原理

由水聲傳播基本理論可知，在理想流體介質(zhì)中，不依賴于時(shí)間變量的穩(wěn)態(tài)聲場(chǎng)的亥姆霍茲（Helmholtz）方程為

在狄利克雷（Dirichlet）邊界條件下，格林函數(shù)為

對(duì)式（2）兩邊進(jìn)行空間二維傅里葉變換，可得：

FFT-NAH的基本思想是由全息面聲壓分布矩陣對(duì)其做空間二維傅里葉變換，將其轉(zhuǎn)化到波數(shù)域，經(jīng)波數(shù)域?yàn)V波后乘以對(duì)應(yīng)的格林函數(shù)矩陣，便可得到空間二維傅里葉變換后的重建面的聲壓分布矩陣，對(duì)其做空間二維傅里葉逆變換即可得到重建面聲壓分布。

2 指數(shù)濾波器

在近場(chǎng)聲全息理論中，當(dāng)全息面聲壓進(jìn)行二維傅里葉變換后將其與對(duì)應(yīng)指數(shù)濾波器相乘，以此來進(jìn)行波數(shù)域?yàn)V波，是提高近場(chǎng)聲全息算法重建精度的一種有效手段。本文在計(jì)算過程中首先對(duì)全息面聲壓數(shù)據(jù)進(jìn)行二維傅里葉變換，然后與對(duì)應(yīng)指數(shù)濾波器相乘進(jìn)行波數(shù)域?yàn)V波，最后進(jìn)行二維傅里葉逆變換得到重建面聲壓分布。

基于空間二維傅里葉變換法的近場(chǎng)聲全息算法中常用的指數(shù)濾波器計(jì)算公式為[11]

由式（5）可以看出，指數(shù)濾波器的設(shè)計(jì)有兩個(gè)關(guān)鍵參數(shù)，分別是窗函數(shù)的陡度系數(shù)α和濾波器的截止波數(shù)kc。已有研究中，通常情況下指數(shù)濾波器窗函數(shù)陡度系數(shù)α的取值范圍為0.1～0.2，指數(shù)濾波器截止波數(shù)kc經(jīng)驗(yàn)公式為：kc=0.6π/Δ，Δ為全息面上的測(cè)量間隔，對(duì)于本文研究工作具有一定的參考意義。本文將針對(duì)窗函數(shù)的陡度系數(shù)α和濾波器的截止波數(shù)kc這兩個(gè)關(guān)鍵參數(shù)的優(yōu)化設(shè)計(jì)，采用最小二乘法分別進(jìn)行仿真分析。

3 仿真分析

3.1 點(diǎn)聲源及全息測(cè)量面模型建立

以頻率為1 kHz、聲源半徑為0.01 m和間距為1.5 m的兩個(gè)模擬點(diǎn)聲源為聲源模型進(jìn)行仿真分析，結(jié)果如圖1所示。同時(shí)，建立規(guī)模為50 m×50 m、測(cè)量間隔分別為 50 m和 255 m、距離聲源中心面1 m的全息測(cè)量面仿真模型。以均勻分布噪聲模擬環(huán)境噪聲等各種噪聲，模擬信噪比為20 dB時(shí)仿真全息測(cè)量面的測(cè)量效果如圖2所示。

圖1 模擬點(diǎn)聲源聲壓傳播圖Fig.1 Sound pressure propagation diagram of simulated point sound source

圖2 模擬全息測(cè)量面測(cè)得聲壓分布圖Fig.2 Sound pressure distribution diagram measured by simulated holographic measuring surface

通常情況下濾波器窗函數(shù)陡度系數(shù)α的取值范圍為0.1～0.2，濾波器截止波數(shù)kc=0.6π/Δ，但通過仿真分析發(fā)現(xiàn)，當(dāng)聲源頻率不同時(shí)，選擇不同的濾波器窗函數(shù)陡度系數(shù)α和濾波器截止波數(shù)kc會(huì)有更優(yōu)的反演效果。例如，當(dāng)聲源頻率為1 kHz，時(shí)若α=0.12，kc=0.6π/Δ，則反演點(diǎn)聲源聲壓傳播圖如圖3所示。當(dāng)聲源頻率為1 kHz，時(shí)若α=0.06，kc=0.3π/Δ，則反演點(diǎn)聲源聲壓傳播圖如圖4所示。

圖3 聲源頻率為1 kHz，α=0.12，kc=0.6π/?，反演點(diǎn)聲源聲壓傳播圖Fig.3 The inversed sound pressure propagation diagram of 1 kHz point source for α=0.12 and kc=0.6π/?

圖4 聲源頻率為1 kHz，α=0.06，kc=0.3π/?，反演點(diǎn)聲源聲壓傳播圖Fig.4 The inversed sound pressure propagation diagram of 1 kHz point source for α=0.062 and kc=0.3π/?

3.2 指數(shù)濾波器參數(shù)設(shè)計(jì)分析

定義p1為模擬點(diǎn)聲源聲壓幅值分布，p2為反演點(diǎn)聲源聲壓幅值分布。定義點(diǎn)聲源聲壓反演結(jié)果的幅值誤差為[12]

3.2.1 濾波器截止波數(shù)設(shè)計(jì)分析

α值為0.12時(shí)，分別取kcΔ/π的值為0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9，對(duì)相應(yīng)的反演聲源面聲壓分布效果對(duì)比分析如表1所示。

表1 α=0.12 時(shí)，kc?/π 與 E 對(duì)應(yīng)關(guān)系Table 1 The correspondence between kcΔ/π and E,when α=0.12

α值為0.15時(shí)，分別取kcΔ/π的值為0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9，對(duì)相應(yīng)的反演聲源面聲壓分布效果對(duì)比分析如表2所示。

表2 α=0.15時(shí)，kc?/π與E對(duì)應(yīng)關(guān)系Table 2 The correspondence between kc?/π and E,when α=0.15

α值為0.18時(shí)，分別取 kcΔ /π的值為0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9，對(duì)相應(yīng)的反演聲源面聲壓分布效果對(duì)比分析如表3所示。

表3 α=0.18時(shí)，kc?/π與E對(duì)應(yīng)關(guān)系Table 3 The correspondence between kc?/π and E,when α=0.18

從表 1、表 2和表 3中不難發(fā)現(xiàn)，當(dāng)α值為0.12、0.15、0.18時(shí)，改變kcΔ/π的值，三種情況下均是當(dāng) kcΔ/π的值為0.3時(shí)，反演效果相對(duì)較好。

3.2.2 窗函數(shù)陡度系數(shù)設(shè)計(jì)分析

當(dāng)kcΔ/π的值為0.3時(shí)，分別取α值為0.02、0.04、0.06、0.08、0.1、0.12、0.14、0.16、0.18、0.2，對(duì)相應(yīng)的反演聲源面聲壓分布效果對(duì)比分析如表4所示。

表4 kc?/π=0.3時(shí)，α與E對(duì)應(yīng)關(guān)系Table 4 The correspondence between α and E,when kc?/π=0.3

kc/Δπ的值為0.2時(shí)，分別取α值為0.02、0.04、0.06、0.08、0.1、0.12、0.14、0.16、0.18、0.2、0.25、0.3、0.35，如表5所示，對(duì)相應(yīng)的反演聲源面聲壓分布效果進(jìn)行對(duì)比分析。

表5 kc?/π=0.2時(shí)，α與E對(duì)應(yīng)關(guān)系Table 5 The of correspondence between α and E,when kc?/π=0.2

kc/Δπ的值為0.4時(shí)，分別取α值為0.02、0.04、0.06、0.08、0.1、0.12、0.14、0.16、0.18、0.2，對(duì)相應(yīng)的反演聲源面聲壓分布效果對(duì)比分析如表6所示。

表6 kc?/π=0.4時(shí)，α與E對(duì)應(yīng)關(guān)系Table 6 The of correspondence between α and E,when kc?/π=0.4

通常情況下α的取值范圍為0.1～0.2，表4和表6中在對(duì)應(yīng)分析條件下，當(dāng)α的取值在0.02～0.2中已找到最優(yōu)取值，而在表5中，α取到了0.35，當(dāng)α=0.3時(shí)E出現(xiàn)最小值，這一結(jié)果說明針對(duì)不同的全息面參數(shù)設(shè)計(jì)、測(cè)量環(huán)境以及不同的聲源頻率，指數(shù)濾波器窗函數(shù)陡度系數(shù)的最優(yōu)取值并不局限于0.1～0.2。

從表1～6中數(shù)據(jù)不難看出，kcΔ/π的值為0.3、α的值為0.06時(shí)，反演效果最好。且通過對(duì)比可以看出，相比較窗函數(shù)陡度系數(shù)α的取值變化，濾波器截止波數(shù)kc的取值變化對(duì)反演效果影響更加明顯。

3.3 指數(shù)濾波器參數(shù)設(shè)計(jì)與聲源頻率對(duì)應(yīng)關(guān)系分析

通過進(jìn)一步仿真分析，可得到在相同點(diǎn)聲源面模型和全息面模型中，對(duì)應(yīng)不同的聲源頻率時(shí)，要取得最優(yōu)重建效果，相對(duì)應(yīng)的濾波器窗函數(shù)陡度系數(shù)α和濾波器截止波數(shù)kc取值如表7所示。

表7 對(duì)應(yīng)不同聲源頻率時(shí)kc?/π值，α與E對(duì)應(yīng)關(guān)系Table 7 The correspondence of kc?/π,α and E for the sound sources of different frequencies

從表 7中數(shù)據(jù)可以看出，聲源頻率在 100～1 500 Hz之間時(shí)，隨著聲源頻率的增加，指數(shù)濾波器窗函數(shù)的陡度系數(shù)α的最優(yōu)設(shè)置逐漸減小，濾波器截止波數(shù)kc的最優(yōu)設(shè)置逐漸增大。

4 結(jié) 論

本文對(duì)平面近場(chǎng)聲全息中指數(shù)濾波器窗函數(shù)設(shè)計(jì)進(jìn)行了分析和優(yōu)化，建立了頻率為1 kHz、聲源半徑為0.01 m和間距為1.5 m的兩個(gè)模擬點(diǎn)聲源模型，以及規(guī)模為50 m×50 m、測(cè)量間隔分別為50 m、255 m和與聲源中心面距離為1 m的全息測(cè)量面仿真模型。同時(shí)以均勻分布噪聲模擬環(huán)境噪聲。通過分析濾波器截止波數(shù)kc和濾波器窗函數(shù)陡度系數(shù)α的取值對(duì)反演誤差的影響，找出了kcΔ/π的值為0.3、α的值為0.06時(shí)，采用指數(shù)濾波器濾波后的反演效果最好。相比較窗函數(shù)陡度系數(shù)α的取值變化，濾波器截止波數(shù)kc的取值變化對(duì)反演效果影響更加明顯。同時(shí)發(fā)現(xiàn)，聲源頻率在100～1 500 Hz之間時(shí)，隨著聲源頻率的增加，指數(shù)濾波器窗函數(shù)的陡度系數(shù)α的最優(yōu)設(shè)置逐漸減小，濾波器截止波數(shù)kc的最優(yōu)設(shè)置逐漸增大?？偨Y(jié)了聲源頻率在100～1 500 Hz之間時(shí)，對(duì)應(yīng)不同聲源頻率，要取得最優(yōu)重建效果，濾波器窗函數(shù)陡度系數(shù)α和濾波器截止波數(shù)kc的最優(yōu)取值，對(duì)平面近場(chǎng)聲全息中指數(shù)濾波器窗函數(shù)的設(shè)計(jì)優(yōu)化具有重要意義。