熊久良，封吉平，韓壯志，何強(qiáng)，徐宏

（1.軍械工程學(xué)院光學(xué)與電子工程系，石家莊050003；2.解放軍94568部隊，鄭州450047）

組網(wǎng)紅外／雷達(dá)協(xié)同間歇式目標(biāo)跟蹤?

熊久良1，封吉平1，韓壯志1，何強(qiáng)1，徐宏2

（1.軍械工程學(xué)院光學(xué)與電子工程系，石家莊050003；2.解放軍94568部隊，鄭州450047）

為進(jìn)一步提高組網(wǎng)火控雷達(dá)間歇跟蹤系統(tǒng)性能，提出了基于組網(wǎng)的紅外／雷達(dá)協(xié)同間歇跟蹤方法。該方法利用紅外／雷達(dá)進(jìn)行協(xié)同間歇跟蹤，通過數(shù)據(jù)擬合對連續(xù)紅外數(shù)據(jù)與間歇雷達(dá)數(shù)據(jù)進(jìn)行時間配準(zhǔn)，并對配準(zhǔn)后的數(shù)據(jù)進(jìn)行壓縮，然后采用集中式序貫濾波融合算法對異類間歇數(shù)據(jù)進(jìn)行融合，形成連續(xù)的目標(biāo)跟蹤航跡。仿真結(jié)果表明，與組網(wǎng)火控雷達(dá)間歇跟蹤相比，該方法能夠得到更高的低截獲性能。

紅外／雷達(dá)；組網(wǎng)火控雷達(dá)；間歇輻射；目標(biāo)跟蹤；數(shù)據(jù)融合

1 引言

隨著科學(xué)技術(shù)的發(fā)展，反輻射導(dǎo)彈性能越來越高，使得戰(zhàn)場火控雷達(dá)的生存能力受到嚴(yán)重威脅。正是基于此背景，文獻(xiàn)［1］將間歇輻射理論應(yīng)用于組網(wǎng)火控雷達(dá)，提出了間歇式目標(biāo)跟蹤方法，該方法在保證火控雷達(dá)基本跟蹤能力的同時有效地提高了雷達(dá)的生存能力。然而，如何進(jìn)一步提高火控雷達(dá)的生存能力，是一個亟待解決的問題。

與雷達(dá)不同，紅外探測器不向空中輻射能量，具有探測精度高、抗干擾能力強(qiáng)等優(yōu)點。但是，紅外不能提供目標(biāo)的距離信息，所以大多數(shù)情況下將紅外與雷達(dá)配合使用，通過信息互補(bǔ)共同完成探測跟蹤任務(wù)。文獻(xiàn)［2，3］通過仿真和分析證明了紅外／雷達(dá)協(xié)同跟蹤可以有效地提高系統(tǒng)的跟蹤性能。文獻(xiàn)［4，5］研究了單紅外／機(jī)載雷達(dá)協(xié)同跟蹤與輻射控制技術(shù)，定義了相關(guān)門限指標(biāo)來控制機(jī)載雷達(dá)開關(guān)機(jī)，減少雷達(dá)輻射時間，以提高雷達(dá)的生存能力，仿真結(jié)果證明了此方法可以在一定程度上提高雷達(dá)的生存能力。結(jié)合文獻(xiàn)［1］的相關(guān)理論，本文提出基于組網(wǎng)的紅外／雷達(dá)協(xié)同間歇式目標(biāo)跟蹤方法，利用紅外與火控雷達(dá)的信息互補(bǔ)來進(jìn)一步減少火控雷達(dá)的輻射時間，提高雷達(dá)的生存能力。

2 紅外／雷達(dá)協(xié)同間歇跟蹤原理分析

考慮紅外／雷達(dá)同平臺配置，每部雷達(dá)配置一部紅外，假定組網(wǎng)系統(tǒng)由M組紅外／雷達(dá)組成，記為Si（i=1，2，…，M）。每組紅外／雷達(dá)工作時，雷達(dá)采用間歇跟蹤，而紅外進(jìn)行固定采樣跟蹤。系統(tǒng)跟蹤過程如圖1所示。

假設(shè)在k時刻紅外／雷達(dá)Si進(jìn)行協(xié)同間歇跟蹤，其余雷達(dá)進(jìn)行靜默跟蹤，紅外／雷達(dá)Si將目標(biāo)實時跟蹤量測數(shù)據(jù)傳送到融合中心進(jìn)行融合處理，處理后的數(shù)據(jù)送到相關(guān)武器系統(tǒng)，供其火力打擊。跟蹤精度評估環(huán)節(jié)和低截獲性能評估環(huán)節(jié)則分別對當(dāng)前的跟蹤質(zhì)量和低截獲性能進(jìn)行計算和評估，并在達(dá)到要求跟蹤精度的前提下，以獲取最小平均截獲概率密度［1］為準(zhǔn)則選擇最優(yōu)的脈沖發(fā)射間隔。此時，間歇輻射控制環(huán)節(jié)根據(jù)該參數(shù)控制雷達(dá)的脈沖發(fā)射間隔；同時，根據(jù)系統(tǒng)截獲概率值的變化，發(fā)送交替選通信號以確定下一時刻進(jìn)行工作的紅外／雷達(dá)組合。不考慮交替輪換工作過程中調(diào)整時間的影響，當(dāng)Si轉(zhuǎn)為靜默跟蹤后，Si+1間歇工作以繼續(xù)跟蹤目標(biāo)。以此循環(huán)往復(fù)，各雷達(dá)在融合控制中心的協(xié)同控制下，共同構(gòu)成對目標(biāo)的閉環(huán)交替間歇跟蹤。

3 基于數(shù)據(jù)擬合的時間配準(zhǔn)算法

考慮紅外／雷達(dá)同平臺配置，同步工作，可以忽略空間誤差。在間歇跟蹤過程中，紅外一直進(jìn)行固定采樣，而雷達(dá)間歇輻射信號，因此，需要對紅外和雷達(dá)的數(shù)據(jù)進(jìn)行時間配準(zhǔn)。本文重點對間歇跟蹤過程中的時間配準(zhǔn)問題進(jìn)行分析。

目前，多傳感器時間配準(zhǔn)算法主要有最小二乘法［6］、內(nèi)插外推法［7］以及數(shù)據(jù)擬合法［8］等。其中，最小二乘法是在假設(shè)兩傳感器的采樣周期之比為整數(shù)的基礎(chǔ)上建立的，對任意采樣周期情況的異類數(shù)據(jù)并不適用；內(nèi)插外推法是在一個選定的時間片內(nèi)進(jìn)行內(nèi)插或外推，而且無法外推到時間片外；數(shù)據(jù)擬合法是基于最小二乘準(zhǔn)則，經(jīng)過數(shù)據(jù)擬合得到一條或幾條平滑曲線，由擬合后的曲線外推得到傳感器在任意時刻的量測值。從本文的研究背景可以知道，間歇數(shù)據(jù)的特點是雷達(dá)采樣間隔自適應(yīng)變化，而紅外采用固定采樣，兩者采樣間隔并不一定成整數(shù)比例，并且往往需要將紅外的數(shù)據(jù)外推到后來的時間點上?？紤]到以上因素，本文采用基于數(shù)據(jù)擬合的時間配準(zhǔn)方法來對異類數(shù)據(jù)進(jìn)行實時配準(zhǔn)。

若己知數(shù)據(jù)集（xi，yi）（i=1，2，…，n），擬合的曲線為P（x），則定義誤差δi為

利用最小二乘法使誤差平方和達(dá)到最小，即

通常采用多項式形式的P（x），即在函數(shù)類φ=｛φ1，φ2，…，φm｝中找一個函數(shù)組合，使得

式中，aj（j=1，2，…，m）為待定系數(shù)；φ1（x），φ2（x），…，φm（x）線性無關(guān)且滿足m＜n。

使誤差平方和最小，即

從而，用最小二乘求擬合的問題，就是求P（x）使‖δ‖2取得最小值。這樣求擬合曲線函數(shù)的過程就轉(zhuǎn)化為求多元函數(shù)極小值的問題。多元函數(shù)如下：

定義

從而式（6）可以表示為

當(dāng)線性無關(guān)時，RTR可逆，方程組（10）有唯一解P（x），即：

4 紅外／雷達(dá)間歇數(shù)據(jù)融合跟蹤算法

4.1 傳感器量測模型

設(shè)目標(biāo)的離散狀態(tài)方程為

紅外可以測量目標(biāo)的方位角θI（k）和俯仰角φI（k），其量測方程可表示為

雷達(dá)可以測量目標(biāo)的距離rR（k）、方位角θR（k）和俯仰角φR（k），其量測方程可表示為

4.2 基于數(shù)據(jù)壓縮的集中式序貫融合算法

從紅外／雷達(dá)協(xié)同間歇跟蹤原理可以看出，整個跟蹤過程中，紅外可以得到連續(xù)等間隔的目標(biāo)量測數(shù)據(jù)，而雷達(dá)由于采用間歇輻射方式，得到的目標(biāo)數(shù)據(jù)是時斷時續(xù)的，并且是非等間隔變化的。系統(tǒng)間歇數(shù)據(jù)模型如圖2所示，圖中長實線表示雷達(dá)量測數(shù)據(jù)，短實線表示紅外量測數(shù)據(jù)，虛線表示融合后的量測數(shù)據(jù)。

從圖2中可以看出，組網(wǎng)內(nèi)的每組紅外／雷達(dá)只有在工作期間才能得到目標(biāo)數(shù)據(jù)，在間歇狀態(tài)時，沒有目標(biāo)量測數(shù)據(jù)。但從整個系統(tǒng)來看，若不考慮工作輪換期間的時間延遲，則整個組網(wǎng)系統(tǒng)將形成連續(xù)的目標(biāo)量測數(shù)據(jù)。本文采用基于數(shù)據(jù)壓縮的集中式序貫濾波融合算法對間歇數(shù)據(jù)進(jìn)行融合。

基于數(shù)據(jù)壓縮的集中式序貫濾波融合算法的中心思想是：首先用雷達(dá)量測進(jìn)行濾波初始化，得到初始狀態(tài)X（k-1／k-1）和初始誤差協(xié)方差矩陣P（k-1／k-1），根據(jù)采樣間隔選擇策略選取下一時刻雷達(dá)的采樣間隔T（k）。實時記錄T（k）時間內(nèi)紅外的采樣數(shù)據(jù)，經(jīng)過數(shù)據(jù)擬合外推到雷達(dá)采樣點上，然后對經(jīng)過時間配準(zhǔn)后的紅外／雷達(dá)量測進(jìn)行數(shù)據(jù)壓縮［9］，將壓縮后的值作為量測值進(jìn)行主動EKF濾波，得到k時刻的估計值X（k／k）和P（k／k）。同時，對雷達(dá)截獲概率進(jìn)行評估，根據(jù)評估結(jié)果選擇下一時刻工作的紅外／雷達(dá)組合。結(jié)合文獻(xiàn)［1］的雷達(dá)輻射控制方法及自適應(yīng)采樣策略，可以得到紅外／雷達(dá)間歇數(shù)據(jù)融合跟蹤步驟如下。

Step 1：利用紅外／雷達(dá)Si量測進(jìn)行濾波初始化，得到初始狀態(tài)X（k-1／k-1）和初始誤差協(xié)方差矩陣P（k-1／k-1）。

Step 2：計算k-1時刻跟蹤精度［1］C（k-1），根據(jù)遞推采樣公式得到下一時刻的采樣間隔T（k）。

Step 3：記錄T（k）時間內(nèi)紅外的量測值，利用二階數(shù)據(jù)擬合將紅外數(shù)據(jù)外推到雷達(dá)量測時刻。

Step 4：對配準(zhǔn)后的數(shù)據(jù)進(jìn)行壓縮。壓縮后的俯仰角量測φeff、方差Rφeff、方位角量測θeff、方差f分別為

式中，φR（k）為雷達(dá)俯仰角量測，RφR為其方差；θR（k）為雷達(dá)方位角量測為其方差；φI（k）為紅外俯仰角量測，RφI為其方差；θI（k）為紅外方位角量測，為其方差。

Step 5：將壓縮后的值作為k時刻的量測值Z（k）進(jìn)行如下主動EKF濾波：

式中，H（k）為k時刻非線性量測函數(shù)線性化Jacobian矩陣。

Step 6：計算截獲概率［1］Pi（t），若Pi（t）＞P0，則進(jìn)行工作輪換，即i=i+1；否則，執(zhí)行Step 7。

Step 7：記錄仿真時長t，若t≥tmax（tmax為跟蹤時間），則執(zhí)行Step 8，否則，令k=k+1，返回Step 2。

Step 8：跟蹤過程結(jié)束。

5 仿真實驗與分析

本文采用自適應(yīng)常加速模型［10］并結(jié)合擴(kuò)展卡爾曼濾波算法對目標(biāo)進(jìn)行跟蹤，在相同條件下對文獻(xiàn)［1］的方法與本文提出的方法進(jìn)行仿真比較。

假定組網(wǎng)系統(tǒng)由3部火控雷達(dá)組成，每部雷達(dá)統(tǒng)一配置一部紅外探測器，不考慮信息傳輸延遲，紅外／雷達(dá)1首先工作，采用輪換工作方式。雷達(dá)脈沖寬度τ=0.5μs，脈沖重復(fù)周期T=1ms，每隔兩個脈沖進(jìn)行一次數(shù)據(jù)處理。采樣周期上下限分別為Tmin=2ms、Tmax=0.4 s，初始采樣周期T0（1）=0.1 s，跟蹤精度門限為5m，調(diào)整因子η=1，各組紅外／雷達(dá)的觀測噪聲矩陣分別為

假定偵察系統(tǒng)頻率搜索周期Tf=0.1 s，瞬時頻率帶寬Δfe=10 MHz，頻率范圍Δf=2 GHz，截獲概率門限值P0=60%。組網(wǎng)間歇跟蹤時，單部雷達(dá)的截獲概率達(dá)到P0的瞬間進(jìn)行工作輪換。當(dāng)3組紅外／雷達(dá)輪換工作循環(huán)一次之后，各組紅外／雷達(dá)采取固定時間長度進(jìn)行工作，工作時間Tchange=10 s。設(shè)目標(biāo)初始位置（x0，y0，z0）為（5 km，5 km，5 km），三維空間運動軌跡為

仿真時間t=50 s，作50次Monte-Carlo仿真。仿真結(jié)果如圖3和圖4所示。為便于對比，文中僅給出了前30 s的仿真結(jié)果。

從圖3（a）和圖4（a）可以看出，兩種方法都能使系統(tǒng)的跟蹤精度在門限值范圍內(nèi)變化，并且由于雷達(dá)輪換工作前后，觀測噪聲、測量精度等條件的不同，使得兩種方法下的系統(tǒng)跟蹤精度都會產(chǎn)生不同程度的波動。

對比圖3（b）、圖4（b）和圖3（c）、圖4（c）可以知道，火控雷達(dá)組網(wǎng)間歇跟蹤情況下，雷達(dá)1截獲概率到達(dá)P0的時間約為11.6 s，工作期間截獲概率密度約為0.16 s-1，而在本文提出的跟蹤方法下時間約為20 s，截獲概率密度約為0.08 s-1?？梢钥闯?，本文提出的方法可以有效地延長雷達(dá)被定位的時間，提高火控雷達(dá)的電子反對抗能力。

6 結(jié)論

本文提出了基于組網(wǎng)的紅外／雷達(dá)協(xié)同間歇跟蹤方法，該方法利用紅外不向外輻射能量的特性，通過紅外與雷達(dá)的信息互補(bǔ)，進(jìn)一步減少了雷達(dá)信號的輻射時間，提高了火控雷達(dá)的低截獲性能。仿真結(jié)果也證明了該方法的有效性。從仿真結(jié)果中也可以看出，紅外／雷達(dá)的不同配置方式將會影響系統(tǒng)的整體效果，其影響機(jī)理及程度如何有待以后的繼續(xù)研究。

［1］熊久良，徐宏，韓壯志，等.一種基于跟蹤精度控制的組網(wǎng)火控雷達(dá)抗ARM方法［J］.電訊技術(shù)，2011，51（4）：64-68.

XIONG Jiu-liang，XU Hong，HAN Zhuang-zhi，et al.A Method for Fire-Control Radar Network against ARM Based on Tracking Accuracy Control［J］.Telecommunication Engineering，2011，51（4）：64-68.（in Chinese）

［2］Yuan Quan，Dong Chaoyang，Wang Qing.An Fuzzy Adaptive Fusion Algorithm for Radar／Infrared Based on Wavelet Analysis［C］／／Proceedings of 2007 IEEE International Conference on Control and Automation.Guangzhou：IEEE，2007：1344-1348.

［3］Anfu Zhu，Zhanrong Jing，Weijun Chen，et al.Data Fusion of Infrared and Radar for Target Tracking［C］／／Proceedings of the 2nd International Symposium on Systems and Control in Aerospace and Astronautics.Shenzhen：IEEE，2008：1-4.

［4］Olivier Baud，Laurent Ratton.Radar Scheduling with Emission Constraints［C］／／Proceedings of the 2002 SPIE Conference on Acquisition Tracking and Pointing XVI.Bellingham，USA：Society of Photo-Optical Instrumentation Engineers，2002：153-156.

［5］Li Anping，Jing Zhongliang，Hu Shiqiang.UKF-based multi-sensor passive tracking with active assistance［J］. Journal of Systems Engineering and Electronics，2006，17（2）：245-250.

［6］BlairW D，Rice TR，AlouaniA T，et al.Asynchronous Data Fusion for Target Tracing with a Multitasking Radar and Optical Sensor［C］／／Proceedings of the 1991 SPIE Conference on Acquisition Tracking and Pointing X.Bellingham，USA：Society of Photo-Optical Instrumentation Engineers，1991：234-245.

［7］王寶樹，李芳社.基于數(shù)據(jù)融合技術(shù)的多目標(biāo)跟蹤算法研究［J］.西安電子科技大學(xué)學(xué)報，1998，24（6）：111-114.

WANG Bao-shu，LIFang-she.The research onmultiple targets tracking based on the data fusion technique［J］.Journal ofXidian University，1998，24（6）：111-114.（in Chinese）

［8］屈耀紅，閆建國.曲線擬合濾波在無人機(jī)導(dǎo)航數(shù)據(jù)處理中的應(yīng)用［J］.系統(tǒng)工程與電子技術(shù)，2004，26（12）：1912-1914.

QU Yao-hong，YAN Jian-guo.Application of curve fitting filtering at data processing of unmanned air vehicle［J］.Systems Engineering and Electronics，2004，26（12）：1912-1914.（in Chinese）

［9］何友，修建娟，張晶煒，等.雷達(dá)數(shù)據(jù)處理及應(yīng)用［M］.北京：電子工業(yè)出版社，2009.

HE You，XIU Jian-juan，ZHANG Jing-wei，etal.Radar data processing with applications［M］.Beijing：Publishing House of Electronics Industry，2009.（in Chinese）

［10］潘泉.現(xiàn)代目標(biāo)跟蹤與信息融合［M］.北京：國防工業(yè)出版社，2009.

PANQuan.Modern Target Tracking and Information Fusion［M］.Beijing：National Defense Industry Press，2009.（in Chinese）

XIONG Jiu-liang was born in Linyi，Shandong Province，in 1986.He is now a graduate student.His research concerns radar signal processing.

Email：xjl3274147＠163.com

封吉平（1961—），男，河北石家莊人，1987年獲碩士學(xué)位，現(xiàn)為副教授、碩士生導(dǎo)師，主要從事導(dǎo)航、制導(dǎo)與控制方面的研究；

FENG Ji-ping was born in Shijiazhuang，Hebei Province，in 1961.He received theM.S.degree in 1987.He isnow an associate professor and also the instructor of graduate students.His research concerns navigation，guidance and control.

韓壯志（1972—），男，河北石家莊人，2001年獲博士學(xué)位，現(xiàn)為副教授、碩士生導(dǎo)師，主要從事雷達(dá)組網(wǎng)方面的研究；

HAN Zhuang-zhiwas born in Shijiazhuang，Hebei Province，in 1972.He received the Ph.D.degree in 2001.He is now an associate professor and also the instructor ofgraduate students.His research concerns radar network.

何強(qiáng)（1972—），男，河北石家莊人，2001年獲博士學(xué)位，現(xiàn)為副教授、碩士生導(dǎo)師，主要從事雷達(dá)信號處理方面的研究；

HEQiangwasborn in Shijiazhuang，HebeiProvince，in 1972. He received the Ph.D.degree in 2001.He isnow an associate professor and also the instructor of graduate students.His research concerns radar signal processing.

徐宏（1981—），男，河南新縣人，2011年獲博士學(xué)位，主要從事雷達(dá)組網(wǎng)方面的研究。

XUHongwas born in Xinxian，Henan Province，in 1981.He received the Ph.D.degree in 2011.His research c oncerns radar network.

Netted Infrared and Radar Synergistic Interm ittent Target Tracking

XIONG Jiu-liang1，F(xiàn)ENG Ji-ping1，HAN Zhuang-zhi1，HEQiang1，XU Hong2
（1.Optical and Electronic Engineering Department，Ordnance Engineering College，Shijiazhuang 050003，China；2.Unit 94568 of PLA，Zhengzhou 450047，China）

In order to improve the performance of intermittent tracking in fire-control radar network，infrared／radar synergistic intermittent trackingmethod is proposed.Themethod uses infrared and radar to track the target，and time registration for the data is done with data fitting，then sequential filtering fusion algorithm based on data compression is used for dissimilar data fusion，and continuous target track is formed.Simulation results show that the new method can achieve higher low interception performance.

infrared／radar；netted fire-control radar；intermittent radiation；target tracking；data fusion

The National Natural Science Foundation of China（No.60672143）

TN95

A

10.3969／j.issn.1001-893x.2011.11.002

熊久良（1986—），男，山東臨沂人，碩士研究生，主要從事雷達(dá)信號處理方面的研究；

1001-893X（2011）11-0005-06

2011-07-11；

2011-08-15

國家自然科學(xué)基金資助項目（60672143）