基于SVD和罰函數(shù)約束的子陣級方向圖控制

欒曉明，肖華飛

(哈爾濱工程大學(xué) 信息與通信工程學(xué)院，黑龍江 哈爾濱150001 )

基于SVD和罰函數(shù)約束的子陣級方向圖控制

欒曉明，肖華飛

(哈爾濱工程大學(xué) 信息與通信工程學(xué)院，黑龍江 哈爾濱150001 )

對于子陣級自適應(yīng)波束形成中旁瓣電平較高的問題，提出一種基于奇異值分解和罰函數(shù)約束(SVD-PFC)的子陣級自適應(yīng)方向圖控制算法。該算法通過對輸入數(shù)據(jù)進行奇異值分解，得到干擾子空間和噪聲子空間，通過干擾子空間修正約束矩陣和約束響應(yīng)矢量，可在小快拍情況下有效抑制干擾，同時結(jié)合罰函數(shù)對自適應(yīng)方向圖約束，使其逼近期望靜態(tài)方向圖。該算法具有良好的副瓣特性，提高了小快拍情況下輸出的信干噪比，仿真實驗證明了該算法的可行性和有效性。

自適應(yīng)方向圖控制；奇異值分解；罰函數(shù)；子陣級

0 引言

子陣級自適應(yīng)波束形成可以有效地抑制干擾[1-3]，但在快拍數(shù)目較少的情況下，主瓣波束容易變形，副瓣電平會明顯升高。自適應(yīng)方向圖控制技術(shù)可以在抑制干擾的同時進行低副瓣控制，一定程度上可抵消自適應(yīng)加權(quán)矢量不能適應(yīng)突變而帶來的性能下降[4-7]。文獻[8-9]提出罰函數(shù)法來增加一個約束，使方向圖逼近靜態(tài)方向圖而降低副瓣。文獻[10-13]提出基于特征向量的線性約束最小方差的自適應(yīng)方向圖控制算法，有效抑制干擾，降低旁瓣。這幾種都是主要針對陣元級。在子陣級上，文獻[13-16]提出基于子陣的劃分方法，抑制旁瓣電平，但會造成一定的輸出信干噪比損失。文獻[17]提出一種基于子空間投影和最優(yōu)失配檢測法，可在沒有干擾的情況下得到靜態(tài)方向圖。文獻[18]提出基于陣列特征空間重構(gòu)的方法，需要求解信號協(xié)方差矩陣并對其進行分解，運算量較大。文獻[19]提出在子陣級噪聲歸一化的基礎(chǔ)上，進行正交投影，抑制旁瓣電平。

上述算法在抑制干擾的同時控制副瓣特性，但存在小快拍數(shù)時，主瓣變形，副瓣起伏較大的問題。此次提出一種基于奇異值分解和罰函數(shù)約束(SVD-PFC)的子陣級自適應(yīng)方向圖控制算法，可在小快拍條件下抑制干擾使其逼近期望靜態(tài)方向圖，主瓣保形，旁瓣降低，并能夠獲得較好的信干噪比性能。

1 子陣信號模型

設(shè)有一個M陣元的等間距直線陣，將其不均勻地劃為L個子陣，波束指向為θ0，設(shè)有K個遠場干擾入射信號，干擾信號的波達方向為θi，其中i=1,2,…,K。陣列從M維變換為L維的變換矩陣為：

T=Φ0T0,

(1)

(2)

式中，T0為M×L維的子陣形成矩陣，此矩陣的列向量兩兩相交，只有第l列中與第l個子陣中陣元序號相對應(yīng)的元素值為1，其余元素值為0。Φ0表示陣元級的移相處理，d為陣元間距，λ表示波長。

在子陣級上，采用LCMV進行波束形成，求得常規(guī)子陣級自適應(yīng)權(quán)矢量：

(3)

(4)

(5)

式中，Csub為子陣級的約束矩陣，fsub為約束響應(yīng)矢量，R為陣元級接收協(xié)方差矩陣。

在子陣級采用LCMV的常規(guī)算法，雖然可以抑制干擾，但在快拍數(shù)有限的情況下，其方向圖主瓣容易發(fā)生變形，旁瓣電平升高。

2 SVD子陣級自適應(yīng)方向圖控制算法

假設(shè)子陣級只有理想噪聲，并無干擾，可得理想噪聲下的子陣級靜態(tài)加權(quán)矢量為：

(6)

考慮陣列接收到快拍數(shù)為I，存在K個干擾信號，得到包含干擾和噪聲的接收數(shù)據(jù)矩陣X，經(jīng)過子陣級變換得到子陣級的接收數(shù)據(jù)Xsub，對其進行SVD分解得到：

(7)

如果按照奇異值大小對U進行分解，有：

(8)

此時，Ui為子陣級的干擾子空間及UP為其子陣級的正交補空間矩陣，通過一些合適的信號源估計方法，可以確定有K個干擾，則在左奇異特征值矢量U中，與K個大奇異值對應(yīng)的可以組成Ui，由剩下的左奇異特征值矢量組成UP。因子陣級干擾導(dǎo)向矢量asub和Ui的列向量張成為同一子空間，也就是子陣級的干擾子空間Ui和自適應(yīng)加權(quán)矢量必須與正交，就可以在干擾方向上形成零陷，即：

(9)

此時子陣級的約束矩陣和約束響應(yīng)矢量可以構(gòu)造如下：

(10)

將式(10)中的約束矩陣和約束矢量帶入式(6)中，得到基于SVD的子陣級自適應(yīng)權(quán)矢量：

(11)

基于SVD的子陣級自適應(yīng)權(quán)矢量，將加權(quán)矢量置于噪聲噪聲子空間內(nèi)，利用噪聲子空間和干擾子空間的正交性抑制干擾，解決了小快拍下噪聲子空間小特征值發(fā)散的問題，因此可以一定程度上抑制干擾，主瓣保形。但由于子陣的不均勻劃分，子陣級的接收數(shù)據(jù)并不一致，旁瓣電平仍然較高。

3 基于SVD和罰函數(shù)約束的自適應(yīng)方向圖 控制算法

在子陣級，考慮加權(quán)靜態(tài)方向圖的向量wsub_q，對子陣級自適應(yīng)進行加權(quán)向量為wsub，則子陣級罰函數(shù)約束定義如下：

(12)

(13)

(14)

式中，asub(θ)=THa(θ)，為子陣級的導(dǎo)向矢量，h(θ)是一個合適的非負加權(quán)函數(shù)，在實際應(yīng)用中常被定義為θ的偶函數(shù)。可以看出罰函數(shù)E是個能量函數(shù)，反映自適應(yīng)方向圖與期望的靜態(tài)方向圖之間的差異。

利用線性約束最小方差準則，對子陣級自適應(yīng)權(quán)矢量進行約束的情況下，使得罰函數(shù)E能量最小，自適應(yīng)權(quán)矢量的最優(yōu)化問題為：

(15)

可以求得自適應(yīng)權(quán)矢量為：

(16)

將式(10)的約束矩陣和約束響應(yīng)矢量，代入式(16)可得SVD-PFC算法的自適應(yīng)加權(quán)矢量。

(17)

與式(11)的基于SVD算法的加權(quán)矢量的波束形成相比，采用式(17)的SVD-PFC算法的子陣級自適應(yīng)波束形成，通過罰函數(shù)進一步約束自適應(yīng)方向圖，在抑制干擾的同時，能夠使得自適應(yīng)方向圖逼近期望的靜態(tài)方向圖，進一步地降低旁瓣電平。

4 仿真結(jié)果分析

本文對子陣級常規(guī)算法、SVD算法和SVD-PFC算法的陣列方向圖和輸出信干噪比進行仿真，對比分析幾種算法的性能。仿真實驗假設(shè)一均勻直線陣，陣元間距為半波長，陣元個數(shù)為48，子陣數(shù)目為10，陣元不重疊，子陣陣元個數(shù)依次為[2 2 2 10 3 5 5 3 10 2 2 2]。

實驗1：假設(shè)波束指向為0o，信噪比為10 dB，干擾入射方向為-22°、12°，信干比均為-30 dB，快拍數(shù)目為25。常規(guī)算法、SVD和SVD-PFC算法的仿真結(jié)果如圖1所示。

(a) 常規(guī)算法和SVD算法

(b) SVD算法和SVD-PFC算法

由圖1可知，在小快拍情況下，常規(guī)算法雖然在干擾方向上形成零陷，但主瓣明顯發(fā)生惡化，出現(xiàn)指向偏移，旁瓣電平升高，而只采用SVD的算法可以很好地保持主瓣方向圖。SVD算法的旁瓣電平保持-8 dB不變，而常規(guī)算法的旁瓣電平為-6 dB，2種算法的零陷深度都達到-70 dB、-75 dB左右。相比常規(guī)算法，SVD算法在快拍變化時，能夠保持主瓣不發(fā)生惡化，但旁瓣電平略高。SVD-PFC算法的旁瓣電平約-13.3 dB，比SVD算法低5 dB左右，兩者的零陷深度也幾乎一樣。

實驗2：期望信號入射方向為0°，信噪比為10 dB，干擾信號入射方向為-15°、20°，信干比均為-30 dB，輸出SINR隨快拍數(shù)變化如圖2所示。

由圖2可知，3種算法的SINR都隨著快拍數(shù)目增大而增大，SVD-PFC算法和SVD算法在快拍數(shù)20左右，SINR就趨于平緩。而常規(guī)算法在快拍數(shù)60以后，才逐漸升高到平緩的范圍。其中在快拍數(shù)目小于40的情況下，SVD算法、SVD-PFC的算法的輸出SINR要比改進前的要好。但是在快拍數(shù)目接近40左右，常規(guī)算法的SINR開始高于SVD算法，而改進后的SVD-PFC算法的SINR始終都高于常規(guī)算法和SVD算法，驗證了結(jié)合SVD和罰函數(shù)的SVD-PFC算法可以有效地提高輸出SINR。

圖2 SINR隨快拍數(shù)目變化

實驗3：設(shè)期望信號入射方向為-40°～40°，信噪比為10 dB。干擾信號入射方向為-15°、20°，信干比均為-30 dB。輸出的SINR如圖3所示。

圖3 輸出SINR隨角度變化

從圖3中可以看到，在同一波束指向處，SVD算法的SINR略低于常規(guī)算法，但是通過罰函數(shù)改進SVD算法得到的SVD-PFC算法則始終略高于常規(guī)算法。

實驗4：期望信號入射方向為0°，干擾入射方向為-15°、20°，信噪比為10 dB，在不同快拍情況下，最大旁瓣電平隨信干比SIR的變化如圖4所示。

從圖4可以看出，在不同快拍情況下，常規(guī)算法的最大旁瓣電平隨著信干比的變化上下起伏波動很大，而SVD-PFC算法的最大旁瓣電平值保持在-13.2 dB左右，表明了SVD-PFC算法克服快拍數(shù)對自適應(yīng)波束形成的影響，具有保持靜態(tài)方向圖副瓣的特性。

圖4 最大旁瓣電平隨信干比SIR變化

5 結(jié)束語

針對小快拍數(shù)情況下，子陣級自適應(yīng)波束形成方向圖的主瓣變形，旁瓣電平升高問題，本文提出的SVD-PFC算法，通過SVD對子陣級接收數(shù)據(jù)分解，得到干擾子空間和噪聲子空間，修正子陣級的約束矩陣和約束響應(yīng)矢量，結(jié)合罰函數(shù)約束的子陣級自適應(yīng)方向圖控制算法。可在小快拍數(shù)時，主瓣波束保形，抑制干擾并保持靜態(tài)方向圖的副瓣特性，同時獲得比較好的輸出SINR。

[1] Nickel U R O. Subarray Configurations for Interference Suppression with Phased Array Radar[C]∥ Communications and Radar. Paris,France:Proc. of Intel. Conf. on radar,1989:82-86.

[2] Ferrier J M. Antenna Subarray Architectures And Anti-jamming Constraints[C]∥ Proc. of Intel. Conf. on Radar,Atlanta,USA,1994:466-469.

[3] Nickel U. Monopulse Estimation with Subarray-adaptive Arrays and Arbitrary Sum and Difference Beams[J]. IEE Proceedings-Radar Sonar and Navigation,1996,143(4):232-238.

[4] Combaud M. Adaptive Processing at the Subarray Level[J]. Aerospace Science & Technology,1999,3(3):93-105.

[5] Lombardo P,Cardinali R,Pastina D,et al. Array Optimization and Adaptive Processing for Sub-array Based Thinned Arrays[C]∥International Conference on Radar. IEEE,2008:197-202.

[6] 李洪濤. 自適應(yīng)數(shù)字波束形成關(guān)鍵技術(shù)研究[D].南京:南京理工大學(xué),2012.

[7] 張宗傲. 陣列雷達自適應(yīng)主瓣干擾抑制算法研究[D].北京:北京理工大學(xué),2015.

[8] Hughes D T,McWhirter J E. Penalty Function Method for Sidelobe Control in Least-squares Adaptive Beam Forming[C]∥SPIE's 1995 International Symposium on Optical Science,Engineering,and Instrumentation. International Society for Optics and Photonics,1995:170-181.

[9] Herbert G M. A New Projection Based Algorithm for Low Sidelobe Pattern Synthesis in Adaptive Arrays[C]∥ Radar. IEEE Xplore,1997:396-400.

[10]李洪濤,賀亞鵬,朱曉華,等. 基于特征向量的線性約束最小方差自適應(yīng)方向圖控制[J]. 南京理工大學(xué)學(xué)報(自然科學(xué)版),2011,35(4):529-533.

[11]蘇保偉，王永良，周良柱. 基于 LCMV 線性約束的自適應(yīng)方向圖控制[J]. 電子與信息學(xué)報，2008，30(2):282-285.

[12]李榮峰，趙士鋒，丁前軍. 一種基于線性約束的自適應(yīng)方向圖控制方法[J]. 雷達科學(xué)與技術(shù)，2009,7(2)：139-142.

[13]Nickel U R O. Subarray Configurations for Digital Beamforming with Low Sidelobes and Adaptive Interference Suppression[C]∥ Radar Conference,1995. Record of the IEEE 1995 International,1995:714-719.

[14]許志勇,保錚,廖桂生. 一種非均勻鄰接子陣結(jié)構(gòu)及其部分自適應(yīng)處理性能分析[J]. 電子學(xué)報,1997(9):20-24.

[15]陳子歡,劉剛,蔣寧. 一種新的子陣結(jié)構(gòu)及其自適應(yīng)性能分析[J]. 電子信息對抗技術(shù),2006,21(3):33-37.

[16]熊子源,徐振海,張亮,等. 陣列雷達最優(yōu)子陣劃分研究[J]. 雷達科學(xué)與技術(shù),2011,09(4):370-377.

[17]Lombardo P,Pastina D. Pattern Control for Adaptive Antenna Processing with Overlapped Sub-arrays[C]∥ Radar Conference,2003. Proceedings of the International. IEEE,2003:188-193.

[18]楊志偉,賀順,廖桂生,等. 子空間重構(gòu)的一類自適應(yīng)波束形成算法[J]. 電子與信息學(xué)報,2012,34(5):1115-1119.

[19]Hu X, Guo L,Li S,et al. Improved Orthogonal Projection Adaptive Beamforming Based on Normalization at Subarray Level[C]∥ Iet International Radar Conference,2015:14-16.

Subarray Level Pattern Control Based on Singular Value Decomposition and Penalty Function Constraint

LUAN Xiao-ming，XIAO Hua-fei

(School of Information and Communication Engineering,Harbin Engineering University,Harbin Heilongjiang 150001,China)

To solve the problem of high sidelobe in adaptive beamforming at subarray level,an adaptive pattern control algorithm based on singular value decomposition and penalty function constraints (SVD-PFC) at subarray level is proposed. The algorithm obtains the interference subspace and noise subspace by the singular value decomposition of the input data,the constraint matrix and the constraint response vector are modified by the interference subspace,and the interference can be effectively suppressed in the condition of small snapshots. At the same time,it combines the penalty function to constrain the adaptive pattern to approximate the expected static pattern. The proposed algorithm has good sidelobe characteristics and improves the signal to interference plus noise ratio (SINR) of the output in the condition of small snapshots. Simulation results confirm the feasibility and effectiveness of the proposed algorithm.

adaptive pattern control;singular value decomposition;penalty function;subarray level

10. 3969/j.issn. 1003-3114. 2017.04.13

欒曉明,肖華飛. 基于SVD和罰函數(shù)約束的子陣級方向圖控制[J].無線電通信技術(shù),2017,43(4):56-59.

[LUAN Xiaoming,XIAO Huafei. Subarray Level Pattern Control Based on Singular Value Decomposition and Penalty Function Constraint [J].Radio Communications Technology,2017,43(4):56-59.]

2017-02-16

欒曉明(1963—)，男，教授，碩士生導(dǎo)師，主要研究方向：寬帶信號檢測、處理與識別。肖華飛(1991— )，男，碩士研究生，主要研究方向：陣列天線、自適應(yīng)信號處理。

TN820.1

A

1003-3114(2017)04-56-4