基于G-S空間平滑處理含期望信號的相干信號波束形成算法*

孫維忠 李成全 王 波

(武漢船舶通信研究所 武漢 430205)

?

基于G-S空間平滑處理含期望信號的相干信號波束形成算法*

孫維忠 李成全 王 波

(武漢船舶通信研究所 武漢 430205)

Gram-Schmidt正交投影算法在自適應(yīng)波束形成中有廣泛應(yīng)用,但是對于連續(xù)波體制雷達(dá),接收信號中存在期望信號,算法會(huì)將期望信號抑制掉,而且作為信號子空間算法,在接收信號相干時(shí)將失效。論文針對陣列接收信號中含期望信號,并且信號相干情況下,提出通過空間平滑處理,對接收信號消除期望信號并解相干,有效抑制干擾,并保證期望信號不被抑制。

自適應(yīng)波束形成; G-S算法; 空間平滑

Class Number TP301.6

1 引言

在自適應(yīng)數(shù)字波束形成中,Gram-Schmidt正交投影算法是一種重要的開環(huán)算法,由于其收斂速度快,算法簡單易行,運(yùn)算量小,因而取得廣泛的應(yīng)用[1～4]。然而,這種算法要求接收數(shù)據(jù)中不含期望信號,對于脈沖體制雷達(dá)而言,通過在雷達(dá)休止期內(nèi)采樣,就可以保證采集到的快拍向量數(shù)據(jù)內(nèi)不含期望信號。而對于連續(xù)波體制雷達(dá)而言,由于不存在休止期,因而采集到的快拍向量數(shù)據(jù)內(nèi)總是包含有期望信號分量,傳統(tǒng)的Gram-Schmidt正交投影算法已不再適用。通過空間平滑對接收信號進(jìn)行預(yù)處理,消除其中期望信號,然后再采用傳統(tǒng)Gram-Schmidt正交化算法。另一方面,當(dāng)信號子空間存在相干或強(qiáng)相關(guān)信號時(shí),傳統(tǒng)的Gram-Schmidt正交投影算法就會(huì)失效[5～7]。采用空間平滑技術(shù),對接收信號進(jìn)行解相干,然后再采用傳統(tǒng)Gram-Schmidt正交化算法[8～9]。然而,在更一般的情況下,接收信號中不僅含有期望信號,而且信號相干或強(qiáng)相關(guān),前述方法都存在無法同時(shí)解決兩方面的問題。本文在前述方法基礎(chǔ)上,結(jié)合兩者優(yōu)點(diǎn),提出基于空間平滑的Gram-Schmidt正交化方法,同時(shí)解決接收信號中含期望信號和相干信號的問題。

2 G-S正交化算法基本原理

2.1 G-S正交化算法基本原理[10]

假設(shè)一均勻線陣,天線陣列由M個(gè)全向陣元組成。若一個(gè)期望信號和P個(gè)窄帶遠(yuǎn)場平面波干擾信號入射到陣上,各陣元噪聲是相互獨(dú)立、功率相等的零均值高斯白噪聲,則陣列接收信號矢量表示為

X(t)=a(θ0)·s0(t)+AJ(t)+N(t)

(1)

其中,X(t)=[x1(t),x2(t),…,xM(t)]T,xm(t)為第m(m=1,2,…,M)個(gè)陣元的接收信號,a(θ0)、s0(t)分別為期望信號的導(dǎo)向矢量和復(fù)包絡(luò),A=[a(θ1),a(θ2),…,a(θP)],a(θP)為第p(p=1,2,…,P)個(gè)干擾信號的導(dǎo)向矢量,J(t)=[J1(t),J2(t),…,JP(t)]T,Jp(t)為第p(p=1,2,…,P)個(gè)干擾信號的復(fù)包絡(luò),N(t)為背景噪聲。L次快拍數(shù)據(jù)得到的自適應(yīng)權(quán)值為w,則第n次快拍自適應(yīng)網(wǎng)絡(luò)輸出表示為

Y(n)=WHX(n)=WHXS(n)+WHXJ(n)

(2)

其中,XS(n)、XJ(n)分別為期望信號和干擾信號(包含噪聲),可見式(2)中前一項(xiàng)為期望輸出,而后一項(xiàng)是干擾(包含噪聲)的輸出,因此需要置零。XJ(n)是P個(gè)干擾的合成,即XJ(n)=X1(n)+X2(n)+…+XP(n),則應(yīng)使

WHXp(n)=0,p=1,2,…,P

(3)

要對干擾進(jìn)行抑制,則需使加權(quán)矢量與每一個(gè)干擾矢量正交,即加權(quán)矢量要與P個(gè)干擾信號組成的干擾子空間正交。在脈沖體制雷達(dá)休止期采集信號(不含期望信號),應(yīng)用Gram-Schmidt正交化方法得到最優(yōu)權(quán)。設(shè)

(4)

2.2 基于G-S空間平滑處理去期望信號

簡化模型,考慮無背景噪聲的情況。將線陣劃分為兩個(gè)子陣,分別包含前M-1個(gè)連續(xù)陣元和后M-1個(gè)連續(xù)陣元。設(shè)n時(shí)刻陣列接收信號矢量

(5a)

(5b)

為消除快拍數(shù)據(jù)中期望信號,對上式作如下處理:

=BX

(6)

可見,X′中只有干擾信號,不含期望信號,再按前述傳統(tǒng)Gram-Schmidt正交化方法,就可以得到在干擾方向上形成零點(diǎn)的自適應(yīng)權(quán)值。

2.3 基于G-S空間平滑處理去相干

(7)

第h(h=1,2,…,H)個(gè)前向子陣接收信號矢量

=AHDi-1S(n)

(8)

第h(h=1,2,…,H)個(gè)后向子陣接收信號矢量的共軛:

(10)

2.4 本文改進(jìn)方法

在更一般的情況下,接收信號中不僅含有期望信號,而且信號相干或強(qiáng)相關(guān),前述方法都將失效,均無法同時(shí)解決兩方面的問題。通過結(jié)合前述兩種方法的優(yōu)勢,本文提出改進(jìn)方法,能有效對接收信號進(jìn)行解相干、抑制干擾信號并保證期望信號不被抑制。具體步驟如下:

1) 由n時(shí)刻陣列接收信號矢量X(n)=[x1(n),x2(n),…,xM(n)]T,作空間平滑處理,構(gòu)造Xs(n),對接收信號進(jìn)行解相干;

2) 對Xs(n)作空間平滑處理,X′(n)=BXs(n),消除Xs(n)中期望信號分量;

3 仿真結(jié)果與分析

假設(shè)陣列由16個(gè)間距d=λ/2的陣元組成的均勻線陣。期望信號方向15°,信噪比30dB,三個(gè)與期望信號互為相干的干擾信號方向分別為-45°、-10°和40°,信干比分別為10dB、20dB和15dB。快拍數(shù)為3,進(jìn)行100次獨(dú)立試驗(yàn)仿真,并且分別與采用僅期望信號消除和僅干擾信號去相干方法仿真結(jié)果進(jìn)行了對比,如圖1、圖2和圖3所示。

圖1 僅期望信號消除法

圖2 僅干擾信號去相干法

圖3 采用本文改進(jìn)方法

可見,在接收信號存在期望信號,并且期望信號與干擾信號相干情況下,僅期望信號消除法能接收期望信號,可是方向圖的零陷已經(jīng)偏離了干擾方向,對于干擾的抑制效果較差,在-10°方向零深僅約為-24dB,而僅干擾信號去相干法已經(jīng)完全失效,在期望信號方向形成約-40dB零陷,無法接收期望信號。本文所采用的方法不僅能保證期望信號不被抑制,同時(shí)在干擾方向形成不高于-38dB的零陷,很好地抑制了干擾。

4 結(jié)語

對于陣列接收信號中含有信號或存在相干信號,傳統(tǒng)的Gram-Schmidt正交投影算法不再適用。本文提出采用空間平滑處理陣列接收信號,消除其中期望信號分量,并對信號分量解相干,再運(yùn)用Gram-Schmidt正交投影算法進(jìn)行波束形成。通過對比仿真分析,在接收信號中含期望信號并且存在信號相干時(shí),相比于將空間平滑處理應(yīng)用于消除期望信號或信號解相干,本文提出的將兩者結(jié)合的方法,能有效抑制干擾,并保證期望信號不被抑制。

[1] Hung E K L, Turner R M. A fast beam-forming algorithm for large arrays[J]. IEEE Trans on AES,1983,19(4):598-607.

[2] GS:陳建春,趙樹杰.基于G-S算法的脈沖追趕式自適應(yīng)DBF實(shí)現(xiàn)技術(shù)[J].西安電子科技大學(xué)學(xué)報(bào),1999,26(5):541-544.

[3] F. Ling. A Recursive Modified Gram Schmidt Algorithm for Least Squares Estimation[J]. IEEE Trans on ASSP,1986,34(8):829-835.

[4] 王俊鳴,齊會(huì)來,張子華,等.基于GS算法的圓柱形共形微帶陣列自適應(yīng)零點(diǎn)形成[J].火力與指揮控制,2011,36(3):118-120.

[5] 王紅哲.數(shù)字波束形成的G-S自適應(yīng)零點(diǎn)形成算法研究[D].鄭州:河南大學(xué),2003:28-49.

[6] 郝亮飛,李勇,李滔,等.改進(jìn)G-S零點(diǎn)控制算法在雷達(dá)中的應(yīng)用[J].探測與控制學(xué)報(bào),2010,32(1):72-76.

[7] 李寧濤,陶海紅,廖桂生.樣本含有期望信號時(shí)的正交投影波束形成算法[J].雷達(dá)科學(xué)與技術(shù),2007,5(4):246-264.

[8] 余繼周.用基于空間平滑的正交化算法分離相干信號[J].電子工程師,2005,31(1):1-4.

[9] Van Trees H L. Optimum array processing, part Ⅳ of detection, estimation and modulation theory[M]. New York: Wiley Press,2002.

[10] 龔耀寰.自適應(yīng)濾波[M].第二版.北京:電子工業(yè)出版社,2003:88-90.

Beam Formation of Coherent Signals with Target Signal Based on G-S Space-Smoothing

SUN Weizhong LI Chengquan WANG Bo

(Wuhan Maritime Communication Research Institute, Wuhan 430205)

Gram-Schmidt orthography projection algorithm is applied in adaptive digital beam forming, but for continuous wave radar, the target signal included in the received signals will be restrained, so in this situation and the traditional algorithm as the signal subspace algorithm will also be disabled if the received signals are coherent with each other. Under the situation that the array signals contain the desired signal, and the signals are coherent, through space smoothing process, the target signal is eliminated from the received signals, and the jamming signal decorrelation is received, so as to effectively suppress interference, and to ensure that the desired signal is not suppressed.

adaptive digital beam forming, G-S algorithm, space-smoothing

2015年2月12日,

2015年3月22日

孫維忠,男,助理工程師,研究方向:相控陣天線設(shè)計(jì)。李成全,男,助理工程師,研究方向:艦載通信天線,抗干擾陣列天線。王波,男,助理工程師,研究方向:相控陣天線設(shè)計(jì)。

TP301.6

10.3969/j.issn1672-9730.2015.08.018