余康林,匡華星,王超宇

(中國船舶集團(tuán)第七二四研究所,南京 211153)

0 引 言

線性調(diào)頻信號(hào)是脈沖壓縮雷達(dá)常用的一種發(fā)射信號(hào)形式，因其具有較大的時(shí)寬帶寬積，使得雷達(dá)能夠同時(shí)擁有較優(yōu)的距離分辨率和速度分辨率。針對(duì)該類型雷達(dá)，Sparrow等[1]提出了多種密集假目標(biāo)干擾，如頻譜彌散(Smeared Spectrum,SMSP)干擾、間歇采樣轉(zhuǎn)發(fā)干擾，其本質(zhì)都是對(duì)干擾機(jī)截獲的雷達(dá)發(fā)射信號(hào)進(jìn)行調(diào)制轉(zhuǎn)發(fā)，因此與雷達(dá)發(fā)射信號(hào)具有較強(qiáng)相干性，脈壓處理后將在多個(gè)距離單元上形成假目標(biāo)，因其與真實(shí)目標(biāo)較為相似，增加了雷達(dá)目標(biāo)檢測的難度。根據(jù)Sparrow等對(duì)SMSP干擾不同參數(shù)下的干擾效果的分析，當(dāng)SMSP干擾的調(diào)頻率取為雷達(dá)發(fā)射信號(hào)的5～7倍時(shí)，干擾效果最優(yōu)。對(duì)于SMSP干擾，考慮自衛(wèi)式干擾場景下，干擾與真實(shí)目標(biāo)回波混疊進(jìn)入雷達(dá)接收機(jī)，由于干擾與真實(shí)目標(biāo)回波時(shí)、頻、空域上的高度混疊，傳統(tǒng)的波束形成[2]、盲源分離[3]等抗干擾措施將會(huì)在干擾抑制的同時(shí)造成嚴(yán)重的信號(hào)功率損失，因此針對(duì)該干擾需要研究相應(yīng)的抗干擾措施。

目前，針對(duì)SMSP干擾的抗干擾措施，主要包括干擾重構(gòu)對(duì)消和干擾剔除。文獻(xiàn)[4-5]根據(jù)SMSP干擾與雷達(dá)發(fā)射信號(hào)的關(guān)系，通過估計(jì)SMSP干擾子脈沖個(gè)數(shù)、時(shí)延、幅度和初始相位重構(gòu)干擾并對(duì)消抑制SMSP干擾，但涉及對(duì)時(shí)延、幅度和初始相位的遍歷估計(jì)，并存在較大誤差。文獻(xiàn)[6]基于壓縮感知理論，通過構(gòu)建干擾基字典及調(diào)頻率匹配，能夠自適應(yīng)重構(gòu)干擾，并通過對(duì)消抑制SMSP干擾，但需較大的計(jì)算量。文獻(xiàn)[7]提出了基于干擾剔除的干擾抑制方法，由于SMSP干擾與真實(shí)目標(biāo)回波的時(shí)頻分布差異，通過剔除SMSP干擾時(shí)頻點(diǎn)抑制SMSP干擾，但該方法會(huì)對(duì)目標(biāo)回波功率造成一定損失。文獻(xiàn)[8]提出了聯(lián)合時(shí)頻重排和雙正交傅里葉變換的SMSP干擾抑制方法，利用Radon變換估計(jì)SMSP干擾子脈沖個(gè)數(shù)，并構(gòu)建參考信號(hào)將SMSP干擾調(diào)制成一個(gè)線性調(diào)頻信號(hào)，通過峰值濾波并解調(diào)重構(gòu)SMSP干擾，最后通過對(duì)消抑制干擾，但是該方法計(jì)算量較大。

本文在現(xiàn)有研究基礎(chǔ)上，根據(jù)干擾重構(gòu)對(duì)消的思想提出了干擾重構(gòu)的另一思路。仿真結(jié)果表明，所提方法計(jì)算量較小，能夠獲得較高的信干噪比增益，在干擾抑制后能夠有效檢測到目標(biāo)。

1 SMSP干擾模型

自衛(wèi)式干擾場景下，雷達(dá)回波模型可假設(shè)為

z(t)=s(t)+j(t)+n(t) 。

(1)

式中：t為時(shí)間變量，s(t)、j(t)、n(t)分別為真實(shí)目標(biāo)回波、干擾和噪聲。假設(shè)雷達(dá)發(fā)射信號(hào)為線性調(diào)頻(Linear Frequency Modulation,LFM)信號(hào)，不考慮載頻，則真實(shí)目標(biāo)回波模型如下：

(2)

根據(jù)SMSP干擾產(chǎn)生原理[4]，其數(shù)學(xué)模型如下：

(3)

式中：Aj、Tj、kj、t1、φj分別為SMSP干擾的幅度、脈寬、調(diào)頻率、時(shí)延和初始相位，Tj、kj與T、k對(duì)應(yīng)關(guān)系如下：

(4)

由式(4)可得，SMSP干擾包含m個(gè)子脈沖，每個(gè)子脈沖的脈寬為雷達(dá)發(fā)射信號(hào)脈寬的1/m，調(diào)頻率為雷達(dá)發(fā)射信號(hào)調(diào)頻率的m倍，m一般取大于1的整數(shù)。

2 SMSP干擾重構(gòu)

由上述對(duì)SMSP干擾的分析可知，SMSP干擾與真實(shí)目標(biāo)回波調(diào)頻率存在較大差異，并且干擾功率一般遠(yuǎn)大于真實(shí)目標(biāo)回波功率，因此較容易辨識(shí)出雷達(dá)回波中是否包含該干擾，若干擾與真實(shí)目標(biāo)回波不存在混疊，則很容易將該干擾剔除，無法對(duì)雷達(dá)造成有效的干擾，因此SMSP干擾一般與真實(shí)目標(biāo)回波存在一定的混疊。

圖1 基于干擾重構(gòu)的SMSP干擾抑制方法流程圖

2.1 SMSP干擾調(diào)頻率估計(jì)

分?jǐn)?shù)階傅里葉變換(Fractional Fourier Transform,FRFT)常用于估計(jì)LFM信號(hào)調(diào)頻率。對(duì)于信號(hào)x(t)，其p階分?jǐn)?shù)階傅里葉變換定義為

(5)

式中：Kp(t,u)定義為

(6)

分?jǐn)?shù)階傅里葉變換對(duì)LFM信號(hào)具有能量聚集性，當(dāng)LFM信號(hào)進(jìn)行旋轉(zhuǎn)角度為α=arccot(-k)的FRFT時(shí)，能量聚集效果最優(yōu)，k為其調(diào)頻率。Renyi熵常用于描述信號(hào)的能量聚集性[9]:信號(hào)的能量聚集性越高，Renyi熵越小;反之，Renyi熵越大。在對(duì)SMSP干擾進(jìn)行最優(yōu)旋轉(zhuǎn)角度搜索時(shí)，可以用Renyi熵評(píng)價(jià)SMSP干擾不同旋轉(zhuǎn)角度下的能量聚集性。由于回波中干擾功率遠(yuǎn)大于真實(shí)目標(biāo)回波功率，因此最小Renyi熵對(duì)應(yīng)的旋轉(zhuǎn)角度即是SMSP干擾對(duì)應(yīng)的最優(yōu)旋轉(zhuǎn)角度，再根據(jù)旋轉(zhuǎn)角度與調(diào)頻率之間的關(guān)系，即可得到SMSP干擾的調(diào)頻率。設(shè)雷達(dá)回波z1,z2,…,zn離散分?jǐn)?shù)階傅里葉變換歸一化幅度譜為X1,X2,…,XN，使其幅度之和為1，則Renyi熵定義為

(7)

式中：q為Renyi熵的階次，一般取正整數(shù)，這里q取4。

由于SMSP干擾的調(diào)頻率是雷達(dá)發(fā)射信號(hào)的m倍，因此其對(duì)應(yīng)最優(yōu)旋轉(zhuǎn)角度一般位于雷達(dá)發(fā)射信號(hào)對(duì)應(yīng)最優(yōu)旋轉(zhuǎn)角度的一側(cè)，在雷達(dá)發(fā)射信號(hào)旋轉(zhuǎn)角度已知的情況下，可以先以一個(gè)較大的角度搜索步長遍歷正半?yún)^(qū)間或負(fù)半?yún)^(qū)間進(jìn)行粗估計(jì)，再以較小角度搜索步長進(jìn)行精估計(jì)。

SMSP干擾調(diào)頻率估計(jì)步驟如下：

(8)

調(diào)頻率估計(jì)相對(duì)誤差定義為

(9)

在估計(jì)得到SMSP干擾的調(diào)頻率后，根據(jù)SMSP干擾調(diào)頻率與雷達(dá)發(fā)射信號(hào)調(diào)頻率對(duì)應(yīng)關(guān)系，SMSP干擾的子脈沖個(gè)數(shù)可由下式計(jì)算得到：

(10)

2.2 SMSP干擾時(shí)延估計(jì)

在估計(jì)SMSP干擾調(diào)頻率和子脈沖個(gè)數(shù)的基礎(chǔ)上，進(jìn)一步估計(jì)該干擾對(duì)應(yīng)的時(shí)延。構(gòu)造參考信號(hào)h(t),其表達(dá)式如下：

(11)

exp{i[πkj(t-t1-lTj-τ)2+φj]}dτ。

(12)

可得匹配濾波器輸出幅度譜為

sinc[πkj(t-t1+(i-l)Tj)(Tj-|t-t1+(i-l)Tj|)]|,

|t-t1+(i-l)Tj|

(13)

對(duì)于雷達(dá)接收機(jī)采樣得到的離散數(shù)據(jù)，可設(shè)為

(14)

(15)

式中：Nj=N/m。

2.3 SMSP干擾幅度、初始相位估計(jì)

(16)

TAjeiφje-i2πft1sinc(πTf)。

(17)

當(dāng)f=0時(shí)，上式取得最大值，幅度為TAj，相位為φj。由于系數(shù)T已知，則SMSP干擾的幅度和初始相位可由下式計(jì)算得到：

(18)

因此在SMSP干擾調(diào)頻率、子脈沖個(gè)數(shù)和時(shí)延估計(jì)誤差較小時(shí)，通過上述過程可以估計(jì)得到SMSP干擾的幅度和初始相位，在估計(jì)得到SMSP干擾的參數(shù)后就能重構(gòu)該干擾，再利用對(duì)消抑制該干擾。重構(gòu)的SMSP干擾表達(dá)式如下：

(19)

(20)

|Xp(u)|=|sinc[πT(ucscα+kt0)]|，

(21)

其輸出表現(xiàn)為一個(gè)sinc函數(shù)，主瓣寬度為2|sinα|/T，此時(shí)分?jǐn)?shù)階傅里葉變換對(duì)真實(shí)目標(biāo)回波具有最佳的能量聚集性，真實(shí)目標(biāo)回波能量聚集在較窄的范圍內(nèi)，通過選擇合適的濾波器，保留真實(shí)目標(biāo)回波的主瓣能量，能夠有效抑制殘余干擾和噪聲。

3 仿真與分析

不失一般性，雷達(dá)發(fā)射信號(hào)的脈寬設(shè)置為10 μs，帶寬設(shè)置為40 MHz，調(diào)頻率為4 MHz/μs，真實(shí)目標(biāo)回波及SMSP干擾參數(shù)設(shè)置如表1所示，干信比(Jamming-to-Signal Ratio,JSR)設(shè)置為10～30 dB，信噪比(Signal-to-Noise Ratio,SNR)設(shè)置為-10～20 dB，每個(gè)信噪比下進(jìn)行10 000次蒙特卡洛實(shí)驗(yàn)。

表1 仿真參數(shù)

3.1 SMSP的干擾效果

令干信比分別為10 dB、20 dB，信噪比為0 dB，真實(shí)目標(biāo)回波與SMSP干擾參數(shù)設(shè)置如表1所示，不同干信比下雷達(dá)回波的脈壓歸一化功率譜及單元平均恒虛警率(Constant False Alarm Rate,CFAR)檢測門限如圖2所示。CFAR門限檢測使用參數(shù)設(shè)置如下：參考單元數(shù)為60，保護(hù)單元數(shù)為8，虛警概率為10-6。由圖2可知，當(dāng)JSR為10 dB時(shí)，脈壓處理后干擾不能壓制住真實(shí)目標(biāo)，經(jīng)CFAR門限檢測仍能有效檢測到目標(biāo)，此時(shí)無法有效干擾雷達(dá)；當(dāng)JSR增加到20 dB時(shí)，脈壓處理后真實(shí)目標(biāo)被干擾完全壓制，經(jīng)CFAR門限檢測無法有效檢測到真實(shí)目標(biāo)。因此對(duì)于SMSP干擾，需要JSR很大時(shí)才能有效干擾雷達(dá)。

(a)JSR為10 dB時(shí)脈壓輸出

3.2 SMSP干擾不同參數(shù)的估計(jì)誤差

令信噪比為-10～20 dB，干信比為10～30 dB，SMSP干擾和真實(shí)目標(biāo)回波參數(shù)不變，分析上述參數(shù)估計(jì)方法估計(jì)SMSP干擾參數(shù)的誤差，并對(duì)比干信比為30 dB時(shí)文獻(xiàn)[5]的參數(shù)估計(jì)誤差。SMSP干擾調(diào)頻率估計(jì)時(shí)，粗估計(jì)搜索步長設(shè)為0.5°，i設(shè)為20，搜索精度設(shè)為0.000 5°。

SMSP干擾的調(diào)頻率估計(jì)誤差如圖3所示。由圖可知，隨著信噪比的升高，SMSP干擾的調(diào)頻率估計(jì)誤差逐漸減小，這表明信噪比越高，噪聲對(duì)SMSP干擾調(diào)頻率估計(jì)造成的影響越?。恍旁氡纫欢〞r(shí)，隨著干信比的升高，SMSP干擾的調(diào)頻率估計(jì)誤差越低，這表明當(dāng)SMSP干擾功率在雷達(dá)回波中占比越高，真實(shí)目標(biāo)回波功率對(duì)SMSP干擾調(diào)頻率估計(jì)造成的影響越小。在估計(jì)得到SMSP干擾的調(diào)頻率后，通過四舍五入取整，在上述仿真參數(shù)下均能準(zhǔn)確估計(jì)得到SMSP干擾的子脈沖個(gè)數(shù)。文獻(xiàn)[5]通過時(shí)延自相關(guān)估計(jì)SMSP干擾的子脈沖個(gè)數(shù)，干信比為30 dB時(shí)能夠準(zhǔn)確估計(jì)得到SMSP干擾的子脈沖個(gè)數(shù)及調(diào)頻率。

圖3 SMSP干擾調(diào)頻率估計(jì)相對(duì)誤差

圖4 SMSP干擾時(shí)延估計(jì)誤差

圖5 SMSP干擾幅度估計(jì)誤差

圖6 SMSP干擾初始相位估計(jì)誤差

3.3 本文方法的抑制效果

令干信比為20 dB，信噪比為0 dB，采用上述干擾抑制方法后，脈壓歸一化功率譜如圖7(a)所示，經(jīng)CFAR門限檢測能準(zhǔn)確檢測到真實(shí)目標(biāo)的位置。不同干信比下，信干比增益隨信噪比變化情況如圖7(b)所示。信干比(Signal-to-jamming Ratio,SJR)增益定義為G=JSR1-JSR2，JSR1為干擾抑制前干信比，JSR2為干擾抑制后干信比。由圖可知，干信比一定時(shí)，信干比增益隨信噪比的提高逐漸提高，這是因?yàn)镾MSP干擾的參數(shù)估計(jì)誤差隨信噪比的提高逐漸減??；信噪比一定時(shí)，信干比增益隨干信比的提高逐漸提高，這表明當(dāng)干擾功率所占比重上升時(shí)，干擾重構(gòu)對(duì)消抑制掉的干擾功率增加；并且上述干擾抑制方法在較低干信比時(shí)也能達(dá)到29 dB以上的信干比增益。干信比為30 dB時(shí)，本文干擾抑制方法與文獻(xiàn)[5]算法的信干比增益如圖7(b)所示，從圖中可以看出相同仿真條件下本文干擾抑制方法干擾抑制效果優(yōu)于文獻(xiàn)[5]算法。

(a)JSR為20 dB時(shí)干擾抑制后脈壓輸出

(a)干擾抑制后回波FRFT幅度譜

4 結(jié) 論

本文提出了一種基于干擾重構(gòu)和分?jǐn)?shù)階濾波的SMSP干擾抑制方法，詳述了SMSP干擾參數(shù)估計(jì)流程，并通過仿真實(shí)驗(yàn)驗(yàn)證了該方法的可行性。仿真結(jié)果表明，本文所提干擾參數(shù)估計(jì)方法具有較好的魯棒性，在低信噪比下也能夠精確估計(jì)得到SMSP干擾的參數(shù)，并且信干比增益仍能達(dá)到29 dB以上；在干擾重構(gòu)對(duì)消的基礎(chǔ)上，可以采用分?jǐn)?shù)階濾波進(jìn)一步抑制殘余干擾和噪聲，以提高雷達(dá)回波信干噪比。