袁西超 何 輝 楊 儉

(中國人民解放軍電子工程學(xué)院 合肥 230037)

基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法*

袁西超 何 輝 楊 儉

(中國人民解放軍電子工程學(xué)院 合肥 230037)

論文設(shè)計并驗證了基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法,通過無人機(jī)在任務(wù)區(qū)域搜索時位置改變引起的信號強(qiáng)度的變化,反向推算出電臺的位置,并且采用了幾種手段有效減弱了噪聲的影響,使得在信噪比較低的情況下仍能達(dá)到戰(zhàn)場通信態(tài)勢的精度需求,最后利用仿真驗證了其有效性。

戰(zhàn)場通信態(tài)勢; 態(tài)勢感知; 無人機(jī); 位置信息融合

(Electronic Engineering Institute of PLA, Hefei 230037)

Class Number TN914

1 引言

本文中的戰(zhàn)場通信態(tài)勢是戰(zhàn)場通信對抗電磁態(tài)勢[1]的子集,主要關(guān)注敵方通信電臺的數(shù)量、工作頻率、位置分布的狀態(tài)以及在空間域、時間域、頻率域上的活動趨勢[2]。戰(zhàn)場通信態(tài)勢感知是指通過某種手段快速搜索戰(zhàn)場,獲取相關(guān)要素信息并處理形成戰(zhàn)場通信態(tài)勢情報[3]。根據(jù)戰(zhàn)場通信態(tài)勢情報,可以輔助分析敵方力量編成、部署以及行動企圖,可以輔助評估作戰(zhàn)效果,可以引導(dǎo)通信偵察干擾行動,具有十分重要的意義[4]。在高技術(shù)戰(zhàn)爭條件下,戰(zhàn)場上通信電臺數(shù)量龐大、信號復(fù)雜、變化快速,傳統(tǒng)的用于搜索截獲、分析識別、測向定位的通信偵察手段難以適應(yīng)未來戰(zhàn)場通信態(tài)勢感知的需求。本文提出的基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法具有頻段寬、速度快、自動化、適應(yīng)能力強(qiáng)的優(yōu)點,充分發(fā)揮電子戰(zhàn)無人機(jī)的優(yōu)勢[5],能夠生成快速更新的戰(zhàn)場通信態(tài)勢情報,可有效對接文獻(xiàn)[6～7]描述的戰(zhàn)場電磁態(tài)勢感知體系。

2 方法概述

2.1 基本思想

無人機(jī)接收敵方電臺的通信信號經(jīng)過地—空傳播路徑,接近自由空間傳播,偵察接收機(jī)輸入端功率為

(1)式(1)中,Pr為偵察接收機(jī)輸入端功率,Pt為發(fā)射功率,Gt發(fā)射天線在接收機(jī)方向增益,Gr接收天線在發(fā)射機(jī)方向增益,λ為信號波長,d為接收機(jī)與發(fā)射機(jī)的距離。從式(1)可以看出,Pr與d的平方成反比。由于無人機(jī)空中飛行,使得無人機(jī)與電臺相對位置不斷變化,導(dǎo)致d和Pr跟隨變化。如果記錄下無人機(jī)飛行位置與對應(yīng)的Pr,則有可能根據(jù)式(1)反推出電臺位置。

2.2 影響因素分析

實際上每個飛行位置對應(yīng)的Pr不僅與d有關(guān),還受到多種因素的影響。根據(jù)式(1)分析如下:Pt不是恒定值,而是受信息源變化影響;Gt與Gr受無人機(jī)與電臺相對方位的變化影響;λ受信號頻率影響;電臺與無人機(jī)間的傳播路徑也不一定嚴(yán)格遵守自由空間傳播,可能受地形地貌的遮擋和多徑衰落影響;Pr的測量值還受到同時接收的背景噪聲影響。

通過準(zhǔn)確分析各項因素的影響,并采用合理手段進(jìn)行消除或限制,是研究基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法的基礎(chǔ)。信息源變化因素,雖然其瞬時值變化較大,但其短時間的平均值變化較小。天線增益因素,假設(shè)天線均采用豎直全向天線,則天線增益不隨方向改變,即使天線采用垂直極化定向天線,在天線主瓣方向內(nèi),也可近似認(rèn)為天線增益不變。信號頻率因素,由于以頻率區(qū)分信號,因此信號存在則說明頻率不變。傳播路徑因素和背景噪聲因素難以預(yù)測,但其一般具有突發(fā)性,可以通過數(shù)學(xué)手段減弱其影響。因此,d以外的其他影響因素都可以在一定程度上有效限制或消除。

2.3 設(shè)備需求

為實現(xiàn)基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法,通信對抗無人機(jī)系統(tǒng)應(yīng)具備以下功能:

1) 偵察天線是豎直全向天線,或垂直極化定向天線。

2) 偵察接收機(jī)能夠?qū)π盘柗糯蟆⒆冾l、采樣變?yōu)樾盘枖?shù)據(jù)。

3) 數(shù)據(jù)鏈設(shè)備能夠?qū)崟r將信號數(shù)據(jù)傳輸至地面控制站。

4) 地面控制站計算機(jī)能夠流暢運行實現(xiàn)本方法的軟件。

2.4 基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知過程

1) 無人機(jī)起飛并前往任務(wù)區(qū)域。

2) 無人機(jī)按照任務(wù)航線規(guī)劃在任務(wù)區(qū)域內(nèi)搜索飛行。

3) 偵察接收機(jī)將任務(wù)頻段內(nèi)的信號變?yōu)樾盘枖?shù)據(jù),并實時傳輸至地面。

4) 地面控制站計算機(jī)運行實現(xiàn)本方法的軟件,實時計算并標(biāo)示電臺的位置和頻率信息。

5) 無人機(jī)將任務(wù)區(qū)域搜索一遍則同步刷新一遍任務(wù)區(qū)域的戰(zhàn)場通信態(tài)勢。

3 方法設(shè)計

3.1 基于無人機(jī)位置和信號頻譜功率的電臺位置求解

信號頻譜功率是指將偵察接收機(jī)采集并傳回地面的信號數(shù)據(jù)進(jìn)行快速傅里葉變換,然后將信號頻譜峰值附近譜線求和,用以表示接收信號功率。顯然,信號頻譜功率是一部分接收信號功率加上信號頻帶內(nèi)的背景噪聲功率。根據(jù)式(1)可得:

(2)

(3)

d2=(xuav-x)2+(yuav-y)2+H2

(4)

式(2)由式(3)、式(4)補(bǔ)充說明,p為信號頻譜功率,σ為占接收信號功率的比例,ρ為傳播路徑因素的影響,根據(jù)影響因素分析,可將除d以外的其他因素歸并為W,由此產(chǎn)生的誤差通過后面的數(shù)學(xué)手段消除,(xuav,yuav)、(x,y)分別是無人機(jī)和電臺的位置(直角坐標(biāo)系表示),H為無人機(jī)相對電臺的飛行高度。

無人機(jī)在任務(wù)區(qū)域搜索時不斷采集電臺信號并形成信號數(shù)據(jù),每次信號數(shù)據(jù)和采集時的無人機(jī)位置都可以列出如式(2)所示的方程。由于只有x、y、W等三個未知數(shù),因此可以任意抽出三個方程組成方程組求解,如式(5)所示。

(5)

式(5)中,三個方程里面的W實際并不相等,但差值不大,在此認(rèn)為相等,因為無法準(zhǔn)確計算它們的值。這樣,就造成了解存在誤差,多個解在一起就會表現(xiàn)出散點,需要進(jìn)一步用數(shù)學(xué)手段找出最優(yōu)解。方程組的解如式(6)所示。

(6)

a=K2+1

(7)

b= 2KB-K[2x1+P(-2x1+2x2)]

-[2y1+P(-2y1+2y2)]

(8)

(9)

(10)

(11)

(12)

式(6)由式(7)～式(12)補(bǔ)充說明。

3.2 經(jīng)緯度坐標(biāo)與直角坐標(biāo)的轉(zhuǎn)化

由于無人機(jī)和電臺的位置一般用經(jīng)緯度表示,而解方程組用直角坐標(biāo)比較方便,因此應(yīng)當(dāng)進(jìn)行經(jīng)緯度坐標(biāo)和平面直角坐標(biāo)的變換。由于需要解方程組的區(qū)域范圍不大,地球曲面可在小范圍內(nèi)近似為平面,因此可選擇附近一點為坐標(biāo)原點,利用經(jīng)緯度距離換算公式求得無人機(jī)位置與原點位置的距離和方位,再通過三角函數(shù)換算為直角坐標(biāo)。同理,最終的電臺位置也可反向換算為經(jīng)緯度坐標(biāo)。

3.3 信號自動判決門限的設(shè)定

實際信號中信號強(qiáng)度和噪聲強(qiáng)度都在不斷變化,判決信號存在與否的門限也應(yīng)當(dāng)跟隨變化,稱為信號自動判決門限。信號自動判決門限由信號數(shù)據(jù)頻譜的均值和方差來表示,如式(13)所示。

J=mean(F)+k*var(F)

(13)

式中,J為判決門限,F為信號數(shù)據(jù)經(jīng)快速傅里葉變換后頻譜的模,其長度為信號數(shù)據(jù)長度的1/2,mean(F)是F的均值,var(F)是F的方差,k為系數(shù),可根據(jù)信號和噪聲情況人工調(diào)整。

3.4 信號頻譜功率曲線的生成

信號頻譜功率曲線是指信號頻譜功率隨采集次數(shù)變化的曲線,表示不同無人機(jī)位置的不同信號強(qiáng)度。生成信號頻譜功率曲線包括四步:

1) 根據(jù)信號自動判決門限,判決信號存在,并求出信號頻帶。

2) 在信號頻帶內(nèi)求出峰值頻率,作為信號中心頻率。

3) 將中心頻率左右Δf內(nèi)的譜線求和,作為信號頻譜功率。

4) 按每次采集并處理得到的信號頻譜功率畫出信號頻譜功率曲線。

其中,Δf可根據(jù)信號帶寬情況人工調(diào)整。

3.5 信號頻譜功率曲線的平滑

由于噪聲的影響,每次采集并處理得到的信號頻譜功率都存在誤差,表現(xiàn)為信號頻譜功率曲線在一定程度(與信噪比相關(guān))內(nèi)忽高忽低,嚴(yán)重影響方程組求解的準(zhǔn)確性,因此應(yīng)利用噪聲的隨機(jī)特性對信號頻譜功率曲線平滑處理,減弱噪聲的干擾。具體方法是將各次采集點的信號頻譜功率變?yōu)槠淝昂髆次采集點信號頻譜功率的均值,m可根據(jù)信號情況和求解情況人工調(diào)整。

3.6 電臺位置解集的生成

求解電臺位置解包括三步:

1) 挑選采集點。在信號功率頻譜曲線中,較大的值能更好反映信號強(qiáng)度,因此挑選的方法是從大到小依次選出3*n個采集點。

2) 組成方程組。將選出的采集點的對應(yīng)的無人機(jī)位置和信號頻譜功率,組成n個方程組。

3) 方程組求解。根據(jù)式(6)求出n個方程組的解。

3.7 電臺鏡像位置的消除方法

從式(6)可以看出,每個方程組都有2套解,其中1套解表示了電臺的位置,另1套解是電臺位置相對無人機(jī)位置的鏡像位置,應(yīng)予以消除。

根據(jù)方程組求解的原理,畫出圖1、圖2、圖3,圖中有向線段表示無人機(jī)航線,實心五角星表示電臺位置,空心五角星表示電臺鏡像位置,黑方塊表示采集點,圓圈的半徑表示方程組求解過程中得到的無人機(jī)與電臺的距離。由圖1看出,無人機(jī)直線飛行時難以分辨鏡像位置。從圖2、圖3看出,無人機(jī)弧形航線飛行時,無論電臺在內(nèi)側(cè)還是外側(cè),鏡像位置的點都比較散,因此可以利用這一特征將鏡像位置識別并消除。

圖1 直線航線

圖2 弧形航線內(nèi)側(cè)

圖3 弧形航線外側(cè)

3.8 電臺位置解集的優(yōu)化

由于噪聲的影響,n個方程組的解不在同一位置,而是出現(xiàn)散點,為進(jìn)一步減弱噪聲的干擾,可去掉較為分散的點,方法包括三步:

1) 計算n個解表示的點之間的距離。

2) 計算每個點到其他點之間距離的總和。

3) 按照每個點的距離總和由小到大排序,舍棄后q%的點。

其中,q可以根據(jù)解的情況人工調(diào)整。

3.9 最優(yōu)電臺位置估計

優(yōu)化后的解集是較為集中的散點,而最終解必須是1個點,因此采用求重心的方法估計出最優(yōu)電臺位置。具體方法是分別求x坐標(biāo)和y坐標(biāo)的均值。

4 仿真驗證

4.1 仿真條件

任務(wù)區(qū)域范圍為200km×200km,設(shè)置兩部電臺分別在(100,40)、(120,80),電臺信號中心頻率分別為70MHz、80MHz,均為BPSK調(diào)制,信息碼速率為100KHz,采樣頻率為200MHz,無人機(jī)沿航線飛行,相對電臺的飛行高度為3000m,仿真如圖4、圖5所示。

圖4 任務(wù)區(qū)域仿真圖

圖4中直線為無人機(jī)航線,航線上的“★”表示無人機(jī)當(dāng)前位置,“”表示電臺位置。

圖5 信號頻譜圖

圖5中可以看出在背景噪聲中70MHz、80MHz處的信號頻譜。

4.2 無噪聲條件下信號頻譜功率曲線與解集

下面仿真驗證式(6)求解是否正確。關(guān)掉噪聲,生成信號頻譜功率曲線如圖6所示,圖中橫坐標(biāo)表示航線上的采樣點數(shù),通過式(6)求解電臺位置如圖7所示,圖7中航線兩側(cè)的點為求解得到的電臺位置,其中1個點與電臺真實位置重合,另1個點在鏡像位置。

圖6 無噪聲條件下信號頻譜功率曲線圖

圖7 無噪聲條件下求解電臺位置

4.3 噪聲條件下信號頻譜功率曲線的平滑

下面檢驗噪聲對求解的影響,在噪聲條件下信號頻譜功率曲線如圖8所示,由其求解得到的電臺位置如圖9所示,可以看出由于噪聲影響,電臺位置變?yōu)榱松Ⅻc。

圖8 無平滑處理的信號頻譜功率曲線圖

圖9 無平滑處理求解電臺位置

下面對圖8所示的信號頻譜功率曲線進(jìn)行平滑處理,結(jié)果如圖10所示,由其求解得到的電臺位置如圖11所示,可以看出散點情況得到了較好收斂。

圖10 平滑處理后信號頻譜功率曲線圖

圖11 平滑處理后求解電臺位置

4.4 直線航線與弧形航線的解集對比

下面驗證通過弧形航線消除鏡像位置。當(dāng)電臺位于弧形外側(cè)時,鏡像位置的散點對比如圖12所示。當(dāng)電臺位于弧形內(nèi)側(cè)時,鏡像位置的散點對比如圖13所示。由圖12、圖13可以看出,鏡像位置的散點均較真實位置的情況嚴(yán)重,電臺位置弧形內(nèi)側(cè)時尤為明顯。

圖12 電臺在弧形外側(cè)的鏡像對比圖

圖13 電臺在弧形內(nèi)側(cè)的鏡像對比圖

4.5 噪聲與弧形航線條件下的最優(yōu)解與真實位置對比

下面驗證采用全部減弱噪聲影響的手段時得到的最優(yōu)解位置與真實位置的偏差。在如圖14所示的信號最強(qiáng)時信噪比條件下,得到如圖15所示最優(yōu)解與真實位置偏差。圖15中“★”號為最優(yōu)解位置,其與真實位置偏差約為800m。該仿真驗證了在較低信噪比條件下基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法得到的最優(yōu)解能夠滿足戰(zhàn)場通信態(tài)勢標(biāo)繪的精度需求。

圖14 信號最強(qiáng)時的頻譜圖

圖15 最優(yōu)解位置與真實位置偏差

5 應(yīng)用前景分析

空中測量電臺位置的手段主要有單站移動測向定位和雙站交叉測向定位,雖然定位精度較高,但均需要多通道接收機(jī)和測向天線陣,并且同時只能對1個信號測向定位,難以滿足戰(zhàn)場通信態(tài)勢快速更新的需求[8～11]。本文介紹的基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法使無人機(jī)在任務(wù)區(qū)域搜索一遍就可以標(biāo)出任務(wù)頻段內(nèi)的所有電臺的位置,具有顯著的搜索效率,并且設(shè)備需求容易達(dá)到,能夠充分發(fā)揮無人機(jī)目標(biāo)小、成本低、無人員傷亡的優(yōu)勢,雖然估計的電臺位置精度不高,但用于引導(dǎo)通信偵察、通信干擾和分析敵方部署和行動已足夠,采用本方法快速更新的戰(zhàn)場通信態(tài)勢情報具有重要意義。

6 結(jié)語

本文從戰(zhàn)場通信態(tài)勢需求出發(fā),設(shè)計并驗證了基于位置信息融合的無人機(jī)戰(zhàn)場通信態(tài)勢感知方法,從基本思想和影響因素分析入手,逐步完成了信號頻譜功率和飛行位置關(guān)系的方程組組建和求解,減弱噪聲影響的幾種手段的設(shè)計,并使用Matlab進(jìn)行了仿真,驗證了方法的有效性。雖然本方法基于全向天線和常規(guī)定頻信號設(shè)計,但對使用定向天線的信號以及直接序列擴(kuò)頻信號、跳頻信號也可能有效,可進(jìn)一步研究。

[1] 周峰,何俊,胡振彪.仿真系統(tǒng)中通信對抗電磁態(tài)勢顯示研究[J].四川兵工學(xué)報,2011,32(4):77-80.

[2] 湯扣林.電磁態(tài)勢的氣象環(huán)境影響分析[J].指揮信息系統(tǒng)與技術(shù),2014,5(5):20-24.

[3] 周倜,王小非,陳煒.海戰(zhàn)場電磁態(tài)勢感知模型[J].火力與指揮控制,2013,38(8):1-5.

[4] 康安康.一體化通信對抗威脅評估與資源分配方法研究[J].信息通信,2015(12):219.

[5] 蔣盤林.無人電子戰(zhàn)作戰(zhàn)平臺系統(tǒng)與技術(shù)發(fā)展研究[J].通信對抗,2008(3):3-8.

[6] 馮德俊,朱江,李方偉.戰(zhàn)場電磁態(tài)勢感知關(guān)鍵技術(shù)研究[J].數(shù)字通信,2013,40(5):20-23.

[7] 薛磊,沈陽,喬亞.戰(zhàn)場電磁態(tài)勢融合感知技術(shù)框架研究[J].指揮控制與仿真,2016,38(2):1-5.

[8] 湯扣林,劉韻,趙春東.單站測向定位技術(shù)研究[J].火力指揮與控制,2009,34(12):112-116.

[9] 韓慧珠.單站測向無源定位技術(shù)研究[D].成都:電子科技大學(xué)碩士學(xué)位論文,2014.

[10] 徐冉.無源測向定位算法研究[D].長春:長春理工大學(xué)碩士學(xué)位論文,2013.

[11] 張秋月.新型無源探測系統(tǒng)中的無源定位技術(shù)研究[D].哈爾濱:哈爾濱工程大學(xué)碩士學(xué)位論文,2005.

A Method of Battlefield Communication Situation Awareness by UAV Based on Position Information Fusion

YUAN Xichao HE Hui YANG Jian

The method of battlefield communication situation awareness by UAV based on position information fusion was designed and validated in the paper. The change of signal intensity was brought on the change of UAV position when UAV was searching in mission area, the method reckoned the broadcasting station position reversely, and weakened the influence of noise effectively by several artifices. The result of the method could reach the precision requirement of communication’s situation on battlefield. Finally, the availability of the method was validated by simulation.

battlefield communication situation, situation awareness, UAV, position information fusion

2016年6月12日,

2016年7月30日

袁西超,男,碩士,講師,研究方向:電子對抗無人機(jī)作戰(zhàn)運用、通信對抗技術(shù)、計算機(jī)仿真技術(shù)。何輝,男,碩士,講師,研究方向:通信對抗技術(shù)、電子對抗情報。楊儉,男,碩士,助教,研究方向:電子對抗無人機(jī)作戰(zhàn)運用、雷達(dá)對抗技術(shù)。

TN914

10.3969/j.issn.1672-9730.2016.12.012