郭忠偉　賴　彬　徐文濤　周兆才

(1.陸軍軍官學(xué)院二系　合肥　230031)(2.陸軍軍官學(xué)院一系　合肥　230031)

?

邊境反擊作戰(zhàn)中電磁影響下的雷達(dá)輻射可視化研究*

郭忠偉1賴彬1徐文濤2周兆才1

(1.陸軍軍官學(xué)院二系合肥230031)(2.陸軍軍官學(xué)院一系合肥230031)

針對(duì)邊境反擊作戰(zhàn)中復(fù)雜電磁環(huán)境對(duì)雷達(dá)的影響，給出了實(shí)現(xiàn)電磁影響下雷達(dá)輻射可視化的方法。首先明確了電磁干擾影響下雷達(dá)輻射的可視化的流程，然后對(duì)電磁干擾影響下的探測(cè)范圍進(jìn)行了建模并用Matlab進(jìn)行了仿真，最后運(yùn)用STK和OpenGL對(duì)電磁干擾影響下的雷達(dá)輻射數(shù)據(jù)進(jìn)行了可視化設(shè)計(jì)與實(shí)現(xiàn)。

邊境反擊作戰(zhàn);雷達(dá)輻射;可視化;數(shù)據(jù)采樣;面繪制

Class NumberTN911.73

1　引言

未來邊境反擊作戰(zhàn)中電磁環(huán)境異常復(fù)雜，電子設(shè)備不僅要受到對(duì)方的干擾，還要受到己方的自擾，雷達(dá)作為戰(zhàn)場(chǎng)上重要的電子裝備深受其影響，若能對(duì)雷達(dá)在戰(zhàn)場(chǎng)電磁干擾影響下的雷達(dá)波束進(jìn)行感知、分析和處理，并以三維等形式對(duì)其探測(cè)范圍進(jìn)行可視化直觀顯現(xiàn)，將不僅有助于指揮員在邊境反擊作戰(zhàn)復(fù)雜的戰(zhàn)場(chǎng)條件下合理運(yùn)用雷達(dá)和正確地實(shí)施指揮，而且有助于指揮人員對(duì)戰(zhàn)場(chǎng)電磁情況進(jìn)行準(zhǔn)確判斷。

2　電磁干擾影響下雷達(dá)輻射可視化流程

根據(jù)Dye提出的數(shù)據(jù)場(chǎng)可視化理論模型，電磁干擾影響下雷達(dá)輻射可視化主要包括三個(gè)階段：電磁干擾影響下雷達(dá)探測(cè)范圍建模、電磁干擾影響下雷達(dá)電磁數(shù)據(jù)采樣、電磁干擾影響下雷達(dá)電磁數(shù)據(jù)可視化繪制。電磁干擾影響下雷達(dá)探測(cè)范圍建模首先依據(jù)雷達(dá)方程計(jì)算雷達(dá)在干擾信號(hào)影響下在水平面上的二維探測(cè)區(qū)域包絡(luò)線，即雷達(dá)在干擾下各個(gè)方位角上的最大探測(cè)距離點(diǎn)在水平面上的投

影相連，構(gòu)成的封閉多邊形，然后計(jì)算由于干擾源造成的在方位角、俯仰角上與雷達(dá)的主瓣方向存在一定的夾角引起探測(cè)范圍的變化；電磁干擾影響下雷達(dá)電磁數(shù)據(jù)采樣則是根據(jù)電磁干擾影響下雷達(dá)探測(cè)范圍模型，將雷達(dá)三維探測(cè)范圍包絡(luò)按照一定的算法進(jìn)行采樣，得到數(shù)據(jù)采樣點(diǎn)，從而計(jì)算出空間各個(gè)點(diǎn)的電磁數(shù)據(jù)值，即生成三維電磁數(shù)據(jù)場(chǎng)；電磁干擾影響下雷達(dá)電磁數(shù)據(jù)可視化繪制則是將離散分布的雷達(dá)三維電磁數(shù)據(jù)場(chǎng)，按照一定的規(guī)則轉(zhuǎn)換為可視元素，生成形象、逼真的場(chǎng)景。

3　電磁干擾影響下雷達(dá)探測(cè)探測(cè)范圍建模

3.1雷達(dá)方程

設(shè)雷達(dá)發(fā)射和接收采用同一天線，則雷達(dá)接收到的回波功率為

(1)

式中，Pt為雷達(dá)發(fā)射峰值功率；Gt雷達(dá)發(fā)射天線增益；R為雷達(dá)與目標(biāo)之間的距離；σ為目標(biāo)的雷達(dá)截面積；λ為雷達(dá)波長(zhǎng)。

雷達(dá)接收到的回波功率Pr必須超過最小可檢測(cè)信號(hào)功率Simin，雷達(dá)才能可靠地發(fā)現(xiàn)目標(biāo)，當(dāng)Pr正好等于Simin時(shí)，就可得到雷達(dá)檢測(cè)該目標(biāo)的最大探測(cè)距離：

(2)

最小可檢測(cè)信號(hào)功率Simin常常用最小輸出信號(hào)噪聲功率比(S/N)omin來表示：

Simin=kTsBnFn(S/N)omin

(3)

式中，k為波爾茲曼常數(shù)1.38×10-23J/K；Ts以絕對(duì)溫度表示的雷達(dá)接收機(jī)噪聲溫度(K)；Bn為接收機(jī)等效噪聲帶寬(Hz)，一般情況下可認(rèn)為Bn≈1/τ，τ為脈沖寬度；Fn為接收機(jī)噪聲系數(shù)(倍)。

于是可得出：

(4)

3.2電磁干擾影響下水平面上雷達(dá)包絡(luò)線的計(jì)算

當(dāng)雷達(dá)未受干擾時(shí)，由于水平面上各個(gè)方位角的最大探測(cè)距離相等，根據(jù)雷達(dá)方程，所以得出二維包絡(luò)線上方位角為θ處的最大探測(cè)距離計(jì)算公式為

(5)

式中，θ為空間某點(diǎn)和雷達(dá)與目標(biāo)連線的夾角(度)，相當(dāng)于方位角；D0為檢測(cè)目標(biāo)信號(hào)所需的最小輸出信噪比，即檢測(cè)因子，D0=(S/N)omin；L為雷達(dá)各部分損耗的損耗因子L(L>1)；CB為帶寬校正因子，表示接收機(jī)帶寬失配所帶來的信噪比損失，匹配時(shí)CB=1，其他時(shí)CB>1。

當(dāng)雷達(dá)受干擾時(shí)(以遠(yuǎn)距離支援干擾為例)，設(shè)有M部干擾機(jī)干擾天線都指向雷達(dá)，第i(1≤i≤M)部干擾機(jī)和目標(biāo)在地面的投影相對(duì)于雷達(dá)的張角為θi(°)，與雷達(dá)距離為Rji，功率為Pji(W)，天線增益為Gji，帶寬為Bji，極化損失為γji(當(dāng)采用圓極化時(shí)取0.5)。把雷達(dá)到干擾機(jī)1的連線作為基線，按照逆時(shí)針方向，基線與雷達(dá)到干擾機(jī)i的連線的夾角分別為ξi(1≤i≤M)(度)，且ξ1=0，如圖1所示。

圖1　多部干擾機(jī)遠(yuǎn)距離支援干擾示意圖

第i部干擾機(jī)進(jìn)入雷達(dá)接收機(jī)輸入端外的干擾信號(hào)功率Prji為

(6)

式中：Pj為干擾機(jī)發(fā)射功率；γj為干擾信號(hào)對(duì)雷達(dá)天線的極化損失(當(dāng)采用圓極化時(shí)取0.5)；Bj為干擾機(jī)帶寬(Hz)。

由雷達(dá)方程，多部干擾機(jī)遠(yuǎn)距離支援干擾條件下雷達(dá)的最大探測(cè)距離為

(7)

將式(6)代入式(7)，得到多部干擾機(jī)遠(yuǎn)距離支援干擾條件下雷達(dá)最大探測(cè)距離為

(8)

在求雷達(dá)的二維包絡(luò)線時(shí)，θi可以采用方位角θ(以雷達(dá)到干擾機(jī)1的連線為基線)為參量代替。將方位角θ按照一定間隔從0°取到360°分別計(jì)算Rmaxj，并連接相應(yīng)的點(diǎn)，形成的曲線就是多部干擾機(jī)遠(yuǎn)距離支援干擾條件下雷達(dá)的二維包絡(luò)線。多部干擾機(jī)遠(yuǎn)距離支援干擾條件下二維包絡(luò)線在方位角θ處的最大探測(cè)距離Rj(θ)為

(9)

M=1時(shí)(特例情況)，即單部干擾機(jī)遠(yuǎn)距離支援干擾情況時(shí)，用Matlab進(jìn)行仿真，方位角θ從0°取到360°，結(jié)果如圖2所示。

圖2　單部干擾機(jī)遠(yuǎn)距離支援干擾時(shí)壓制區(qū)仿真圖

3.3電磁干擾影響下垂直面上雷達(dá)探測(cè)距離的計(jì)算

當(dāng)雷達(dá)未受到干擾時(shí)，雷達(dá)在方位角θ和俯仰角φ處的最大探測(cè)距離，等于該方位角θ處的最大探測(cè)距離與天線垂直面方向性函數(shù)fy(φ)的平方根的乘積：

式中，fy(φ)為雷達(dá)天線垂直面方向性函數(shù)；R(θ)為二維包絡(luò)線在方位角θ處的最大探測(cè)距離。

當(dāng)雷達(dá)受到干擾時(shí)，在方位角上干擾機(jī)與雷達(dá)的主瓣方向存在一定的夾角，造成進(jìn)入雷達(dá)的能量變化而引起探測(cè)距離的變化，同時(shí)在俯仰角上干擾機(jī)與雷達(dá)天線的主瓣方向不同的夾角也會(huì)引起雷達(dá)探測(cè)距離的變化。為能夠真實(shí)反映干擾在垂直面上的影響，需對(duì)天線垂直面方向性函數(shù)進(jìn)行修正：

(10)

式中，δ0.5為雷達(dá)的垂直波瓣寬度；δ為干擾機(jī)與雷達(dá)在垂直面上的夾角；φz為波束指向角。K常數(shù)，對(duì)于余割平方方向性函數(shù)仿真情況取0.04～0.06，對(duì)于辛格方向性函數(shù)和高斯方向性函數(shù)，取0.07～0.10。

在雷達(dá)受到干擾情況下，經(jīng)過對(duì)天線垂直面方向性函數(shù)的修正，雷達(dá)在方位角θ和俯仰角φ處的最大探測(cè)距離：

(11)

4　電磁干擾影響下雷達(dá)輻射數(shù)據(jù)采樣

電磁波可視化要求數(shù)據(jù)必須是離散化的，所以要將電磁干擾影響下雷達(dá)電磁輻射包絡(luò)按照一定的算法進(jìn)行采樣，得到數(shù)據(jù)采樣點(diǎn)。

以求點(diǎn)Aφn(θm)的坐標(biāo)為例，數(shù)據(jù)采樣算法步驟為：

步驟1：根據(jù)方位角θm、干擾源數(shù)量、干擾方式計(jì)算最大探測(cè)距離Rj(θm)；

步驟3：受氣象因素的影響對(duì)雷達(dá)站至采樣點(diǎn)Aφn(θm)的距離進(jìn)行修正：

exp(-0.115δj(φn)Raj(φn))

其中δj為在最大探測(cè)距離處大氣影響因素的衰減系數(shù)之和，Raj為電波傳播途徑上大氣衰減存在的距離；

步驟4：求采樣點(diǎn)Aφn(θm)的大地坐標(biāo)：

M_CoordTransform.RD_RAzElToBLH(Rjφn(θm),θm,φn,BA,LA,HA)(BA,LA,AA)為A點(diǎn)的雷達(dá)地理坐標(biāo)；

步驟5：求采樣點(diǎn)Aφn(θm)的地心直角坐標(biāo)

M_CoordTransform.BLHToXYZ (BA,LA,HA,XA,YA,ZA)。

按照一定角度間隔，將φn從0°取到90°，θm從0°取到360°，按照上述五步計(jì)算，得到每個(gè)采樣點(diǎn)的地心直角坐標(biāo)，從而建立采樣點(diǎn)集合AN×M。

5　電磁干擾影響下雷達(dá)探測(cè)電磁波可視化實(shí)現(xiàn)

圖3　電磁干擾影響條件下雷達(dá)探測(cè)范圍基本曲面繪制示意圖

電磁干擾影響下雷達(dá)探測(cè)電磁波的可視化實(shí)現(xiàn)需要將采樣點(diǎn)數(shù)據(jù)集AN×M與OpenGL圖元建模的方法結(jié)合起來。OpenGL繪制復(fù)雜曲面時(shí)，是按照一定次序繪制由復(fù)雜曲面分解出的一系列基本曲面的方法進(jìn)行的。由于雷達(dá)三維探測(cè)范圍包絡(luò)曲面非常復(fù)雜，因此要將其分解為一系列基本曲面，而后再進(jìn)行面繪制。如圖3所示，陰影部分是由集合AN×M中的三個(gè)相鄰的點(diǎn)Aφi(θm)、Aφi(θn)、Aφj(θn)圍成的基本曲面。以這個(gè)基本曲面的繪制為例，研究電磁干擾影響下雷達(dá)探測(cè)電磁波的可視化方法。

利用OpenGL的glBegin(GL_TRIANGLES)函數(shù)，結(jié)合三個(gè)頂點(diǎn)坐標(biāo)，就可以繪制基本曲面，偽代碼如下：

glPolygonMode(GL_FRONT_AND_BACK,GL_LINE);//線繪制模式

glBegin(GL_TRIANGLES));//繪制基本曲面

glVertex3d(Xφi(θm),Yφi(θm),Zφi(θm));

glVertex3d(Xφi(θn),Yφi(θn),Zφi(θn));

glVertex3d(Xφj(θn),Yφj(θn),Zφj(θn));

glEnd();//結(jié)束繪制

電磁干擾影響條件下雷達(dá)三維探測(cè)范圍包絡(luò)就是由許多個(gè)這樣的基本曲面組成的。從集合AN×M中取俯仰角φn及其相鄰的φn+1進(jìn)行采樣的兩組采樣數(shù)據(jù)點(diǎn)，分別為

Aφn={Aφn(θm)|φn}(1≤m≤M,1≤n≤N-1)

Aφn+1={Aφn+1(θm)|φn+1}(1≤m≤M,1≤n≤N-1)

將這兩組點(diǎn)集按照?qǐng)D3的繪制方法，方位角遍歷方位角集θM，就可以繪制俯仰角φn、φn+1之間的曲面。然后，改變?chǔ)課的取值，遍歷俯仰角集φN，按照前述繪制方法，即可繪制出整個(gè)曲面，從而實(shí)現(xiàn)電磁干擾影響下雷達(dá)輻射的三維可視化。

設(shè)計(jì)的雷達(dá)輻射可視化系統(tǒng)綜合運(yùn)用了STK和OpenGL。該系統(tǒng)主要包括基于STK/Terrain的地理信息系統(tǒng)、雷達(dá)電磁波展現(xiàn)子系統(tǒng)、采樣數(shù)據(jù)庫(kù)、想定數(shù)據(jù)庫(kù)以及系統(tǒng)展現(xiàn)過濾與控制模塊等構(gòu)成。圖4所示為某雷達(dá)在電磁干擾影響下電磁波可視化效果圖。由圖可以看出，在雷達(dá)受到電磁干擾的方向上，雷達(dá)探測(cè)區(qū)域被“撕開”一個(gè)缺口，受到干擾后，雷達(dá)的探測(cè)性能明顯降低。

圖4　電磁干擾影響下的雷達(dá)電磁波效果圖

6　結(jié)語(yǔ)

針對(duì)邊境反擊作戰(zhàn)中電磁環(huán)境對(duì)雷達(dá)的影響，給出了電磁干擾影響下雷達(dá)輻射建模和基于面繪制的可視化方法，并進(jìn)行了設(shè)計(jì)與實(shí)現(xiàn)。該技術(shù)可作為支撐技術(shù)嵌入雷達(dá)指揮系統(tǒng)，對(duì)雷達(dá)在復(fù)雜電磁環(huán)境下的探測(cè)范圍進(jìn)行可視化顯現(xiàn)，為雷達(dá)的運(yùn)用提供支持。

[1]陳淦濤,許稼.有源壓制干擾下雷達(dá)探測(cè)距離分析與計(jì)算[J].雷達(dá)科學(xué)與技術(shù)，2011，9(1)：13-17.

[2]李鋒，張宗佩.戰(zhàn)場(chǎng)電磁態(tài)勢(shì)3維可視化技術(shù)[J].測(cè)繪科學(xué)技術(shù)學(xué)報(bào)，2012，29(3)：223-225.

[3]張濤,曹婉,陳振宇.戰(zhàn)場(chǎng)環(huán)境與可視化技術(shù)[M].北京:軍事科學(xué)出版社,2008:290-303.

[4]齊鋒,劉雅奇,婁寧.一種簡(jiǎn)化的雷達(dá)三維探測(cè)數(shù)據(jù)生成方法研究[J].計(jì)算機(jī)仿真,2010，27(3)：12-15.

[5]鄧晶,張明智,李志強(qiáng).基于STK的信息作戰(zhàn)態(tài)勢(shì)三維可視化表現(xiàn)方法研究[J].系統(tǒng)仿真學(xué)報(bào)，2010,22(11)：2654-2659.

[6]陳淦濤,許稼.有源壓制干擾下雷達(dá)探測(cè)距離分析與計(jì)算[J].雷達(dá)科學(xué)與技術(shù)，2011，9(1):13-17.

[7]黃潔,黨同心,趙擁軍.VC和STK集成的途徑及其在仿真中的應(yīng)用[J].計(jì)算機(jī)仿真,2007,24(1):291-294.

[8]郭忠偉,楊烽.戰(zhàn)場(chǎng)電磁環(huán)境建模與仿真研究[J].炮兵學(xué)院學(xué)報(bào)，2011,31(1)：36-38.

[9]錢江林,陳國(guó)華,周雯雯.虛擬戰(zhàn)場(chǎng)環(huán)境三維可視化關(guān)鍵技術(shù)研究[J].裝備學(xué)院學(xué)報(bào)，2013,24(4)：81-84.

[10]王春雨,田斌,馮超.基于APM模型的電磁波三維可視化研究與實(shí)現(xiàn)[J].艦船電子工程，2014,34(4)：122-124.

[11]湯大榮,韓慧,袁可亮.戰(zhàn)場(chǎng)電磁環(huán)境本質(zhì)及可視化描述方法研究[J].兵工自動(dòng)化，2014,33(11)：40-43.

Visualization Research of Radar Radiation under the Influence of Electromagnetic in Border Counterattack

GUO Zhongwei1LAI bin1XU Wentao2ZHOU Zhaocai1

(1.2nd Department,Officer Academy of PLA,Hefei230031)(2.1st Department,Officer Academy of PLA,Hefei230031)

Aiming at the influence of the complex electromagnetic environment to radar,a method of visualization for radar radiation under the influence of electronic jamming is presented.Firstly,the flow of visualization is presented for radar radiation.Secondly,the radar detection range under the influence of electromagnetic jamming is modeled and the model is simulated with Matlab.Finally,the 3-D visualization scene of radar radiation under the influence of electromagnetic jamming is gotten by surface rendering the radiation data with STK and OpenGL.

border counter attack,radar radiation,visualization,data acquisition,surface rendering

2016年3月14日,

2016年4月29日

郭忠偉，男，博士，副教授，研究方向：有作戰(zhàn)指揮學(xué)、智能信息處理。賴彬，男，碩士研究生，研究方向：作戰(zhàn)指揮學(xué)。徐文濤，男，碩士，講師，研究方向：作戰(zhàn)指揮學(xué)、軍事地形學(xué)。周兆才，男，碩士研究生，研究方向：作戰(zhàn)指揮學(xué)。

TN911.73DOI：10.3969/j.issn.1672-9730.2016.09.033