信雜噪比約束下基于互信息的雷達(dá)波形設(shè)計(jì)

肖 宇, 鄧正宏

(1. 西北工業(yè)大學(xué)自動(dòng)化學(xué)院, 陜西 西安 710072; 2. 空軍工程大學(xué), 陜西 西安 710051)

0 引 言

雷達(dá)通過(guò)發(fā)射信號(hào)波形來(lái)探測(cè)目標(biāo)和環(huán)境,隨著目標(biāo)的運(yùn)動(dòng),探測(cè)環(huán)境在發(fā)生變化,根據(jù)目標(biāo)和環(huán)境之間的相對(duì)關(guān)系設(shè)計(jì)發(fā)射波形是認(rèn)知雷達(dá)研究[1]的關(guān)鍵問(wèn)題。因此,發(fā)射信號(hào)波形設(shè)計(jì)成為影響雷達(dá)探測(cè)性能的關(guān)鍵因素之一。

波形設(shè)計(jì)的本質(zhì)是發(fā)射信號(hào)能量的分配,現(xiàn)有的波形設(shè)計(jì)方法主要以提高雷達(dá)目標(biāo)參數(shù)估計(jì)和目標(biāo)檢測(cè)性能等[2]為目的,通過(guò)優(yōu)化準(zhǔn)則建立波形優(yōu)化目標(biāo)函數(shù),結(jié)合波形約束條件,提出優(yōu)化算法等[3-4]措施實(shí)現(xiàn)波形設(shè)計(jì)。常用的優(yōu)化準(zhǔn)則包括最小均方誤差[5](minimum mean squared error,MMSE)、互信息(mutual information, MI)[6]、信雜噪比(signal to interference noise ratio, SINR)[7]等。在頻域波形設(shè)計(jì)方面,Bell[8]較早地將信息論引入雷達(dá)波形設(shè)計(jì),基于最大信噪比(signal to noise ratio,SNR)設(shè)計(jì)波形提升目標(biāo)檢測(cè)性能,基于最大MI設(shè)計(jì)波形提升目標(biāo)信息提取能力,文獻(xiàn)[9]針對(duì)Bell[8]提出的water-filing波形中的拉格朗日乘子區(qū)間不同情況進(jìn)行了討論,本質(zhì)上是能量約束的問(wèn)題。Romero[10]在Bell[8]的基礎(chǔ)上,較系統(tǒng)地給出了基于SNR和MI針對(duì)確定目標(biāo)、隨機(jī)目標(biāo)在雜波和噪聲情況下的波形設(shè)計(jì)方法。文獻(xiàn)[11]針對(duì)目標(biāo)先驗(yàn)信息不確定情況下,提出了基于最大SNR與最大MI的魯棒波形設(shè)計(jì)方法,主要特點(diǎn)是針對(duì)目標(biāo)檢測(cè)和參數(shù)估計(jì)在能量有限約束情況下的波形設(shè)計(jì)。隨著電子對(duì)抗技術(shù)的發(fā)展,波形設(shè)計(jì)需要考慮到雷達(dá)干擾技術(shù)的問(wèn)題,文獻(xiàn)[12-13]以Bell[8]提出的波形設(shè)計(jì)目標(biāo)函數(shù)模型為基礎(chǔ),將信號(hào)相關(guān)干擾作為優(yōu)化目標(biāo)函數(shù),提出了相應(yīng)的干擾波形優(yōu)化方法,并給出了魯棒的干擾波形。文獻(xiàn)[14]從博弈論的角度給出了干擾波形的優(yōu)化方法,上述研究都是在某種優(yōu)化準(zhǔn)則下,給出了定能量約束的優(yōu)化波形設(shè)計(jì)方法,在SINR較低或者SINR約束情況下,并未給出具體的波形優(yōu)化方法。而文獻(xiàn)[15]從認(rèn)知雷達(dá)波形設(shè)計(jì)方面,給出了SNR約束條件下的波形設(shè)計(jì)方法,討論了SNR與MI之間的均衡關(guān)系,但并未明確SNR約束條件的來(lái)由,也未從理論上給出SNR的作用閾,同時(shí)并未考慮存在信號(hào)相關(guān)雜波情況下的波形設(shè)計(jì)。

本文擬從二元假設(shè)檢驗(yàn)?zāi)Ｐ偷南鄬?duì)熵入手,結(jié)合相對(duì)熵的非負(fù)性條件,挖掘MI與SINR的約束關(guān)系,以及SINR對(duì)最大化MI的影響區(qū)間,并基于最大化MI準(zhǔn)則,結(jié)合SINR的約束條件,提出一種新的波形設(shè)計(jì)方法,以期為認(rèn)知雷達(dá)的波形設(shè)計(jì)提供理論支撐。

1 二元假設(shè)檢驗(yàn)?zāi)Ｐ?/h2>

當(dāng)存在信號(hào)相關(guān)雜波干擾c(t)的情況下,雷達(dá)發(fā)射波形x(t)來(lái)探測(cè)具有脈沖響應(yīng)h(t)的目標(biāo),接收信號(hào)y(t)[10]可以表示為

y(t)=h(t)*x(t)+c(t)*x(t)+n(t)

(1)

對(duì)于目標(biāo)檢測(cè)問(wèn)題,在二元假設(shè)檢驗(yàn)?zāi)Ｐ椭?接收信號(hào)可以產(chǎn)生兩種不同的輸出結(jié)果,當(dāng)接收信號(hào)中僅有雜波和噪聲信號(hào)時(shí),表示為假設(shè)H0,當(dāng)接收信號(hào)中同時(shí)存在目標(biāo)、雜波和噪聲信號(hào)時(shí),表示為假設(shè)H1,即

(2)

(3)

2 MI與SINR的約束關(guān)系

信息論中通常使用相對(duì)熵(Kullback-Leibler divergence,KLD)來(lái)描述兩個(gè)概率密度函數(shù)(probability density function,PDF)間的“距離”度量[16]。令假設(shè)H0和假設(shè)H1的相對(duì)熵表示為D(p0(Y)‖p1(Y)),其中p0(Y)、p1(Y)分別表示輸出信號(hào)Y(f)在假設(shè)H0和H1下的PDF,為了簡(jiǎn)化算式表達(dá),將D(p0(Y)‖p1(Y))簡(jiǎn)寫(xiě)為D(p0‖p1)。同理,假設(shè)H1和假設(shè)H0的相對(duì)熵可以表示為D(p1‖p0)。

當(dāng)雷達(dá)接收端對(duì)接收信號(hào)的采樣數(shù)據(jù)量N趨于無(wú)窮大時(shí),假設(shè)H0、假設(shè)H1之間的相對(duì)熵[17]可以表示為

(4)

(5)

將式(3)模型代入相對(duì)熵式(4)和式(5),并進(jìn)行變換,可得

(6)

(7)

由于任意兩個(gè)PDF的KLD滿(mǎn)足條件KLD≥0,因而假設(shè)H0和H1之間的KLDD(p0‖p1)≥0、D(p1‖p0)≥0恒成立,聯(lián)合變換式(6)和式(7),可得

(8)

從式(8)中可以看出目標(biāo)響應(yīng)與接收信號(hào)的互信息MI與SINR之間的約束關(guān)系,若令

(9)

(10)

式(8)可表示為

XINR≤MI≤SINR

(11)

其中,

(12)

(13)

式(12)所體現(xiàn)的SINR(f)與MI之間的關(guān)系,與文獻(xiàn)[10]相一致,式(13)顯示了XINR隨著SINR(f)的變化關(guān)系。從式(11)可知,在給定波形|X(f)|2的情況下,MI總被限定于某個(gè)區(qū)間,并在每個(gè)頻率時(shí)刻受到SINR(f)約束的影響。如果要從最大化MI的角度,對(duì)信號(hào)波形進(jìn)行設(shè)計(jì),則必須要綜合考慮以SINR(f)作為約束條件。在許多以MI為準(zhǔn)則設(shè)計(jì)雷達(dá)波形的文獻(xiàn)中,缺少了SINR(f)的約束,即使在文獻(xiàn)[15]中有SNR的約束,但是并未對(duì)SNR的作用閾進(jìn)行限定討論。

從式(10)可見(jiàn),SINR(f)≥0,若以SINR(f)作為坐標(biāo)橫軸,對(duì)MI、XINR的核函數(shù)繪制示意圖,如圖1所示。

圖1 MI、SINR、XINR核函數(shù)關(guān)系示意圖Fig.1 Relationship among MI, SINR and XINR kernel functions

3 雜波環(huán)境下基于MI的波形設(shè)計(jì)

基于MI準(zhǔn)則的波形優(yōu)化方法,主要通過(guò)優(yōu)化發(fā)射波形使接收信號(hào)中目標(biāo)信息量最大化。在雜波環(huán)境下,給定發(fā)射信號(hào),最大化接收信號(hào)與目標(biāo)脈沖響應(yīng)之間的MI[10]可表示為

MI(|X(f)|2)=

(14)

式中:Ty表示接收信號(hào)的持續(xù)時(shí)間。

3.1 波形設(shè)計(jì)的目標(biāo)函數(shù)

在文獻(xiàn)[8-10]中,基于MI最大化波形設(shè)計(jì),主要受到能量的限制。因此,根據(jù)第2節(jié)的討論結(jié)果,加上SINR作為約束條件的發(fā)射波形優(yōu)化目標(biāo)函數(shù)模型可表示為

SINR≤SINR0

(15)

結(jié)合第2節(jié)的分析,對(duì)于式(15),在SINR約束情況下,基于最大化MI獲得的最優(yōu)波形將不同于無(wú)約束情況的最優(yōu)波形;當(dāng)SINR0較小時(shí),相較于無(wú)SINR約束情況,此時(shí)目標(biāo)函數(shù)MI將無(wú)法達(dá)到最大值。

根據(jù)SINR的定義,可知

(16)

其中

(17)

同理,可得

(18)

(19)

式中:h=[h[1],h[2], …,h[Lh]]T,C=[c[1],c[2], …,c[Lh]]T,則

(20)

(21)

因此,對(duì)于SINR較低的情況,也就是SINR≤SINR0時(shí),可得

(22)

根據(jù)維納-辛欽定理,發(fā)射信號(hào)x(t)的功率譜密度|X(f)|2可表示為自相關(guān)函數(shù)Rxx(t)的傅里葉變換,經(jīng)離散采樣后,可得

(23)

其中,

Jf=[1,2cos(2πf),2cos(4πf),…,cos(2π(Lx-1)f)]T

(24)

令Sxx=[|X(f1)|2, |X(f2)|2, …, |X(fN)|2]T,J=[Jf1,Jf2, …,JfN]T,那么Sxx=JRxx。此時(shí)通過(guò)矩陣廣義逆方法,可得Rxx=J?Sxx。因此,SINR的約束條件即可表示為

(25)

根據(jù)前面的推導(dǎo),當(dāng)存在SINR約束條件時(shí),基于最大化MI的發(fā)射波形優(yōu)化目標(biāo)函數(shù)模型可表示為

Sxx=[|X(f1)|2, |X(f2)|2, …, |X(fN)|2]T

(26)

從式(26)來(lái)看,該目標(biāo)函數(shù)為非凸函數(shù),因此可采取拉格朗日乘子法對(duì)波形PSD進(jìn)行優(yōu)化求解。

3.2 SINR的作用區(qū)間范圍

從式(11)可知,MI的取值介于XINR與SINR之間,同時(shí)由于MI受到SINR的約束。因此,在SINR約束條件下,基于最大化MI優(yōu)化雷達(dá)發(fā)射波形,需要對(duì)SINR0的取值范圍進(jìn)行討論。也就是對(duì)于式(26)而言,需要考慮當(dāng)SINR0處于什么區(qū)間范圍時(shí),會(huì)對(duì)最大化MI的波形設(shè)計(jì)有影響。本節(jié)討論的出發(fā)點(diǎn)主要基于MI、XINR、SINR三者的相互關(guān)系,進(jìn)而為波形設(shè)計(jì)提供三者間的均衡關(guān)系。從式(12)和式(13)可知二者的導(dǎo)數(shù)為

MI(f)′=ln[1+SINR(f)]

(27)

(28)

對(duì)于頻率點(diǎn)k=1,2,…,N,存在SINR(fk)≥0,根據(jù)式(27)和式(28)的結(jié)果,可以推導(dǎo)出MI(f)′≥0、XINR(f)′≥0,可見(jiàn)MI(f)、SINR(f)、XINR(f)三者均為單調(diào)遞增函數(shù),由于MI(f)、XINR(f)與SINR(f)之間存在相互約束關(guān)系,因此三者間呈同步單調(diào)遞增趨勢(shì)。當(dāng)SINR0趨于無(wú)窮大時(shí),相當(dāng)于無(wú)SINR約束,式(26)的目標(biāo)函數(shù)并不能實(shí)現(xiàn)MI的無(wú)窮大,因?yàn)榇藭r(shí)存在有限能量的約束條件。對(duì)于給定頻域帶寬W的情況,無(wú)SINR約束時(shí),式(26)的目標(biāo)函數(shù)優(yōu)化波形轉(zhuǎn)變?yōu)閮H有能量約束的基于最大化MI的波形設(shè)計(jì),從而可得波形[10]為

|XMI(f)|2=max {0,B(f)[A-D(f)]}

(29)

式中:A為常數(shù);

(30)

(31)

根據(jù)前面的討論,可以得到SINR0最大值SINR(|XMI(f)|2),同理,從式(11)可以看出MI總是大于XINR,即使最小的MI也將大于最大的XINR。所以,在無(wú)約束情況下,采取拉格朗日乘子法獲得最大化XINR的優(yōu)化波形為

(32)

此時(shí)得到SINR0最小值SINR(|XXINR(f)|2),因此SINR的作用閾可表示為

SINR(|XXINR(f)|2)≤SINR0≤SINR(|XMI(f)|2)

(33)

在式(33)區(qū)間之外的SINR約束將不會(huì)對(duì)最大化MI的波形設(shè)計(jì)有任何影響。

4 仿真及結(jié)果分析

本節(jié)主要對(duì)基于最大化MI的波形設(shè)計(jì)進(jìn)行驗(yàn)證,分別對(duì)雷達(dá)探測(cè)環(huán)境中存在雜波干擾、不存在雜波干擾等兩種情況,對(duì)目標(biāo)函數(shù)中存在SNR/SINR約束和不存在SNR/SINR約束等兩種條件,進(jìn)行綜合仿真對(duì)比實(shí)驗(yàn)。

4.1 SNR/SINR作用閾與MI的約束關(guān)系

當(dāng)不存在雜波干擾時(shí),在給定目標(biāo)脈沖響應(yīng)情況下,此時(shí)式(26)轉(zhuǎn)化為

Sxx=[|X(f1)|2, |X(f2)|2, …, |X(fN)|2]T

(34)

從式(34)的目標(biāo)模型來(lái)看,可以采取凸優(yōu)化的方式進(jìn)行計(jì)算,根據(jù)凸優(yōu)化工具箱,可計(jì)算當(dāng)SNR0取任意值時(shí),最大化MI的取值情況,仿真結(jié)果如圖2所示。

圖2 無(wú)雜波時(shí)SNR與MI的約束關(guān)系Fig.2 Constraint relationship between SNR and MI without clutter

從圖2可以看出,當(dāng)SNR<-11 dB或SNR>12 dB時(shí),無(wú)論SNR再小或再大,都不會(huì)影響最大化MI的取值;當(dāng)SNR介于-11 dB到12 dB區(qū)間時(shí),隨著SNR的增大,MI在逐漸增大,二者呈單調(diào)遞增關(guān)系,這與第3.2節(jié)理論推導(dǎo)及式(33)是相符的;同時(shí)圖2曲線(xiàn)的上下邊界,也就是MI的最大值與最小值之間的關(guān)系,與式(11)相吻合。

同理,當(dāng)存在雜波干擾時(shí),需要針對(duì)SINR0取不同值的情況,對(duì)式(26)進(jìn)行優(yōu)化,從式(26)的目標(biāo)函數(shù)來(lái)看,難以使用凸優(yōu)化的方式進(jìn)行計(jì)算,本文采取拉格朗日乘子法對(duì)目標(biāo)函數(shù)進(jìn)行優(yōu)化,獲得SINR與MI的約束關(guān)系如圖3所示。

圖3 有雜波時(shí)SINR與MI的約束關(guān)系Fig.3 Constraint relationship between SNR and MI with clutter

圖3與圖2所示結(jié)果具有相似的情況,符合第3.2節(jié)的理論推導(dǎo),當(dāng)SINR介于-15 dB到8 dB區(qū)間時(shí),會(huì)對(duì)基于最大化MI設(shè)計(jì)的波形有影響,該區(qū)間不同于圖2所示區(qū)間的變化范圍,SINR的影響區(qū)間總體發(fā)生了左移,主要原因在于雜波引起了約束的變化。由于雜波的加入,使得最大化MI值的區(qū)間范圍有所擴(kuò)展,并且相較于圖2的MI值減小了近10倍,這個(gè)減小范圍將根據(jù)雜波PSD的不同而不同。

4.2 有SINR約束和無(wú)SINR約束的波形對(duì)比

考慮到對(duì)比分析,本節(jié)分別針對(duì)無(wú)雜波、有雜波兩種情況,對(duì)SINR約束下基于最大化MI的波形進(jìn)行仿真。從SNR/SINR的作用閾中選取1 dB SINR作為約束閾值,并使用Sxx-NoSNR/SINR表示無(wú)SNR/SINR約束時(shí)的優(yōu)化波形,使用Sxx-SNR/SINR表示有SNR/SINR約束時(shí)的優(yōu)化波形。在無(wú)雜波影響情況下,使用式(34),采用凸優(yōu)化方法,對(duì)波形進(jìn)行仿真,并進(jìn)行歸一化處理,如圖4所示。

圖4 無(wú)雜波時(shí)有、無(wú)SNR約束的波形設(shè)計(jì)對(duì)比Fig.4 Comparison of waveform designs with and without SNR constraints when there is no clutter

在能量有限情況下,圖4所示基于最大化MI的波形信號(hào)能量,分布于具有目標(biāo)脈沖響應(yīng)的所有頻率點(diǎn),在目標(biāo)PSD越強(qiáng)的地方,分配的波形信號(hào)能量越多,波形隨頻率的變化關(guān)系與目標(biāo)PSD隨頻率的變化關(guān)系相一致。當(dāng)在目標(biāo)函數(shù)中加入SNR約束進(jìn)行波形優(yōu)化時(shí),波形的總體變化趨勢(shì)不變,不同之處在于不同的頻率點(diǎn)分布能量會(huì)出現(xiàn)有限的波動(dòng),也就是在不同頻率點(diǎn),相較于無(wú)SNR約束情況,獲取的MI存在有限的波動(dòng),主要原因在于不存在雜波的情況下,SNR約束與能量約束存在等效約束關(guān)系。

在有雜波影響情況下,使用式(26),采用拉格朗日乘子法對(duì)波形進(jìn)行優(yōu)化仿真,并進(jìn)行歸一化處理,如圖5所示。

圖5 有雜波時(shí)有、無(wú)SINR約束的波形設(shè)計(jì)對(duì)比Fig.5 Comparison of waveform designs with and without SINRconstraints when there is clutter

由于受到雜波的影響,在無(wú)SINR約束時(shí),波形能量主要分布在目標(biāo)PSD強(qiáng)、雜波PSD弱的地方,而在雜波PSD強(qiáng)、目標(biāo)PSD弱的地方分配較少的能量,這與文獻(xiàn)[8,10,15]所示結(jié)論是一致的,因?yàn)椴ㄐ蝺?yōu)化方法都是基于water-filing方法。當(dāng)存在SINR約束時(shí),波形能量分配發(fā)生了變化,傾向于在有限能量下,更好地利用能量來(lái)提取目標(biāo)信號(hào),即在目標(biāo)PSD強(qiáng)、雜波PSD弱的頻率點(diǎn),波形分配能量有所下降,因?yàn)榇藭r(shí)分配的能量足夠提取目標(biāo)信號(hào);而在存在目標(biāo)、且雜波PSD弱的頻率點(diǎn),波形分配能量有所增加,因?yàn)樵诖瞬糠诸l率點(diǎn)分配能量,能夠更多地提取目標(biāo)信號(hào)。總之,基于最大化MI設(shè)計(jì)的波形需要在不同頻率點(diǎn)獲得目標(biāo)的相關(guān)信息,當(dāng)存在SINR約束時(shí),波形能量分配更傾向于有目標(biāo)響應(yīng)且無(wú)雜波或雜波信號(hào)弱的頻率點(diǎn),以便增大MI,提取目標(biāo)信號(hào)。

5 結(jié) 論

在雜波環(huán)境下,為提升目標(biāo)估計(jì)性能,本文綜合考慮SINR準(zhǔn)則與MI準(zhǔn)則間的均衡關(guān)系,提出了一種SINR約束下基于最大化MI的波形設(shè)計(jì)方法。同時(shí),基于KLD的約束條件,從理論上推導(dǎo)了SINR的作用閾,進(jìn)一步限定了基于最大化MI的波形設(shè)計(jì)。仿真結(jié)果表明,在SINR的約束下,能夠充分發(fā)揮最大化MI的特性增大目標(biāo)信息提取量,相較于無(wú)SINR約束的情況,波形信號(hào)能量會(huì)在存在目標(biāo)信號(hào)、且無(wú)雜波信號(hào)或雜波信號(hào)弱的頻率點(diǎn)有所增加。根據(jù)MI與SINR之間的約束均衡關(guān)系,有助于深入開(kāi)展認(rèn)知雷達(dá)波形設(shè)計(jì)方面的研究工作。