国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

基于水下高速小目標(biāo)末彈道測(cè)量的實(shí)孔徑成像聲吶波束形成研究

2019-11-25 14:25張思宇何心怡
艦船科學(xué)技術(shù) 2019年10期
關(guān)鍵詞:旁瓣聲吶波束

張思宇,何心怡,陳 菁,王 磊

(1.海軍研究院,北京 100161;2.中國(guó)科學(xué)院聲學(xué)研究所,北京 100097)

0 引 言

在聲吶系統(tǒng)的設(shè)計(jì)過(guò)程中,波束形成是一個(gè)核心基礎(chǔ)環(huán)節(jié),利用一定幾何形狀的基陣和信號(hào)處理技術(shù),使發(fā)射和接收的信號(hào)形成指向性波束。在發(fā)射信號(hào)時(shí),通過(guò)指向性波束提高聲吶的作用距離;在接收信號(hào)時(shí),在空間域中抑制環(huán)境噪聲和混響干擾,提高聲吶接收機(jī)輸入端的信噪比,獲得更好的聲吶圖像和目標(biāo)的分辨率,其直接表征為旁瓣級(jí),旁瓣級(jí)越低,系統(tǒng)抑制旁瓣區(qū)域的干擾能力越強(qiáng)。

在沒(méi)有水下彈道測(cè)量系統(tǒng)的條件下,設(shè)計(jì)測(cè)量水下高速小目標(biāo)末彈道的聲吶成像系統(tǒng)的主要考慮:1)目標(biāo)接近成像系統(tǒng)的方向未知,可能從任一方向接近;2)目標(biāo)運(yùn)動(dòng)速度快,回波多普勒偏移量大,可能達(dá)到寬帶信號(hào)的標(biāo)準(zhǔn);3)目標(biāo)距離近,成像系統(tǒng)受到混響的干擾強(qiáng)等因素的限制[1]。在聲吶系統(tǒng)的陣型設(shè)計(jì)中,只有使用均勻排布的圓形陣列,才具有 0°~360°的方位角范圍內(nèi),形成均勻指向性波束的能力,避免目標(biāo)的方位模糊,圓陣波束形成數(shù)學(xué)模型簡(jiǎn)單,設(shè)計(jì)的工程實(shí)現(xiàn)性強(qiáng)。但是圓陣的常規(guī)波束形成的旁瓣級(jí)高,圖像的目標(biāo)分辨率低,彈道測(cè)量誤差大,本文以一種對(duì)均勻圓陣進(jìn)行Dolph-Chebyshev 加權(quán)的波束形成方法,并給出在圓陣?yán)么朔úㄊ纬傻募s束條件,并通過(guò)仿真實(shí)驗(yàn)驗(yàn)證方法的可行性。

1 常規(guī)波束形成模型

由于聲信號(hào)到達(dá)基陣各個(gè)陣元的距離不同,導(dǎo)致各個(gè)陣元輸出的信號(hào)也不同,需要通過(guò)延時(shí)等信號(hào)處理方法補(bǔ)償差異,這是常規(guī)波束形成的基本思想[2-3]。假設(shè)入射信號(hào)是平面波,基陣是N 個(gè)陣元組成。以其中某個(gè)陣元為參考點(diǎn), 設(shè)第 i個(gè)陣元接收的信號(hào)為s[t +τi(θ0)], 其中 θ0為預(yù)設(shè)的信號(hào)入射角,對(duì)各通道信號(hào)延時(shí)補(bǔ)償 τi(θ0) , 使各路輸出信號(hào)均為 s(t), 再對(duì)N 路信號(hào)求和、平方和積分就得到,其中為信號(hào)功率。此時(shí),信號(hào)如從 θ角入射,則第 i個(gè)陣元輸出信號(hào)為 s[t +τi(θ)?τi(θ0)],陣列輸出為:

圖 1 常規(guī)波束形成模型Fig.1 Conventional beamforming model

2 圓陣波束形成的信號(hào)模型

目前,陣元發(fā)射信號(hào)受到換能器諧振頻點(diǎn)的限制,一般只能發(fā)射窄帶信號(hào)。接收信號(hào)則根據(jù)信號(hào)的不同分為窄帶信號(hào)模型和寬帶信號(hào)模型,通常當(dāng)目標(biāo)回波信號(hào)帶寬大于信號(hào)中心頻率的1/10 時(shí),被認(rèn)為是寬帶信號(hào),由于發(fā)射信號(hào)和接收信號(hào)的波束形成原理相同[3],文中以接收的信號(hào)模型為例說(shuō)明。

被測(cè)目標(biāo)到成像系統(tǒng)的距離一般是遠(yuǎn)大于成像系統(tǒng)的圓陣直徑,陣元接收的信號(hào)屬于遠(yuǎn)場(chǎng)窄帶信號(hào)。假設(shè)回波信號(hào)為 s(t) , 圓陣的陣元數(shù)為 N,圓陣接收的信號(hào)可以表示為:

式中, f0為 信號(hào)的中心頻率, a (t) 為信號(hào)幅度, φ (t)為信號(hào)相位,其解析信號(hào)可以表示為:

假設(shè)接收信號(hào)的時(shí)延量為τ,圓陣接收信號(hào)為:

但是,窄帶信號(hào)的幅度是慢變化的,則a (t?τ)=a(t)??梢?jiàn),窄帶信號(hào)時(shí)延只會(huì)引起信號(hào)復(fù)包絡(luò)的相移變化,對(duì)信號(hào)的幅值的影響基本可以忽略。假設(shè)圓陣的參考陣元的接收信號(hào)為第i個(gè)陣元接收信號(hào)相對(duì)參考陣元接收信號(hào)的時(shí)延量為 τi,其信號(hào)可表示為:

圓陣在第n 次快拍接收到的信號(hào)模型可用矩陣表示為:

3 圓陣波束形成的數(shù)理模型

3.1 均勻加權(quán)波束形成

假設(shè)圓陣半徑為 R, 在其圓周上等間距排布 N 個(gè)陣元,以原點(diǎn)作為參考點(diǎn),依次確定信號(hào)到達(dá)各陣元的時(shí)延量,如圖2 所示。

圖 2 均勻圓陣陣元分布圖Fig.2 Distribution of uniform circular array elements

其中, ?R 為參考原點(diǎn)和第i 個(gè)陣元的波程差。則第i 個(gè)陣元和參考原點(diǎn)之間的相位差為:

整個(gè)均勻圓陣的波束響應(yīng)可以表示為:

式中, si(θ,φ)為 第 i 個(gè) 陣元的接收信號(hào), ωi為 第 i 個(gè)陣元接收信號(hào)的加權(quán)值。

在均勻圓陣上,采用均勻加權(quán)法進(jìn)行波束形成,則有:

整個(gè)均勻圓陣的波束響應(yīng)可以表示為:

3.2 旁瓣抑制加權(quán)波束形成

通過(guò)Dolph-Chebyshev 加權(quán)法抑制旁瓣高度,改善圓陣波束形成的性能,但是以上加權(quán)法只適用于線性陣列,故需要通過(guò)相位模式激勵(lì)法將圓陣等效成線列陣,再進(jìn)行后續(xù)的波束形成[4-5]。為了方便討論波束性能,令則方位角為 φ,此時(shí)參考陣元 i的接收信號(hào)響應(yīng)為:

定義矩陣 J 和 F分別為:

圖 3 圓陣響應(yīng)模式轉(zhuǎn)換模型圖Fig.3 Model diagram of response mode transformation of circular array

經(jīng)過(guò)轉(zhuǎn)化后,陣元信號(hào)響應(yīng)矩陣為:

形成了線列陣的響應(yīng)模式,具有范德蒙行列式的結(jié)構(gòu)形式,在此基礎(chǔ)上,可以按照Dolph-Chebyshev 加權(quán)法對(duì)輸出信號(hào) yi(n),i = ?h,···,h進(jìn)行波束形成,假設(shè)波束指向 φ1方 向,設(shè)置權(quán)向量 W:

則整個(gè)整理的波束響應(yīng)為:

B(φ)

即,波束響應(yīng) 可以轉(zhuǎn)化為:

令 β為主瓣幅值和旁瓣幅值的比值,將式(22)表示成根值 z ∈[?1,1]且 還有 2h項(xiàng)的Chebyshev 多項(xiàng)式:

式 中, z = z0cos(u(φ)),z ∈[?z0,z0] 。 其 中, z0有T2h(z0)= β >1決定,可以進(jìn)一步得到:

圓陣響應(yīng)模式可表示為:

此時(shí),圓陣波束響應(yīng)模式已經(jīng)等效為均勻線列陣Dolph-Chebyshev 多項(xiàng)式的響應(yīng)模式。因?yàn)?u(φ) 和 φ1成線性關(guān)系,對(duì)于任一取值的 φ1, 圓陣波束響應(yīng) B(φ)主瓣寬度都相同,即任一波束指向方向的圓陣波束響應(yīng)的主瓣寬度一致。對(duì)于任一方向接近的目標(biāo),其波束響應(yīng)的強(qiáng)度相同,有利于均衡圖像的背景,提高目標(biāo)的圖像分辨率。

值得注意的是上述陣列模式等效轉(zhuǎn)換需要滿(mǎn)足矩陣 BJF可逆的條件,即圓陣的陣元數(shù)圓陣的陣元間距必須小于回波信號(hào)的半波長(zhǎng)。下面通過(guò)仿真實(shí)驗(yàn)經(jīng)過(guò)Dolph-Chebyshev 加權(quán)后的旁瓣抑制效果。

4 圓陣波束形成仿真實(shí)驗(yàn)及性能分析

仿真實(shí)驗(yàn)1:假設(shè)均勻圓陣共有33 個(gè)各向同性的陣元,圓陣半徑為 R, 接收信號(hào)的波長(zhǎng)為 λ。要求波束指向 (0°,0°)方向,經(jīng)過(guò)Matlab 仿真得到圓陣的波束指向性圖,如圖4 所示。

由圖4 可知,圓陣三維波束形成圖可以獲得均勻加權(quán)圓陣的俯仰角和方位角信息,反映了均勻加權(quán)波束形成性能。但是,實(shí)孔徑聲吶成像系統(tǒng)和被測(cè)目標(biāo)基本在同一深度,故方位角信息更值得關(guān)注。為了方便分析圓陣波束形成水平方向的指向性特征及討論不同陣元數(shù)和半徑的波束指向性能。

圖 4 均勻加權(quán)圓陣三維波束形成圖Fig.4 Uniform weighted circular array three-dimensional beamforming

圖 5 均勻加權(quán)圓陣二維波束形成圖Fig.5 Two-dimensional beamforming of uniformly weighted circular array

由圖5 可知,均勻加權(quán)圓陣水平方向的半功率波束寬度1 1.5°, 最大旁瓣高度約為 ? 8 dB,遠(yuǎn)離預(yù)成方向的旁瓣呈現(xiàn)遞減趨勢(shì)。

仿真實(shí)驗(yàn)3:其他實(shí)驗(yàn)假設(shè)條件不變,保持 R/λ=0.66 不變,改變陣元數(shù) N,得到圖6。

圖 6 信號(hào)頻率一定不同陣元數(shù)圓陣二維波束形成圖Fig.6 Two-dimensional beamforming of circular arrays with different signal frequencies

由圖6 可知,由于圓陣的非線性特點(diǎn),在調(diào)整圓陣半徑和接收信號(hào)波長(zhǎng)比值不變的情況下,無(wú)論怎樣改變圓陣陣元的數(shù)量,均勻加權(quán)圓陣的第一旁瓣高度總是約為-8 dB,明顯高于線列陣的第一旁瓣-13.5 dB的高度。高旁瓣會(huì)引起高的虛警概率,造成大的彈道測(cè)量誤差,應(yīng)當(dāng)在設(shè)計(jì)中采用適當(dāng)?shù)姆椒ㄒ种婆园辏WC測(cè)量精度。

仿真實(shí)驗(yàn)4:實(shí)驗(yàn)1 的其他假設(shè)條件不變,改變圓陣接收信號(hào)頻率,即改變圓陣的半徑和信號(hào)波長(zhǎng)比,得到圖7。

圖 7 圓陣半徑一定不同頻率信號(hào)二維波束形成圖Fig.7 Two-dimensional beamforming of signals with certain circular array radius and different frequencies

由圖7 所示,在圓陣陣形固定的條件下,提高發(fā)射信號(hào)頻率可以減小波束形成的主瓣寬度,改善波束指向性能,但是最大旁瓣高度仍為-8 dB 左右,虛警概率和軌跡測(cè)量依然較大,且提高發(fā)射信號(hào)頻率會(huì)導(dǎo)致信號(hào)吸收損失的增強(qiáng),限制了聲吶系統(tǒng)的作用距離,需要慎重選擇發(fā)射頻率。

仿真實(shí)驗(yàn)5:設(shè)圓陣具有32 個(gè)均勻排布的陣元且陣元各向同性,半徑波長(zhǎng)比為0.75,主瓣幅值和旁瓣幅值比為 β =10( ? 20 dB ),令波束指向 ( 0°,0°)方向,俯仰角,最大相位模式 h = 5,則圓陣chebyshev 加權(quán)部分可以等效為11 個(gè)陣元均勻排布的線列陣,利用Matlab 仿真,分別得到均勻加權(quán)圓陣波束響應(yīng)圖和Dolph-Chebyshev 加權(quán)圓陣波束響應(yīng)圖,如圖8 和圖9 所示。

5 結(jié) 語(yǔ)

圖 8 均勻加權(quán)圓陣二維波束形成圖Fig.8 Two-dimensional beamforming diagram of uniform weighted circular array

圖 9 Dolph-Chebyshev 加權(quán)圓陣二維波束形成圖Fig.9 Dolph-Chebyshev weighted circular array 2-D beamforming

設(shè)計(jì)基于水下高速小目標(biāo)末彈道測(cè)量的實(shí)孔徑成像聲吶應(yīng)該優(yōu)先考慮均勻排布的圓陣,保證沒(méi)有觀測(cè)盲區(qū),提高觀測(cè)效率。在圓陣做波束形成時(shí),Dolph-Chebyshev 加權(quán)法能夠有效抑制旁瓣高度,克服均勻加權(quán)圓陣波束形成旁瓣高的缺陷,同時(shí)能夠保證在各個(gè)波束方向上,波束響應(yīng)的幅值相同,為后續(xù)圖像信號(hào)處理提供高質(zhì)量的輸入信號(hào)。但是,圓陣陣元數(shù)必須滿(mǎn)足才能等效成線列陣的響應(yīng)模式進(jìn)行Dolph-Chebyshev 加權(quán),否則常規(guī)波束形成的方法旁瓣抑制效果均不理想,故為了獲得良好的目標(biāo)回波圖像采用Dolph-Chebyshev 加權(quán)時(shí),必須合理設(shè)計(jì)圓陣半徑、發(fā)射信號(hào)頻率和陣元數(shù)量的關(guān)系以滿(mǎn)足 BJF矩陣可逆。

猜你喜歡
旁瓣聲吶波束
約束優(yōu)化的空間變跡算法的旁瓣抑制應(yīng)用
基于圓柱陣通信系統(tǒng)的廣義旁瓣對(duì)消算法
空投聲吶浮標(biāo)空中運(yùn)動(dòng)軌跡研究
基于共形超表面的波束聚焦研究
一種基于方位估計(jì)的雙站被動(dòng)聲吶航跡關(guān)聯(lián)方法
強(qiáng)干擾背景下水中微弱信號(hào)的垂直陣波束形成研究
60 GHz無(wú)線通信系統(tǒng)中臨近波束搜索算法研究
一種基于線性規(guī)劃的頻率編碼旁瓣抑制方法
超波束技術(shù)在岸基光纖陣中的應(yīng)用
航母中程反潛區(qū)聲吶艦對(duì)潛警戒能力分析