原　慧, 王春陽, 安　磊, 李　欣

(1. 空軍工程大學(xué)防空反導(dǎo)學(xué)院, 陜西 西安 710051; 2. 中國人民解放軍94221部隊(duì), 山東 日照 276824)

0　引　言

目前,間歇采樣轉(zhuǎn)發(fā)干擾(interrupted-sampling repeater jamming, ISRJ)主要包括固定周期間歇采樣直接轉(zhuǎn)發(fā)干擾[1]、間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾[2]、參差周期間歇采樣干擾[3]、間歇采樣逐次循環(huán)轉(zhuǎn)發(fā)干擾[4]、間歇采樣非均勻重復(fù)轉(zhuǎn)發(fā)干擾[5]等,都是基于“存儲(chǔ)-轉(zhuǎn)發(fā)-存儲(chǔ)-轉(zhuǎn)發(fā)”干擾思想提出來的,以數(shù)字射頻存儲(chǔ)(digital radio frequency memory,DRFM)技術(shù)為實(shí)現(xiàn)基礎(chǔ),針對(duì)大時(shí)寬帶寬積線性調(diào)頻(linear frequency modulation, LFM)信號(hào)的新型靈巧干擾樣式。ISRJ的特點(diǎn)主要體現(xiàn)在天線系統(tǒng)和信號(hào)采樣存儲(chǔ)階段,即ISRJ既可用于收發(fā)系統(tǒng)隔離、同時(shí)工作的干擾設(shè)備,也可用于天線收發(fā)共用、分時(shí)工作的干擾設(shè)備,且無需對(duì)整段雷達(dá)脈沖信號(hào)進(jìn)行無失真的全采樣和存儲(chǔ)。因此,ISRJ可以克服對(duì)收發(fā)天線同時(shí)工作帶來的隔離度要求高、對(duì)寬帶雷達(dá)信號(hào)的高速采樣等難題,在彈載等干擾機(jī)平臺(tái)較小的應(yīng)用場(chǎng)合具有很大的實(shí)際應(yīng)用價(jià)值。

ISRJ與文獻(xiàn)[6]中提出的切片轉(zhuǎn)發(fā)(chopping and interleaving, C&I)干擾相比,具有一定的相似性,但也存在本質(zhì)區(qū)別:①對(duì)于收發(fā)共用天線,C&I干擾只能在下一個(gè)或幾個(gè)雷達(dá)脈沖重復(fù)周期內(nèi)完成干擾,而ISRJ可以在當(dāng)前雷達(dá)脈沖重復(fù)周期內(nèi)完成干擾;②C&I干擾信號(hào)在雷達(dá)的一個(gè)脈沖持續(xù)時(shí)間內(nèi)是連續(xù)的,而ISRJ是間斷的。由于ISRJ可以在當(dāng)前雷達(dá)脈沖重復(fù)周期內(nèi)完成干擾,因此,通過在脈沖重復(fù)周期間改變或調(diào)整發(fā)射信號(hào)或其參數(shù)的方法,例如相位擾動(dòng)[6]、調(diào)頻斜率調(diào)整[7]、脈沖分集[8-9]就失去了對(duì)抗效果。

理論上,可以通過在時(shí)頻面內(nèi)提取雷達(dá)信號(hào)和干擾的參數(shù),實(shí)現(xiàn)目標(biāo)參數(shù)提取和干擾抑制。但是,線性時(shí)頻變換方法受限于不確定性原理,時(shí)頻聚集性差,而非線性時(shí)頻變換又存在交叉項(xiàng)問題,并且通常當(dāng)干擾存在時(shí),干擾的功率都強(qiáng)于雷達(dá)信號(hào)的功率,這就導(dǎo)致雷達(dá)信號(hào)在時(shí)頻面被淹沒。因此,對(duì)于較小的采樣脈沖和較大的干擾信號(hào)比(jamming-to-signal ratio,JSR),通過在時(shí)頻面內(nèi)實(shí)現(xiàn)目標(biāo)參數(shù)提取和干擾抑制是十分困難的。此外,由于干擾與目標(biāo)回波信號(hào)的調(diào)頻斜率相同,所以利用分?jǐn)?shù)階傅里葉變換(fractional Fourier transform,FRFT)進(jìn)行掃頻濾波[10]的方法也失效。文獻(xiàn)[11]根據(jù)ISRJ信號(hào)的時(shí)域不連續(xù)性,以時(shí)頻分析的方法檢測(cè)只有回波信號(hào)存在的時(shí)間單元,根據(jù)回波單元所在的時(shí)頻特性構(gòu)建帶通濾波器,從而實(shí)現(xiàn)對(duì)干擾的濾除,但是其并沒有給出只有回波信號(hào)存在的時(shí)間單元長度的確定方法,并且在一定的干擾參數(shù)下,所構(gòu)建的帶通濾波器并不一定能將干擾完全濾除。

受此方法的啟發(fā),并結(jié)合壓縮感知(compressed sensing,CS)信號(hào)重構(gòu)思想,本文提出了一種新的ISRJ抑制方法。該方法同樣利用ISRJ信號(hào)的時(shí)域不連續(xù)性,但不需要進(jìn)行短時(shí)傅里葉變換等時(shí)頻分析,只需要對(duì)雷達(dá)接收信號(hào)進(jìn)行簡(jiǎn)單的能量函數(shù)運(yùn)算,然后通過與設(shè)定的閾值比較,即可得到未受干擾影響的目標(biāo)回波信號(hào)段;然后對(duì)提取到的未受干擾影響的目標(biāo)回波信號(hào)段進(jìn)行解線調(diào)處理,以保證解調(diào)后的信號(hào)在頻域是稀疏的(線性調(diào)頻信號(hào)經(jīng)解調(diào)頻后變?yōu)閱晤l信號(hào));最后,將提取到的未受干擾影響的目標(biāo)回波信號(hào)看作是目標(biāo)回波信號(hào)的壓縮數(shù)據(jù),根據(jù)其與目標(biāo)回波信號(hào)稀疏頻域之間的線性關(guān)系,建立CS最小問題求解模型并利用正交匹配追蹤算法(orthogonal matching pursuit,OMP)重構(gòu)目標(biāo)信號(hào),實(shí)現(xiàn)對(duì)ISRJ的抑制。

1　ISRJ信號(hào)特性分析

不失一般性,本文用間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的信號(hào)模型進(jìn)行具體分析。圖1給出了間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的原理示意圖及其與目標(biāo)回波信號(hào)的時(shí)頻分布示意圖,其中T和B分別為雷達(dá)發(fā)射LFM信號(hào)的脈寬和帶寬;τ為采樣切片的長度;Ts為切片采樣周期;M重復(fù)轉(zhuǎn)發(fā)當(dāng)前采樣信號(hào)的次數(shù)。

圖1　間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾原理示意圖及其與目標(biāo)回波信號(hào)的 時(shí)頻分布示意圖

由圖1可以看出,干擾信號(hào)和雷達(dá)信號(hào)具有不同的時(shí)頻變化特性,兩者在任意時(shí)刻的瞬時(shí)頻率都是不同的;此外ISRJ干擾是不連續(xù)的,即在干擾機(jī)的采樣存儲(chǔ)階段不發(fā)射干擾信號(hào)。

雷達(dá)發(fā)射的幅度歸一化LFM信號(hào)的數(shù)學(xué)模型可以表示為

(1)

(2)

ISRJ多用于突防環(huán)境下的自衛(wèi)式密集假目標(biāo)干擾,根據(jù)圖1,雷達(dá)接收到的ISRJ信號(hào)可表示為

(3)

式中,τd表示干擾機(jī)的系統(tǒng)轉(zhuǎn)發(fā)延遲。根據(jù)干擾方程可得干擾信號(hào)幅度Aj與干擾機(jī)的發(fā)射峰值功率、雷達(dá)天線主瓣增益、雷達(dá)發(fā)射信號(hào)波長以及干擾機(jī)與雷達(dá)的距離有關(guān)。

綜合式(2)和式(3)可得雷達(dá)接收信號(hào)為

z(t)=sr(t)+sj(t)+n(t)

(4)

將信號(hào)s(t)的能量函數(shù)定義為其模值的平方,在一個(gè)脈沖持續(xù)時(shí)間內(nèi)可以不考慮目標(biāo)起伏的影響,即可以認(rèn)為在一個(gè)脈沖持續(xù)時(shí)間內(nèi),目標(biāo)回波信號(hào)的幅度是恒定的。那么,目標(biāo)回波信號(hào)sr(t)、干擾信號(hào)sj(t)、以及不考慮噪聲時(shí)干擾條件下的雷達(dá)接收信號(hào)x(t)=sr(t)+sj(t)的能量函數(shù)分別為

(5)

Ej(t)=|sj(t)|2=

(6)

(7)

式中

由于|cos[2πk(mτ+τd)t+φ]|≤1,所以在干擾轉(zhuǎn)發(fā)時(shí)間段有(Aj-Ar)2≤Ex(t)≤(Aj+Ar)2。圖2給出了目標(biāo)回波信號(hào)、干擾信號(hào)以及不考慮噪聲時(shí)干擾條件下的雷達(dá)接收信號(hào)的能量函數(shù)隨時(shí)間的變化曲線。

圖2　能量函數(shù)隨時(shí)間的變化曲線

由圖2可以看出,在間歇采樣時(shí)間段,由于干擾機(jī)不發(fā)射干擾信號(hào),雷達(dá)接收到的只有目標(biāo)回波信號(hào),信號(hào)能量平穩(wěn)且較小,在干擾機(jī)發(fā)射干擾信號(hào)時(shí)間段,信號(hào)能量較大且起伏劇烈。

2　基于CS信號(hào)重構(gòu)的抗ISRJ方法

2.1　目標(biāo)回波信號(hào)壓縮數(shù)據(jù)的提取

γ=mean(envminp(Ex(t)))

(8)

式中,envminp (Ex(t))表示取能量函數(shù)的極小值包絡(luò);mean(·)表示取均值。

因?yàn)?/p>

(9)

所以有

(10)

當(dāng)噪聲不可以忽略時(shí),能量函數(shù)Ex(t)的包絡(luò)不再是恒定的,而是有一定的起伏,且當(dāng)干信比不是很大時(shí),再利用上述方法求取的閾值來提取未受干擾影響的目標(biāo)回波信號(hào)段時(shí),可能會(huì)存在很大誤差,即可能所得到的目標(biāo)回波壓縮數(shù)據(jù)并不“純正”,很可能還包含有干擾信號(hào)數(shù)據(jù)。為了最大限度地去除干擾信號(hào)的影響,而又不損失太多有用的目標(biāo)回波數(shù)據(jù),受經(jīng)驗(yàn)?zāi)B(tài)分解思想的啟發(fā),提出一種改進(jìn)的提取未受干擾影響的目標(biāo)回波信號(hào)段方法,具體步驟如下:

步驟1求雷達(dá)接收信號(hào)z(t)的能量函數(shù)Ez(t)的極大值包絡(luò)envmaxp(Ez(t))和極小值包絡(luò)envminp(Ez(t))。

步驟2求極大值包絡(luò)和極小值包絡(luò)的均值包絡(luò),得

(11)

顯然envmoyp(Ez(t))相對(duì)于能量函數(shù)Ez(t)更加平滑,具有一定的抑制噪聲影響能力。

步驟3取閾值為

γ=ρ×mean(envminp(Ez(t)))

(12)

式中,0<ρ≤1為閾值修正因子。因?yàn)楦蓴_信號(hào)間歇采樣的占空比小于等于50%且通常有Aj>2Ar,所以由式(8)得到的閾值偏大,加上噪聲的影響,以此為閾值得到的信號(hào)中可能仍包含干擾信號(hào)。理論和仿真分析可得:在一定的信噪比(signal-to-noise ratio,SNR)和占空比條件下或在一定的SNR和JSR條件下,隨著JSR或占空比的減小,ρ取逐漸減小的值可獲得較純凈的目標(biāo)回波信號(hào)數(shù)據(jù);在一定的SNR條件下,大占空比和大JSR對(duì)提取未受干擾影響的目標(biāo)回波信號(hào)數(shù)據(jù)都是有利的條件,此時(shí)可以取ρ=1。

(13)

2.2　基于CS的目標(biāo)回波信號(hào)重構(gòu)

(14)

(15)

圖3　基于CS信號(hào)重構(gòu)的抗ISRJ干擾方法流程圖

3　仿真與分析

3.1　對(duì)中帶信號(hào)的仿真分析

雷達(dá)發(fā)射帶寬B=10 MHz,脈沖寬度T=32 μs的LFM信號(hào),采樣頻率fs=20 MHz,則采樣點(diǎn)數(shù)N=640。假設(shè)間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的采樣周期Ts=4 μs,不考慮干擾機(jī)的系統(tǒng)轉(zhuǎn)發(fā)延遲,即τd=0。根據(jù)文獻(xiàn)[13]的分析和式(16),當(dāng)τ取不同值時(shí)可以形成以下4種主要的干擾效果:①Ts和τ滿足式(16)所示關(guān)系時(shí),所產(chǎn)生的假目標(biāo)都是獨(dú)立的;②Ts不滿足式(16),并且τ不為Δt(Δt=T/(TsB)=0.8 μs為各主假目標(biāo)間的分布間距)的整數(shù)倍時(shí),主假目標(biāo)群和非主假目標(biāo)群將可能發(fā)生混疊,使其中有的假目標(biāo)幅度呈現(xiàn)起伏變化;③Ts不滿足式(16),但τ為Δt的整數(shù)倍時(shí),由于相鄰幾次轉(zhuǎn)發(fā)的主次假目標(biāo)間相互重疊,假目標(biāo)幅度不再是起伏的;④當(dāng)Ts和τ不滿足式(16)時(shí),假目標(biāo)之間發(fā)生相互混疊而出現(xiàn)了多個(gè)分布假目標(biāo),而且主假目標(biāo)形成的分布假目標(biāo)幅度最高。

(16)

式中,T、Ts、τ的單位均為μs;B的單位為MHz。

因此,根據(jù)分析,選取τ=2 μs、1.2 μs、0.8 μs、0.1 μs分別對(duì)應(yīng)4種情況進(jìn)行仿真分析。仿真中,定義干擾抑制前,即輸入SNR和輸入JSR的計(jì)算公式分別為

(17)

表1　τ=2 μs時(shí)本文方法的成功抗干擾概率Panti

表2　τ=1.2 μs時(shí)本文方法的成功抗干擾概率Panti

表3　τ=0.8 μs時(shí)本文方法的成功抗干擾概率Panti

表4　τ=0.1 μs時(shí)本文方法的成功抗干擾概率Panti

此外,為了與文獻(xiàn)[14]提出的帶通濾波方法(簡(jiǎn)記為BF,本文所提方法簡(jiǎn)記為CS)進(jìn)行比較,定義干擾抑制后輸出信號(hào)的信干噪比(signal-to-jamming-and-noise ratio,SJNR)為

(18)

圖4～圖7分別給出了τ分別為2 μs、1.2 μs、0.8 μs、0.1 μs,對(duì)應(yīng)的JSR分別為5 dB、10 dB、15 dB、20 dB,ρ分別取0.9、0.7、0.5、1時(shí),干擾抑制前的相關(guān)處理結(jié)果、與BF方法相比的干擾抑制后的脈壓結(jié)果以及CS信號(hào)重構(gòu)效果。

圖4　τ=2 μs時(shí)的仿真結(jié)果(JSR=5 dB,ρ=0.9)

圖5　τ=1.2 μs時(shí)的仿真結(jié)果(JSR=10 dB,ρ=0.7)

圖6　τ=0.8 μs時(shí)的仿真結(jié)果(JSR=15 dB,ρ=0.5)

圖7　τ=0.1 μs時(shí)的仿真結(jié)果(JSR=20 dB,ρ=1)

由圖4～圖7可以看出,當(dāng)τ不是非常小時(shí),本文所提方法除了可以完全抑制干擾外,還基本可以相位不失真地重構(gòu)目標(biāo)回波信號(hào);但當(dāng)τ非常小(τ=0.1 μs)時(shí),雖然仍然可以完全抑制干擾,但可能存在一定的相位失真。而BF方法不能完全抑制干擾,因?yàn)閹V波器的帶寬不能無限小,且?guī)V波器的旁瓣可能非常高。

圖8給出了τ分別為2 μs、1.2 μs、0.8 μs、0.1 μs時(shí),BF方法和CS方法的干擾抑制性能曲線。其中,干擾抑制性能曲線的橫坐標(biāo)表示輸入信號(hào)的信干比(signal-to-jamming ratio of input,SJR),顯然,SJR=-JSR,縱坐標(biāo)為SJNR,當(dāng)輸入SJR相同時(shí),輸出SJNR越大,表明干擾的抑制效果越好。干擾抑制性能曲線都是蒙特卡羅仿真100次的平均結(jié)果。

圖8　中等帶寬LFM信號(hào)條件下,不同τ值條件下的干擾抑制性能曲線

3.2　對(duì)寬帶信號(hào)的仿真分析

寬帶LFM信號(hào)是雷達(dá)實(shí)現(xiàn)高分辨率成像的一種常用信號(hào)。對(duì)這種寬帶信號(hào),首先進(jìn)行去斜處理以顯著降低系統(tǒng)的中頻處理帶寬[14]。文獻(xiàn)[15]分析了ISRJ信號(hào)經(jīng)雷達(dá)去斜處理后形成假目標(biāo)的基本原理,研究了假目標(biāo)的時(shí)域、頻域、幅度、空間分布特性以及真假目標(biāo)的空間分布關(guān)系,并指出:和針對(duì)利用匹配濾波處理技術(shù)的雷達(dá)的干擾相比,兩者在形成假目標(biāo)的幅度特性及各階假目標(biāo)之間的距離方面沒有差別,從形成假目標(biāo)的空間分布來看,兩者存在較大差別,匹配濾波雷達(dá)中形成的主假目標(biāo)必然滯后于真目標(biāo)回波,而在去斜體制雷達(dá)中,可通過時(shí)延控制實(shí)現(xiàn)主假目標(biāo)的超前干擾,且有多組時(shí)延參數(shù)可供選擇。

此外,去斜處理不會(huì)影響能量函數(shù)的分布情況[16],因此也不會(huì)對(duì)未受干擾影響目標(biāo)回波信號(hào)段的提取產(chǎn)生影響。

可以看出,理論上本文所提方法對(duì)去斜體制雷達(dá)對(duì)抗ISRJ依然有效。下面進(jìn)行仿真驗(yàn)證。其中仿真參數(shù)如下:雷達(dá)發(fā)射LFM信號(hào)帶寬B=1 000 MHz,脈沖寬度為T=128 μs,去斜處理后的中頻帶寬為10 MHz。假設(shè)目標(biāo)為一單散射體目標(biāo),間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的采樣周期Ts=10 μs,根據(jù)文獻(xiàn)[16]的分析以及式(16),取3組τ和τd進(jìn)行仿真(對(duì)應(yīng)于式(16)的第4種情況要求τ<2 ns,這對(duì)干擾機(jī)要求較高,因此這里不考慮這一情況)。①τ=5 μs,τd=2.8 μs;②τ=2 μs,τd=1.7 μs;③τ=0.256 μs,τd=1.5 μs。

τ=5 μs、τ=2 μs和τ=0.256 μs對(duì)應(yīng)的干擾轉(zhuǎn)發(fā)占空比分別為50%,20%和2.56%,與第3.1節(jié)中的τ=2 μs、τ=0.8 μs和τ=0.1 μs對(duì)應(yīng)的干擾轉(zhuǎn)發(fā)占空比(分別為50%,20%和2.5%)相同或相近,因此寬帶條件下,以上3種τ值條件下,不同JSR和閾值修正因子ρ條件下,本文所提方法的成功抗干擾概率也類似于表1、表3和表4,在此不再重復(fù)給出。圖9給出了寬帶條件下,τ分別取5 μs、2 μs和0.256 μs時(shí),BF方法和CS方法的干擾抑制性能曲線。

圖9　寬帶LFM信號(hào)條件下,不同τ值的干擾抑制性能曲線

由圖9可以看出,對(duì)于寬帶LFM去斜體制雷達(dá),在適當(dāng)?shù)拈撝敌拚蜃訔l件下,CS方法的干擾抑制性能都基本能達(dá)到40 dB以上,遠(yuǎn)優(yōu)于BF方法,說明CS方法也適用于寬帶LFM去斜體制雷達(dá)對(duì)抗ISRJ。

4　結(jié)　論

本文利用間歇采樣轉(zhuǎn)發(fā)干擾信號(hào)時(shí)域的不連續(xù)特性,通過分析干擾與目標(biāo)回波信號(hào)的能量函數(shù)特征,提出了一種有效可行地提取未受干擾影響的目標(biāo)回波信號(hào)段的方法,經(jīng)過解線調(diào)處理后作為壓縮數(shù)據(jù),利用解線調(diào)后的目標(biāo)回波信號(hào)頻域上的稀疏性,構(gòu)建其之間的線性關(guān)系,得到CS最小問題求解模型,通過OMP算法實(shí)現(xiàn)了對(duì)目標(biāo)回波信號(hào)的重構(gòu)和對(duì)ISRJ干擾的抑制。仿真結(jié)果表明,通過設(shè)置合適的閾值,本文所提方法不僅適用于中帶LFM匹配濾波雷達(dá)對(duì)抗ISRJ,而且適用于寬帶LFM去斜體制雷達(dá)對(duì)抗ISRJ。而且,相比于文獻(xiàn)[14]的方法,本文所提方法不僅提出了有效可行的未受干擾影響的目標(biāo)回波信號(hào)段提取方法,而且干擾抑制性能也更好。

參考文獻(xiàn)：

[1] WANG X S, LIU J C, ZHANG W M, et al. Mathematics principles of interrupted-sampling repeater jamming[J]. Science in China Series F: Information Sciences, 2006, 50(1): 113-123.

[2] 劉忠,王雪松,劉建成,等.基于數(shù)字射頻存儲(chǔ)器的間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾[J].兵工學(xué)報(bào),2008, 29 (4): 405-410.

LIU Z, WANG X S, LIU J C, et al. Jamming technique of interrupted-sampling and periodic repeater based on digital radio frequency memory[J].Acta Armamentarii,2008,29(4):405-410.

[3] 周政,唐宏,張永順,等.基于時(shí)域采樣的靈巧噪聲干擾研究[J].現(xiàn)代雷達(dá),2010,32(5): 53-55.

ZHOU Z, TANG H, ZHANG Y S, et al. A study on smart noise jamming based on time domain sampling[J]. Modern Radar, 2010, 32(5): 53-55.

[4] 劉建成,王雪松,劉忠.對(duì)線性調(diào)頻脈壓雷達(dá)的導(dǎo)前假目標(biāo)群干擾[J].電子與信息學(xué)報(bào),2008, 30(6):1350-1353.

LIU J C, WANG X S, LIU Z. Preceded false target groups jamming against LFM pulse compression radars[J]. Journal of Electronic and Information Technology,2008,30(6):1350-1353.

[5] 張養(yǎng)瑞,李云杰,李曼玲,等.間歇采樣非均勻重復(fù)轉(zhuǎn)發(fā)實(shí)現(xiàn)多假目標(biāo)壓制干擾[J].電子學(xué)報(bào), 2016, 44(1): 46-53.

ZHANG Y R, LI Y J, LI M L, et al. Suppress jamming technique of multiple false targets on interrupted-sampling and non-uniform periodic repeater[J].Acta Electronica Sinica,2016,44(1):46-53.

[6] SOUMEKH M. SAR-ECCM using phase-perturbed LFM chirp signals and DRFM repeat jammer penalization[J]. IEEE Trans.on Aerospace and Electronic Systems, 2006, 42(1): 191-205.

[7] 馮祥芝,許小劍.隨機(jī)線性調(diào)頻斜率SAR抗欺騙干擾方法研究[J].系統(tǒng)工程與電子技術(shù),2009, 31(1): 69-73.

FENG X Z, XU X J. Study of countermeasures to deceptive jamming using random linear modulation frequency ratio SAR[J]. Systems Engineering and Electronics, 2009, 31(1): 69-73.

[8] AKHTAR J. An ECCM signaling approach for deep fading of jamming reflectors[C]∥Proc.of the IET International Conference on Radar Systems, 2007:1-5.

[9] SCHUERGER J, GARMATYUK D. Performance of random OFDM radar signals in deception jamming scenarios[C]∥Proc.of the IEEE Radar Conference, 2009:1-6.

[10] SUN M H, TANG B. Suppression of smeared spectrum ECM signal[J]. Journal of the Chinese Institute of Engineers, 2009, 32(3): 407-413.

[11] GONG S X, WEI X Z, LI X. ECCM scheme against interrupted sampling repeater jammer based on time-frequency analysis[J]. Journal of Systerms Engineering and Electronics, 2014, 25(6): 996-1003.

[12] TROPP J, GILBERT A. Signal recovery from random measurement via orthogonal matching pursuit[J]. Transactions on Information Theory, 2007, 53(12):4655-4666.

[13] 王雪松,肖順平,馮德軍,等.現(xiàn)代雷達(dá)電子戰(zhàn)系統(tǒng)建模與仿真[M]. 北京:電子工業(yè)出版社, 2010: 167-169.

WANG X S, XIAO S P, FENG D J, et al. Modeling and simulation of modern radar and electronic warfare systems[M]. Beijing: Publishing House of Electronics Industry, 2010: 167-169.

[14] WILLAM J C. Strech: a time-transformation technique[J].IEEE Trans.on Aerospace and Electronic Systems,1971,7(2): 269-278.

[15] 馮德軍, 陶華敏, 楊勇, 等.對(duì)去斜體制雷達(dá)的間歇采樣轉(zhuǎn)發(fā)干擾[J]. 中國科學(xué):信息科學(xué),2012,42(2): 184-193.

FENG D J, TAO H M, YANG Y, et al. Jamming de-chirping radar using interrupted-sampling repeater[J]. Science in China: Information Sciences, 2012, 42(2): 184-193.

[16] YUAN H, WANG C Y, AN L, et al. A method against interrupted-sampling repeater jamming based on energy function detection and band-pass filtering[J]. International Journal of Antennas and Propagation, 2017, 2017 (1): 1-9.