基于STFrFT的間歇采樣轉(zhuǎn)發(fā)干擾抑制

摘 要：干擾機(jī)對(duì)雷達(dá)脈沖快速切片、轉(zhuǎn)發(fā)形成間歇采樣轉(zhuǎn)發(fā)干擾，若從雷達(dá)主瓣進(jìn)入，將對(duì)雷達(dá)目標(biāo)檢測(cè)形成嚴(yán)重威脅。從波形設(shè)計(jì)和時(shí)頻分析的角度出發(fā)，提出一種基于短時(shí)分?jǐn)?shù)階傅里葉變換（short-time fractional Fourier transform， STFrFT）的主瓣間歇采樣轉(zhuǎn)發(fā)干擾抑制方法。首先設(shè)計(jì)脈內(nèi)捷變頻信號(hào)，提升干擾與目標(biāo)信號(hào)的差異;接著采用STFrFT進(jìn)行時(shí)頻分析，相比傳統(tǒng)短時(shí)傅里葉變換大大提升了時(shí)頻分辨率，相比分?jǐn)?shù)階傅里葉變換類方法沒(méi)有對(duì)信號(hào)參數(shù)的限制;最后結(jié)合圖像學(xué)方法對(duì)干擾進(jìn)行剔除，對(duì)從主瓣進(jìn)入的高干信比干擾也可形成有效抑制。仿真結(jié)果表明，所提方法可在多種環(huán)境下有效對(duì)抗間歇采樣轉(zhuǎn)發(fā)干擾。

關(guān)鍵詞： 間歇采樣轉(zhuǎn)發(fā)干擾; 短時(shí)分?jǐn)?shù)階傅里葉變換; 圖像處理; 頻率捷變

中圖分類號(hào)： TN 974 文獻(xiàn)標(biāo)志碼： A""" DOI：10.12305/j.issn.1001-506X.2024.10.09

Interrupted-sampling repeater jamming suppression based on STFrFT

LI Jinjie CAO Yunhe ZHANG Yulin WANG Meng²

（1. National Laboratory of Radar Signal Processing， Xidian University， Xi’an 71007 China;

2. Beijing Electro-Mechanical Engineering Institute， Beijing 100074， China）

Abstract： Jammers can rapidly slice and forward radar pulse to form interrupted-sampling repeater jamming. If it enters from the main lobe of radar， it will pose a serious threat to radar target detection. Starting from waveform design and time-frequency analysis， a main lobe interrupted-sampling repeater jamming suppression method is proposed based on short-time fractional Fourier transform （STFrFT）. Firstly， an intra-pulse frequency agility signal is harnessed to improve the distinguishability between jamming and target signal.Moreover， the STFrFT is applied to time-frequency analysis， which significantly improves the time-frequency analysis accuracy compared with traditional short-time Fourier transform methods and removes the restrictions on signal parameters compared with fractional Fourier transform methods. Finally， through eliminating jamming by using image processing methods， the interference of high jamming-to-signal ratio from the main lobe of radar is suppressed effectively. The simulation results indicate that the proposed method can effectively suppress interrupted-sampling repeater jamming in various environments.

Keywords： interrupted-sampling repeater jamming; short-time fractional Fourier transform （STFrFT）; image processing; frequency agility

0 引 言

隨著電子技術(shù)的發(fā)展，數(shù)字射頻存儲(chǔ)器（digital radio frequency memory， DRFM）在電子對(duì)抗中的應(yīng)用愈加廣泛^［13］，基于DRFM的間歇采樣轉(zhuǎn)發(fā)干擾^［47］對(duì)雷達(dá)信號(hào)進(jìn)行多次采樣并轉(zhuǎn)發(fā)，形成與雷達(dá)信號(hào)十分接近的相干假目標(biāo)串，在雷達(dá)接收端形成大量假目標(biāo)。由于只需要對(duì)雷達(dá)脈沖進(jìn)行部分采樣和轉(zhuǎn)發(fā)^［8］，DRFM具有較高的時(shí)效性，又由于干擾脈沖與雷達(dá)脈沖相參，可以獲得雷達(dá)脈壓增益，因此DRFM對(duì)干擾機(jī)功率要求較低，從而大大增加了雷達(dá)抗干擾的難度。

常規(guī)的雷達(dá)抗干擾方法諸如脈間頻率捷變^［910］、脈間射頻掩護(hù)^［11］、波形設(shè)計(jì)^［1215］等對(duì)間歇采樣轉(zhuǎn)發(fā)干擾難以有效抑制，而當(dāng)干擾在雷達(dá)波束的主瓣位置時(shí)，旁瓣對(duì)消等空域?yàn)V波方法也失去效果。針對(duì)間歇采樣轉(zhuǎn)發(fā)干擾，文獻(xiàn)［16］ 以短時(shí)傅里葉變換（short-time Fourier transform， STFT）進(jìn)行時(shí)頻分析，采用脈內(nèi)步進(jìn)頻線性調(diào)頻（linear frequency modulation， LFM）信號(hào)，在時(shí)頻域進(jìn)行干擾抑制;文獻(xiàn)［17］設(shè)計(jì)脈內(nèi)頻率捷變信號(hào)，對(duì)分段脈壓結(jié)果進(jìn)行分析后剔除被干擾子脈沖，但脈壓損失較大。文獻(xiàn)［18］在脈內(nèi)隨機(jī)頻率編碼的基礎(chǔ)上，利用分?jǐn)?shù)階傅里葉變換（fractional Fourier transform， FrFT）在分?jǐn)?shù)階域剔除干擾，但該方法需準(zhǔn)確估計(jì)干擾參數(shù)，實(shí)現(xiàn)難度大。在此基礎(chǔ)上，文獻(xiàn)［19］設(shè)計(jì)了脈內(nèi)雙參數(shù)頻率捷變信號(hào)，使調(diào)頻斜率正負(fù)隨機(jī)切換，進(jìn)一步增強(qiáng)了目標(biāo)信號(hào)與干擾信號(hào)的差異，但干擾參數(shù)預(yù)估的問(wèn)題仍沒(méi)有解決。文獻(xiàn)［20］提出一種應(yīng)用于寬帶多輸入多輸出（multiple-input multiple-output， MIMO）雷達(dá)的波形設(shè)計(jì)方案以對(duì)抗間歇采樣轉(zhuǎn)發(fā)干擾，但該方法會(huì)導(dǎo)致信噪比（signal-to-noise ratio， SNR）的損失。在時(shí)頻分析方面，常見的思路有信號(hào)重構(gòu)和構(gòu)造時(shí)頻濾波器?；谛盘?hào)重構(gòu)思想，文獻(xiàn)［21］通過(guò)估計(jì)干擾信號(hào)采樣、轉(zhuǎn)發(fā)次數(shù)及采樣寬度等參數(shù)，重構(gòu)干擾信號(hào)并與回波信號(hào)相消，從而抑制干擾，但這種方法對(duì)參數(shù)估計(jì)精確度要求高，計(jì)算也比較復(fù)雜;文獻(xiàn)［22］利用時(shí)頻分析估計(jì)干擾信號(hào)參數(shù)，針對(duì)干擾設(shè)計(jì)脈內(nèi)頻率編碼信號(hào)，設(shè)計(jì)濾波器對(duì)脈壓結(jié)果進(jìn)行濾波，但該方法對(duì)多目標(biāo)或多干擾的場(chǎng)景無(wú)效?；跁r(shí)頻域帶通濾波器濾除干擾的思想，文獻(xiàn)［23］基于間歇采樣轉(zhuǎn)發(fā)干擾的不連續(xù)性構(gòu)造帶通濾波器在脈沖壓縮域?qū)Ω蓴_進(jìn)行濾除;文獻(xiàn)［24］通過(guò)構(gòu)造最大時(shí)頻函數(shù)反映最大信號(hào)分量振幅隨時(shí)間的變化情況，進(jìn)而自動(dòng)提取無(wú)干擾信號(hào)片段以構(gòu)造濾波器消除干擾，對(duì)無(wú)干擾信號(hào)段的提取更加精確。文獻(xiàn)［25］避開計(jì)算復(fù)雜度較高的時(shí)頻分析，采用稀疏分解對(duì)回波信號(hào)進(jìn)行參數(shù)估計(jì)，然后利用聚類方法將干擾與目標(biāo)區(qū)分開。另外，文獻(xiàn)［26］利用目標(biāo)回波和干擾在差分特征域的區(qū)別，設(shè)計(jì)判決量成功辨識(shí)和抑制間歇采樣轉(zhuǎn)發(fā)干擾;基于深度學(xué)習(xí)，文獻(xiàn)［27］提出了基于疊式雙向門控循環(huán)單元（stacked bidirectional gated recurrent unit， SBiGRU）網(wǎng)絡(luò)的方法，準(zhǔn)確提取無(wú)干擾片段以構(gòu)造帶通濾波器，對(duì)干擾進(jìn)行抑制。

本文結(jié)合波形設(shè)計(jì)和時(shí)頻分析方法，首先構(gòu)造脈內(nèi)捷變頻信號(hào)以增加目標(biāo)信號(hào)與干擾信號(hào)的分離度，然后利用短時(shí)分?jǐn)?shù)階傅里葉變換（short-time fractional Fourier transform， STFrFT）進(jìn)行時(shí)頻分析，與傳統(tǒng)的STFT相比，大大提升了時(shí)頻分辨率^［28］，最后結(jié)合圖像學(xué)方法對(duì)干擾進(jìn)行剔除。分辨率的提升，極大地減少了在使用圖像學(xué)方法的過(guò)程中目標(biāo)信號(hào)的損失，與STFT圖像學(xué)方法的一般應(yīng)用產(chǎn)生了本質(zhì)的區(qū)別;該方法與文獻(xiàn)［19］采用的FrFT方法相比，無(wú)需提前感知干擾參數(shù)，設(shè)計(jì)發(fā)射信號(hào)使得子脈沖脈寬與干擾采樣時(shí)間一致，具有更強(qiáng)的適用性。當(dāng)干擾從雷達(dá)主瓣中進(jìn)入時(shí)，由于主瓣帶來(lái)的高接收增益往往具有較高功率，而本文方法在高干信比（jamming-to-signal ratio， JSR）情況下仍能夠有效抑制干擾。實(shí)驗(yàn)結(jié)果表明，該方法在有效抑制干擾的同時(shí)，極大地保留了目標(biāo)信號(hào)，脈壓積累損失較小，且可以較好地適應(yīng)干擾采樣延遲、采樣時(shí)長(zhǎng)與子脈沖脈寬不匹配等現(xiàn)實(shí)問(wèn)題，具有良好的抗干擾效果。

1 信號(hào)模型

1.1 間歇采樣轉(zhuǎn)發(fā)干擾

間歇采樣轉(zhuǎn)發(fā)干擾由干擾機(jī)對(duì)雷達(dá)信號(hào)的切片、存儲(chǔ)和轉(zhuǎn)發(fā)形成，在目標(biāo)信號(hào)周圍形成功率高、密集度大的假目標(biāo)，主要可以分為間歇采樣直接轉(zhuǎn)發(fā)干擾、重復(fù)轉(zhuǎn)發(fā)干擾和循環(huán)轉(zhuǎn)發(fā)干擾，本文重點(diǎn)以直接轉(zhuǎn)發(fā)干擾和重復(fù)轉(zhuǎn)發(fā)干擾為例進(jìn)行分析。

間歇采樣直接轉(zhuǎn)發(fā)干擾對(duì)截獲的雷達(dá)信號(hào)進(jìn)行切片轉(zhuǎn)發(fā)，如圖1（a）所示。設(shè)s（t^，t_n）為干擾機(jī)截獲雷達(dá)所發(fā)射的第n個(gè)脈沖信號(hào)，其中t^代表快時(shí)間，t_n=（n－1）T_r代表慢時(shí)間，T_r為脈沖重復(fù)周期，則間歇采樣直接轉(zhuǎn)發(fā)干擾可以表示為

j₁（t^，t_n）=A_j∑Pp=1rectt^－（p－1）T_sj－τ_j－τ_cτ_js（t^－τ_j，t_n）（1）

式中：recttT= |t|lt;T_p2

0， 其他 表示窗函數(shù)；A_j表示干擾信號(hào)幅度；τ_j表示干擾采樣寬度；τ_c表示干擾采樣延遲；T_sj表示干擾采樣周期;P=［T_p/T_sj］表示當(dāng)前脈沖內(nèi)干擾采樣次數(shù)，其中［·］表示取整運(yùn)算，T_p表示發(fā)射信號(hào)脈寬。

間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾對(duì)截獲的雷達(dá)信號(hào)進(jìn)行切片，并重復(fù)轉(zhuǎn)發(fā)數(shù)次，如圖1（b）所示。同樣設(shè)s（t^，t_n）為干擾機(jī)截獲的雷達(dá)發(fā)射的第n個(gè)脈沖信號(hào)，那么間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾可以表示為

j₂（t^，t_n）=A_j∑Pp=1∑Qq=1rectt^－（p－1）T_sj－qτ_j－τ_cτ_j·

s（t^－τ_j，t_n）（2）

式中：Q=［T_sj/τ_j］－1表示重復(fù)轉(zhuǎn)發(fā)次數(shù)。

此時(shí)，雷達(dá)收到的第n個(gè)接收信號(hào)可以表示為

x（t^，t_n）=s_r（t^，t_n）+j（t^，t_n）+n（t^，t_n）（3）

式中：s_r（t^，t_n）表示目標(biāo)回波信號(hào)；j（t^，t_n）表示間歇采樣轉(zhuǎn)發(fā)干擾；n（t^，t_n）表示回波中的噪聲。

1.2 脈內(nèi)捷變頻信號(hào)模型

脈間捷變頻雷達(dá)可以在脈沖間進(jìn)行發(fā)射信號(hào)載頻的快速跳變，使干擾機(jī)難以跟蹤、截獲雷達(dá)信號(hào)，因此具有良好的抗干擾能力。但間歇采樣轉(zhuǎn)發(fā)干擾可以在一個(gè)脈沖時(shí)間內(nèi)對(duì)雷達(dá)信號(hào)進(jìn)行切片轉(zhuǎn)發(fā)，使得脈間捷變頻雷達(dá)抗干擾性能下降。結(jié)合間歇采樣干擾由采樣轉(zhuǎn)發(fā)模式導(dǎo)致的時(shí)域不連續(xù)的特點(diǎn)，在信號(hào)脈間捷變頻的基礎(chǔ)上加入脈內(nèi)捷變頻信號(hào)，利用脈內(nèi)捷變頻信號(hào)進(jìn)行交替掩護(hù)，使轉(zhuǎn)發(fā)的干擾信號(hào)與被干擾信號(hào)載頻不同，從而增強(qiáng)干擾信號(hào)與目標(biāo)信號(hào)的隔離度。同時(shí)，脈內(nèi)捷變頻信號(hào)使得相鄰子脈沖脈壓后的干擾殘余能夠均勻分布，防止了干擾的積累，從而增加了信干噪比（signal-to-interference-plus-noise ratio， SINR）^［29］。

脈內(nèi)捷變頻信號(hào)模型如圖2（a）所示，假設(shè)雷達(dá)發(fā)射脈沖數(shù)為N，單個(gè)脈沖內(nèi)子脈沖數(shù)為M，子脈沖調(diào)制為線性調(diào)頻，則第n個(gè)發(fā)射脈沖、第m個(gè)脈內(nèi)子脈沖可以表示為

s_sub（t^，n，m）=rectt^－（m－1）T_{p_sub}T_{p_sub}·

exp{j2πf_mt^}·exp{jπμ（t^－（m－1）T_{p_sub}）²}（4）

式中：μ=B_sub/T_{p_sub}表示子脈沖的調(diào)頻斜率;B_sub=B/M表示子脈沖帶寬;B為單個(gè)脈沖帶寬;T_{p_sub}=T_p/M表示子脈沖脈寬;f_m=f_n+c（m）Δf表示第m個(gè)子脈沖的載頻，Δf為脈內(nèi)跳頻間隔，c（m）∈［0，1，…，M－1］為脈內(nèi)隨機(jī)跳頻序列，f_n表示第n個(gè)脈沖對(duì)應(yīng)的載頻。

在得到單個(gè)子脈沖的表達(dá)方式后，將m個(gè)子脈沖組合起來(lái)，可以得到第n個(gè)脈沖的表達(dá)方式：

s_T（t^，t_n）=exp{j2πf_n（t^+t_n）}·∑Mm=1s_sub（t^，n，m）=

exp{j2πf_n（t^+t_n）}·

∑Mm=1rectt^－（m－1）T_{p_sub}T_{p_sub}·

exp{j2πf_mt^}·

exp{jπμ（t^－（m－1）T_{p_sub}）²}（5）

將n個(gè)脈沖結(jié)合起來(lái)，得到雷達(dá)發(fā)射脈沖串信號(hào)，表示為

s（t）=∑Nn=1s_T（t^，t_n）=∑Nn=1exp{j2πf_n（t^+t_n）}·

∑Mm=1s_sub（t^，n，m）=

∑Nn=1exp{j2πf_n（t^+t_n）}·

∑Mm=1rectt^－（m－1）T_{p_sub}T_{p_sub}·

exp{j2πf_mt^}·

exp{jπμ（t^－（m－1）T_{p_sub}）²}（6）

假設(shè)發(fā)射信號(hào)初始載頻為f₀，則第n個(gè)脈沖載頻為f_n=f₀+d（n）ΔF，d（n）∈［0，1，…，N－1］為脈間隨機(jī)調(diào)頻序列，ΔF≥B表示脈間跳頻間隔。假設(shè)子脈沖脈寬T_{p_sub}=5 μs，子脈沖個(gè)數(shù)M=8，子脈沖帶寬與子脈沖調(diào)頻間隔相等，B_sub=Δf=5 MHz，圖2（b）所示為單個(gè)脈沖的時(shí)頻分布仿真圖。

2 干擾抑制方法

2.1 基于STFrFT的回波變換

時(shí)頻分析^［30］可以展示信號(hào)頻率隨時(shí)間變化的動(dòng)態(tài)過(guò)程，通常分為兩類：① 線性時(shí)頻分析，如STFT、小波變換（wavelet transform， WT）等;② 二次型時(shí)頻分析，如Cohen類時(shí)頻分布、維格納維利分布（Wigner-Ville distribution， WVD）等。線性時(shí)頻分析應(yīng)用廣泛，但時(shí)頻分辨率較低^{［28，31］};而雙線性時(shí)頻分析存在自相關(guān)函數(shù)交叉項(xiàng)的干擾^［32］。

FrFT是對(duì)傅里葉變換的推廣，沒(méi)有交叉項(xiàng)的困擾，且對(duì)LFM信號(hào)具有良好的能量聚焦性。因此，本文選擇在FrFT基礎(chǔ)上產(chǎn)生的STFrFT^{［28，31，3334］}作為信號(hào)時(shí)頻分析的工具，在此基礎(chǔ)上區(qū)分目標(biāo)信號(hào)和干擾信號(hào)，并對(duì)干擾信號(hào)進(jìn)行抑制。

如圖3所示，一次傅里葉變換旋轉(zhuǎn)角度為90°，可以看作信號(hào)從時(shí)域到頻域的變換，即F{f（t）}=f（w），F(xiàn){·}表示一次傅里葉變換。而對(duì)信號(hào)進(jìn)行任意角度的旋轉(zhuǎn)，即可實(shí)現(xiàn)信號(hào)從時(shí)域到分?jǐn)?shù)階域的變換，完成FrFT，對(duì)信號(hào)x（t）的FrFT可以表示為

X_p（u）=∫^+∞_－∞K_p（u，t）x（t）dt（7）

式中：p表示FrFT的旋轉(zhuǎn)階數(shù)。定義旋轉(zhuǎn)角度α=pπ/2，則FrFT的核函數(shù)K_p（u，t）可以表示為

K_p（u，t）=A_αexp［jπ（u²cot α－2utcsc α+t²cot α）］， α≠nπ

δ（u－t）， α=2nπ

δ（u+t）， α=（2n±1）π（8）

式中：A_α=exp［－jπsgn（sin α）/4+jα/2］|sin α|¹²。當(dāng)p=1即α=π/2時(shí)，核函數(shù)K_p（u，t）=e^－2jπut 為傅里葉變換核，此時(shí)FrFT等同于傅里葉變換。

對(duì)于p階FrFT，可以通過(guò)－p階FrFT求取逆變換，表達(dá)如下：

x（t）=∫^+∞_－∞K_－p（t，u）X_p（u）du（9）

對(duì)于形如式（10）所示的LFM信號(hào)，F(xiàn)rFT在某一合適的旋轉(zhuǎn)角度具有良好的能量聚焦性，LFM信號(hào)可以形成譜峰，即

s（t）=exp（2jπf₀t+jπμt²）（10）

將式（10）代入式（7），可得LFM信號(hào)的FrFT表達(dá)式為

S_p（u）=∫^+∞_－∞K_p（u，t）s（t）dt=

∫^+∞_－∞K_p（u，t）exp（2jπf₀t+jπμt²）dt=

∫^+∞_－∞A_αexp［jπ（u²cot α－2utcsc α+t²cot α）］·

exp（2jπf₀t+jπμt²）dt（11）

式（11）可化簡(jiǎn)為

S_p（u）=

∫^+∞_－∞A_αexp［jπ（u²cot α）+2t（f₀－ucsc α）+

t²（μ+cot α）］dt（12）

當(dāng)cot α=－μ且f₀=ucsc α?xí)r，LFM信號(hào)的FrFT結(jié)果形成峰值。因此，α_opt=arccot（－μ）為最佳旋轉(zhuǎn)角度，此時(shí)分?jǐn)?shù)階域峰值位置u=f₀/csc α_opt。LFM信號(hào)的FrFT如圖4所示。

結(jié)合STFT，將FrFT推廣到STFrFT，STFT可以表示為

STFT_x（t，f）=∫x（τ）g^*（τ－t）e^－jfτdτ（13）

式中：g（t）表示窗函數(shù)。結(jié)合式（7）和式（13）可得STFrFT為

STFrFT_{（x，p）}（t，u）=∫^+∞_－∞x（τ）g（τ－t）K_p（u，τ）dτ（14）

根據(jù)文獻(xiàn)［35］，高斯窗函數(shù)在時(shí)頻域均具有良好的聚焦性，其表達(dá)式為

g（t）=12πδ_texp-t²2δ²_t（15）

式中：δ_t表示高斯信號(hào)的時(shí)域分辨率，即3 dB寬度。本文選擇高斯窗函數(shù)為變換所需的窗函數(shù)。

在圖2（b）所示的雷達(dá)發(fā)射信號(hào)回波中加入干信比JSR為20 dB的間歇采樣直接轉(zhuǎn)發(fā)干擾，干擾采樣周期為10 μs，采樣寬度為5 μs，采樣延遲τ_c=0。圖5（a）所示為回波信號(hào)的STFT時(shí)頻圖，圖5（b）所示為回波信號(hào)的STFrFT時(shí)頻圖，圖5（a）和圖5（b）中的亮帶均代表干擾信號(hào)，暗帶代表目標(biāo)回波信號(hào)。

從圖5（a）可以看到，間歇采樣轉(zhuǎn)發(fā)干擾信號(hào)經(jīng)STFT后在頻率維存在較高的旁瓣，易與目標(biāo)信號(hào)發(fā)生混疊，在將亮帶剔除時(shí)還會(huì)留下大量幅度較高的干擾殘留，造成SINR下降;而由圖5（b）可以看到，間歇采樣轉(zhuǎn)發(fā)干擾信號(hào)經(jīng)STFrFT后與目標(biāo)回波呈平行狀態(tài)，旁瓣較低，與目標(biāo)之間的交疊極小，易于進(jìn)行亮帶剔除等信號(hào)處理。

2.2 基于圖像處理的干擾提取

對(duì)于回波信號(hào)的STFrFT二維圖像，若按照信號(hào)的幅值進(jìn)行干擾亮帶剔除，將仍然殘留部分與信號(hào)幅度接近或低于信號(hào)幅度的干擾。因此，本文采用圖像處理的方式將信號(hào)與干擾區(qū)分開，并確定干擾的位置，進(jìn)行整體剔除，以最大程度地減小干擾信號(hào)的剩余，同時(shí)保留每一個(gè)目標(biāo)回波子脈沖。與一般的直接剔除被干擾子脈沖的方法相比，目標(biāo)信號(hào)可以得到更好的積累。對(duì)回波信號(hào)的圖像處理分為二值化處理、形態(tài)學(xué)處理和連通域分析，本節(jié)將作出仿真。為突出圖像處理的作用，設(shè)子脈沖脈寬T_{p_sub}=1 μs，干擾采樣周期為2 μs，采樣寬度為1 μs，其他參數(shù)與圖5一致，具體分析如下。

（1） 二值化處理。二值化處理可以消除回波信號(hào)中噪聲的影響，突出信號(hào)與干擾的分布，因此二值化門限應(yīng)略大于噪聲幅值，以在不保留噪聲的同時(shí)盡量保留干擾。同時(shí)，當(dāng)JSR過(guò)大時(shí)，干擾旁瓣仍會(huì)影響二值化處理，因此可以設(shè)二值化門限為TH₁=A_n+η₁A_j，其中A_n為STFrFT二維圖像中的噪聲幅值。由于噪聲在STFrFT中不會(huì)得到能量聚焦而產(chǎn)生峰值，因此可以通過(guò)SNR估算得出A_n;A_j為STFrFT二維圖像中的干擾幅值，可以用STFrFT二維圖像的最大值代替;η1為比例因子，可以選取大于0而遠(yuǎn)小于1的數(shù)字以模擬干擾旁瓣的幅值，在實(shí)際中可以結(jié)合實(shí)際數(shù)據(jù)進(jìn)行修正。圖6（a）所示為回波信號(hào)STFrFT時(shí)頻圖的二值化結(jié)果，圖像保留了全部干擾信息而剔除了噪聲。

（2） 形態(tài)學(xué)處理。數(shù)學(xué)形態(tài)學(xué)可以對(duì)圖像進(jìn)行簡(jiǎn)化，突顯出圖像的主要特征，其基本運(yùn)算有4個(gè)：腐蝕、膨脹、開啟和閉合。腐蝕使物體邊界向內(nèi)部收縮，消除噪點(diǎn);膨脹使物體邊界外擴(kuò)，將兩個(gè)物體連接在一起;開啟運(yùn)算為先腐蝕后膨脹，可以平滑物體的輪廓，剔除突起;閉合運(yùn)算為先膨脹后腐蝕，可以彌合窄縫，填充物體內(nèi)部的空洞。同時(shí)，合適的結(jié)構(gòu)窗口設(shè)計(jì)也對(duì)處理效果十分重要^［3637］，針對(duì)STFrFT后的細(xì)長(zhǎng)干擾塊，本文采用碟形結(jié)構(gòu)進(jìn)行形態(tài)學(xué)處理。二值化后，對(duì)圖像進(jìn)行開閉運(yùn)算，如圖6（b）所示，進(jìn)一步消除時(shí)頻圖像上剩余的噪點(diǎn)、毛刺和粘連等。

（3） 連通域分析。經(jīng)過(guò)二值化和形態(tài)學(xué)處理后的STFrFT圖像存在多個(gè)連通區(qū)域，將這些區(qū)域進(jìn)行識(shí)別并記錄其位置，結(jié)合區(qū)域的最大幅值，即可確定干擾的位置并將其剔除。相鄰像素一般可定義為4領(lǐng)域或8領(lǐng)域，本文采用8領(lǐng)域準(zhǔn)則，連通域?yàn)閳D像中像素值相同且相鄰的區(qū)域，對(duì)形態(tài)學(xué)處理后的圖像進(jìn)行連通域分析，如圖6（c）所示，不同顏色的區(qū)域代表不同的連通域，從圖中可以確定出10個(gè)連通域。

進(jìn)行連通域分析后，對(duì)每個(gè)區(qū)域的幅值進(jìn)行判斷，若某區(qū)域幅值遠(yuǎn)大于信號(hào)幅值，則判斷該連通域?yàn)楦蓴_域，對(duì)該區(qū)域進(jìn)行整體剔除以消除干擾。設(shè)干擾判斷門限TH₂=η₂A_j，其中A_j為STFrFT二維圖像中的干擾幅值，可以用STFrFT二維圖像的最大值代替，η₂∈（10^－JSR/20，1）為比例因子，應(yīng)選取略小于1的數(shù)字以保證干擾被選中而信號(hào)不被選中。由于JSR通常較大，因此η₂的取值范圍較大，JSR估值對(duì)門限取值影響較小。設(shè)STFrFT二維圖像中某一點(diǎn)的幅值為A（t，u），某連通區(qū)域所包含點(diǎn)的幅值集合為A_p={A（t，u），A（t+1，u+1），…，A（t+i，u+j）}，則A_p中存在幅值A(chǔ)（t+m，u+n）gt;TH₂，其中0≤m≤j，0≤n≤j，令A(yù)_p=0。

由此將干擾塊剔除后，對(duì)STFrFT二維圖像進(jìn)行逆STFrFT，可以得到干擾抑制后的信號(hào)。

2.3 分段脈壓

由于脈內(nèi)捷變頻信號(hào)的每個(gè)子脈沖載頻不同，無(wú)法用一個(gè)完整脈沖信號(hào)作為匹配濾波器，傳統(tǒng)的脈沖壓縮方法不再適用。因此，本文采用分段脈壓的方法，當(dāng)子脈沖個(gè)數(shù)為M時(shí)，利用子脈沖構(gòu)造M個(gè)匹配濾波器對(duì)回波信號(hào)分別進(jìn)行脈壓處理后將結(jié)果累積，完成對(duì)回波信號(hào)的脈沖壓縮。

設(shè)干擾抑制后的回波信號(hào)為x～（t^，t_n），利用每個(gè)子脈沖發(fā)射信號(hào)構(gòu)造M個(gè)匹配濾波器，對(duì)回波信號(hào)進(jìn)行分段脈沖壓縮處理，第m個(gè)子脈沖的脈壓的具體表達(dá)式為

y（t^，t_n，m）=x～（t^，t_n）s_sub^*（-t^，n，m）=

x～（t^，t_n）*rect-t^T_{p＼_sub}·exp{j2πf_mt^}·exp{-jπμ（t^）²}（16）

式中：s_sub^*（-t^，n，m）為第n個(gè)脈沖中第m個(gè)子脈沖的匹配濾波器;^*表示卷積運(yùn)算。

將所有子脈沖的分段脈壓結(jié)果積累起來(lái)，得到回波信號(hào)的分段脈壓結(jié)果，具體表示為

y（t^，t_n）=∑M－1m=0y（t^－mT_{p_sub}，t_n，m）（17）

綜上所述，基于STFrFT的間歇采樣干擾抑制流程如圖7所示。

3 仿真實(shí)驗(yàn)

為驗(yàn)證上述基于STFrFT的抗間歇采樣轉(zhuǎn)發(fā)干擾方法的有效性及穩(wěn)健性，本節(jié)設(shè)計(jì)仿真實(shí)驗(yàn)分別分析算法對(duì)間歇采樣直接轉(zhuǎn)發(fā)干擾、重復(fù)轉(zhuǎn)發(fā)干擾的抑制作用，并通過(guò)固定其他參數(shù)，分別分析JSR、SNR、干擾采樣延遲和子脈沖脈寬對(duì)抗干擾性能的影響。具體發(fā)射波形及目標(biāo)參數(shù)如表1所示。

3.1 直接轉(zhuǎn)發(fā)干擾抑制

設(shè)干擾樣式為間歇采樣直接轉(zhuǎn)發(fā)干擾，對(duì)雷達(dá)發(fā)射信號(hào)進(jìn)行同步采樣，即采樣延遲τ_c=0 μs?？紤]一般情況，設(shè)干擾采樣時(shí)長(zhǎng)與子脈沖脈寬不同，干擾采樣周期T_sj=6 μs，干擾采樣時(shí)長(zhǎng)τ_j=3 μs，即采用“切一轉(zhuǎn)一”的策略，干擾前置300 m，目標(biāo)SNR=0 dB，JSR=25 dB，圖像處理二值化門限比例因子η₁=0.05，干擾判斷門限比例因子η₂=0.1。仿真結(jié)果如圖8所示。

對(duì)回波信號(hào)進(jìn)行STFT得到的時(shí)頻分布如圖8（a）所示，由于干擾采樣時(shí)長(zhǎng)與子脈沖脈寬不同，采樣切片分布十分不規(guī)則，而干擾的前置與較高的JSR進(jìn)一步增強(qiáng)了目標(biāo)識(shí)別的難度。脈內(nèi)頻率捷變提升了目標(biāo)回波信號(hào)與干擾信號(hào)的隔離度，但仍存在干擾旁瓣過(guò)高、與目標(biāo)回波信號(hào)混疊嚴(yán)重的問(wèn)題，在時(shí)頻域上直接對(duì)干擾進(jìn)行濾除會(huì)導(dǎo)致干擾旁瓣的殘留和目標(biāo)回波信號(hào)的損失。對(duì)回波信號(hào)進(jìn)行STFrFT得到時(shí)頻分布，如圖8（b）所示，可見經(jīng)過(guò)STFrFT后目標(biāo)回波與干擾的隔離度大大增強(qiáng)，干擾旁瓣也較低，對(duì)該時(shí)頻圖進(jìn)行圖像學(xué)處理后剔除干擾信號(hào)塊，可以在完整剔除干擾信號(hào)的同時(shí)，最大限度地保留目標(biāo)回波信號(hào)。圖8（c）所示為干擾抑制前的分段脈壓結(jié)果，目標(biāo)附近存在著嚴(yán)重的假目標(biāo)干擾，將極大地影響真實(shí)目標(biāo)的探測(cè)。圖8（d）所示為按前文所述算法進(jìn)行干擾抑制后的分段脈壓結(jié)果，可見干擾假目標(biāo)大部分被消除，僅剩一點(diǎn)殘余，目標(biāo)信號(hào)可以被探測(cè)出來(lái)。

圖9（a）所示為SNR相同時(shí)，本文所提干擾抑制方法和無(wú)干擾情況的脈壓結(jié)果的對(duì)比圖;圖9（b）所示為本文所提方法、文獻(xiàn)［16］所提方法、文獻(xiàn)［17］所提方法以及直接剔除被干擾片段方法的脈壓結(jié)果對(duì)比圖。文獻(xiàn)［16］采用脈內(nèi)步進(jìn)LFM信號(hào)，利用STFT進(jìn)行時(shí)頻分析，以時(shí)間維投影方差為門限進(jìn)行時(shí)間維干擾抑制，進(jìn)而在時(shí)頻域?yàn)V除干擾亮帶。文獻(xiàn)［17］和直接剔除被干擾片段方法不進(jìn)行時(shí)頻分析，均采用脈內(nèi)頻率捷變波形，根據(jù)各個(gè)子脈沖的脈壓結(jié)果判斷子脈沖是否被干擾，從而提取合適的子脈沖進(jìn)行積累。不同之處是，文獻(xiàn)［17］采用最大類間方差法根據(jù)各子脈沖脈壓結(jié)果方差值自適應(yīng)計(jì)算門限，可適當(dāng)保留部分受干擾不嚴(yán)重的子脈沖，而直接剔除被干擾片段方法以各子脈沖脈壓結(jié)果能量的最小值為標(biāo)準(zhǔn)，能量近似的子脈沖被保留，能量高于此標(biāo)準(zhǔn)的子脈沖則被直接剔除，包括受干擾不嚴(yán)重的子脈沖。

由圖9（b）可見，本文方法和文獻(xiàn)［16］方法可有效抑制干擾，文獻(xiàn)［16］在時(shí)頻域?qū)Ω蓴_進(jìn)行抑制，可以將干擾徹底剔除，但由于時(shí)頻分辨率較低，在剔除干擾的同時(shí)會(huì)使得目標(biāo)信號(hào)被一起剔除，造成嚴(yán)重的脈壓損失，并且留有較高的干擾旁瓣。由于干擾前置與目標(biāo)回波信號(hào)有部分重疊，子脈沖被干擾情況比較嚴(yán)重且脈壓后干擾幅值不同，導(dǎo)致文獻(xiàn)［17］方法閾值選取不準(zhǔn)確，部分被干擾子脈沖也被積累，因此干擾剔除不徹底，使得目標(biāo)未被凸顯出來(lái)。另外，由于最大類間方差法按閾值對(duì)部分被干擾片段進(jìn)行剔除，導(dǎo)致部分目標(biāo)回波信號(hào)與干擾片段一同被剔除，脈壓積累的幅值受到一定損失。直接對(duì)被干擾片段進(jìn)行剔除可以徹底清除干擾，但當(dāng)被干擾子脈沖數(shù)過(guò)多時(shí)，會(huì)導(dǎo)致大量目標(biāo)回波子脈沖被剔除，使得脈壓結(jié)果損失較大。本文所提方法在STFrFT時(shí)頻圖上結(jié)合圖像學(xué)方法對(duì)干擾進(jìn)行剔除，在完整剔除干擾的同時(shí)，可以最大化保留所有目標(biāo)回波信號(hào)，因此脈壓損失極小，與無(wú)干擾情況下的脈壓結(jié)果接近，同時(shí)干擾殘留較少。

3.2 重復(fù)轉(zhuǎn)發(fā)干擾抑制

設(shè)干擾樣式為間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾，對(duì)雷達(dá)發(fā)射信號(hào)進(jìn)行同步采樣，干擾采樣周期T_sj=12 μs，干擾采樣時(shí)長(zhǎng)τ_j=3 μs，即采用“切一轉(zhuǎn)三”的策略，圖像處理二值化門限比例因子η₁=0.067，干擾判斷門限比例因子η₂=0.1，其他仿真條件與第3.1節(jié)相同，仿真結(jié)果如圖10所示。

圖10（a）所示為回波信號(hào)進(jìn)行STFT后的時(shí)頻分析圖，脈沖信號(hào)被采樣了4次，每個(gè)采樣片段被重復(fù)轉(zhuǎn)發(fā)了3次，大量子脈沖被干擾，頻率維存在幅值較高的干擾旁瓣，對(duì)目標(biāo)信號(hào)造成了影響。對(duì)脈沖回波進(jìn)行STFrFT的結(jié)果如圖10（b）所示，可見存在著多個(gè)干擾“亮帶”，但干擾“亮帶”與目標(biāo)“暗帶”區(qū)分度相較STFT時(shí)頻圖更大，且干擾旁瓣較低。圖10（c）所示為干擾抑制前的回波信號(hào)的脈壓結(jié)果，目標(biāo)信號(hào)周圍存在幅值較大的假目標(biāo)群，難以區(qū)分，對(duì)目標(biāo)檢測(cè)造成了影響。圖10（d）所示為經(jīng)本文算法進(jìn)行干擾抑制后的脈壓結(jié)果，可見干擾被剔除，真實(shí)目標(biāo)得到積累，可以被檢測(cè)出來(lái)。

圖11（a）所示為無(wú)干擾情況和經(jīng)本文所提方法進(jìn)行干擾抑制后的分段脈壓結(jié)果對(duì)比。圖11（b）為本文方法、文獻(xiàn)［16］方法、文獻(xiàn)［17］方法以及直接剔除被干擾片段方法抗間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的對(duì)比圖。由于在間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的情況下被干擾子脈沖較多，文獻(xiàn)［16］方法在時(shí)頻域直接剔除干擾，在分辨率不足的情況下造成目標(biāo)信號(hào)嚴(yán)重?fù)p失。文獻(xiàn)［17］方法所用的最大類間方差法使部分被干擾子脈沖被剔除，部分被干擾子脈沖被保留，所得脈壓結(jié)果存在較大的脈壓增益損失。而在被干擾子脈沖較多的情況下，直接剔除被干擾子脈沖的方法則極大地造成了脈壓增益損失，給目標(biāo)檢測(cè)造成了一定的影響。本文方法在剔除干擾的同時(shí)完整保留了所有子脈沖回波，脈壓積累結(jié)果幅值接近無(wú)干擾情況，且脈壓結(jié)果寬度更窄，因此具有更好的抗干擾效果與目標(biāo)檢測(cè)效果。

3.3 性能分析

為對(duì)本文算法性能進(jìn)行進(jìn)一步測(cè)試分析，分別設(shè)計(jì)實(shí)驗(yàn)驗(yàn)證脈內(nèi)捷變頻波形相對(duì)傳統(tǒng)波形的優(yōu)勢(shì)、干擾參數(shù)感知的作用以及算法在不同JSR、SNR、采樣延遲和子脈沖個(gè)數(shù)條件下的抗干擾效果。

為顯示脈內(nèi)捷變頻波形對(duì)干擾抑制的效果，分別采用傳統(tǒng)LFM信號(hào)和脈內(nèi)捷變頻信號(hào)作為發(fā)射信號(hào)，設(shè)LFM信號(hào)脈寬T_P=50 μs，帶寬B=50 MHz，脈內(nèi)捷變頻信號(hào)和間歇采樣干擾與第3.1節(jié)和第3.2節(jié)中的參數(shù)一致，對(duì)比LFM信號(hào)與脈內(nèi)捷變頻信號(hào)的抗干擾效果，對(duì)比結(jié)果如圖12所示。由于缺乏跳頻產(chǎn)生的頻譜隔離度，LFM信號(hào)干擾抑制不完全，存在大量干擾殘留，且目標(biāo)信號(hào)損失也比脈內(nèi)捷變頻信號(hào)嚴(yán)重。因此，脈內(nèi)捷變頻信號(hào)的選擇可以大大提升干擾抑制效果。

本文所提方法不需要提前獲取干擾采樣時(shí)間、采樣周期等參數(shù)，但對(duì)干擾有關(guān)參數(shù)的感知可以使抗干擾效果更好。結(jié)合雷達(dá)對(duì)干擾的感知識(shí)別，可以通過(guò)估計(jì)間歇采樣轉(zhuǎn)發(fā)干擾的采樣時(shí)間、采樣周期及采樣延遲等，自適應(yīng)優(yōu)化調(diào)整發(fā)射信號(hào)的子脈沖脈寬，使干擾信號(hào)與目標(biāo)信號(hào)的區(qū)分度更高，從而形成主動(dòng)抗干擾，提升抗干擾效果。當(dāng)子脈沖脈寬與間歇采樣轉(zhuǎn)發(fā)干擾的采樣時(shí)間相同時(shí)，目標(biāo)信號(hào)不易被干擾覆蓋，區(qū)分度較大，實(shí)驗(yàn)假設(shè)干擾參數(shù)已知。

設(shè)雷達(dá)發(fā)射信號(hào)子脈沖數(shù)為8，子脈沖脈寬T_{p_sub}=6 μs，其他條件與第3.1節(jié)和第3.2節(jié)相同。通過(guò)干擾認(rèn)知獲取間歇采樣轉(zhuǎn)發(fā)干擾采樣時(shí)長(zhǎng)τ_j=3 μs，對(duì)發(fā)射信號(hào)進(jìn)行參數(shù)優(yōu)化，將子脈沖脈寬調(diào)整為T_{p_sub}=3 μs，子脈沖數(shù)為16，對(duì)比優(yōu)化前后信號(hào)的抗干擾效果，對(duì)比結(jié)果如圖13所示。由于優(yōu)化前子脈沖脈寬與干擾采樣時(shí)間不一致，使得干擾切片與目標(biāo)信號(hào)混疊嚴(yán)重，干擾抑制效果較差，干擾殘留和目標(biāo)損失較多;優(yōu)化后的子脈沖脈寬與干擾采樣時(shí)間一致，增強(qiáng)了目標(biāo)信號(hào)與干擾的隔離度，使抗干擾效果大大增強(qiáng)。

對(duì)本文算法性能進(jìn)行進(jìn)一步測(cè)試分析，由于SINR是衡量抗干擾性能的常用度量，可以同時(shí)表征目標(biāo)積累程度與干擾抑制程度，本文采用脈壓輸出SINR為評(píng)估指標(biāo)，設(shè)計(jì)蒙特卡羅仿真實(shí)驗(yàn)，分別驗(yàn)證算法在不同JSR、SNR、采樣延遲和子脈沖個(gè)數(shù)條件下的干擾抑制性能。

假設(shè)干擾前置300 m，雷達(dá)接收回波脈沖脈壓前JSR在5～40 dB內(nèi)變化，設(shè)置SNR為0 dB，干擾機(jī)對(duì)雷達(dá)發(fā)射信號(hào)進(jìn)行同步采樣，用本文算法和文獻(xiàn)［17］算法分別對(duì)間歇采樣直接轉(zhuǎn)發(fā)干擾和重復(fù)轉(zhuǎn)發(fā)干擾兩種干擾樣式進(jìn)行干擾抑制。蒙特卡羅實(shí)驗(yàn)次數(shù)為100，得到脈壓后SINR隨JSR變化的曲線如圖14所示。相對(duì)于文獻(xiàn)［17］方法，本文算法對(duì)于直接轉(zhuǎn)發(fā)干擾有至少10 dB的SINR提升，重復(fù)轉(zhuǎn)發(fā)干擾有至少5 dB的SINR提升。在SNR相同的情況下，本文算法在無(wú)干擾時(shí)的分段脈壓SNR為41.25 dB，由于間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的干擾轉(zhuǎn)發(fā)次數(shù)更多，受干擾子脈沖數(shù)更多，導(dǎo)致間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的抗干擾性能略低于間歇采樣直接轉(zhuǎn)發(fā)干擾。在低JSR情況下，干擾抑制后的SINR接近于無(wú)干擾情況，而在JSR逐漸增大的過(guò)程中，由于STFrFT時(shí)頻圖仍然存在一定旁瓣，使得抗干擾性能有所下降，但在高JSR條件下依然有不錯(cuò)的抗干擾性能。

假設(shè)雷達(dá)接收回波脈沖脈壓前SNR在-10～10 dB變化，設(shè)置脈壓前JSR為25 dB，其他條件不變，用本文算法和文獻(xiàn)［17］算法分別對(duì)間歇采樣直接轉(zhuǎn)發(fā)干擾和重復(fù)轉(zhuǎn)發(fā)干擾兩種干擾樣式進(jìn)行干擾抑制。蒙特卡羅實(shí)驗(yàn)次數(shù)為100，得到SINR隨SNR的變化曲線如圖15所示。相對(duì)于文獻(xiàn)［17］方法，本文算法對(duì)于直接轉(zhuǎn)發(fā)干擾有至少20 dB的SINR提升，對(duì)于重復(fù)轉(zhuǎn)發(fā)干擾有大約10 dB的SINR提升，可以看出本文算法SINR在SNR在-10 dB到10 dB變化的過(guò)程中始終高于27.5 dB，證明了本文算法在不同的SNR下對(duì)兩種干擾均具有良好的抗干擾能力。而文獻(xiàn)［17］算法在SNR變化過(guò)程中較為穩(wěn)定，但始終低于20 dB;隨著SNR的增大，算法抗干擾性能也同步上升。當(dāng)SNR≥－8 dB后，直接轉(zhuǎn)發(fā)干擾抑制后脈壓SINR達(dá)到30 dB以上;當(dāng)SNR≥－4 dB后，重復(fù)轉(zhuǎn)發(fā)干擾抑制后脈壓SINR穩(wěn)定在29 dB以上，具有良好的干擾抑制效果。

上述仿真均假設(shè)干擾機(jī)對(duì)雷達(dá)脈沖信號(hào)進(jìn)行同步采樣，而現(xiàn)實(shí)中干擾機(jī)對(duì)雷達(dá)脈沖信號(hào)的采樣通常是非同步的，與干擾采樣周期和子脈沖脈寬不同，干擾位置前置的前提條件加在一起使得干擾信號(hào)的排布更加雜亂，多數(shù)子脈沖都會(huì)受到干擾的影響，目標(biāo)信號(hào)與干擾信號(hào)更加難以區(qū)分，令多數(shù)抗干擾方法失效。因此，假設(shè)干擾機(jī)采樣延遲τ_c在0～5 μs變化，設(shè)置SNR為0 dB，JSR為25 dB，其他條件不變，分別對(duì)間歇采樣直接轉(zhuǎn)發(fā)干擾和重復(fù)轉(zhuǎn)發(fā)干擾兩種干擾樣式進(jìn)行100次蒙特卡羅實(shí)驗(yàn)，得到SINR隨τ_c變化的曲線，如圖16所示。可見受到采樣延遲的影響，SINR低于同步采樣時(shí)的數(shù)值，在一定范圍內(nèi)有所波動(dòng)，這是因?yàn)椴蓸友舆t導(dǎo)致一些脈寬較短的干擾片段存在，在對(duì)這些片段進(jìn)行STFrFT時(shí)，能量無(wú)法像長(zhǎng)片段一樣良好聚焦，因此形成了一些幅值較低的干擾塊，可能低于預(yù)設(shè)的干擾檢測(cè)門限，使得干擾有所殘留，并使SINR有所下降。但這些干擾脈寬較小，對(duì)脈壓結(jié)果的影響不大，因此本文方法可以有效地抵抗由干擾采樣延遲造成的性能下降。

脈內(nèi)捷變頻信號(hào)增強(qiáng)了目標(biāo)信號(hào)與干擾信號(hào)之間的隔離度，但子脈沖數(shù)越多，雷達(dá)波形設(shè)計(jì)和產(chǎn)生也越復(fù)雜，因此設(shè)計(jì)實(shí)驗(yàn)驗(yàn)證子脈沖個(gè)數(shù)與干擾抑制效果的關(guān)系。假設(shè)信號(hào)脈沖總脈寬固定為50 μs，取子脈沖個(gè)數(shù)M分別為4，10，20，即子脈沖脈寬分別為12.5 μs，5 μs，2.5 μs，設(shè)置SNR為0 dB，脈壓前JSR在5～40 dB變化，其他條件不變。干擾機(jī)對(duì)雷達(dá)發(fā)射信號(hào)進(jìn)行同步采樣，分別對(duì)間歇采樣直接轉(zhuǎn)發(fā)干擾和重復(fù)轉(zhuǎn)發(fā)干擾兩種干擾樣式進(jìn)行100次蒙特卡羅實(shí)驗(yàn)，得到不同子脈沖個(gè)數(shù)下脈壓后SINR隨JSR變化的曲線，如圖17（a）和圖17（b）所示。可以看出，當(dāng)JSR較低時(shí)，M=10時(shí)，脈壓后SINR較高，干擾抑制效果較好;當(dāng)JSR逐漸增高，M=4時(shí)，脈壓后SINR較M=10和M=20時(shí)更低;在JSRgt;25 dB后，M=20的脈壓后SINR最高，抗干擾效果最佳，因此在高JSR環(huán)境下干擾抑制效果隨子脈沖個(gè)數(shù)的增加而提升。綜上所述，當(dāng)干擾幅度較低時(shí)，可以適當(dāng)減少子脈沖個(gè)數(shù)以降低雷達(dá)信號(hào)設(shè)計(jì)復(fù)雜度，當(dāng)干擾幅度較高時(shí)，因采用子脈沖個(gè)數(shù)較多的發(fā)射信號(hào)，提升了干擾抑制效果。

4 結(jié) 論

本文針對(duì)間歇采樣轉(zhuǎn)發(fā)干擾，設(shè)計(jì)了脈內(nèi)捷變頻脈沖信號(hào)，從波形設(shè)計(jì)的角度提升了目標(biāo)信號(hào)與干擾信號(hào)之間的隔離度，使得雷達(dá)發(fā)射信號(hào)具有主動(dòng)抗干擾特性。之后，采用STFrFT，將目標(biāo)信號(hào)與干擾信號(hào)在STFrFT時(shí)頻圖上進(jìn)行分離，采用圖像學(xué)處理的方法，精確定位干擾塊位置，從而最大化地剔除干擾，同時(shí)保留目標(biāo)信號(hào)。通過(guò)仿真分析，驗(yàn)證了算法在不同JSR、SNR、采樣延遲和子脈沖個(gè)數(shù)條件下的干擾抑制性能。本文算法在高JSR、低SNR環(huán)境下依然有良好的抗干擾性能，同時(shí)良好的信號(hào)干擾區(qū)分度和對(duì)干擾的精確定位使得干擾采樣延遲對(duì)算法的影響很小，說(shuō)明本文算法具有較高的適用性。

參考文獻(xiàn)

［1］ ZHANG J D， ZHU D Y， ZHANG G. New antivelocity deception jamming technique using pulses with adaptive initial phases［J］. IEEE Trans.on Aerospace and Electronic Systems， 2013， 49（2）： 12901300.

［2］ GRECO M， GINI F， FARINA A， et al. Radar detection and classification of jamming signals belonging to a cone class［J］. IEEE Trans.on Signal Processing， 2008， 56（5）： 19841993.

［3］ ROOMW S J. Digital radio frequency memory［J］. Electronics Communications Engineering Journal， 1990， 2（4）： 147153.

［4］ LI C Z， SU W M， GU H， et al. Improved interrupted sampling repeater jamming based on DRFM［C］∥Proc.of the IEEE International Conference on Signal Processing， Communications and Computing， 2014： 254257.

［5］ WANG X S， LIU J C， ZHANG W M， et al. Mathematic principles of interrupted-sampling repeater jamming （ISRJ）［J］. Science China （Information Sciences）， 2007， 50（1）： 113123.

［6］ ZHOU C， LIU Q H， CHEN X L. Parameter estimation and suppression for DRFM-based interrupted sampling repeater jammer［J］. IET Radar， Sonar amp; Navigation， 2018， 12（1）： 5663.

［7］ WEI Z H， LIU Z， PENG B， et al. ECCM scheme against interrupted sampling repeater jammer based on parameter-adjusted waveform design［J］. Sensors， 2018， 18（4）： 11411156.

［8］ YUAN H， WANG C Y， LI X， 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（2）： 6759169.

［9］ LIU Z X， QUAN Y H， WU Y J， et al. Range and Doppler reconstruction for sparse frequency agile linear frequency modulation-orthogonal frequency division multiplexing radar［J］. IET Radar， Sonar amp; Navigation， 202 16（6）： 10141025.

［10］ CAROTENUTO V， AUBRY A， MAIO A， et al. Assessing agile spectrum management for cognitive radar on measured data［J］. IEEE Aerospace and Electronic Systems Magazine， 2020， 36（6）： 2032.

［11］ DING S， SHANG C X， HAN Z Z， et al. Continuous wave radio frequency-screen signal［J］. Modern Defense Technology， 2018， 46（1）： 190194.

［12］ ZHOU C， LIU F F， LIU Q H. An adaptive transmitting scheme for interrupted sampling repeater jamming suppression［J］. Sensors， 2017， 17（11）： 24802495.

［13］ YUAN R F， GAN R B， TANG G F， et al. Range-Doppler and anti-interference performance of cognitive radar detection waveform［C］∥Proc.of the IEEE 12th International Conference on Electronic Measurement amp; Instruments， 2015： 607612.

［14］ HANBALI S， KASTANTIN R. Technique to counter active echo cancellation of self-protection ISRJ［J］. Electronics Letters， 2017， 53（10）： 680681.

［15］ 張彥， 葉春茂， 魯耀兵， 等. 雙曲調(diào)頻波形抑制間歇采樣轉(zhuǎn)發(fā)干擾效果分析［J］. 系統(tǒng)工程與電子技術(shù)， 202 43（11）： 31693176.

ZHANG Y， YE C M， LU Y B， et al. Hyperbolic frequency modulated waveform for interrupted sampling repeater jamming suppression［J］. Systems Engineering and Electronics， 202 43（11）： 31693176.

［16］ 張建中， 穆賀強(qiáng)， 文樹梁， 等. 基于脈內(nèi)步進(jìn)LFM時(shí)頻分析的抗間歇采樣干擾方法［J］. 北京理工大學(xué)學(xué)報(bào)， 2020， 40（5）： 543551.

ZHANG J Z， MU H Q， WEN S L， et al. Anti-intermittent sampling repeater jamming method based on stepped LFM joint time-frequency analysis［J］. Transactions of Beijing Institute of Technology， 2020， 40（5）： 543551.

［17］ 董淑仙， 全英匯， 沙明輝， 等. 捷變頻雷達(dá)聯(lián)合脈內(nèi)頻率編碼抗間歇采樣干擾［J］. 系統(tǒng)工程與電子技術(shù)， 202 44（11）： 33713379.

DONG S X， QUAN Y H， SHA M H， et al. Frequency agile radar combined with intra-pulse frequency coding to resist intermittent sampling jamming［J］. Systems Engineering and Electronics， 202 44（11）： 33713379.

［18］ 劉智星， 杜思予， 吳耀君， 等. 脈間脈內(nèi)捷變頻雷達(dá)抗間歇采樣干擾方法［J］. 雷達(dá)學(xué)報(bào)， 202 11（2）： 301312.

LIU Z X， DU S Y， WU Y J， et al. Anti-interrupted sampling repeater jamming method for interpulse and intrapulse frequency-agile radar［J］. Journal of Radars， 202 11（2）： 301312.

［19］ LIU Z X， QUAN Y H， DU S Y， et al. A novel ECCM scheme against interrupted-sampling repeater jamming using intra-pulse dual-parameter agile waveform［J］. Digital Signal Processing， 202 129： 103652.

［20］ ZHENG H， JIU B， LIU H W. Waveform design based ECCM scheme against interrupted sampling repeater jamming for wideband MIMO radar in multiple targets scenario［J］. IEEE Sensors Journal， 202 22（2）： 16521669.

［21］ ZHOU C， LIU Q H， CHEN X L. Parameter estimation and suppression for DRFM-based interrupted sampling repeater jammer［J］. IET Radar， Sonar amp; Navigation， 2018， 12（1）： 5663.

［22］ GONG S X， WEI X Z， LI X. ECCM scheme against interrupted sampling repeater jammer based on time-frequency analysis［J］. Journal of Systems Engineering and Electronics， 2014， 25（6）： 9961003.

［23］ 杜思予， 劉智星， 吳耀君， 等. 頻率捷變波形聯(lián)合時(shí)頻濾波器抗間歇采樣轉(zhuǎn)發(fā)干擾［J］. 系統(tǒng)工程與電子技術(shù)， 2023， 45（12）： 38193827.

DU S Y， LIU Z X， WU Y J， et al. Frequency agility waveform combined with time-frequency filter to suppress interrupted-sampling repeater jamming［J］. Systems Engineering and Electronics， 2023， 45（12）： 38193827.

［24］ CHEN J， WU W Z， XU S Y， et al. Band pass filter design against interrupted-sampling repeater jamming based on time-frequency analysis［J］. IET Radar， Sonar amp; Navigation， 2019， 13（10）： 16461654.

［25］ DUAN J， ZHANG L， WU Y F， et al. Interrupted-sampling repeater jamming suppression with one-dimensional semi-parametric signal decomposition［J］. Digital Signal Processing， 202 127： 103546.

［26］ 蓋季妤， 姜維， 張凱翔， 等. 基于差分特征的間歇采樣轉(zhuǎn)發(fā)干擾辨識(shí)與抑制方法［J］. 雷達(dá)學(xué)報(bào)， 202 11（4）： 111.

GAI J Y， JIANG W， ZHANG K X， et al. A method for interrupted-sampling repeater jamming identification and suppression based on differential features［J］. Journal of Radars， 202 11（4）： 111.

［27］ CHEN J， XU S， ZOU J W， et al. Interrupted-sampling repeater jamming suppression based on stacked bidirectional gated recurrent unit network and infinite training［J］. IEEE Access， 2019， 7： 107428107437.

［28］ CAPUS C， BROWN K. Short-time fractional Fourier methods for the time-frequency representation of chirp signals［J］. Journal of the Acoustical Society of America， 2003， 113（6）： 32533263.

［29］ 張建中， 穆賀強(qiáng)， 文樹梁， 等. 基于脈內(nèi)LFM-Costas頻率步進(jìn)的抗間歇采樣干擾方法［J］. 系統(tǒng)工程與電子技術(shù)， 2019， 41（10）： 21702177.

ZHANG J Z， MU H Q， WEN S L， et al. Anti-intermittent sampling jamming method based on intra-pulse LFM-Costas frequency stepping［J］. Systems Engineering and Electronics， 2019， 41（10）： 21702177.

［30］ BOASHASH B. Time-frequency signal analysis and processing： a comprehensive reference［M］. 2nd ed. Amsterdam： Elsevier Ltd， 2016.

［31］ TAO R， LI Y L， WANG Y. Short-time fractional Fourier transform and its applications［J］. IEEE Trans.on Signal Processing， 2010， 58（5）： 25682580.

［32］ LINNET L， MORRISON S， NICHOLSON P. The analysis of signals containing mixtures of linear chirps［J］. Proceedings of Symposium on Underwater Bio-sonar and Bioacoustic Systems， 200 23（4）： 5562.

［33］ DENG B， JIN D， LUAN J B. Adaptive short-time fractional Fourier transform based on minimum information entropy［J］. Journal of Beijing Institute of Technology， 202 30（3）： 265273.

［34］ ZAKRZEWSKI M， RAITTINEN H， VANHALA J. Comparison of center estimation algorithms for heart and respiration monitoring with microwave Doppler radar［J］. IEEE Sensors Journal， 201 12（3）： 627634.

［35］ TAO R， DENG B， WANG Y. Research progress of the fractional Fourier transform in signal processing［J］. Science in China. Series F， Information Sciences， 2006， 49（1）： 125.

［36］ DENG L Z， ZHANG J K， XU G X， et al. Infrared small target detection via adaptive M-estimator ring top-hat transformation［J］. Pattern Recognition， 202 112（1）： 107729.

［37］ BAI X Z， ZHOU F G. Analysis of new top-hat transformation and the application for infrared dim small target detection［J］. Pattern Recognition， 2010， 43（6）： 21452156.

作者簡(jiǎn)介

李晉杰（1998—），男，碩士研究生，主要研究方向?yàn)槔走_(dá)智能抗干擾及波形優(yōu)化。

曹運(yùn)合（1978—），男，教授，博士，主要研究方向?yàn)槔走_(dá)信號(hào)處理、陣列信號(hào)處理、高速實(shí)時(shí)信號(hào)處理、DSP與FPGA開發(fā)。

張鈺林（2001—），男，博士研究生，主要研究方向?yàn)榻葑冾l信號(hào)處理及抗干擾。

王 蒙（1987—），女，工程師，碩士，主要研究方向?yàn)槔走_(dá)探測(cè)處理、分布式協(xié)同探測(cè)處理。