頻譜共存下面向多目標(biāo)跟蹤的組網(wǎng)雷達功率時間聯(lián)合優(yōu)化算法

時晨光 董 璟 周建江

(南京航空航天大學(xué)雷達成像與微波光子技術(shù)教育部重點實驗室 南京 210016)

1 引言

當(dāng)前，空戰(zhàn)逐漸由信息化向智能化過渡，通過多雷達協(xié)同能夠提高空天目標(biāo)的預(yù)警探測能力與情報保障能力，特別是能夠提高復(fù)雜戰(zhàn)場條件下對非合作特種軍用目標(biāo)的連續(xù)探測跟蹤能力，已成為國內(nèi)外研究人員的共識。與此同時，隨著軍事裝備技術(shù)的飛速發(fā)展，多雷達協(xié)同的作戰(zhàn)環(huán)境也變得越來越復(fù)雜，電磁頻譜環(huán)境的復(fù)雜性、時變性以及難以預(yù)測性，對雷達系統(tǒng)提出了嚴峻挑戰(zhàn)。如何讓雷達系統(tǒng)在與通信基站共享同一工作頻段的條件下，通過優(yōu)化自身的射頻輻射資源配置以提升其探測跟蹤性能，已經(jīng)成為一個熱點問題[1–7]。

對于共享頻譜的雷達和通信系統(tǒng)，需要設(shè)計有效的管理方案來減少兩者間的相互干擾，從而保證彼此都能正常工作[8–18]。針對該問題，Zheng等人[9]提出了一種非合作雷達/通信共存的自適應(yīng)干擾消除方法，優(yōu)先考慮保護通信系統(tǒng)性能，并通過求解凸問題對雷達參數(shù)和通信解調(diào)誤差進行估計，仿真結(jié)果驗證了算法的有效性與優(yōu)越性。魯彥希[10]以雷達對通信系統(tǒng)的干擾能量為約束條件，對組網(wǎng)雷達節(jié)點選擇與功率參數(shù)進行自適應(yīng)優(yōu)化設(shè)計。文獻[11]提出了用于多雷達系統(tǒng)中多目標(biāo)跟蹤的通信感知資源調(diào)度策略，考慮在最小化雷達對通信干擾的同時，對雷達發(fā)射功率和采樣間隔進行優(yōu)化分配，仿真結(jié)果表明，該算法能夠有效提高系統(tǒng)多目標(biāo)跟蹤性能。文獻[12]針對異構(gòu)雷達與通信系統(tǒng)之間的頻譜共存問題，通過分配有限的發(fā)射功率、駐留時間和共享帶寬資源，來提高組網(wǎng)雷達的跟蹤性能，同時保持通信下行鏈路的吞吐量水平。文獻[13]對通信系統(tǒng)傳輸協(xié)方差矩陣以及基于稀疏感知和矩陣補全的多輸入多輸出(Multi-Input Multi-Output,MIMO)雷達采樣方案進行聯(lián)合設(shè)計，旨在保持通信系統(tǒng)的平均容量和發(fā)射功率的同時，最大限度地降低雷達接收機處的有效干擾功率。Hessar等人[14]采用具有低復(fù)雜度的空頻域隔離機制，并結(jié)合具體的雷達通信系統(tǒng)操作策略解決頻譜共存問題。韓凱峰等人[15]提出了一種雷達通信頻譜共存場景下基于塊對角化的通信雷達波束設(shè)計方案，該方案以通信對雷達無干擾及通信用戶間無干擾為約束，在滿足功率資源約束的條件下，最大化通信系統(tǒng)和速率。仿真結(jié)果表明，所提算法適用于多種場景，且與現(xiàn)有算法相比，具有更低的計算復(fù)雜度和更好的通信性能。文獻[16]針對多雷達通信一體化系統(tǒng)相互干擾的問題，提出了基于正交頻分復(fù)用(Orthogonal Frequency Division Multiplexing,OFDM)子載波分配的雷達通信一體化系統(tǒng)，利用子載波間的互相正交性，消除雷達與通信用戶間干擾以及通信用戶之間的相互干擾。文獻[17]提出了頻譜共存下基于Stackelberg博弈的多基地分布式雷達功率分配算法，在滿足預(yù)先設(shè)定的目標(biāo)探測信干噪比閾值的情況下，最小化各雷達輻射功率，從而起到保護通信基站不受雷達射頻輻射干擾。

組網(wǎng)雷達具有多雷達探測資源協(xié)同運用與信息融合緊密結(jié)合的技術(shù)體制特點，且在目標(biāo)檢測、跟蹤、識別等諸多場景中均具有顯著優(yōu)勢[19]。射頻資源管理是組網(wǎng)雷達目標(biāo)跟蹤任務(wù)中重要的一環(huán)，通過對輻射功率、駐留時間、信號帶寬和輻射采樣間隔等射頻參數(shù)進行優(yōu)化，能夠有效提升雷達系統(tǒng)性能[20–35]。Zhang等人[20]對雜波環(huán)境下分布式MIMO雷達多目標(biāo)跟蹤的資源分配問題進行了研究，通過對雷達節(jié)點和功率資源進行自適應(yīng)優(yōu)化配置，最小化多目標(biāo)跟蹤誤差，并提出了次優(yōu)解排序方法對該優(yōu)化模型進行求解。Lu等人[21]提出了一種機載雷達路徑規(guī)劃與輻射資源聯(lián)合優(yōu)化算法，該算法將機載雷達路徑規(guī)劃和輻射資源作為優(yōu)化參數(shù)，并采用多目標(biāo)跟蹤誤差效用函數(shù)作為優(yōu)化目標(biāo)，采用多步分解法和梯度投影法對該問題進行求解，仿真結(jié)果表明，該算法有效提升了雷達系統(tǒng)多目標(biāo)跟蹤性能。文獻[22]提出了面向低截獲概率的機載雷達功率資源分配算法，該算法以包含雷達多目標(biāo)跟蹤誤差和截獲概率的加權(quán)為目標(biāo)函數(shù)，并采用兩步分解法對模型進行求解。仿真結(jié)果表明，該算法能夠有效降低跟蹤過程中雷達被敵方無源探測系統(tǒng)截獲的概率。文獻[23]提出了一種多目標(biāo)跟蹤場景下的相控陣雷達最優(yōu)資源分配算法，該算法以多目標(biāo)跟蹤性能作為約束條件，以最小化雷達資源消耗為優(yōu)化目標(biāo)，對雷達的功率和帶寬參數(shù)進行聯(lián)合優(yōu)化分配。Cheng等人[24]提出了同時多波束集中式MIMO雷達自適應(yīng)時空資源與波形聯(lián)合優(yōu)化算法，該算法采用智能優(yōu)化方法對雷達采樣周期、發(fā)射能量、子陣列數(shù)和多波束方向矢量等多個參數(shù)進行聯(lián)合優(yōu)化設(shè)計，能夠自適應(yīng)選擇工作模式，從而在最大化目標(biāo)跟蹤精度的同時，達到最小化雷達資源消耗的目的。文獻[25]以最小化雷達駐留時間和輻射功率資源加權(quán)和為優(yōu)化目標(biāo)，結(jié)合內(nèi)點法和匈牙利算法對雷達節(jié)點分配方式、輻射功率和駐留時間進行優(yōu)化，仿真結(jié)果表明，相較于資源均勻分配算法，該算法能夠有效提升雷達系統(tǒng)的射頻隱身性能。文獻[26]對多基地雷達之間的博弈優(yōu)化模型進行分析，將預(yù)先設(shè)定的信干噪比閾值作為約束條件，結(jié)合凸優(yōu)化方法和博弈思想，對多基地雷達功率分配和波束形成進行聯(lián)合優(yōu)化，仿真結(jié)果驗證了該算法對雷達間的干擾抑制能力。文獻[27]提出了基于Stackelberg博弈的組網(wǎng)雷達功率分配算法，該算法旨在最小化功率消耗的同時最優(yōu)化目標(biāo)跟蹤精度，仿真結(jié)果表明，該算法能夠有效提高雷達系統(tǒng)的資源利用率。嚴俊坤等人[28]提出了基于機會約束的集中式MIMO雷達功率分配算法，該算法以MIMO雷達高概率滿足多目標(biāo)跟蹤精度為約束條件，以最小化MIMO雷達的發(fā)射功率為優(yōu)化目標(biāo)，并采用凸優(yōu)化方法對優(yōu)化模型進行求解。文獻[29]提出了分布式MIMO雷達多目標(biāo)跟蹤接收波束資源分配算法，該算法能夠基于跟蹤周期中的反饋信息，采用內(nèi)點法實現(xiàn)接收波束與目標(biāo)之間的最優(yōu)分配。

綜上所述，上述研究成果為優(yōu)化組網(wǎng)雷達目標(biāo)跟蹤性能和解決雷達通信頻譜共存問題奠定了堅實的基礎(chǔ)。然而，針對多目標(biāo)跟蹤場景，文獻[10,11]等已有研究并未充分考慮雷達通信頻譜共存環(huán)境下組網(wǎng)雷達功率時間資源的聯(lián)合優(yōu)化配置，解決該問題對于提升我方雷達探測系統(tǒng)的目標(biāo)跟蹤精度、射頻資源利用率及對敵作戰(zhàn)效能具有重要意義。因此，本文圍繞多目標(biāo)跟蹤場景下雷達通信頻譜共存和組網(wǎng)雷達輻射資源管理領(lǐng)域，研究了頻譜共存下面向多目標(biāo)跟蹤的組網(wǎng)雷達射頻資源聯(lián)合優(yōu)化分配問題。首先，由于各優(yōu)化參數(shù)的約束不同，采用發(fā)射能量作為優(yōu)化參數(shù)的算法復(fù)雜度較高，因此針對頻譜共存環(huán)境下多目標(biāo)跟蹤問題，建立了組網(wǎng)雷達功率時間資源聯(lián)合優(yōu)化模型，以最小化多目標(biāo)跟蹤貝葉斯克拉默-拉奧下界(Bayesian Cramér-Rao Lower Bound,BCRLB)為優(yōu)化目標(biāo)，以給定的組網(wǎng)雷達射頻資源和預(yù)先設(shè)定的通信基站最大可容忍干擾能量閾值為約束條件，通過聯(lián)合優(yōu)化雷達節(jié)點選擇、發(fā)射功率和駐留時間等射頻輻射參數(shù)，提升組網(wǎng)雷達的多目標(biāo)跟蹤精度。仿真結(jié)果驗證了本文所提算法的可行性和穩(wěn)健性。

2 系統(tǒng)模型

考慮一個由N部雷達組成的組網(wǎng)雷達系統(tǒng)，各雷達節(jié)點分散部署于二維直角坐標(biāo)系中，并保持時間、空間、頻率同步。第n部雷達的位置坐標(biāo)可以表示為xR,n=(xR,n,yR,n),n=1,2,...,N。另外，假設(shè)組網(wǎng)雷達系統(tǒng)的探測范圍內(nèi)存在Q個運動目標(biāo)和C個通信基站。

2.1 目標(biāo)運動模型

式中，rq表示過程噪聲強度。

2.2 雷達量測模型

2.3 干擾模型

通信基站在一定頻帶范圍內(nèi)持續(xù)性向外輻射電磁信號，會使雷達系統(tǒng)的探測跟蹤性能下降。另外，在頻譜共存環(huán)境下，雷達在執(zhí)行目標(biāo)探測跟蹤任務(wù)時也會對通信基站的正常工作造成負面影響。因此，本節(jié)對雷達和通信基站兩者的相互干擾進行分析。首先，考慮通信基站對雷達的干擾，此類干擾可以看作多種隨機調(diào)制信號的組合。根據(jù)中心極限定理，假設(shè)信號形式可近似表示為占據(jù)一定帶寬的循環(huán)對稱帶限的復(fù)高斯序列，該信號的功率譜密度(Power Spectral Density,PSD)在雷達的工作通帶范圍內(nèi)表現(xiàn)出了均勻特性[36]，即

式中，Nm,n表示信號PSD強度；fU,m,n和fL,m,n分別表示第m個通信基站和第n部雷達共同覆蓋頻帶的上界和下界。定義第m個通信基站對第n部雷達產(chǎn)生的干擾信號的時域協(xié)方差矩陣Ξm,n可以表示為[10]

式中，Ξm,n的第u行第v列的元素為對應(yīng)信號延時的自相關(guān)函數(shù)，同時也是通信基站信號PSD的離散時間傅里葉逆變換；Ts表示信號采樣周期?？紤]到距離衰減等因素影響，第m個通信基站對第n部雷達施加的干擾噪聲為

式中，χm,n表示傳輸信道的強度系數(shù)；表示第m個通信基站在第n部雷達頻帶范圍產(chǎn)生的干擾序列；表示通信基站m相對于雷達n的方位角。假設(shè)各通信基站的發(fā)射信號相互正交，則所有通信基站對雷達n的干擾時域協(xié)方差矩陣可以表示為

式中，xC,m表示第m個通信基站的位置。由式(8)可知，(fU,m,n-fL,m,n)值越大，則通信基站對雷達的干擾能量越強；如果第m個通信基站和第n部雷達在頻譜上沒有覆蓋，則fU,m,n-fL,m,n=0，即可認為對應(yīng)的干擾強度為0。

雷達在跟蹤運動目標(biāo)時對通信基站的干擾是時變的，且與目標(biāo)的運動狀態(tài)有關(guān)。本文使用雷達發(fā)射信號能量的聯(lián)合空頻域分布來表征干擾強度，可以計算為

式中，sn=[sn(1),sn(2),...,sn(L)]T表示具有有限間隔的雷達波形序列；表示雷達發(fā)射信號能量的空間分布[10]，即

3 頻譜共存下面向多目標(biāo)跟蹤的組網(wǎng)雷達功率時間聯(lián)合優(yōu)化算法

針對頻譜共存環(huán)境下多目標(biāo)跟蹤資源分配問題，提出一種面向多目標(biāo)跟蹤的組網(wǎng)雷達功率時間聯(lián)合優(yōu)化分配算法，旨在滿足給定雷達系統(tǒng)射頻資源和預(yù)先設(shè)定的通信基站最大可容忍干擾能量閾值等約束的條件下，對雷達節(jié)點選擇、發(fā)射功率和駐留時間等射頻參數(shù)進行自適應(yīng)優(yōu)化配置，從而在保證通信基站正常工作的同時，有效提升組網(wǎng)雷達的多目標(biāo)跟蹤精度。

3.1 優(yōu)化模型建立

在建立優(yōu)化模型之前，需要推導(dǎo)表征目標(biāo)跟蹤精度衡量指標(biāo)的解析表達式。文獻[37]中指出，在參數(shù)無偏估計的條件下，BCRLB為運動目標(biāo)跟蹤的均方誤差提供了下界，可以用來表征組網(wǎng)雷達跟蹤運動目標(biāo)的性能衡量指標(biāo)。首先，推導(dǎo)出k時刻目標(biāo)q的預(yù)測貝葉斯信息矩陣表達式為

本文提出了一種頻譜共存下面向多目標(biāo)跟蹤的組網(wǎng)雷達功率時間聯(lián)合優(yōu)化算法，通過聯(lián)合優(yōu)化分配雷達節(jié)點選擇、發(fā)射功率和駐留時間等射頻輻射參數(shù)，在滿足給定的雷達射頻資源和預(yù)先設(shè)定的通信基站最大可容忍干擾能量閾值的條件下，最大限度地降低組網(wǎng)雷達系統(tǒng)的多目標(biāo)跟蹤誤差，可構(gòu)建如下數(shù)學(xué)優(yōu)化模型：

式中，Tmax和Tmin分別表示各雷達駐留時間的上界和下界；Pmax和Pmin分別表示各雷達發(fā)射功率的上界和下界；Ttotal和Ptotal分別表示照射各目標(biāo)的所有雷達駐留時間和發(fā)射功率之和；Emax表示通信基站所能容忍的最大雷達干擾能量閾值；表示k時刻每部雷達最多可跟蹤一個目標(biāo)；?max表示各時刻組網(wǎng)雷達跟蹤目標(biāo)q可選擇的最多雷達節(jié)點數(shù)目。

3.2 優(yōu)化模型求解

同樣地，設(shè)定發(fā)射功率和駐留時間兩個參數(shù)的優(yōu)先級相同，采用SDP算法對優(yōu)化模型(19)進行求解，得到所選擇雷達節(jié)點的發(fā)射功率和駐留時間分配值，進而可以得到(?max-1)個射頻資源配置方案及其對應(yīng)的多目標(biāo)跟蹤誤差。

步驟3 選取步驟2所有方案中最小的多目標(biāo)跟蹤誤差與其對應(yīng)的射頻資源優(yōu)化分配方案作為備選方案，并將該資源優(yōu)化分配結(jié)果作為下一次循環(huán)的初始值。跳轉(zhuǎn)至步驟1，直到連續(xù)兩次得到的備選方案多目標(biāo)跟蹤誤差差值小于某一設(shè)定的閾值時，跳出循環(huán)，并將得到的優(yōu)化分配方案作為k時刻跟蹤目標(biāo)q的最優(yōu)雷達節(jié)點選擇方案和功率時間優(yōu)化分配方案。

步驟4 確定下一個跟蹤的目標(biāo)，移除步驟3已選取的雷達節(jié)點，跳轉(zhuǎn)至步驟1，直到對所有目標(biāo)的跟蹤方案都完成優(yōu)化，即可得到k時刻組網(wǎng)雷達跟蹤多目標(biāo)時的雷達節(jié)點選擇方案和功率時間聯(lián)合優(yōu)化分配方案。上述算法流程如算法1所示。

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

為了驗證頻譜共存下面向多目標(biāo)跟蹤的組網(wǎng)雷達功率時間聯(lián)合優(yōu)化算法的可行性與有效性，本節(jié)進行了如下仿真：假設(shè)組網(wǎng)雷達系統(tǒng)由N=6部位置固定且已知的雷達組成，各部雷達的系統(tǒng)參數(shù)均相同。組網(wǎng)雷達系統(tǒng)在k時刻需要同時跟蹤Q=2個目標(biāo)，目標(biāo)1的初始位置為[-70,0] km，以速度[900,400] m/s勻速飛行；目標(biāo)2的初始位置為[70,80] km，以速度[-900,-400] m/s勻速飛行。在組網(wǎng)雷達探測區(qū)域中存在C=2個通信基站，如圖1所示，兩個通信基站位置分別為[50,50] km和[–40,50] km。雷達采樣間隔為ΔT=3 s，跟蹤持續(xù)過程時間為 150 s，共50幀。k時刻跟蹤某一目標(biāo)的雷達節(jié)點數(shù)目最大值?max=3，通信基站最大可容忍干擾能量閾值為Emax=4.2J。其中，組網(wǎng)雷達仿真參數(shù)設(shè)置如表1所示。本文針對RCS不變、RCS變化和Emax閾值變化3種仿真場景分別進行仿真。

表1 仿真參數(shù)設(shè)置Tab.1 Simulation parameter settings

圖1 組網(wǎng)雷達布陣及多目標(biāo)運動軌跡圖Fig.1 Deployment of radar networks and trajectories of multiple targets

4.1 RCS不變

在該仿真場景中，假設(shè)雷達節(jié)點觀測運動目標(biāo)的RCS均為1 m2，組網(wǎng)雷達布陣及多目標(biāo)運動軌跡如圖1所示。

圖2給出了組網(wǎng)雷達跟蹤目標(biāo)1和目標(biāo)2的節(jié)點選擇和功率時間資源優(yōu)化分配結(jié)果，黑色部分表示該雷達節(jié)點在相應(yīng)幀未被賦予目標(biāo)跟蹤任務(wù)，即=0；反之，則表示該雷達節(jié)點在相應(yīng)幀被賦予目標(biāo)跟蹤任務(wù)，即=1。從圖2可以看出，組網(wǎng)雷達系統(tǒng)會優(yōu)先選擇距離目標(biāo)較近的雷達節(jié)點進行跟蹤，且各部雷達的發(fā)射功率和駐留時間分配情況均隨著目標(biāo)運動狀態(tài)變化自適應(yīng)進行調(diào)整。例如，前28幀組網(wǎng)雷達系統(tǒng)均選擇了距離目標(biāo)1更近的雷達1和雷達2進行跟蹤，隨著目標(biāo)1逐漸遠離雷達1而靠近雷達2，分配給雷達1的射頻資源也在不斷增加；對于目標(biāo)2，前8幀組網(wǎng)雷達系統(tǒng)選擇了距離目標(biāo)2更近的雷達5和雷達6進行跟蹤，并向與目標(biāo)2距離變大的雷達6分配更多的射頻資源。

圖2 雷達節(jié)點選擇與功率時間資源優(yōu)化分配結(jié)果Fig.2 Radar node selection and allocation results of power and dwell time resources

同時，組網(wǎng)雷達系統(tǒng)的節(jié)點選擇結(jié)果也受到通信基站最大可容忍干擾能量閾值的影響。例如，從第40幀開始，目標(biāo)1與通信基站距離小，雷達對通信基站的干擾能量接近通信基站最大可容忍干擾能量閾值，導(dǎo)致每一幀都會選擇不同的雷達節(jié)點對目標(biāo)1進行跟蹤。同樣地，對于目標(biāo)2，從第9幀開始，組網(wǎng)雷達系統(tǒng)選擇對通信基站干擾較弱的雷達3替代雷達5對目標(biāo)2進行跟蹤。

為了進一步驗證本文所提算法的優(yōu)越性，將所提算法的多目標(biāo)跟蹤誤差與其他4種對比算法進行比較。

(1) 無通信干擾算法。該算法不考慮通信基站與組網(wǎng)雷達之間的相互干擾，采用本文所提兩步分解求解方法對組網(wǎng)雷達節(jié)點選擇、發(fā)射功率和駐留時間等射頻輻射參數(shù)進行優(yōu)化配置。

(2) 隨機節(jié)點選擇算法。該算法對各時刻跟蹤目標(biāo)的雷達節(jié)點進行隨機選擇，并采用本文所提求解方法對選擇的雷達進行發(fā)射功率和駐留時間優(yōu)化配置。

(3) 均勻時間分配算法。該算法固定各雷達駐留時間，僅優(yōu)化雷達節(jié)點選擇和發(fā)射功率分配。

(4) 均勻功率時間分配算法。該算法固定各雷達的發(fā)射功率和駐留時間，僅優(yōu)化雷達節(jié)點選擇。

定義k時刻目標(biāo)跟蹤的平均均方根誤差(Average Root Mean Square Error,ARMSE)為

式中，NMC為蒙特卡羅實驗次數(shù)；為第n次蒙特卡羅實驗得到的目標(biāo)位置估計，本文設(shè)定NMC=100。本文所提算法和上述4種算法的ARMSE對比結(jié)果如圖3所示。從圖3可以看出，本文所提算法相較于有固定資源分配算法和隨機節(jié)點選擇算法具有更低的多目標(biāo)跟蹤誤差，這說明雷達節(jié)點選擇和射頻資源優(yōu)化分配對提升多目標(biāo)跟蹤精度有較大影響，而且優(yōu)化的射頻資源越多，各目標(biāo)跟蹤精度也越好。由于無通信干擾算法中不存在通信基站，即不考慮雷達與通信基站間的相互干擾，其所得多目標(biāo)跟蹤誤差低于本文所提算法，這也進一步說明了通信基站的存在對組網(wǎng)雷達多目標(biāo)跟蹤性能有一定影響?？偟膩碚f，本文所提算法能夠在給定組網(wǎng)雷達射頻資源條件和預(yù)先設(shè)定的通信基站最大可容忍干擾能量閾值下，通過聯(lián)合優(yōu)化雷達節(jié)點選擇以及輻射功率與駐留時間等射頻資源配置，獲得相較于其他對比算法更優(yōu)的多目標(biāo)跟蹤精度。

圖3 RCS不變場景下ARMSE對比結(jié)果Fig.3 Comparison results of ARMSE with same RCS

4.2 RCS變化

在該仿真場景中，主要研究目標(biāo)RCS的變化對組網(wǎng)雷達射頻資源配置與多目標(biāo)跟蹤精度的影響。目標(biāo)1和目標(biāo)2相對于雷達3和雷達4的RCS模型為Swerling I型，其中，各目標(biāo)RCS變化情況如圖4所示，目標(biāo)相對于其他雷達節(jié)點的RCS值仍固定為1 m2。

圖4 各目標(biāo)RCS數(shù)值Fig.4 RCS values of each target

圖5給出了RCS變化情況下組網(wǎng)雷達跟蹤目標(biāo)1和目標(biāo)2的雷達節(jié)點選擇和功率時間資源優(yōu)化分配結(jié)果。從圖5可以看出，由于目標(biāo)回波信噪比與其RCS值呈正比關(guān)系，組網(wǎng)雷達系統(tǒng)傾向于選擇相對目標(biāo)RCS值較大的雷達節(jié)點來完成跟蹤任務(wù)。例如，各目標(biāo)相對于雷達3和雷達4的RCS值在大多數(shù)跟蹤時刻均高于其他雷達，相較于RCS不變條件下的節(jié)點選擇方案，由于目標(biāo)1和目標(biāo)2相對于雷達3和雷達4的RCS值在多數(shù)時刻均高于RCS不變場景下的RCS設(shè)定值，因此組網(wǎng)雷達選擇雷達3和雷達4來跟蹤目標(biāo)1和目標(biāo)2的時刻顯著增多；同時，為了獲得更好的目標(biāo)跟蹤精度，本文所提算法更多選擇3部雷達對目標(biāo)進行跟蹤。

圖5 RCS變化場景下雷達節(jié)點選擇與功率時間資源優(yōu)化分配結(jié)果Fig.5 Radar node selection and allocation results of power and dwell time resources with various RCS

同樣地，RCS變化情況下本文所提算法與其他對比算法的ARMSE對比結(jié)果如圖6所示。由于雷達節(jié)點選擇受到目標(biāo)RCS變化的影響，組網(wǎng)雷達的跟蹤性能略差于RCS不變情況下的跟蹤性能，但仍優(yōu)于同一場景下的其他對比算法，體現(xiàn)了本文所提算法的穩(wěn)健性。

圖6 RCS變化場景下ARMSE對比結(jié)果Fig.6 Comparison results of ARMSE with various RCS

4.3 Emax變化

在該仿真場景中，主要研究通信基站最大可容忍干擾能量閾值，即Emax變化對組網(wǎng)雷達多目標(biāo)跟蹤精度的影響。圖7給出了Emax=3 J時組網(wǎng)雷達跟蹤目標(biāo)1和目標(biāo)2的節(jié)點選擇和功率時間資源優(yōu)化分配結(jié)果。從圖7可以看出，相較于第1種仿真場景，由于Emax值變小，通信基站對雷達干擾能量的可容忍度降低，此時，組網(wǎng)雷達選擇了滿足約束條件下的最優(yōu)雷達射頻資源分配方案；類似地，組網(wǎng)雷達更多地選擇3部雷達對多目標(biāo)進行跟蹤。

圖7 Emax=3 J時雷達節(jié)點選擇與功率時間資源優(yōu)化分配結(jié)果Fig.7 Radar node selection and allocation results of power and dwell time resources with Emax=3 J

同樣地，圖8顯示不同Emax條件下本文所提算法與其他4種算法的ARMSE對比結(jié)果。從圖中可以看出，隨著Emax值不斷增加，通信基站對雷達干擾能量的可容忍度不斷提升，使得可供組網(wǎng)雷達配置的射頻輻射資源增多，于是，系統(tǒng)能夠得到更低的多目標(biāo)跟蹤誤差，所得ARMSE值也更逼近無通信干擾條件下的跟蹤誤差。

圖8 不同Emax條件下ARMSE對比結(jié)果Fig.8 Comparison results of ARMSE under different Emax conditions

5 結(jié)語

本文考慮了頻譜共存環(huán)境下通信基站對組網(wǎng)雷達多目標(biāo)跟蹤性能的影響，提出了一種面向多目標(biāo)跟蹤的組網(wǎng)雷達功率時間聯(lián)合優(yōu)化算法，以滿足組網(wǎng)雷達射頻資源和預(yù)先設(shè)定的通信基站最大可容忍干擾能量閾值為約束條件，以最小化多目標(biāo)跟蹤BCRLB為優(yōu)化目標(biāo)，通過聯(lián)合優(yōu)化分配雷達節(jié)點選擇、發(fā)射功率和駐留時間等參數(shù)，提高了組網(wǎng)雷達的多目標(biāo)跟蹤精度。仿真結(jié)果表明，多目標(biāo)跟蹤精度不僅與組網(wǎng)雷達和各目標(biāo)的空間位置關(guān)系及目標(biāo)散射特性有關(guān)，而且還受通信基站最大可容忍干擾能量閾值的影響；另外，與現(xiàn)有算法相比，本文所提算法能夠在保證通信基站正常工作的條件下，有效提升多目標(biāo)跟蹤性能。下一步將圍繞頻譜共存環(huán)境下組網(wǎng)雷達射頻輻射參數(shù)與飛行路徑聯(lián)合優(yōu)化問題進行研究。