一種分級子陣聚焦mvdr波束形成方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于聲納數(shù)字信號處理領(lǐng)域，特別涉及一種分級子陣聚焦MVDR(MinimumV arianceDistortionlessResponse，最小方差無失真響應(yīng)算法）波束形成方法，可以應(yīng)用于 各種線列陣聲納系統(tǒng)，實現(xiàn)快速高精度的波束形成。
【背景技術(shù)】
[0002] 受瑞利限的限制，對于小孔徑基陣，采用傳統(tǒng)的波達估計算法難以分辨方位上接 近的 2 個或多個目標。1969 年 Capon 在文獻"High-resolution frequency-wavenumbersp ectrumanalysis"中提出的MVDR算法是一種在理論上具有高分辨性能的波達估計算法，但 是運算過程中會涉及到高維矩陣求逆，復(fù)雜度很高，而且算法對信噪比要求比較高，穩(wěn)定性 也較差。為 了解決這些問題，Swingler. D. N 在文獻"Alow-complexityMVDRbeamformerforu sewithshortobservationtimes"中提出一種連續(xù)不重疊不等長子陣波束形成方法，田彪等 在文獻"多子陣高分辨實時波達估計算法研究"中提出了一種類似的基于均勻鄰接子陣列 處理方法的MVDR算法，這些方法在降低計算量的同時還能保持MVDR算法的高分辨性能，實 時性較好。另外，通過子陣處理提高了單個陣元的信噪比，對信噪比的寬容性也更好。
[0003] 雖然子陣聚焦MVDR算法在一定程度上有效降低了運算復(fù)雜度，但是，距離一般聲 納系統(tǒng)所要求的實時性還有一定差距。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明目的在于，克服常規(guī)波束形成面臨的分辨率有限，無法實現(xiàn)對目標的高精 度測向，以及MVDR等自適應(yīng)高分辨波束形成方法面臨的運算復(fù)雜度高、穩(wěn)定性差的問題， 從而提出了一種分級子陣聚焦MVDR波束形成方法，以同時滿足聲納系統(tǒng)波束形成在實時 性、高分辨能力和寬容性方面的需求。
[0005] 本發(fā)明考慮普通波束頻域?qū)拵Рㄊ纬煞椒ㄔ趯崟r性和寬容性方面的巨大優(yōu)勢， 為了解決小孔徑水平基陣高精度被動測向的問題，本發(fā)明提出了一種分級波束形成的思路 是：將普通快速頻域?qū)拵Рㄊ纬膳c子陣聚焦MVDR算法結(jié)合起來，由快速頻域?qū)拵Рㄊ?成算法完成對目標方位粗測，再由子陣聚焦MVDR算法完成對目標方位的精測，從而提高整 個聲納定向系統(tǒng)的性能，同時滿足實時性、高分辨和寬容性的需要。一種線列陣聲納裝置， 可以是一條拖曳陣，也可以是舷側(cè)陣，所述線列陣由多個水聽器組成。這里設(shè)陣元數(shù)目N，陣 元間距d，目標入射方向Θ，陣元接收信號表示為 x(t)，信號頻帶范圍為fmin~fmax;聲速為 c，數(shù)據(jù)快拍長度為L。
[0006] 本發(fā)明的方法詳細流程如附圖1所示，首先對線列陣接收到的數(shù)據(jù)通過一次時域 FFT變換和空域相位補償，完成快速寬帶頻域波束形成，由于普通波束形成分辨力有限，這 一步只是完成對目標方位的初步估計，可以將波束間隔設(shè)置的比較寬。然后，對普通波束形 成曲線進行峰值檢測，發(fā)現(xiàn)目標后，對頻域陣列信號進行子陣劃分，只對目標附近的方位做 子陣聚焦MVDR波束形成，然后進行二次方位檢測，得到更精確的目標方位。
[0007] 具體步驟描述如下：
[0008] 11)線列陣數(shù)據(jù)輸入
[0009] 用線陣接收空間信號，得到M個陣元的時域信號，取L采樣點長度的數(shù)據(jù)形成一個 數(shù)據(jù)快拍，作為數(shù)據(jù)輸入。
[0010] X (t) = Lx^k) X2 (k) ··· xM(k)]T
[0011] 其中，xm(k) = [xm(l) xm(2)…Xm(L)]
[0012] 12)快速寬帶頻域波束形成
[0013] 對數(shù)據(jù)快拍X(t)在頻帶范圍fmin~fmax內(nèi)做快速頻域?qū)拵Рㄊ纬?，具體過程如 下。
[0014] 首先對陣列信號x(t)的每一個陣元信號做時域FFT，截取fmin~f max范圍的頻域 信號表不如下：
[0015] X(D = Lx^f) x2(f) ··· xM(f)]T [0016]其中，Xni (f) = [Xni U …Xni (f；J ]
[0017] 對不同頻率分量&構(gòu)造導(dǎo)引矢量：
【主權(quán)項】
1. 一種分級子陣聚焦MVDR波束形成方法，該方法將普通快速頻域?qū)拵Рㄊ纬伤惴?與子陣聚焦MVDR算法結(jié)合起來，首先由快速頻域?qū)拵Рㄊ纬赏瓿蓪δ繕朔轿淮譁y，再劃 分子陣，最終由子陣聚焦MVDR算法完成對目標方位的精測，實現(xiàn)快速高精度的波束形成。
2. 根據(jù)權(quán)利要求1所述的分級子陣聚焦MVDR波束形成方法，其特征在于，所述方法包 括以下步驟： 1) 首先，對線列陣聲納裝置接收到的數(shù)據(jù)通過一次時域FFT變換和空域相位補償，完 成快速寬帶頻域波束形成，完成對目標方位的初步估計； 2) 然后，對普通波束形成曲線進行峰值檢測，發(fā)現(xiàn)目標后，對頻域陣列信號進行子陣劃 分，并且只對目標附近的方位做子陣聚焦MVDR波束形成，進行二次方位檢測，得到更精確 的目標方位。
3. 根據(jù)權(quán)利要求2所述的分級子陣聚焦MVDR波束形成方法，其特征在于，所述線列陣 聲納裝置由多個水聽器陣元組成，是一條拖曳陣或是舷側(cè)陣；設(shè)水聽器陣元數(shù)目N，陣元間 距d，目標入射方向Θ，陣元接收信號表示為X (t)，信號頻帶范圍為fmin~fmax ;聲速為c，數(shù) 據(jù)快拍長度為L。
4. 根據(jù)權(quán)利要求3所述的分級子陣聚焦MVDR波束形成方法，其特征在于，所述方法的 具體步驟如下： 11) 線列陣數(shù)據(jù)輸入； 線列陣聲納裝置接收空間信號，得到M個陣元的時域信號，取L采樣點長度的數(shù)據(jù)形成 一個數(shù)據(jù)快拍，作為數(shù)據(jù)輸入： X(t) = [X1 (k) x2(k) ··· xM(k)]T； 其中，xm(k) = [xm(l) xm(2)…xm(L)]; 12) 快速寬帶頻域波束形成； 對數(shù)據(jù)快拍X(t)在頻帶范圍fmin~fmax內(nèi)做快速頻域?qū)拵Рㄊ纬桑? 首先，對陣列信號X (t)的每一個陣元信號做時域FFT，截取fmin~fmax范圍的頻域信號 表示如下： X(f) = [xi(f) x2(f) *·· xM(f)]T； 其中，Xm(f) - [Xm (fmin)…Xm (fmax)]; 對不同頻率分量fi構(gòu)造導(dǎo)引矢量： 然后，進行相位補償，得到各個頻率分量的波束輸出： y (f, Q1)= w(fi，Θ x)X(f)； 最后，將不同頻率分量的各個波束輸出進行平方相加得到波束圖： PaA^) = y(fA)-ynf：A)
21) -次方位檢測； 對于波束曲線Pcbf(Q1) (1 = 1，2，···，180° )進行方位檢測，首先，對波束曲線進行峰 值檢測，對每一個峰值計算檢測系數(shù)DT :
其中，PraF(U表示目標峰值，EtPc^ed]表示波束曲線的均值，δ [Pcbf(Q1)]表示 波束曲線的起伏值； 假設(shè)檢測到了 Q個目標，方位分別表示為eq(k= 1，2, ···，〇); 22)子陣聚焦MVDR波束形成； 對檢測到的每個目標Θ q(k = 1，2,…，Q)，對其附近的方位做子陣聚焦MVDR波束形成： 將線列陣劃分為N個子陣，每個子陣含有Ms個陣元，則有 ? 最后，一個子陣長度為 定義X' n(t)表示第η個子陣的時域陣列信號，則第η個子陣接收信號的頻域矢量X' η 表示為： X " = +1(,) ; X 2(/);…；Χ?η_ηΜ，Μ;(/)]. 其中，Xm(f) = [XJf1) Xm(f2)…Xm(fK)]，K表示頻率分量的數(shù)目； 對觀測方向Θ，構(gòu)造頻點fk對應(yīng)的子陣方向矢量為： ?(θ f) - [l e-J^fkd^ieVce-J27rJk (A/a-I}·^Cos(^)Zc jT 對各子陣頻域陣列信號作波束方向相位補償，得到子陣輸出波束為： Υη(θ) = [YnAf1) Yn(0，f2)…Yn(0，fK)]; 其中，Yn(9，fk) =aH(0，fk)X'n(fk); 將Υη( θ )視為第η個虛擬陣元的陣元信號，則虛擬陣的陣列信號表示為： Υ(Θ) = [Υ"Θ) ; Υ2(θ);…；γΝ(θ)]; 按照平面波假設(shè)，以第一個子陣為基準，對每一個子陣進行延時聚焦： > | I | 備 _ c }2π I1(U-I)M Jcos(O) c "· (2π fk(n-\)M J ζο&(θ) cj ^ 這樣，得到子陣聚焦對齊后虛擬陣的協(xié)方差矩陣估計為： R(的=士細)_々〃(的； M 其中，9(0=[夂(的；十2(0;…；Vw(0二 * 則Θ方向的輸出功率譜估計為：
23)對子陣聚焦MVDR得到的波束曲線θ )進行二次方位檢測，得到更精確的目 標方位Θ,'（k= 1，2, ···，〇)并輸出。
5. 根據(jù)權(quán)利要求4所述的分級子陣聚焦MVDR波束形成方法，其特征在于，所述步驟 22)中，一般將線列陣劃分成3~5個子陣。
6. 根據(jù)權(quán)利要求4所述的分級子陣聚焦MVDR波束形成方法，其特征在于，所述步驟 21) 中，一般取檢測系數(shù)DT=6dB，這時檢測概率為95. 4%。
7. 根據(jù)權(quán)利要求4所述的分級子陣聚焦MVDR波束形成方法，其特征在于，所述步驟 22) 中，對檢測到的每個目標0q(k=l，2，~，Q)附近的方位做子陣聚焦MVDR波束形成，所 述的附近的方位取值為Θ Θ~Θ q+A Θ，Λ Θ取值范圍為5~10。。
8. 根據(jù)權(quán)利要求4所述的分級子陣聚焦MVDR波束形成方法，其特征在于，所述步驟 23) 中的二次方位檢測同一次方位檢測的方法完全一致。
【專利摘要】本發(fā)明涉及一種分級子陣聚焦MVDR波束形成方法，該方法將普通快速頻域?qū)拵Рㄊ纬伤惴ㄅc子陣聚焦MVDR算法結(jié)合起來，首先由快速頻域?qū)拵Рㄊ纬赏瓿蓪δ繕朔轿淮譁y，再劃分子陣，最終由子陣聚焦MVDR算法完成對目標方位的精測，實現(xiàn)快速高精度的波束形成；包括：1）首先，對線列陣聲納裝置接收到的數(shù)據(jù)通過一次時域FFT變換和空域相位補償，完成快速寬帶頻域波束形成，完成對目標方位的初步估計；2）然后，對普通波束形成曲線進行峰值檢測，發(fā)現(xiàn)目標后，對頻域陣列信號進行子陣劃分，并且只對目標附近的方位做子陣聚焦MVDR波束形成，進行二次方位檢測，得到更精確的目標方位。本發(fā)明提高了整個聲納定向系統(tǒng)的性能，同時滿足實時性、高分辨和寬容性的需要。
【IPC分類】G01S7-539
【公開號】CN104730513
【申請?zhí)枴緾N201310706004
【發(fā)明人】李崢, 黃海寧
【申請人】中國科學(xué)院聲學(xué)研究所
【公開日】2015年6月24日
【申請日】2013年12月19日