陳文東,湯 斌

(1.中國航天科工集團(tuán)第二研究院,北京100854;2.北京遙感設(shè)備研究所,北京100854)

復(fù)雜電磁環(huán)境下雷達(dá)抗干擾性能評(píng)估方法

陳文東1,2,湯 斌2

(1.中國航天科工集團(tuán)第二研究院,北京100854;2.北京遙感設(shè)備研究所,北京100854)

根據(jù)現(xiàn)有的雷達(dá)抗干擾性能評(píng)估方法存在的問題,提出了一種基于高維量測特征空間幾何分析的雷達(dá)抗干擾性能評(píng)估方法。該方法將抗干擾試驗(yàn)中雷達(dá)得到的量測數(shù)據(jù)映射到高維量測特征空間中,利用量測數(shù)據(jù)的多元密度分布的幾何特征定義量測誤差影響因子(MEIF)、量測穩(wěn)定性影響因子(MSIF)和干擾影響因子(JIF),定量分析干擾對(duì)雷達(dá)的影響以及雷達(dá)的綜合抗干擾性能。計(jì)算機(jī)仿真表明,新的評(píng)估指標(biāo)可以準(zhǔn)確反映雷達(dá)的抗干擾性能,實(shí)驗(yàn)結(jié)果與理論分析相符。

復(fù)雜電磁環(huán)境;雷達(dá);抗干擾;性能評(píng)估

0 引言

現(xiàn)代戰(zhàn)場上,雷達(dá)通常需要在復(fù)雜電磁環(huán)境中工作,在面對(duì)各式各樣干擾的情況下,同樣需要具備精確可靠探測的能力,即抗干擾能力[1]。如何對(duì)雷達(dá)的抗干擾性能進(jìn)行綜合的、定量的、有針對(duì)性的描述成為抗干擾研究中一個(gè)復(fù)雜而且極其重要的問題。近年來,研究人員在雷達(dá)抗干擾效果以及性能評(píng)估方面做了大量工作[2-4],并且針對(duì)不同體制雷達(dá)、不同干擾類型提出了一系列方法和指標(biāo),但是至今并未形成一個(gè)統(tǒng)一的評(píng)估方法和度量標(biāo)準(zhǔn)。

雷達(dá)抗干擾性能評(píng)估的意義在于科學(xué)地描述雷達(dá)在電磁環(huán)境中的工作能力,發(fā)現(xiàn)雷達(dá)在對(duì)抗各種干擾時(shí)的薄弱環(huán)節(jié),有針對(duì)性地提升雷達(dá)的綜合抗干擾能力。因此,一個(gè)科學(xué)、有效的雷達(dá)抗干擾性能評(píng)估方法,必須能夠綜合分析干擾對(duì)雷達(dá)影響效果的同時(shí),發(fā)現(xiàn)雷達(dá)抗干擾能力的瓶頸,以及針對(duì)不同的抗干擾方法、不同的干擾方式,給出相應(yīng)的抗干擾性能分析,為雷達(dá)的抗干擾能力建設(shè)提供強(qiáng)有力的數(shù)據(jù)支撐和理論依據(jù)。

本文從雷達(dá)抗干擾能力建設(shè)的實(shí)際需求出發(fā),提出了一種基于高維量測特征空間的統(tǒng)計(jì)建模方法,利用多元密度函數(shù)分析干擾效果以及雷達(dá)的綜合抗干擾能力,同時(shí)兼具了單項(xiàng)抗干擾性能分析,考慮了不同抗干擾手段、不同干擾方式的性能評(píng)估問題。

1 雷達(dá)抗干擾性能評(píng)估的一般理論與主要問題

1.1 評(píng)估準(zhǔn)則

評(píng)估雷達(dá)抗干擾性能的指標(biāo)一般都遵循如下三個(gè)準(zhǔn)則[5]:

功率準(zhǔn)則:通過干擾前后功率變化描述抗干擾能力,比如抗干擾改善因子(EIF)[6]等。一般適用于壓制式干擾下的抗干擾性能評(píng)估。

概率準(zhǔn)則:通過干擾成功概率或者雷達(dá)識(shí)別干擾成功概率來描述抗干擾能力。

信息準(zhǔn)則:通過干擾前后雷達(dá)信息的損失程度描述雷達(dá)的抗干擾能力。

1.2 評(píng)估方案

雷達(dá)抗干擾性能評(píng)估方案[7]包括:

試驗(yàn)法:通過試驗(yàn)獲取的實(shí)際雷達(dá)數(shù)據(jù)的分析,對(duì)雷達(dá)的抗干擾能力進(jìn)行評(píng)估。

仿真法:通過半實(shí)物仿真實(shí)驗(yàn),或者全數(shù)字仿真實(shí)驗(yàn),模擬真實(shí)電子對(duì)抗環(huán)境,對(duì)雷達(dá)的抗干擾能力進(jìn)行評(píng)估。

數(shù)學(xué)解析法:通過分析雷達(dá)和干擾機(jī)的性能指標(biāo)之間的關(guān)系,建立評(píng)估數(shù)學(xué)模型對(duì)雷達(dá)系統(tǒng)的抗干擾性能進(jìn)行評(píng)估。

1.3 面臨的主要問題

1)試驗(yàn)法與解析計(jì)算法融合的問題

試驗(yàn)法的優(yōu)點(diǎn)在于數(shù)據(jù)來源可靠,但考慮到成本問題,數(shù)據(jù)量有限。而數(shù)學(xué)解析法則沒有成本的限制,但是其數(shù)據(jù)來源于雷達(dá)和干擾機(jī)的性能指標(biāo),實(shí)時(shí)性較差?，F(xiàn)有的方法很少能夠?qū)崿F(xiàn)兩者優(yōu)勢的共享與劣勢的互補(bǔ)。

2)綜合能力評(píng)估和單項(xiàng)性能分析的問題

一般的抗干擾性能評(píng)估指標(biāo)都是有針對(duì)性的評(píng)估指標(biāo),不能對(duì)綜合能力與單項(xiàng)性能統(tǒng)一進(jìn)行評(píng)估,從而很難直接反映出雷達(dá)抗干擾能力的薄弱環(huán)節(jié),不利于抗干擾能力的建設(shè)和提升。

3)單項(xiàng)干擾與復(fù)合干擾統(tǒng)一評(píng)估的問題

在實(shí)戰(zhàn)中,不同類型雷達(dá)面臨的干擾不同,并且干擾的部署靈活可控,可以采用單一干擾,也可以采用復(fù)合干擾。一般評(píng)估指標(biāo)的最大問題在于評(píng)價(jià)依據(jù)與干擾方式相關(guān),缺少綜合評(píng)估的方法和指標(biāo)。

2 雷達(dá)抗干擾性能評(píng)估方法

2.1 基本方法

本文針對(duì)抗干擾評(píng)估面臨的四個(gè)問題,提出了一種新的評(píng)估理論,具體方法如圖1所示。

圖1 雷達(dá)抗干擾性能評(píng)估方法

首先,在實(shí)驗(yàn)環(huán)境中,干擾信號(hào)經(jīng)過雷達(dá)接收處理之后,輸出為量測數(shù)據(jù)。利用量測數(shù)據(jù)作為性能評(píng)估模型的輸入,可以屏蔽單項(xiàng)干擾與復(fù)合干擾對(duì)于評(píng)估指標(biāo)的影響。其次,該方法結(jié)合試驗(yàn)法與數(shù)學(xué)解析法,通過對(duì)量測數(shù)據(jù)多元概率分布的統(tǒng)計(jì)分析(如圖2所示),可以利用有限樣本數(shù)據(jù)分析總體性能,并且通過向具體量測特征軸上的映射,可以對(duì)雷達(dá)單項(xiàng)量測數(shù)據(jù)進(jìn)行分析,找出性能瓶頸。在實(shí)際應(yīng)用中,量測數(shù)據(jù)一般符合正態(tài)分布模型,尤其是當(dāng)樣本數(shù)據(jù)充足的時(shí)候,根據(jù)中心極限定理,可以利用正態(tài)模型進(jìn)行分析。但是在某些干擾環(huán)境中,很難用數(shù)學(xué)解析方法確定量測數(shù)據(jù)在特征空間中滿足的概率分布。因此,對(duì)于輸入的量測數(shù)據(jù),首先需要檢驗(yàn)其正態(tài)性。對(duì)于非正態(tài)數(shù)據(jù),可以通過正態(tài)轉(zhuǎn)換方法對(duì)數(shù)據(jù)進(jìn)行預(yù)處理,對(duì)于經(jīng)過正態(tài)轉(zhuǎn)換之后依然不滿足正態(tài)分布的數(shù)據(jù),可以通過概率密度估計(jì)的方法進(jìn)行建模。

2.2 理論分析

1)量測數(shù)據(jù)與真實(shí)目標(biāo)特性

定義試驗(yàn)獲取的N組雷達(dá)量測數(shù)據(jù)Xm為

式中,量測向量x i=(x1i,x2i,…,x Mi)T,i=1,2,…,N。且各量測向量之間滿足獨(dú)立同分布假設(shè),定義x i的期望和協(xié)方差矩陣為μm和Σm,即

圖2 量測數(shù)據(jù)解析流程

定義真實(shí)目標(biāo)信息為x=(x1,x2,…,x M)T,且x的期望和協(xié)方差矩陣為μ和Σ,即

要求Σ和Σm為正定矩陣,當(dāng)各量測量之間不相關(guān)時(shí),Σ和Σm為對(duì)角陣。

2)多元正態(tài)性檢驗(yàn)與數(shù)據(jù)預(yù)處理

由于雷達(dá)的量測向量的統(tǒng)計(jì)特性為多元分布函數(shù),因此需要選擇多元正態(tài)性檢驗(yàn)方法對(duì)數(shù)據(jù)進(jìn)行檢驗(yàn)。常用的多元正態(tài)性檢驗(yàn)方法包括卡方圖檢驗(yàn)法和主成分檢驗(yàn)法[8]?？ǚ綀D檢驗(yàn)法相對(duì)簡單,但是精度較低,而且在小數(shù)據(jù)情況下檢驗(yàn)性能較差。主成分檢驗(yàn)法步驟如下:

步驟1 對(duì)于量測數(shù)據(jù){Xm}M×N,根據(jù)式(6)、式(7)計(jì)算均值和協(xié)方差矩陣

式中,Λ=diag(λ1,λ2,…,λM),且λi＞0,i=1,2,…,M;Φ=(a1,a2,…,a M)為單位正交矩陣。

步驟3 令

步驟4 依次對(duì)各個(gè)主成分分量作一元正態(tài)性檢驗(yàn)[9]。當(dāng)且僅當(dāng)所有均滿足正態(tài)性檢驗(yàn)時(shí),認(rèn)為量測數(shù)據(jù)符合多元正態(tài)分布。對(duì)于不滿足正態(tài)性檢驗(yàn)的分量,需要進(jìn)行正態(tài)轉(zhuǎn)換預(yù)處理。一般的轉(zhuǎn)換方法有冪次轉(zhuǎn)換、Box-Cox轉(zhuǎn)換[10]和Johnson轉(zhuǎn)換[11]等,可以在一定程度上將非正態(tài)分布數(shù)據(jù)轉(zhuǎn)換為近似正態(tài)分布的數(shù)據(jù)。

3)多元正態(tài)分布建模

正態(tài)性檢驗(yàn)之后經(jīng)過正態(tài)轉(zhuǎn)換的數(shù)據(jù),可以通過正態(tài)或高斯密度函數(shù)進(jìn)行建模。M元正態(tài)密度的形式如下:

同理,對(duì)于真實(shí)目標(biāo)特性的建模如下:

式中,μ由目標(biāo)靜態(tài)位置決定,Σ由雷達(dá)性能決定,|Σ|為矩陣Σ的行列式,Σ-1為矩陣Σ的逆矩陣,記p(x)～N(μ,Σ)。

4)性能評(píng)估指標(biāo)

定義量測誤差影響因子(Measurement Error Impact Factor,MEIF):

圖3 多元正態(tài)分布和一元投影變換

定義量測穩(wěn)定性影響因子(Measurement Stability Impact Factor,MSIF):

式中,V表示距離μ為r的超橢球體積,Vm表示距離為r的超橢球體積,V M為單位超橢球體積。量測特征空間中超橢球體積大小Vm決定了量測數(shù)據(jù)的離散程度,而V代表了無干擾時(shí)雷達(dá)量測量的離散程度。因此,MSIF體現(xiàn)了干擾對(duì)雷達(dá)量測穩(wěn)定性的影響程度。

定義干擾影響因子(Jamming Impact Factor,JIF):

從定義可以看出,干擾影響因子(JIF)綜合體現(xiàn)了干擾對(duì)量測誤差和量測穩(wěn)定性的影響,因此可以體現(xiàn)出雷達(dá)的綜合抗干擾性能。

定義投影向量e i為

式中,除了第i個(gè)元素為1外,其余元素均為0。

那么,令

定義為真實(shí)目標(biāo)特征中心到實(shí)際量測數(shù)據(jù)中心在方向x i上投影的馬氏距離。那么,MEIF,MSIF和JIF在x i方向上定義如下:

3 仿真分析

MEIF,MSIF和JIF不僅可以用于定量驗(yàn)證抗干擾方法的實(shí)際效果,還可以評(píng)價(jià)復(fù)雜電磁環(huán)境中干擾對(duì)于雷達(dá)量測性能的影響。在仿真中,首先針對(duì)搜索雷達(dá),分析了無源沖淡箔條干擾對(duì)于雷達(dá)搜索目標(biāo)的影響。然后對(duì)于跟蹤雷達(dá)常見的距離-速度聯(lián)合拖引干擾,驗(yàn)證了雷達(dá)采用抗干擾手段之后的效果。雷達(dá)參數(shù)、目標(biāo)參數(shù)與干擾參數(shù)的設(shè)置詳見表1。

表1 仿真參數(shù)設(shè)定

3.1 干擾對(duì)雷達(dá)量測性能的影響評(píng)估

搜索雷達(dá)以發(fā)現(xiàn)與定位目標(biāo)為主要目的,其量測值為搜索坐標(biāo)系下目標(biāo)坐標(biāo),包括距離向坐標(biāo)和方位向坐標(biāo)。在箔條干擾環(huán)境中,量測數(shù)據(jù)及其分布如圖4所示。從量測數(shù)據(jù)仿真結(jié)果可以看出:雷達(dá)受到干擾,將目標(biāo)鎖定為干擾中心點(diǎn)附近。且量測數(shù)據(jù)經(jīng)過Anderson-Darling正態(tài)分布檢驗(yàn)后,均判定為正態(tài)分布。量測特征空間正態(tài)分布估計(jì)情況如圖5所示。

干擾對(duì)雷達(dá)量測性能的影響評(píng)估詳見表2。在這種情況下,MEIF與JIF接近于0,可見在這種場景下無源沖淡箔條干擾會(huì)對(duì)搜索類雷達(dá)產(chǎn)生較大的影響。

圖4 有干擾時(shí)搜索量測數(shù)據(jù)及其分布情況

圖5 量測特征空間正態(tài)分布估計(jì)情況

表2 干擾對(duì)雷達(dá)量測性能的影響評(píng)估

3.2 雷達(dá)抗干擾方法效能評(píng)估

以距離-速度聯(lián)合拖引干擾為例,對(duì)跟蹤雷達(dá)采用抗干擾方法前后的抗干擾性能進(jìn)行仿真。量測數(shù)據(jù)包括距離、方位角和目標(biāo)速度。

場景1:雷達(dá)無抗干擾方法,受到距離-速度聯(lián)合拖引干擾影響。

場景2:雷達(dá)通過增加速度門限,對(duì)抗距離-速度聯(lián)合拖引干擾。

場景1時(shí),量測數(shù)據(jù)及其分布如圖6所示。

從量測數(shù)據(jù)仿真結(jié)果可以看出:雷達(dá)在距離上被拖引;距離分布為均勻分布并且無法通過預(yù)處理轉(zhuǎn)變?yōu)檎龖B(tài)數(shù)據(jù)。因此將量測分布向方位角速度聯(lián)合特征空間進(jìn)行映射,映射后量測特征空間正態(tài)分布估計(jì)情況如圖7(a)所示。然后將量測分布分別投影到方位角和速度量測軸上,分布情況如圖7(b)、(c)所示。從圖中可以看出,距離 速度聯(lián)合拖引干擾對(duì)于雷達(dá)的影響主要體現(xiàn)在速度量測值上。因此,可以通過在速度域上加門限的手段,限制距離-速度聯(lián)合拖引干擾。

場景2時(shí),增加速度門限,量測數(shù)據(jù)及其分布如圖8所示。從量測數(shù)據(jù)仿真結(jié)果可以看出:加入速度門限之后,雷達(dá)在距離上沒有受到拖引干擾影響。特征空間正態(tài)分布估計(jì)情況如圖9所示。

圖6 距離拖引干擾時(shí)跟蹤量測數(shù)據(jù)及其分布情況

圖7 量測特征空間正態(tài)分布估計(jì)情況

圖8 增加對(duì)抗方法后量測數(shù)據(jù)及其分布情況

圖9 量測特征空間正態(tài)分布估計(jì)情況

雷達(dá)抗干擾方法效能評(píng)估結(jié)果詳見表3。從表3可以看出,增加速度門限后,MEIF與JIF顯著提升,說明雷達(dá)對(duì)于距離-速度聯(lián)合拖引干擾的抵抗能力顯著增強(qiáng)。

4 結(jié)束語

本文針對(duì)一般抗干擾性能評(píng)估理論中存在的問題,提出了一種基于高維量測特征空間統(tǒng)計(jì)模型分析的雷達(dá)在復(fù)雜電磁環(huán)境下的抗干擾性能定量評(píng)估方法。該方法將試驗(yàn)數(shù)據(jù)映射到高維量測特征空間中,通過統(tǒng)計(jì)理論分析量測值與真實(shí)目標(biāo)特性的差異,提出了MEIF、MSIF和JIF三種指標(biāo),既可以描述干擾對(duì)于雷達(dá)的影響,也可以反映雷達(dá)的抗干擾能力,并且對(duì)于單一干擾、復(fù)合干擾不敏感,對(duì)復(fù)雜電磁干擾環(huán)境具有良好的適用性。通過仿真,可以看出該方法理論可行,對(duì)雷達(dá)的抗干擾能力建設(shè)具有指導(dǎo)性意義。

表3 雷達(dá)抗干擾方法效能評(píng)估

[1]李浩,邱超凡,趙小亮.雷達(dá)抗干擾效能的多層次模糊評(píng)估方法[J].雷達(dá)科學(xué)與技術(shù),2012,10(2):143-149,155.

[2]徐沙,張潔.一種基于層次分析法的雷達(dá)抗干擾效能評(píng)估方法[J].航天電子對(duì)抗,2016,32(3):13-15.

[3]劉松濤,姜寧,劉振興.艦載電子對(duì)抗干擾效果在線評(píng)估方法[J].電子信息對(duì)抗技術(shù),2016,31(4):63-67.

[4]周振宇,賀志毅,劉建欣.末制導(dǎo)雷達(dá)抗干擾性能評(píng)估系統(tǒng)工程實(shí)現(xiàn)[J].現(xiàn)代防御技術(shù),2011,39(2):138-141.

[5]劉江波,舒敬環(huán).復(fù)雜電磁環(huán)境下的雷達(dá)抗干擾效能評(píng)估指標(biāo)體系構(gòu)建[J].信息化研究,2015,41(5):1-5.

[6]李亞南,韓壯志.雷達(dá)抗干擾性能評(píng)價(jià)指標(biāo)體系及發(fā)展趨勢[J].火力與指揮控制,2015,40(4):1-5.

[7]聶紅霞.雷達(dá)抗干擾效能評(píng)估方法探討[J].艦船電子工程,2013,33(4):75-77.

[8]朱寧,趙肖肖,唐慶華.多元正態(tài)性檢驗(yàn)三種方法的比較及SAS程序設(shè)計(jì)[J].蘇州大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,28(3):20-25.

[9]梁小筠.我國正在制定“正態(tài)性檢驗(yàn)”的新標(biāo)準(zhǔn)[J].應(yīng)用概率統(tǒng)計(jì),2002,18(3):269-276.

[10]BOX G E P,COX D R.An Analysis of Transformations[J].Journal of the Royal Statistical Society,1964,26(2):211-252.

[11]李長江,鄧文平,曹元元,等.基于Box-Cox變換與Johnson變換非正態(tài)過程能力分析[J].齊齊哈爾大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,31(1):66-70.

A Radar ECCM Performance Evaluation Method in Complex Electromagnetic Environment

CHEN Wendong1,2,TANG Bin2
(1.The Second Academy of China Aerospace Science and Industry Corporation,Beijing100854,China;2.Beijing Institute of Remote Sensing Equipment,Beijing100854,China)

According to the problems in radar ECCM performance evaluation,a statistic method based on geometry analysis in high-dimensional measurement feature space is proposed.The method projects the measurement information into the high-dimensional measurement space and utilizes geometry feature of the multivariate density distribution to define the measurement error impact factor(MEIF),the measurement stability impact factor(MSIF)and the jamming impact factor(JIF)to evaluate the ECM effect as well as the comprehensive ECCM performance.It is proved by computer simulation that the method is able to describe the radar ECCM performance exactly and the results are consistent with theoretic analysis.

complex electromagnetic environment;radar;electronic counter-countermeasures(ECCM);performance evaluation

TN973.3

A

1672-2337(2017)01-0073-08

10.3969/j.issn.1672-2337.2017.01.013

2016-06-05;

2016-09-01

陳文東男,1987年出生于山西永濟(jì),在讀博士研究生,主要研究方向?yàn)槔走_(dá)總體設(shè)計(jì)與雷達(dá)電子反對(duì)抗技術(shù)。E-mail;chen_eecs@sina.com

湯斌男,1977年出生,高級(jí)工程師,主要研究方向?yàn)槔走_(dá)總體設(shè)計(jì)與雷達(dá)電子反對(duì)抗技術(shù)。