顧 兵

(中國船舶重工集團公司第七二三研究所，江蘇 揚州 225101)

0 引 言

現(xiàn)代戰(zhàn)爭中雷達面臨日趨復(fù)雜的電磁環(huán)境威脅，各種軍用有意干擾、民用無意干擾大量涌入雷達接收機，導(dǎo)致雷達探測性能下降[1]。尤其是多部雷達之間的同頻干擾[2]和射頻采樣存儲技術(shù)[3]的發(fā)展，采用全脈沖復(fù)制、切片復(fù)制、部分復(fù)制形式密集存儲轉(zhuǎn)發(fā)，有相參處理增益，干擾效率大大提高[4]。多相編碼雷達信號，其信號形式復(fù)雜，具備低截獲優(yōu)勢，同樣長度的多相編碼信號比二相編碼信號擁有更高的匹配濾波主副比[4]。

本文利用多相編碼信號復(fù)雜波形帶來的低截獲優(yōu)勢，正交波形的互相關(guān)低峰值特性，設(shè)計一種有效對抗轉(zhuǎn)發(fā)式欺騙干擾的低截獲雷達波形。多相編碼信號結(jié)構(gòu)形式復(fù)雜，能夠有效降低電子戰(zhàn)設(shè)備的截獲概率或者延長截獲時間，用遺傳算法優(yōu)化出多組正交多相編碼子脈沖信號，將這些子脈沖信號按照一定的形式拼接成雷達發(fā)射波形。雷達在處理時進行分段分時處理，當(dāng)前處理信號段大概率與轉(zhuǎn)發(fā)干擾信號正交，具備較好的抗轉(zhuǎn)發(fā)式欺騙干擾能力。

1 遺傳算法

遺傳算法是模擬自然界生物在繁殖過程中遺傳和進化過程而發(fā)展出的一種高效、并行、全局搜索的數(shù)學(xué)算法，在函數(shù)優(yōu)化領(lǐng)域有較為廣泛的應(yīng)用，尤其是一些多模型、非線性的優(yōu)化問題。

遺傳算法是一個典型的迭代優(yōu)化過程，其主要步驟有：

(1) 對參數(shù)集合域選擇適當(dāng)?shù)牟呗赃M行編碼，使其轉(zhuǎn)化為對應(yīng)的位串結(jié)構(gòu)空間；

(2) 設(shè)定衡量各個樣本個體優(yōu)劣尺度的適應(yīng)度函數(shù)；

(3) 確定遺傳策略及參數(shù)，包括種群大小、進化次數(shù)、交叉和變異概率、初始種群數(shù)、復(fù)制、交叉、變異方法等；

(4) 隨機產(chǎn)生規(guī)定量的初始種群；

(5) 計算出初始種群中所有個體的適應(yīng)度評價值；

(6) 根據(jù)確定的遺傳策略，對初始種群進行復(fù)制、交叉、變異等操作進而形成下一代種群；

(7) 判斷進化后的群體是否滿足判斷條件，或者已經(jīng)達到所設(shè)定的進化代數(shù)，若滿足則停止進化，若不滿足則重新進行進化并計算適應(yīng)度評價值。

其算法流程如圖1所示。

圖1 遺傳算法流程圖

2 正交多相碼

二相編碼信號具有理想的相關(guān)特性，但是由于可選擇的序列有限，使得二相編碼信號形式相對簡單，難以滿足雷達低截獲波形的設(shè)計要求，構(gòu)造出具有理想相關(guān)特性的多相編碼雷達信號，突破了二相編碼的局限性[5-6]。

設(shè)一個正交多相碼序列中有L個信號，其中每個信號中包含K個子脈沖，那么這個正交多相編碼信號模型如下：

{fl(x)=a(x)ejφl(k),k=1,2,…,K},l=1,2,…,L

(1)

(2)

式中：τ為子脈沖時寬；φl(k)(0≤φl(k)≤2π)為第l個信號中第k個子脈沖的相位。

如果信號采用N相編碼，則子脈沖的編碼相位從下面集合中選取：

φl(k)∈{φ1，φ2，…，φN}=

(3)

對于一個編碼信號個數(shù)為L，編碼長度為K的多相位編碼集合可用K×L相位矩陣來表示：

(4)

L個多相位編碼信號的自相關(guān)表達式如下：

P(φl,m)=

(5)

式中：l=1,2,…,L。

L個多相位編碼信號的互相關(guān)表達式如下：

Q(φa,φb,m)=

(6)

式中：a≠b,a,b=1,2,…,L。

3 波形設(shè)計

在雷達的時寬持續(xù)時間內(nèi)利用基于遺傳算法優(yōu)化出的正交多相編碼子脈沖進行拼接，設(shè)計出具備低截獲、抗轉(zhuǎn)發(fā)式干擾能力的雷達發(fā)射波形，其結(jié)構(gòu)框架設(shè)計示意圖如圖2所示。

圖2 波形結(jié)構(gòu)示意圖

波形的具體優(yōu)化方法如圖3所示。

(1) 根據(jù)雷達的實際需要確定雷達信號的整體參數(shù)、重頻和脈寬；

(2) 確定各子脈沖的數(shù)量和持續(xù)時間，子脈沖持續(xù)時間之和為雷達信號的時寬，持續(xù)時間可選取相同值，也可選取不同值，具體可結(jié)合預(yù)計的截獲時間和轉(zhuǎn)發(fā)時間進行確定；

(3) 根據(jù)需要的信號復(fù)雜度確定多相編碼的相位數(shù)量；

(4) 設(shè)置遺傳算法的初始種群數(shù)量、適應(yīng)度函數(shù)、交叉概率、變異概率、進化代數(shù)等參數(shù)；

(5) 利用遺傳算法對多個多相位編碼信號進行正交優(yōu)化，優(yōu)化選擇出1組正交多相位編碼序列；

(6) 利用優(yōu)化出的正交多相編碼序列生成相應(yīng)的信號；

(7) 對各子脈沖信號進行自相關(guān)處理和互相關(guān)處理，可進行加權(quán)處理，并提取自相關(guān)處理和自相關(guān)處理絕對峰值、主副比、各子脈沖隔離度等數(shù)據(jù)；

(8) 根據(jù)實際需要對生成的信號是否滿足抗干擾要求進行判決，若滿足要求則輸出信號梳理，若不滿足要求返回第2步重新進行優(yōu)化，直至優(yōu)化出所需要的滿足要求的波形并將信號數(shù)據(jù)輸出。

4 仿真分析

計算機Matlab仿真選取5個等長子脈沖，采用四相編碼，其相位值為0，π/4, π/2, 3π/4之一，自相關(guān)加權(quán)采用自適應(yīng)低副瓣優(yōu)化權(quán)值。

仿真參數(shù)設(shè)定：信號帶寬B=3 MHz，信號時寬T=50 μs，采樣率fs=4 MHz，信號波長為0.03 m，子脈沖數(shù)量為5，子脈沖碼長等長L=150，相位數(shù)量為4，交叉概率為0.7，變異概率為0.12，最大進化代數(shù)為500，初始種群數(shù)為80。

圖4 遺傳算法適應(yīng)度進化曲線

仿真結(jié)果：遺傳算法適應(yīng)度進化曲線如圖4所示。利用遺傳算法優(yōu)化出5個子脈沖的正交四相編碼相位變化如圖5所示。圖6～圖15為利用優(yōu)化出的5組正交四相編碼序列產(chǎn)生的信號的自相關(guān)結(jié)果和利用自適應(yīng)低副瓣加權(quán)進行適配濾波的結(jié)果。

表1給出了5個子脈沖的自相關(guān)處理和自適應(yīng)加權(quán)后處理后的自相關(guān)峰值、自相關(guān)主副比、失配后的主副比數(shù)據(jù)結(jié)果。表2給出了5個子脈沖相關(guān)處理后的峰值數(shù)據(jù)結(jié)果。表3給出了5個子脈沖相關(guān)處理后的主副比數(shù)據(jù)結(jié)果。表4給出了5個子脈沖的隔離度數(shù)據(jù)結(jié)果。

信號段1和5、2和5、3和5、4和5互相關(guān)仿真結(jié)果如圖16～圖19所示。

圖5 優(yōu)化的正交四相編碼相位變化

圖6 信號段1自相關(guān)處理結(jié)果

圖7 信號段1失配濾波處理結(jié)果

表1 自相關(guān)與失配濾波處理結(jié)果

表2 相關(guān)處理峰值數(shù)據(jù)結(jié)果

圖8 信號段2自相關(guān)處理結(jié)果

圖9 信號段2失配濾波處理結(jié)果

圖10 信號段3自相關(guān)處理結(jié)果

圖11 信號段3失配濾波處理結(jié)果

圖12 信號段4自相關(guān)處理結(jié)果

圖13 信號段4失配濾波處理結(jié)果

圖14 信號段5自相關(guān)處理結(jié)果

圖15 信號段5失配濾波處理結(jié)果

圖16 信號段1和5互相關(guān)結(jié)果

圖17 信號段2和段5互相關(guān)結(jié)果

圖18 信號段3和5互相關(guān)結(jié)果

圖19 信號段4和5互相關(guān)結(jié)果

由以上仿真數(shù)據(jù)結(jié)果可知，采用多子脈沖拼接的雷達波形，每個子脈沖自相關(guān)有良好的主副比性能，進行自適應(yīng)加權(quán)處理后主副比進一步降低。各子脈沖之間隔離度好，互相關(guān)無明顯峰值。

5 結(jié)束語

本文提出了一種有效抗轉(zhuǎn)發(fā)式干擾的雷達波形優(yōu)化設(shè)計方法，利用多相編碼的低截獲優(yōu)勢，降低雷達信號被截獲的概率或增加截獲時間。利用遺傳算法優(yōu)化出正交性能良好的子波形，電子戰(zhàn)設(shè)備截獲并轉(zhuǎn)發(fā)欺騙式干擾需要一定時間，雷達采用分段分時處理策略，轉(zhuǎn)發(fā)的波形有很大概率與待處理段波形正交，從而使雷達具備較好的抗轉(zhuǎn)發(fā)式干擾能力。

表3 相關(guān)處理主副比數(shù)據(jù)結(jié)果

表4 信號段1～5的隔離度