障板表面圓環(huán)陣頻率不變響應(yīng)波束形成方法研究

張英豪， 汪勇， 楊益新

(西北工業(yè)大學(xué) 航海學(xué)院，陜西 西安 710072)

在聲吶、雷達(dá)、通信以及語音工程等領(lǐng)域常常需要利用傳感器陣列接收寬帶目標(biāo)信號(hào)，并利用波束形成技術(shù)以達(dá)到抑制噪聲和干擾、提高輸出信噪比、改善方位估計(jì)精度等目的。在應(yīng)用常規(guī)波束形成方法處理接收到的寬帶信號(hào)時(shí)，常常會(huì)造成信號(hào)畸變，且指向性較小。針對(duì)這個(gè)問題，超指向性頻率不變響應(yīng)波束形成方法能夠在寬帶范圍提高指向性[1-4]，且可保證在一定頻率范圍內(nèi)保證信號(hào)無失真。在實(shí)際應(yīng)用中，常常把圓環(huán)陣安裝固定在某些水下平臺(tái)上，而平臺(tái)產(chǎn)生的散射聲場(chǎng)會(huì)使陣元處接收聲場(chǎng)發(fā)生畸變，進(jìn)而影響波束方法的性能。

近年來，很多學(xué)者研究了障板條件下基陣波束形成方法和性能分析。針對(duì)剛性圓球表面的球形傳感器陣列，Elko等[5]提出了相位模態(tài)域波束形成方法，利用陣型的對(duì)稱性，將入射聲場(chǎng)進(jìn)行正交分解，得到正交的球諧波分量，然后乘以相應(yīng)的系數(shù)，獲得所需要的指向性。鄢社鋒等[6-7]提出了多約束優(yōu)化方法，可以在不同的性能指標(biāo)之間合理折中，進(jìn)而獲得想要的波束。Rafael[8]提出的相位模態(tài)域波束形成方法也可以在性能指標(biāo)間進(jìn)行折中。此外，對(duì)于彈性圓柱表面圓環(huán)陣；Zou等[9]利用相位模態(tài)理論研究了此時(shí)圓環(huán)形矢量水聽器陣列的方位估計(jì)方法；楊德森等[10]也在相位模態(tài)域?qū)χ螐椥哉习灞砻鎴A環(huán)形矢量水聽器陣列的方位估計(jì)性能進(jìn)行了研究，并給出了實(shí)驗(yàn)結(jié)果；汪勇等[11]利用噪聲協(xié)方差矩陣的循環(huán)矩陣特性，最小化期望波束和合成波束的均分誤差，得到了合成寬帶頻率不變響應(yīng)波束各個(gè)頻率下的權(quán)值向量，進(jìn)而獲得合成頻率不變響應(yīng)波束。

考慮障板對(duì)波束方法性能的影響，Daigle等[12]對(duì)剛性和阻抗性球體及無限長(zhǎng)柱體表面圓環(huán)陣的常規(guī)方法進(jìn)行了深入研究，結(jié)果表明散射界面可以改善常規(guī)方法的性能；Bertilone等[13]求解了無限長(zhǎng)剛性柱形障板影響下幾種類型的噪聲互譜矩陣，并將其應(yīng)用于MVDR波束形成器中，同時(shí)對(duì)其性能進(jìn)行了分析。然而這些工作只針對(duì)常規(guī)和 MVDR波束形成器，沒有給出障板影響下圓環(huán)陣的穩(wěn)健超指向性方法下的性能變化。

本文主要研究圓環(huán)陣安裝在不同圓柱表面時(shí)頻率不變響應(yīng)波束方法的性能變化。首先考慮聲透明和剛性圓柱障板下的頻率不變響應(yīng)波束形成，并對(duì)獲得的仿真結(jié)果進(jìn)行分析。實(shí)際中，只有當(dāng)界面材料的聲特性阻抗遠(yuǎn)大于介質(zhì)聲特性阻抗時(shí)，剛性假設(shè)才能近似成立，而這雖在空氣中容易做到，但在水下則不易滿足。因此本文構(gòu)建了彈性雙層圓柱障板聲場(chǎng)模型，并給出相應(yīng)的散射聲場(chǎng)仿真結(jié)果，進(jìn)而獲得圓柱障板在不同材質(zhì)和厚度時(shí)的頻率不變響應(yīng)波束，并分析不同情況下的性能變化。

1 信號(hào)模型

考慮如圖1所示半徑為a的無限長(zhǎng)聲透明圓柱表面M元圓環(huán)陣，陣元均為無指向性傳感器，且θm=90°。單位幅度平面波從方向(θ,φ)入射，基于聲散射理論，圓環(huán)陣第m號(hào)傳感器接收到的聲壓：

pm(θ,φ)=eika sin θ cos(φ-φm)

(1)

圖1 圓環(huán)陣示意Fig.1 Circular array

式(1)可展開為：

(2)

式中：φm=mβ,β=2π/M,bn=inJn(krm)；rm是m號(hào)陣元到圓心間距；Jn是第n階柱形Bessel函數(shù)。

當(dāng)圓環(huán)陣位于無限長(zhǎng)剛性圓柱表面時(shí)有：

(3)

波束響應(yīng)是指波束形成器對(duì)某方位入射單位功率平面波信號(hào)的響應(yīng)，表達(dá)式為[14]：

B(θ,φ)=wH(θ0,φ0)P(θ,φ)=
ωH(θ0,φ0)E(θ,φ)

(4)

式中：w(θ,φ)是陣元域權(quán)值向量;E是變換后的新陣列流形向量且E=CP，ω(θ,φ)是變換后的新權(quán)值向量，有：

w=CHω，

式中：

C=[v0…vm…vM-1]H

vm=M-1/2[1eimβ…ei(M-1)mβ]T

指向性因子(DF)可用來描述陣列對(duì)空間各向同性均勻噪聲的抑制能力，可表示為：

(5)

且指向性指數(shù)DI=10lg(DF)，其中R是二維各向同性均勻噪聲場(chǎng)的歸一化噪聲互譜矩陣，其元素為：

ρmm′=J0(2πdmm′/λ)

(6)

Λ=CRCH=diag(λ0,λ1,…λM-1)

(7)

白噪聲增益(WNG)可以用來衡量波束形成方法的穩(wěn)健性，其定義為：

WNG=‖w‖-2=‖ω‖-2

(8)

2 合成波束權(quán)值向量

在二維各向同性均勻噪聲場(chǎng)時(shí)，僅考慮水平范圍內(nèi)的性能水平，俯仰角為90°，此時(shí)期望波束可以表示為：

(9)

式中：

Pd=[e-iN(φ-φ0)…1…eiN(φ-φ0)]T

a=[a-N…a0…aN]T

N是選取的最大階數(shù)，為正整數(shù)。

通過最小化合成波束與期望波束之間的均方誤差，可推導(dǎo)獲得合成波束所需權(quán)值向量元素，合理選取最大階數(shù)N，當(dāng)N

(10)

不顯示的值取零。

3 彈性雙層介質(zhì)圓柱聲場(chǎng)模型

在實(shí)際應(yīng)用中，圓環(huán)陣經(jīng)常安裝在金屬制球體或柱體表面。在空氣中時(shí)，障板剛性假設(shè)容易滿足(金屬的聲特性阻抗一般遠(yuǎn)大于空氣特性阻抗)，然而在水下，剛性假設(shè)很難成立。而且當(dāng)傳感器固定在障板外時(shí)，其位置很容易發(fā)生改變，因此可考慮如圖2所示的彈性柱形障板結(jié)構(gòu)。

圖2 具有覆蓋層的柱體表面圓環(huán)陣俯視圖Fig.2 Circular array on the surface of cylinder with covering layer

該彈性柱體由內(nèi)外2個(gè)覆蓋層組成。外層材料是粘彈性的，一般是橡膠之類，可以固定傳感器；內(nèi)層材料是彈性的，一般是金屬質(zhì)地；內(nèi)腔中介質(zhì)為空氣。此時(shí)對(duì)于傳感器距圓心的距離大小近似等于圓柱的外徑長(zhǎng)度Rc。

考慮平面波入射在圖2所示的雙層介質(zhì)圓柱體上的情況。圓柱體垂直放置在圓柱坐標(biāo)系的原點(diǎn)，整體沒入水中，用ρ和c分別表示水的密度和水中的聲速。半徑為Rc的外覆蓋層由密度為ρc的粘彈性材料制成，其彈性和粘性Lamé常數(shù)分別為λe2、μe2、λv2、μv2。圓柱框架由半徑為Rb密度為ρb彈性Lamé常數(shù)為λe3、μe3的彈性材料制成。半徑為Ra的圓筒內(nèi)腔充滿了密度為ρa(bǔ)ir和聲速為cair的空氣。

由彈性體散射的相關(guān)理論可得，柱體表面圓環(huán)陣的第m號(hào)陣元接收到的平面波信號(hào)仍如式(2)所示，此時(shí)的bn表達(dá)式發(fā)生了改變，有[15]：

(11)

式(11)中Fn,1是向量的第1個(gè)元素：

F=[Fn,1Fn,2…Fn,10]T

(12)

且有：

GF=A

(13)

式中G和A的元素由文獻(xiàn)[15]給出

聲場(chǎng)中任意點(diǎn)位置散射聲壓為：

(14)

彈性圓柱參數(shù)如表2所示，材質(zhì)選取為不銹鋼，有：Rc=0.5 m,Rb=0.40 m,Ra=0.35 m。

圓柱障板周圍散射聲場(chǎng)如圖3所示。圓柱半徑有a=0.5 m，4幅圖分別表示不同障板和頻率條件下的散射聲場(chǎng)，圖(a)、(b)為ka=2時(shí)，剛性圓柱和彈性圓柱周圍散射場(chǎng)，圖(c)、(d)為ka=4時(shí)，剛性圓柱和彈性圓柱周圍散射場(chǎng)。由圖3(a)和(b) 可得，低頻時(shí)，剛性障板和彈性障板周圍散射場(chǎng)大致相同。由圖4(c)和(d) 可得，在高頻時(shí)，剛性障板和彈性障板周圍散射場(chǎng)是有著出入的，水下雙層彈性圓柱不能近似看作是剛性的。

圖3 圓柱表面散射場(chǎng)Fig.3 Plot of the scattering field for cylinder

4 不同圓柱障板表面圓環(huán)陣頻率不變響應(yīng)波束形成性能分析

以12元均勻圓環(huán)陣為例進(jìn)行仿真，考慮聲透明、剛性圓柱和彈性圓柱障板3種情況。

期望權(quán)值向量和期望波束如表1和圖4所示。

表1 望波束的權(quán)值向量Table 1 Weight vector of desired beampatterns

對(duì)于圓柱障板表面圓環(huán)陣，歸一化噪聲互譜矩陣的元素為[12]：

(15)

本文仿真所研究的彈性雙層圓柱界面材料參數(shù)如表2所示，界面尺寸半徑與第3小節(jié)彈性圓柱相同，且材質(zhì)為不銹鋼。由式(10)可得聲透明圓柱、剛性圓柱和彈性圓柱表面圓環(huán)陣合成波束權(quán)值向量，進(jìn)而獲得寬帶合成波束，如圖5、6、7所示。

表2 彈性圓柱相關(guān)參數(shù)Table 2 Relevant parameters of elastic cylinder

圖4 期望波束Fig.4 Desired beampattern

圖5 聲透明圓柱表面圓環(huán)陣寬帶合成波束Fig.5 Broadband synthetic beampattern of acoustic transparent circular array

圖6 剛性圓柱表面圓環(huán)陣寬帶合成波束Fig.6 Broadband synthetic beampattern of circular array on rigid cylindrical surface

圖7 彈性圓柱表面圓環(huán)陣寬帶合成波束Fig.7 Broadband synthetic beampattern of circular array on elastic cylindrical surface

由圖5～8和圖9可知，不同圓柱障板表面圓環(huán)陣所得到的波束均有較好的頻率不變性。由圖8可知，頻率越低，穩(wěn)健性越差，而聲透明圓環(huán)陣在較高頻率范圍內(nèi)，白噪聲增益具有很多零點(diǎn)，當(dāng)剛性圓柱障板或者彈性圓柱障板存在時(shí)可以解決零點(diǎn)存在問題。相較聲透明圓環(huán)陣，彈性雙層障板與剛性障板在給定頻率范圍內(nèi)均可提升穩(wěn)健性。

圖8 白噪聲增益Fig.8 White noise gain

圖9 指向性指數(shù)Fig.9 Directivity index

5 彈性圓柱內(nèi)覆蓋層不同材質(zhì)時(shí)頻率不變響應(yīng)波束形成性能分析

本文仿真研究的彈性雙層圓柱界面材料參數(shù)如表2所示，選取不銹鋼，鋁合金和鈦合金作為內(nèi)覆蓋層的材料，進(jìn)行頻率不變響應(yīng)波束形成，并觀察分析性能變化。此時(shí)頻率不變響應(yīng)波束形成方法的指向性指數(shù)在頻率范圍內(nèi)均穩(wěn)定在7.2dB左右。

白噪聲增益，指向性指數(shù)隨頻率變化如圖10所示。界面尺寸為：

Ra=0.25 m,Rb=0.35 m,Rc=0.50 m

由圖10(a)可得，當(dāng)內(nèi)覆蓋層材質(zhì)為鋁合金時(shí)，白噪聲增益的極小值最小，下降趨勢(shì)最快，鈦合金其次，相比之下材質(zhì)為不銹鋼時(shí)，白噪聲增益隨頻率變化相對(duì)更穩(wěn)定。且白噪聲增益極小值點(diǎn)頻率向更高的頻率偏移。內(nèi)覆蓋層材質(zhì)選取為不銹鋼時(shí)，可以獲得略好的白噪聲增益，更貼近剛性障板情況下的白噪聲增益。由圖10(b)可知，當(dāng)材質(zhì)不同時(shí)，頻率不變響應(yīng)波束形成方法在給定頻率內(nèi)保持著良好的頻率不變性，指向性指數(shù)穩(wěn)定在7.2 dB。

圖10 內(nèi)覆蓋層不同材質(zhì)時(shí)的白噪聲增益和指向性指數(shù)Fig.10 WNGs and DIs of inner cover with different materials

6 彈性圓柱覆蓋層不同厚度時(shí)頻率不變響應(yīng)波束形成性能分析

本節(jié)彈性雙層圓柱障板外覆蓋層半徑有Rc=0.5 m，改變內(nèi)覆蓋層半徑Rb和內(nèi)腔半徑Ra，材料參數(shù)如表2所示。分析外覆蓋層和內(nèi)覆蓋層厚度對(duì)頻率不變響應(yīng)波束形成方法性能的影響。內(nèi)覆蓋層材質(zhì)選取為不銹鋼，外覆蓋層材質(zhì)為橡膠。

只改變內(nèi)覆蓋層半徑Rb，此時(shí)內(nèi)腔和外徑有Ra=0.3 m,Rc=0.5 m，此時(shí)彈性圓柱中不銹鋼和橡膠的厚度在發(fā)生變化，白噪聲增益變化如圖11(a)所示。隨著不銹鋼層厚度的增加和橡膠層厚度的減小，白噪聲增益在頻率范圍的極小值點(diǎn)向更高的頻率偏移，且白噪聲增益幅度變化越來越大。在低頻段，白噪聲增益與剛性圓柱的基本相同，表明在低頻段，彈性圓柱材料參數(shù)的改變對(duì)頻率不變響應(yīng)波束方法穩(wěn)健性影響較小。內(nèi)腔大小確定時(shí)，橡膠厚度越大，穩(wěn)健性越好。由圖11(b)可知改變彈性圓柱的內(nèi)外覆蓋層厚度大小，頻率不變響應(yīng)波束形成方法都有很好的頻率不變性。

圖11 內(nèi)腔大小固定時(shí)的白噪聲增益和指向性指數(shù)Fig.11 WNGs and DIs of fixed size intracavity

只改變內(nèi)腔半徑Ra，彈性圓柱的不銹鋼層厚度發(fā)生變化，此時(shí)Rb=0.42 m,Rc=0.50 m，橡膠厚度固定為0.08 m。白噪聲增益如圖12(a)所示。當(dāng)不銹鋼層厚度逐漸增加時(shí)，白噪聲增益的變化趨勢(shì)逐漸平緩，高頻處極小值點(diǎn)向高頻偏移且數(shù)值變大。當(dāng)不銹鋼層厚度越小時(shí)，在極小值點(diǎn)頻率向右偏移，白噪聲增益變化幅度越大，入射波頻率對(duì)頻率不變響應(yīng)波束形成方法穩(wěn)健性影響越大。當(dāng)橡膠厚度固定時(shí)，不銹鋼厚度的適當(dāng)增加可以使適用頻率范圍內(nèi)的白噪聲增益更加穩(wěn)定。由圖12(b)可知改變彈性圓柱的內(nèi)腔和內(nèi)覆蓋層厚度大小，頻率不變響應(yīng)波束形成方法都有著很好的頻率不變性。

圖12 橡膠厚度固定時(shí)的白噪聲增益和指向性指數(shù)Fig.12 WNGs and DIs of fixed rubber thickness

維持不銹鋼層厚度為0.1 m，改變內(nèi)腔半徑Ra值，Rb=Ra+0.1 m,Rc=0.5 m，白噪聲增益隨頻率變化如圖13(a)所示。由圖13(a)可知，當(dāng)不銹鋼厚度固定時(shí)，橡膠越薄，使得白噪聲增益在頻率范圍內(nèi)的極小值越小，白噪聲增益對(duì)頻率的變化越敏感，在一定頻率范圍內(nèi)使得頻率不變響應(yīng)波束方法穩(wěn)健性變差。由圖13(b)可得，改變彈性圓柱的外覆蓋層和內(nèi)腔尺寸大小，頻率不變響應(yīng)波束形成方法都有很好的頻率不變性。

圖13 不銹鋼層厚度固定時(shí)的白噪聲增益和指向性指數(shù)Fig.13 WNGs and DIs of fixed steel thickness

從本節(jié)白噪聲增益隨頻率變化的3幅圖中，可以看到剛性障板條件下和彈性圓柱障板條件下的白噪聲增益在一定頻率范圍內(nèi)是相吻合的。整體上剛性障板可以看作是理想情況下的彈性障板條件，但實(shí)際在較高頻，彈性圓柱障板的白噪聲增益受圓柱參數(shù)影響較大，很難達(dá)到剛性障板時(shí)的白噪聲增益，因此研究分析彈性障板條件下傳感器陣列頻率不變響應(yīng)波束形成方法的性能是有意義的。

7 結(jié)論

1)在給定頻率范圍內(nèi)，彈性雙層障板與剛性障板均可保證頻率不變響應(yīng)波束形成方法良好的頻率不變性，且彈性雙層障板覆蓋層總的尺寸固定時(shí)厚度和材料的改變對(duì)頻率不變性影響不大。且整體上剛性障板可以看作是理想的彈性障板條件，但在水下彈性圓柱障板很難達(dá)成剛性假設(shè)下的性能。

2)在較低頻率范圍內(nèi)，彈性雙層障板與剛性障板均可改善頻率不變響應(yīng)波束形成方法的穩(wěn)健性，且彈性雙層障板的參數(shù)對(duì)性能的影響可以忽略，此時(shí)信號(hào)頻率越高，圓環(huán)陣頻率不變響應(yīng)波束方法穩(wěn)健性越好，且白噪聲增益隨頻率增加迅速提升。

3)在較窄的高頻范圍內(nèi)，不銹鋼層和橡膠層的加厚，均有利于改善頻率不變響應(yīng)波束形成方法的整體穩(wěn)健性，減小了白噪聲增益極小值和該點(diǎn)頻率左右變化幅度。當(dāng)不銹鋼層和橡膠層總的厚度確定時(shí)，此時(shí)橡膠層的加厚更能提升頻率不變響應(yīng)波束形成方法的穩(wěn)健性，減小白噪聲增益極小值和該點(diǎn)頻率左右變化幅度。

總之，當(dāng)使用本文中尺寸確定的彈性障板表面?zhèn)鞲衅麝嚵袝r(shí)，在允許的范圍內(nèi)，內(nèi)腔應(yīng)盡可能小，橡膠層應(yīng)盡可能厚，不銹鋼層厚度應(yīng)盡可能薄，但其強(qiáng)度需足以支撐和固定整體框架。