傅黎黎,彭巧樂,周逍宙,周 敏，黃文濤

(1. 上海無線電設(shè)備研究所, 上海，200090；2. 上海機電工程研究所，上海，201109)

0 引 言

信號分選作為星載電子偵察系統(tǒng)的關(guān)鍵環(huán)節(jié)，直接影響著電子偵察系統(tǒng)的性能，在現(xiàn)代高技術(shù)戰(zhàn)爭中起著至關(guān)重要的作用。常用的信號分選方法在星載電子偵察系統(tǒng)工作環(huán)境中暴露出各種缺點和不足[1]，因此有必要研究信號分選技術(shù)，以適應(yīng)復(fù)雜的信號形式和電磁環(huán)境。常用的信號分選技術(shù)需要對同一區(qū)域的信號進行長時間偵收，并且會對接收的大量信號全部進行信號分選以及輻射源參數(shù)統(tǒng)計，不適用于星載電子偵察系統(tǒng)中衛(wèi)星快速實時運行的情況[2]。

為解決該問題，本文以星載電子偵察系統(tǒng)為應(yīng)用對象，把基礎(chǔ)、核心和關(guān)鍵技術(shù)信號分選作為研究對象，探索和研究具有自主性、快速性的算法，以提高目前星載電子偵察系統(tǒng)的信號分選能力，提高信號分選快速性[3]。

1 算法原理

聚類分選是根據(jù)微波前端接收的信號脈沖流中雷達參數(shù)的特征相似程度，將信號劃分成不同類的過程，同一類中的信號有較大的相似性，不同類的信號之間有較大的相異性[4]。

對于測量接收到的雷達特征參數(shù)，我們稱之為脈沖描述字(Pulse Description Word，PDW)。假設(shè)接收到n個信號，{pdw1,pdw2,…,pdwn}為接收到的信號脈沖流的脈沖描述字集合，若用p個雷達參數(shù)表示每個信號的屬性，那么參數(shù)矩陣可以寫為

(1)

同樣的，可以用相異度矩陣表示這n個信號彼此之間的差異

(2)

其中，d(i,j)表示信號i和信號j之間的差異(不相似程度)，通常為一個非負數(shù)，數(shù)值越大，表示信號i和信號j越不相似。D一般通過歐式距離計算得到，因此也可以稱為距離矩陣。

根據(jù)距離矩陣，聚類分選可以將接收到的信號集合分類成m個子集C1,C2,…,Cm，并且滿足以下三個條件：

a)Ci≠? (i=1,2,…m)，即各個子集非空，不會分選出未出現(xiàn)信號的目標雷達；

c)Ci∩Cj=? (i,j=1,2,…,m;i≠j)，即任意兩個子集的交集為空，就是說一個信號不可能屬于兩個目標雷達，保證了分選的準確性。

為完成一個聚類分析任務(wù)，一般需要包含以下四項內(nèi)容：

a) 特征選擇：選擇合適的特征，盡可能多地包含任務(wù)相關(guān)信息；

b) 聚類算法：選擇合適的聚類算法對數(shù)據(jù)集進行聚類，揭示數(shù)據(jù)集的內(nèi)在結(jié)構(gòu)；

c) 聚類有效性：通過運行聚類算法得到聚類結(jié)果后，對結(jié)果進行有效性驗證；

d) 結(jié)果解釋：相關(guān)應(yīng)用領(lǐng)域的專家必須用其他實驗證據(jù)來分析和理解聚類結(jié)果，最后得出正確的結(jié)論[5]。

2 算法步驟

根據(jù)接收到的信號的參數(shù)矩陣PDW，通過初始參數(shù)設(shè)置、距離矩陣計算、聚類中心分裂及合并、迭代等步驟，可以完成聚類分選。詳細的聚類分選算法步驟如下：

a) 設(shè)置初始聚類參數(shù)

需要設(shè)置的初始聚類參數(shù)包括：最大迭代次數(shù)、初始聚類中心、合并判別閾值和分裂判別標準差。

b) 歸一化聚類中心

根據(jù)設(shè)置的c個初始聚類中心，形成聚類中心矩陣，其形式與PDW的形式一致。

c) 計算距離矩陣

計算每個樣本(接收到的n個信號)到各個聚類中心的歐氏距離，構(gòu)成距離矩陣D。

d) 樣本歸類

將所有樣本歸入距離最近的聚類中心，并修正各類中心的值：即將每類的所有樣本的平均值作為新的中心，形成新的聚類中心矩陣。

e) 計算平均距離

計算各類中諸樣本到聚類中心的平均距離以及所有類平均距離的總平均距離。

f) 分裂、合并及迭代的判別

根據(jù)上次迭代時聚類中心的數(shù)量和當前聚類中心的數(shù)量進行判別，若當前聚類中心數(shù)量小于上次聚類中心數(shù)量的一半，或者進行偶數(shù)次迭代后當前聚類中心數(shù)量小于兩倍的上次聚類中心數(shù)量時，進行合并處理，進入第g)步，否則進行分裂處理，進入第h)步。如果迭代次數(shù)已達到最大迭代次數(shù)，則執(zhí)行第i)步，運算結(jié)束。

g) 合并處理

計算全部聚類中心間的距離，若某兩個聚類中心的距離小于合并判別閾值，則進行合并處理，當前聚類中心數(shù)量減少1，并得出合并后的新的中心值。

h) 分裂處理

計算所有的聚類中心的標準差向量。找出所有類中標準差的最大值，當標準差最大值大于設(shè)定的分裂判別標準差時，則選取此類中距離最大的兩個樣本作為新的中心，當前聚類中心數(shù)量增加1，并且給出分裂后的新的中心值。

i) 如果是最后一次迭代，則結(jié)束循環(huán)。否則轉(zhuǎn)入第b)步，迭代次數(shù)加1。

根據(jù)如上的聚類分選步驟，可以給出聚類分選的流程圖，如圖1所示。

圖1 聚類分選流程圖Fig.1 Clustering and sorting flow chart

3 快速聚類分選

在星載電子偵察中，衛(wèi)星相對地球表面高速運動，對同一區(qū)域的信號偵收時間有限，并且星載設(shè)備需要滿足體積小、功耗低的要求，傳統(tǒng)的通用接收機加聚類分選算法的組合已經(jīng)不能滿足要求。為此，提出一種快速聚類分選算法，優(yōu)化和改進接收機及搜索策略，減少分選的迭代次數(shù)，實現(xiàn)滿足星載電子偵察的快速聚類分選。

對于傳統(tǒng)的通用搜索策略來說，會在偵察頻段內(nèi)設(shè)置多個窄帶引導(dǎo)接收機，并將幾部接收機中的一部專門用于搜索。搜索接收機的帶寬通常大于引導(dǎo)接收機的帶寬，而且以最大的實際速率進行掃描。它把被測信號的頻率和其他信息送入處理器，處理器將各個信號分派給引導(dǎo)接收機，引導(dǎo)接收機通過較長時間的脈沖積累，采用聚類分選算法對脈沖進行處理，提取所需的詳細信息，具體偵察分選工作流程如圖2所示。

圖2 常規(guī)星載分選偵察流程Fig.2 Conventional satellite-borne sorting reconnaissance flow

在小型化星載偵察應(yīng)用背景下，為減小體積重量，偵察接收機則采用窄帶頻率搜索策略，只設(shè)置窄帶搜索接收機，該接收機既實現(xiàn)頻率搜索的功能，又實現(xiàn)偵察和處理的功能。通過對聚類分選算法距離矩陣的改進，減少迭代步驟、加快聚類分選算法的收斂過程，從而降低分析算法對脈沖樣本量的需求，提高星載偵察的快速性。具體工作流程如圖3所示。

圖3 快速星載分選偵察流程Fig.3 Fast satellite-borne sorting reconnaissance flow

a) 距離矩陣的改進

加權(quán)矩陣的取值利用熵值確定，某一特征的熵值越小，說明其包含的確定信息就越多，對分選結(jié)果的影響就越大，應(yīng)當設(shè)置大權(quán)值，反之則應(yīng)設(shè)置小權(quán)值。加權(quán)矩陣取值的具體步驟如下：

(a) 確定特征參數(shù)的聚類分選區(qū)間，為每個聚類中心的特征雷達參數(shù)設(shè)置范圍區(qū)間。

(b) 計算權(quán)值

按熵值原理計算待分選脈沖每個特征參數(shù)的熵值H(pdwij)。由熵值的性質(zhì)可知，熵值H(pdwij)越小，各聚類中心在特征參數(shù)pdwij上的差別就越大，可分離性也就越大，在分選中所占的權(quán)重也應(yīng)當越大。綜合p個特征參數(shù)的熵值，就可以給出每個特征參數(shù)的權(quán)值，從而得到加權(quán)矩陣。

當然，在采用窄帶頻率搜索策略，簡化了接收機之后，實際應(yīng)用中會遇到輻射源頻率范圍大于接收機帶寬、接收到輻射源信號時間不定，信號積累時間不確定的問題。針對這些問題，快速聚類分選算法進一步優(yōu)化搜索策略，充分利用接收到的信號信息，并且在分選過程中，改進了距離矩陣的計算，減少了算法的迭代次數(shù)，實現(xiàn)快速聚類分選。

b) 輻射源頻率范圍大于接收機帶寬

在實際應(yīng)用中，一些雷達輻射源的頻率范圍可能大于星載電子偵察接收機的帶寬，這樣會出現(xiàn)同一輻射源的信號在多個搜索頻段中出現(xiàn)的情況，如圖4所示。

圖4 輻射源頻率范圍大于接收機帶寬Fig.4 Radiation source frequency range is larger than the receiver bandwidth

圖4中輻射源信號出現(xiàn)在中心頻率分別為f1、f2、f3的搜索頻段中，如果每個搜索頻段進行聚類分選是獨立的，這樣就誤把一個輻射源判斷為三個輻射源。為此，采用如圖5所示的分選策略，在進行若干個單頻段分選之后，將相鄰多個搜索頻段的聚類分選結(jié)果再進行一次聚類，從而避免了誤判的情況。

f0f1f2f3f4單頻段分選單頻段分選單頻段分選單頻段分選單頻段分選多頻段分選多頻段分選

圖5 多頻段分選

Fig.5 Multiband sorting

c) 接收輻射源信號時間不定

受衛(wèi)星的運動、頻率搜索中心頻率的變化以及輻射源開機時間不確定等因素的影響，電子偵察接收機開始接收到輻射源信號的時間也是不定的，信號積累時間也是不確定的。在某一搜索頻段輻射源信號積累時間足夠長時，可以分選出此輻射源；信號積累時間不夠，則不能分選出此輻射源信號。

圖6 輻射源信號出現(xiàn)在兩次頻率搜索中Fig.6 Radiation source signal exists in two frequency search process

對于圖6中所示的情況，輻射源信號在兩次頻率搜索中出現(xiàn)。輻射源信號分別出現(xiàn)在t3和t7時段，每次的積累時間分別為τ1和τ2。若單獨對這兩次頻率搜索中接收的輻射源信號進行聚類分選，由于積累時間τ1和τ2較短，不能夠分選出輻射源。為此，需要將兩次頻率搜索中接收的輻射源信號結(jié)合到一塊使用，也就是說將t7時段接受到的輻射源信號與t3時段接收到的該頻段的輻射源信號與合并使用，此時輻射源信號的積累時間可認為是τ1+τ2，足夠分選出輻射源。圖7給出了上述情況的分選策略示意，將同一頻段在不同時段接收的信號合并使用，避免了因積累時間過短而無法分選的情況。

t1t2t3t4t5t6t7t8t1頻段分選t2頻段分選t3頻段分選t4頻段分選t1+t5頻段分選t2+t6頻段分選t3+t7頻段分選t4+t8頻段分選

圖7 多次頻率搜索組合分選

Fig.7 Multiple frequency search combination sorting

4 性能仿真

為了評估快速星載分選偵察流程結(jié)合快速聚類分選算法與常規(guī)星載偵察流程與分析算法的性能，設(shè)定比較典型的觀測場景對算法的分選性能進行仿真。選擇輻射源信號的載頻頻率、脈沖寬度和到達角作為快速聚類分選的特征參數(shù)。假設(shè)快速聚類分選的載頻頻率范圍為1～10 GHz，脈沖寬度范圍為0.5～10 μs，到達角范圍為0°～10°。

首先設(shè)定場景一，產(chǎn)生三個輻射源的簡單脈沖信號，其特征參數(shù)分別為[1 GHz,1 μs,1°]、[1 GHz,1 μs,2°]、[1 GHz,1 μs,3°]。對接收到的信號進行聚類分選，屬于同一輻射源的信號用同一種圖標標明，不同的輻射源的信號用不同的圖標標明。當脈沖積累時間為100 ms時，快速分選算法和常規(guī)分選算法的分選結(jié)果如圖8所示。仿真結(jié)果可見，脈沖積累時間比較長的條件下，兩種算法都能夠正確分選。

(a) 脈沖積累時間100 ms時常規(guī)算法分選結(jié)果(a) Conventional sorting results of 100 ms pluse accumulation

(b) 脈沖積累時間100 ms時快速聚類分選結(jié)果(b) Fast clustering sorting results of 100ms pluse accumulation圖8 脈沖積累時間100 ms兩種分選算法效果Fig.8 Two sorting algorithm results of 100 ms pluse accumulation

在該場景下，將脈沖積累時間縮短至60 ms，分別使用普通算法及快速分選算法進行分選，結(jié)果如圖9所示。從圖中可以看到，常規(guī)分選算法分選結(jié)果出現(xiàn)錯誤，而快速分選算法在脈沖樣本減少的情況下，仍然能夠正確分選出三個輻射源。

再設(shè)定場景二，產(chǎn)生三個輻射源的簡單脈沖信號，其特征參數(shù)分別為[1 GHz,1 μs,1°]、[1 GHz,2 μs,2°]、[1 GHz,1 μs,3°]。同樣，對接收到的信號進行聚類分選，屬于同一輻射源的信號用同一種圖標標明，不同的輻射源的信號用不同的圖標標明。當脈沖積累時間為100 ms時，快速分選算法和常規(guī)分選算法的分選結(jié)果如圖10所示；當脈沖積累時間為60 ms時，快速分選算法和普通分選算法的分選結(jié)果如圖11所示。

(a) 脈沖積累時間60 ms時常規(guī)算法分選結(jié)果(a) Conventional sorting results of 60 ms pluse accumulation

(b) 脈沖積累時間60 ms時快速聚類分選結(jié)果(b) Fast clustering and sorting results of 60 ms pulse accumulation圖9 脈沖積累時間60 ms兩種分選算法效果Fig.9 Two sorting algorithm results of 60 ms pulse accumulation

(a) 脈沖積累時間100 ms時常規(guī)算法分選結(jié)果(a) Conventional sorting results of 100 ms pluse accumulation

(b) 脈沖積累時間100 ms時快速聚類分選結(jié)果(b) Fast clustering sorting results of 100 ms pluse accumulation圖10 脈沖積累時間100 ms兩種分選算法效果Fig.10 Two sorting algorithm results of 100 ms pluse accumulation

(a) 脈沖積累時間60 ms時常規(guī)算法分選結(jié)果(a) Conventional sorting results of 60 ms pluse accumulation

(b) 脈沖積累時間60 ms時快速聚類分選結(jié)果(b) Fast clustering sorting results of 60 ms pluse accumulation圖11 脈沖積累時間60 ms兩種分選算法效果Fig.11 Two sorting algorithm results of 60 ms pulse accumulation

從場景二的分選仿真結(jié)果可以看到，在脈沖積累時間為100 ms時兩種分選算法都能正確分選三個輻射源，而當脈沖積累時間減少至60 ms時，只有快速分選算法正確地完成了三個輻射源的分選。

5 結(jié)束語

本文針對星載電子偵察對體積重量有比較大的約束、對實時性要求比較高的應(yīng)用需求，在傳統(tǒng)常規(guī)偵察接收機和聚類信號分選算法的基礎(chǔ)上，對聚類分選算法進行改進，提出了一種快速聚類分選算法，使用加權(quán)歐式距離代替?zhèn)鹘y(tǒng)歐式距離來表示信號之間的差異，提高聚類算法的收斂速度，從而提高分選處理的快速性。在典型的目標態(tài)勢下，常規(guī)聚類分選算法需要積累脈沖時間為100 ms量級，才能夠?qū)崿F(xiàn)較好的分選效果。通過仿真分析，采用改進的快速分析算法，脈沖積累時間只需60 ms量級就能夠?qū)崿F(xiàn)通道的輻射源分離效果，大幅地改進提高了分選的快速性。通過快速算法的應(yīng)用，能夠?qū)π禽d偵察接收機進行簡化，使其更易應(yīng)用于工程實踐中。

[1] 吳技．國外航天偵察系統(tǒng)的現(xiàn)狀與發(fā)展[J]．電訊技術(shù),2009,49(5),108-114．

[2] 余渝生，王志誠.雷達點預(yù)案聚類算法性能比較與分析[J].制導(dǎo)與引信,2015,36(1):29-33.

[3] Mardia H K．New techniques for the de-interleaving of repetitive sequences. IEEPROCEEDINGS[J]．August 1989,Vol.136,No.4:149-154．

[4] Milojevic D J, Popovic B M．Improved algorithm for the de-interleaving of radar pulses[J]．IEEPROCEEDINGS-F．February 1992,Vol.139,No.1:98-104．

[5] Craig A. Hanna．The Associative Comparator: Adds New Capabilities to ESM Signal Processing Defense Electronics[J]．February.1984：51-54．

[6] Hanna C A, Cirucci J A, Prunty D R．VLSI Application to ESM Processing[J]．IBM Federal Systems Division, Owego．1987,Vol.13,No.1：32-36．