基于隨機規劃的投擲式通信干擾兵力部署建模研究

沈楠，何俊，齊鋒，田志峰

（解放軍電子工程學院，安徽 合肥230037）

1 引 言

投擲式通信干擾機干擾距離近，能以較小功率獲得較高的干信比，同時具有“數量多”的優點，能近距離、廣空域、寬頻段、大數量地壓制無線戰術通信的有效運行，已成為通信對抗裝備發展的一種趨勢。

在典型作戰背景下，投擲式通信干擾機壓制無線戰術通信的兵力部署問題是高效發揮其作戰效能、實現科學輔助決策的關鍵。目前，國內外文獻主要在分析投擲式通信干擾機作戰效能[1，2，3]的基礎上，研究了對任務區域內的指揮、報知通信實施投擲式干擾的兵力部署問題[3，4]，建立了矩形陣形和蜂窩狀陣形的部署模型，依據效費比指標論證了蜂窩狀陣形優于矩形陣形，但蜂窩狀陣形在發射機與任務區域距離較近時，具有干擾機需求量過大的不足；對協同通信實施投擲式干擾的兵力部署問題鮮有研究。

本文依據對無線戰術通信實施投擲式干擾的作戰任務，建立了基于隨機機會約束規劃的投擲式通信干擾兵力部署模型，并采用隨機模型、遺傳算法相結合的混合智能算法求解模型，實現了對無線戰術通信干擾壓制概率滿足一定置信水平的約束下，干擾機需求量最少的部署規劃。

2 問題描述與模型建立

2.1 作戰任務描述

對無線戰術通信實施投擲式干擾的作戰任務為[5，6]：在作戰的關鍵時節，對主要作戰方向或作戰地域上重要頻段內的無線戰術通信實施多批次、大數量、寬頻段壓制性干擾，使任務區域內的通信接收機受壓制干擾而不能正常工作。對無線戰術指揮、報知通信的投擲式干擾如圖1所示，對無線戰術協同通信的投擲式干擾如圖2所示，圖1、圖2中T表示通信發射機，R表示通信接收機，Ji表示投擲式通信干擾機。

通過以上分析，投擲式通信干擾兵力部署的目標為：規劃干擾機的數量和排布，實現對無線戰術通信（指揮、報知、協同通信）的干擾壓制概率滿足一定置信水平的條件下最小化干擾機的需求量。

2.2 隨機機會約束規劃簡介

隨機機會約束規劃[7，8]主要用于約束條件中含有隨機變量，且必須在觀測到隨機變量的實現之前做出決策的優化問題。考慮到所做決策在一些比較極端的情況下可能不滿足約束條件，而這些情況出現的概率很低，為避免由此引起的優化方案過于保守，隨機機會約束規劃方法允許所做的決策在一定程度上不滿足約束條件，但該決策應使約束條件成立的概率不小于某一置信水平。隨機機會約束規劃的一種常見形式為：

式（1）中，x為n維決策向量；ξ為已知概率密度函數Φ（ξ）的隨機向量；f（x）為目標函數；gj（x，ξ）（j=1，2，…，k）為隨機約束函數；Pr{?}表示{?}中的事件成立的概率；α為給定約束條件的置信水平。

2.3 基于隨機機會約束規劃的投擲式通信干擾兵力部署模型

式（2）中，N是預設干擾機數量的上限；Pi是通信接收機輸入端接收第i個干擾機的干擾功率；Ps是通信接收機輸入端接收發射機的信號功率；C是常數，且：

式（3）中，Pj為干擾機發射功率，Gj為干擾天線主瓣增益，γj為極化損失，Brj為干擾信號進入接收機的比例，εr為相對介電常數，Hj為干擾天線高度，Pt為發射機發射功率，Gt為發射天線主瓣增益，Ht為發射天線高度。

本文以干擾機需求量最小化為目標，允許部署方案在某些比較極端的情況下不滿足對通信接收機的干擾壓制，但這種情況發生的概率必須小于某一置信水平，在此前提下建立基于隨機機會約束規劃的投擲式通信干擾兵力部署模型：

式（4）中，nu是實際干擾機需求量，即干擾機部署方案m中坐標不為（—1，—1）的干擾機數量；N是預設干擾機數量的上限α表示干擾機部署方案m滿足對通信接收機干擾壓制概率達到置信水平α；length是通信干擾任務區的長度；width是通信干擾任務區的寬度。

如果約束條件的可行集為空集，那么干擾任務就無法完成，此時定義懲罰函數如下：

3 基于混合智能算法的模型求解

3.1 兵力部署模型中機會約束的檢驗

對給定的決策變量m即干擾兵力部署方案，采用隨機模擬方法檢驗機會約束是否滿足的方法如下：

步驟1：給定允許的隨機模擬次數M，置隨機模擬計數器t=0，置機會約束成立計數器N′=0；

步驟2：產生一個[0，length]內的均勻分布的偽隨機數xrt，產生一個[0，width]內的均勻分布的偽隨機數yrt，（xrt，yrt）為隨機變量通信接收機位置的實現值；

3.2 采用混合智能算法求解兵力部署模型

采用基于隨機模擬、遺傳算法的混合智能算法求解式（4）描述的機會約束干擾兵力部署模型，用懲罰函數方法處理約束，即目標函數和懲罰的約束項一起形成適應度函數，主要的求解步驟如下：

步驟1：輸入原始數據；

步驟2：輸入遺傳算法中要求的染色體（候選干擾兵力部署方案）個數，以及交叉和變異概率；

步驟3：采用隨機方法產生一組初始部署方案，作為遺傳算法的初始種群。對于種群中的每一個染色體采用隨機模擬方法檢驗其是否滿足機會約束，即根據干擾兵力部署方案計算通信接收機是否被干擾壓制；

步驟4：計算所有染色體對應的目標函數值，即所有干擾兵力部署方案的干擾機需求量；

步驟5：以步驟4中求得的目標函數值為基礎，對違反約束的染色體采用懲罰值計算各染色體的適應度；對滿足約束條件的染色體，以步驟4中求得的目標函數值作為適應度；

步驟6：采用輪盤賭方法對種群中的染色體進行選擇操作；

步驟7：對種群中的染色體進行交叉和變異操作，得到新一代的染色體，之后采用隨機模擬方法檢驗是否滿足機會約束；

步驟8：重復步驟4～步驟7，達到給定的最大允許迭代次數為止；

步驟9：以求解過程中發現的最好的染色體作為最后的干擾兵力部署方案。

3.3 混合智能算法中遺傳算法的規則設計

3.3.1 解的編碼

目前，遺傳算法中主要采用兩種編碼方式：一種是基于二進制，另一種是基于實數的編碼方案。針對投擲式通信干擾兵力部署問題，引入基于干擾機坐標的實數編碼方案。干擾機設置方案的基因編碼映射如圖3所示。

設種群的個體數為T，每個個體的染色體由N個基因組構成，第i個基因組中的兩個基因分別表達第i個干擾機在x軸和y軸方向的坐標，若第i個基因組的表達式為（—1，—1），表示第i個干擾機不存在。

3.3.2 種群的初始化

在作戰任務區域中，預設干擾機數量上限為N，初始化種群個數為T，則初始種群為：

式（6）中，gij=（xij，yij）（1≤i≤T，1≤j≤N）代表第i個個體中第j個基因組所代表的干擾機坐標。

3.3.3 交叉操作和變異操作

為了使子代較好地保留父代的個體基因，本文采取一種循環交叉的方法實現交叉操作。設父代中選取K個個體進行繁殖：（P1，P2，…，PK），則第一次選取（P1，P2）作為雙親，第二次選取（P2，P3）作為雙親，第i次選取（Pi，Pi+1）作為雙親，以此類推共產生（K—1）個后代。

在交叉操作中，以交叉概率pc在父代中選取K個個體進行繁殖。設父代雙親為（Pi，Pi+1），如果雙親的等位基因組gik和g（i+1）k所代表的干擾機都存在，那么后代個體的等位基因組Chdik所代表干擾機的坐標為：

否則，子代個體的等位基因組Chdik所代表干擾機的坐標位置以概率p4設為（xPik，yPik），以概率1—p4設為（xP（i+1）k，yP（i+1）k），p4～U（0，1）。

變異操作能夠為種群提供新的內容。在所有后代個體構成的基因組集合中按概率pm選取基因組進行變異操作。如果原基因組所代表的干擾機不存在，則該基因組以概率p5變異為（v5，v6），以概率1—p5保持不變，即干擾機坐標為（—1，—1）；如果原基因組所代表的干擾機存在且坐標為（v1，v2），則該基因組以概率p6變異為干擾機不存在，以1—p6變異為（v1+ξ3，v2+ξ4）。其中ξ3，ξ4～N（0，σ2）。

3.3.4 評估和選擇操作

設父代種群個體數為T1，新繁殖的個體數為T2，則采樣空間的個體總數為T1+T2。設采樣空間中個體k目標函數值為fk，則目標函數適值集合為：（f1，f2，…，fT1，…，fT1+T2）。在此基礎上計算采用空間中每一個個體的生存概率，個體k的生存概率pk為：

由此生存概率構造出轉輪，新種群的產生首先通過轉輪法選出一對個體，然后將每對中的目標函數適值高（即方案較優）的個體加入新種群，直到新種群恢復原種群數目為止。

4 實例分析

4.1 作戰背景設置

落地式通信干擾機的作戰任務為對矩形區域內的無線戰術通信實施壓制干擾[10]，如圖4所示，矩形區域四個頂點的坐標為A（0Km，8Km），B（0Km，0Km），C（15Km，0Km），D（15Km，8Km），通信發射機的位置（xt，yt）已知，通信接收機的位置（xr，yr）在矩形區域中隨機變化，持續壓制時間為20分鐘，干擾頻帶為48.5～49.5MHz，壓制系數為2。

敵超短波通信電臺[11]的發射功率為10W，工作頻帶為48.998MHz～49.002MHz，全向天線，天線高度4m，采用語音調頻工作樣式，接收機帶寬為10KHz。

落地式通信干擾機的干擾功率為8W，持續干擾時間為50min，干擾頻帶為48～50MHz[10]。

為了研究在多種情況下落地式通信干擾機的部署結果，設定以下兩個方案進行仿真模擬：

仿真方案1：對指揮、后方、報知通信實施投擲式干擾，此時發射電臺在任務區域外部，如圖4所示。

仿真方案2：對協同通信實施投擲式干擾，此時發射電臺在任務區域內部，如圖6所示。

4.2 結果分析

在仿真方案1、2中，預設干擾機數目上限為500，交叉概率為0.9，種群個體數為50，變異概率為0.3，δ為400m，σ為400m，可變的因素包括通信發射機位置（xt，yt）和置信水平α。

在仿真方案1中，當通信發射機位置為（10Km，20Km），壓制概率的置信水平分別為1.0時，落地式通信干擾機的部署結果如圖4所示；當通信發射機的位置分別為（10Km，10Km）、（10Km，11Km）、（10Km，15Km）、（10Km，20Km）、（10Km，25Km），干擾壓制概率的置信水平為1.0時，采用隨機機會約束規劃和蜂窩狀陣形的干擾機需求量如圖5所示。

在仿真方案2中，當通信發射機位置為（3.75Km，6Km）時，對協同通信干擾壓制概率的置信水平為0.5時，落地式通信干擾機部署結果如圖6所示；當通信發射機位置分別為（0.5Km，7.5Km）、（3.75Km，6Km）、（7.5Km，4Km）時，對協同通信干擾壓制概率的置信水平依次為0.4、0.45、0.5、0.6、0.7、0.8、0.9、1時，落地式通信干擾機的需求量圖7所示。

通過對仿真結果的分析，可以得出以下結論：

（1）從圖4和圖6可以看出，基于隨機模擬、遺傳算法的混合智能算法能有效地求解投擲式通信干擾兵力部署的隨機機會約束規劃模型，由于該隨機機會約束規劃的目標為最小化干擾機需求量，因此，混合智能算法的解算結果為在偵察到通信接收機位置之前，滿足對通信接收機干擾壓制概率不小于一定置信水平的條件下，干擾機需求量最小的部署方案。

（2）從圖5可以看出，當對指揮、報知通信實施投擲式干擾時：①在完成相同干擾任務的條件下，隨機機會約束規劃的干擾機需求量較蜂窩狀陣形少；②隨著通信發射機與任務區域中最近點距離d的減小，隨機機會約束規劃中干擾機需求量雖然有所增加，但增加的速率較蜂窩狀陣形低，克服了蜂窩狀陣形中d較小時干擾機需求量急劇增加的不足；③在作戰任務區域內，干擾機的投擲密度隨著與發射機距離的增大而減小。這是因為與發射機距離較近時，接收機接收通信信號的功率較大，在滿足干擾壓制概率置信水平的條件下，干擾機需求量較大，因此投擲密度大；反之，投擲密度小。

（3）從圖7可以看出，當對協同通信實施投擲式干擾時：①在相同干擾壓制概率置信水平的條件下，干擾機需求量隨著發射機與任務區域中心位置距離的增加而減小；②干擾壓制概率置信水平在某個區間變化時，干擾機需求量保持不變；干擾壓制概率的置信水平越過某個數值時，干擾機需求量增加。這是因為在干擾機數量一定的條件下，遺傳算法通過優化干擾機的部署，使干擾壓制概率的置信水平在某個區間內變化；但遺傳算法的優化能力有限，在超出其優化能力的極限時，必須增加干擾機數量來提高干擾壓制概率的置信水平；③干擾壓制概率的置信水平為1時，干擾機需求量無窮大，即不存在一個干擾機部署方案對協同通信干擾壓制概率的置信水平為1。這是因為在距離通信發射機位置很近的區域內，通信信號的功率很大，不論干擾機數量多大、如何規劃干擾機的部署都不能使通信接收機輸入端的干信比大于壓制系數。

5 結束語

本文引入隨機機會約束規劃數學方法建立了投擲式干擾機壓制無線戰術通信的兵力部署模型，給出了基于隨機模擬和遺傳算法的求解方法，設計了遺傳算法中種群初始化、交叉、變異、評估和選擇的規則，實例分析表明，所建立的模型和求解算法可以恰當處理干擾無線戰術通信時接收機位置的不確定性，在約束條件的處理上也更加靈活，為實現對無線戰術通信實施投擲式干擾的輔助決策提供了參考。

1 邵國培，曹志耀，何俊，等.電子對抗作戰效能分析[M].北京：解放軍出版社，1998.

2 李新星，楊正，丁鯤.炮射投擲式通信干擾機作戰效能評估[J].電子對抗，2008，4（4）：35—40.

3 唐克，杜濤，王存威.無線電通信干擾彈作戰效能研究[J].兵工自動化，2009，7（7）：63—65.

4 李新星，趙偉，胡源.炮射投擲式通信干擾機的空間排布及風飄修正計算[J].通信對抗，2008，2（101）：44—48.

5 葛海龍，宋穎風.分布式通信干擾設備的應用及設計[J].無線電工程，2004，4（34）：29—31.

6 曾慶紅.分布式通信對抗發展設想[J].通信對抗，2008，2（2）：13—15.

7 楊寧，文福拴.基于機會約束規劃的輸電系統規劃方法[J].電力系統自動化，2004，14（28）：23—27.

8 劉寶碇，趙瑞清，王綱.不確定規劃及應用[M].北京：清華大學出版社，2006.

9 雷英杰，張善文，李續武.Matlab遺傳算法工具箱及應用[M].西安：電子科技大學出版社，2004.

10 熊群力，陳潤生，楊小牛.綜合電子戰[M].北京：國防工業出版社，2007.

11 陳太一.無線戰術電臺[M].北京：軍事科學出版社，2000.