基于噪声卷积调制的SAR虚假信号生成新方法黄大通①      邢世其*①      刘业民②      李永祯①      肖顺平①①(国防科技大学 电子信息系统复杂电磁环境效应国家重点实验室   长沙   410073)②(中国人民解放军第32579部队   桂林   541001)摘   要：针对传统噪声卷积调制的合成孔径雷达(SAR)虚假信号方法存在距离向位置滞后、方位向压制范围不可控的缺陷，该文提出了一种改进的虚假信号生成方法。该方法首先对截获信号作快时间域移频调制，以控制掩护面的距离向位置；接着将其与经过了慢时间域滤波处理的噪声模板卷积，以控制掩护面积。理论分析与仿真结果表明，相比于传统噪声卷积调制，该文所提的方法可有效控制掩护面的距离向位置和面积，即使在较大侦察误差下仍能对局部场景实施掩护，提高了相同条件下的干扰能量利用率，对实际工程应用具有一定的参考价值。关键词：合成孔径雷达；噪声卷积调制；移频调制；慢时间滤波；相参压制中图分类号：TN95 文献标识码：A 文章编号：2095-283X(2020)05-0898-10DOI: 10.12000/JR20094引用格式：黄大通, 邢世其, 刘业民, 等. 基于噪声卷积调制的SAR虚假信号生成新方法[J]. 雷达学报, 2020, 9(5):898–907. doi: 10.12000/JR20094.Reference format: HUANG Datong, XING Shiqi, LIU Yemin, et al. Fake SAR signal generation method basedon noise convolution modulation[J]. Journal of Radars, 2020, 9(5): 898–907. doi: 10.12000/JR20094.Fake SAR Signal Generation Method Based on NoiseConvolution ModulationHUANG Datong①      XING Shiqi*①      LIU Yemin②      LI Yongzhen①      XIAO Shunping①①(State Key Laboratory of Complex Electromagnetic Environment Effects on Electronics and InformationSystem, National University of Defense Technology, Changsha 410073, China)②(Unit 32579 of the PLA, Guilin 541001, China)Abstract: Suppression position of the traditional noise convolution modulation Synthetic Aperture Radar (SAR)jamming lags behind in range and suppression area in azimuth is uncontrollable. Considering this defect, anenhanced jamming method is proposed herein. First, the frequency of the intercepted signal is shifted in fast-time to control the suppression position in range. Then, the convolution with the noise is implemented, whichhas been filtered in slow-time, to control the suppression area in azimuth. Theoretical analysis and simulationresults demonstrate that the proposed jamming method can efficiently control the jamming position in rangeand suppression area when compared with the traditional noise convolution modulation jamming. Even if somereconnaissance errors exist, the local scenario can still be shielded effectively. Furthermore, the utilizationefficiency of jamming energy is also improved under the same condition, which will provide some referencevalues and inputs for engineering applications.Key words: Synthetic Aperture Radar (SAR); Noise convolution modulation; Frequency shift modulation; Slow-time filtering; Coherent suppression  收稿日期：2020-07-07；改回日期：2020-09-17；网络出版：2020-10-07*通信作者： 邢世其 xingshiqi_paper@163.com        *Corresponding Author: XING Shiqi, xingshiqi_paper@163.com基金项目：国家自然科学基金(61971429, 61901499)Foundation Items: The National Natural Science Foundation of China (61971429, 61901499)责任主编：周峰     Corresponding Editor: ZHOU Feng第9卷第5期 雷  达  学  报 Vol. 9No. 52020年10月 Journal of Radars Oct. 20201    引言合成孔径雷达[1](Synthetic Aperture Radar,SAR)作为一种高分辨率成像雷达，具有全天时、全天候和透视性等优良特性，被广泛应用于情报搜集、战场监控以及导弹末端图像匹配制导等众多军事领域，已成为现代战场“察打一体化”武器系统的核心传感器[2,3]，对我军重要的军事设施和部署构成了极大威胁。鉴于此，展开针对SAR干扰技术的研究，对提升我国的国土防御能力和战时生存能力都具有重要军事价值。SAR干扰技术可分为压制干扰[4]和欺骗干扰[5]，其中普通噪声干扰与SAR信号间不具备相干性，因此对干扰机发射功率需求较大，易受到敌方的反辐射打击；而欺骗干扰则需较多的高精度侦察参数，算法复杂度较高，干扰实时性较差，不利于工程实现。1999年，国际著名电子战专家Schleher[6]在《Electronic Warfare in the Information Age》一书中提出了“灵巧噪声干扰”的概念，由于该类干扰具备自动瞄准信号频率，可获得雷达信号处理增益以及调制方式灵活、干扰效果多样等优点，因此一经提出便受到了科研工作者的广泛关注[7–16]。文献[11]通过对灵巧噪声干扰输出要求的逆向研究，推导出了其卷积调制的数学模型。文献[12]分析了灵巧噪声干扰对SAR系统所具有的特点和优势，提出了利用数字射频存储器来实现干扰调制的方法。文献[13]提出了随机脉冲卷积的干扰方法，通过控制随机脉冲的时延和范围，以此控制干扰带宽度。文献[14]则通过选取不同的噪声时长以及噪声筛选门限值来产生密集假目标干扰，兼具了欺骗与压制效果。文献[15]研究了相干和非相干情况下的噪声卷积的效果，指出可通过改变干扰调制函数产生欺骗和压制的效果。文献[16]研究了脉冲卷积调制干扰对超宽带SAR的效果。然而，上述的传统噪声卷积调制干扰却在距离向始终滞后于干扰机，无法掩护其前方场景，且干扰分散于整个方位向，其能量利用率较低。针对以上问题，本文提出了一种对传统噪声卷积调制干扰的改进方法，其创新之处在于：通过快时间移频调制，实现了对干扰距离向位置的灵活控制，解决了传统噪声卷积干扰滞后的问题；通过对噪声调制模板的慢时间域滤波，实现了对干扰方位向压制范围的精准控制，提高了能量利用率；由于本文所提的方法是在离线模式下对噪声调制模板作预处理，因此不会影响干扰实时性，便于工程实现。本文的主要内容分为3个部分，第1部分首先建立了干扰的信号模型，进而推导分析了其对SAR的干扰效果；第2部分对所提干扰方法的性能作了理论分析；第3部分通过仿真实验对所提方法的干扰效果作了验证。2    基于噪声卷积调制的SAR干扰改进方法2.1  对SAR的干扰信号模型hvxta =0 (vta;0;h) (xj;yj;0)在图1所示的干扰场景几何模型中，SAR平台在高度 的航迹线上以速度 沿 轴正向飞行，记慢时间 的地面投影为坐标原点O，平台位置为。场景中部署了一部干扰机，其坐标位置为 ，则经泰勒公式[1]展开后干扰机到SAR平台的瞬时斜距为R(ta)=√(vta xj)2+y2j+h2R0 xjvR0ta+v22R0t2a(1) R0 =√x2j+y2j+h2其中， 表示干扰机到SAR平台的初始斜距。雷达发射信号为线性调频(Linear FrequencyModulation, LFM)信号s(tr)=rect(trTp)exp(jπkrt2r)exp(j2πfctr) (2)trTp krfcrect(tr)={1; jtrj  1/20; 其他式中， 表示快时间， 为信号时宽， 表示调频斜率， 为信号载频， 。为破坏SAR系统对重要军事设施和部署的侦察，干扰机在被雷达波束照射的合成孔径时间内，对其信号进行截获，并作下变频处理。si(tr;ta) =s(tr)[tr  (ta)]exp( j2πfctr)=rect(tr  (ta)Tp)exp{jπkr[(tr  (ta))]2}
expf j2πfc(ta)g (3) si(tr;ta)(ta)=R(ta)c c式中， 表示干扰机的截获信号， 表示快时间上的卷积运算， 表示雷达到干扰机单程斜距的时延， 为光速。 xy成像区域Ota=0 v干扰机(xj,  yj, 0)R(ta)SAR平台h 图 1 干扰场景的几何模型Fig. 1  The geometry model of jamming scenario第5期 黄大通等：基于噪声卷积调制的SAR虚假信号生成新方法 899根据LFM信号的时频耦合特性[17–20]：信号频率的偏差将造成其时延的改变，故干扰机对截获信号作式(4)所示的移频调制，以此实现对干扰距离向位置的控制。m1(tr;ta) =si(tr;ta)expfj2π∆ftrg=rect(tr  (ta)Tp)exp{jπkr[(tr  (ta))]2}
expfj2π∆ftrgexpf j2πfc(ta)g (4) ∆f式中， 为干扰机的移频调制量。m1(tr;ta)接着，干扰机将调制信号 与预先载入的噪声作卷积处理。m2(tr;ta) =n(tr;ta)m1(tr;ta) (5) n(tr;ta) TnTsTs式中， 是快时间宽度为 ，慢时间宽度为的二维高斯白噪声调制模板， 为雷达的合成孔径时间，噪声模板在快、慢时间上均表现为非相关。m2(tr;ta) sj(tr;ta)最后，干扰机将调制信号 经上变频处理后，形成干扰信号 。sj(tr;ta) =m2(tr;ta)[t j] expfj2fctrg (6) j其中， 为干扰机的固定转发时延。由此，雷达接收到的基带干扰信号为sj0(tr;ta) =sj(tr;ta)[t (ta)]exp( j2πfctr)=n(tr;ta)s′j0(tr;ta)(7)s′j0(tr;ta)其中， 是干扰机处经移频调制后的基带回波s′j0(tr;ta) =rect[tr  2(ta) jTp]rect[ta xj/vTs]
exp{jπkr[tr  2(ta) j]2}
expfj2π∆f[tr  (ta) j]g
expf j4πfc(ta)g (8)2.2  对SAR的干扰结果分析下面以经典的R-D成像算法[1]为例，对SAR的干扰原理作推导分析。SAR系统对接收到的基带干扰作距离向压缩和距离迁徙矫正后，即有jrc(tr;ta) =sj0(tr;ta)hr(tr)=n(tr;ta)[s′j0(tr;ta)hr(tr)]=n(tr;ta)psf(tr)rect[ta xj/vTs]expf j4πfc(ta)g(9)hr(tr)=rect(tr/Tp)exp{ jπkrt2r}psf(tr)=sinc{Br[tr  20 j +∆fkr]}Br0 =R0/c其中， 为距离向压缩函数，为雷达信号的点扩展函数， 为信号带宽， 。可以看出，噪声卷积调制干扰的距离向压缩是噪声模板与点扩展函数的卷积，故干扰的距离向增益主要来自于雷达信号的点扩展函数，这也是该类干扰能够降低干扰机发射功率的理论依据[21]。sinc(
) (
)鉴于SAR信号带宽较大， 函数可近似为函数[22]，故式(9)所示的干扰距离向压缩结果可近似为jrc(tr;ta)n(tr  20 j +∆fkr;ta)
rect[ta xj/vTs]expf j4πfc(ta)g(10)在方位向上，为便于作进一步推导分析，只考虑式(10)中的慢时间项[23]，并代入式(1)所示的瞬时斜距，则慢时间上的干扰信号可近似为ja(ta)na(ta)rect[ta xj/vTs]exp{j4πvxjR0ta}
exp{ jπkat2a}(11) na(ta)ka =2v2R0 式中， 仅为慢时间上的噪声， 是方位向调频斜率， 为信号波长。na(ta)显然，式(11)即是对LFM信号的噪声乘积调制干扰模型[24]，噪声 可看作是分布在频谱上的无数个单频信号的合成，如式(12)所示(忽略幅度的影响)。na(ta)=Bn ∑fi= Bnexpfj2πfitag (12)Bn其中， 为噪声模板的慢时间单边带宽。则式(11)可转化为ja(ta)=Bn ∑fi= Bnrect[ta xj/vTs]exp{j4πvxjR0ta}
expfj2πfitagexp{ jπkat2a}(13)由此可得干扰的方位向压缩结果Ia(ta)=ja(ta)ha(ta)=Bn ∑fi= Bnsinc[Ba(ta xjv fika)]Bn ∑fi= Bn(ta xjv fika)(14)ha(ta)=rect[ta/Ts]exp{jπkat2a}Ba Bn其中， 为方位向压缩函数， 为方位向带宽。干扰的方位向成像结果是无数个冲击函数的合成，即形成了压制效果，且鉴于调制模板为双边谱，故干扰压制范围是以干扰机为中心向左右两边扩散，其扩散长度由噪声调制模板的慢时间带宽 决定。900 雷  达  学  报 第9卷20+j  ∆f/kr Tn ∆f Tn ∆f=0综上所述，在距离向上，由式(10)可知，干扰的成像结果等效于将噪声模板作平移，噪声的起始时刻为 ，压制时长为 。故，干扰方可通过设置移频调制量 和噪声模板的快时间长度 来实现对干扰的距离向压制位置和范围的灵活控制，具体计算方法分别如式(15)和式(16)所示。然而，移频调制却是以干扰功率的失配损耗为代价[25–28]，因此移频量往往不易过大。反之，对于传统噪声卷积调制干扰，由于移频量 ，故噪声将始终滞后于干扰机，无法掩护其前方场景。∆f= 2kr∆Rshiftc(15) Tn =2∆Rjc(16)∆Rshift ∆Rshift >0 ∆Rshift <0 ∆Rj式中， 为干扰方设定的距离向偏移量，，向后移动， ，向前移动；为干扰方设定的距离向压制长度。Bn在方位向上，由式(14)可知，干扰的压制范围取决于噪声模板的慢时间单边带宽 ，因此本文所提的干扰方法通过在离线状态下预先对噪声调制模板作慢时间域的低通滤波，限制其慢时间带宽，以控制干扰的方位向压制范围，具体计算方法如式(17)所示。反之，对于传统噪声卷积调制干扰[15,22]，由于噪声模板未经过滤波处理，因此其干扰能量将扩散于整个方位向。fstop =Bn =∆XjvR0(17)∆Xjfstop式中， 为干扰方设定的方位向压制长度，为低通滤波器的截止频率。3    对SAR干扰性能的理论分析3.1  侦察误差的影响分析kr vR0 ∆kr ∆v ∆ ∆R0"kr=∆kr/kr"v =∆v/v "= ∆/ "R0=∆R0/R0由于本文所提干扰为相参干扰，干扰的执行必然离不开对SAR平台相关参数的侦察，因此有必要分析侦察误差对干扰效果的影响情况。从上一小节可知，调制干扰所需的侦察参数包括信号的调频斜率 ，SAR平台飞行速度，信号波长 以及干扰机与雷达间的初始斜距 ，设其对应的绝对侦察误差分别为 ,  ,  以及 ，则可将相应的相对误差分别记作 ,  , 以及 。根据式(15)和式(17)，可得侦察误差下的移频调制量和滤波器截止频率分别为∆f′=(1 +"kr)∆f (18)f′stop=(1+"v)(1+")(1 +"R0)fstop(19)结合式(10)和式(14)，可得侦察误差下的干扰距离向起始位置和方位向压制范围分别为R′=c0+cj2 c∆f2kr(1+"kr) (20) ∆X′=(1+"v)(1+")(1 +"R0)∆X (21)"kr"v""R0可以看出，侦察误差并不会影响干扰的输出形式，只会造成实际干扰结果与预定效果的偏差：干扰的距离向起始位置与 成正比，而干扰的方位向压制范围则与 成正比，与 和 成反比。因此，即使在侦察设备不能提供准确雷达参数的情况下，我方依然可以利用先验信息(如典型SAR场景的参数范围)，采用直观估值的方法估计出来，完成对场景的掩护，故本文所提的干扰方法对侦察设备的要求较低，战场生存能力强。3.2  功率增益分析对常规噪声干扰，由于其与雷达信号间不具备相参性，无法获得信号处理增益，故对干扰功率的需求较大[29]；对传统噪声卷积调制干扰，其仅能获得距离向信号处理增益，无法控制干扰的方位向能量分布。然而，本文所提的干扰方法却可获得雷达的二维信号处理增益，降低对干扰功率的需求。Tn+Tp1/BrTn+1/Br在距离向上，与传统噪声卷积调制干扰的分析方法相同，这里不考虑移频调制的影响。从式(5)可知，压缩前的干扰信号时宽为 ；从式(9)可知，压缩后的雷达信号时宽为 ，噪声调制模板不变，故压缩后的干扰信号时宽为 。由于压缩网络是无源，故压缩前后的干扰能量保持不变，即有Ji(Tn+Tp)=Jo(Tn+ 1/Br) (22)JiJo式中， 和 分别为距离向压缩前后的干扰功率。Gr从而推得干扰的距离向增益 为Gr =JoJi=1 +Tp 1/BrTn+1/Br(23)Tn ≫1/BrTn ≪1/Br Gr TpBr可以看出，在SAR信号的时宽带宽一定时，干扰的距离向增益直接取决于噪声调制模板的快时间宽度：当 时，干扰增益近似为1，即未获得任何信号处理增益，干扰效果与常规噪声干扰相同，而由式(16)可知，此时干扰的距离向压制范围却较大；当 时，干扰增益近似为 ，即干扰所获增益与目标回波所获增益基本相等，干扰所需功率将极大降低，但此时的距离向压制范围却很小。在方位向上，与文献[24]对噪声乘积调制干扰的分析方法相同，由式(11)和(14)所示的干扰模型第5期 黄大通等：基于噪声卷积调制的SAR虚假信号生成新方法 901Ts2BnTsBa Ga和压缩结果可知，压缩前后的干扰信号时宽分别为和 。同样，根据能量守恒定律可得干扰的方位向增益 为Ga =Ba2Bn(24)Bn可以看出，干扰的方位向增益与噪声调制模板的慢时间带宽成反比： 越小，干扰增益越大，而由式(17)可知，其方位向压制范围却将越小；反之亦然。G故，本文所提干扰获得的二维总增益 为G=Gr Ga(25)3.3  算法复杂度分析o(Nslog2Ns) o(
) Nso(Ns) o(Nslog2Ns)通过第2.1节的推导分析可知，本文所提干扰方法主要包括了离线状态下的干扰调制和实时调制两个阶段。在离线状态下，干扰机首先产生满足距离向压制范围要求的噪声模板，并将其变换到慢时间频域；接着，结合前期的电子侦察所获参数，按照方位向的压制范围要求，对模板作慢时间频域的低通滤波处理，并将其反变换到慢时间域。在该阶段中，噪声模板的慢时间FFT变换对应的计算复杂度[30]为 ，其中 表示“在···量级”或“与···呈正比”， 为慢时间FFT长度；噪声模板的慢时间频域低通滤波为向量间的点对点乘积，其相应的计算复杂度为 ；通过逆FFT将调制好的噪声模板变换到慢时间域的计算复杂度为。由于以上过程均在离线状态下进行，因此干扰机具有较为充裕的运算时间。在实时调制阶段，干扰机首先对截获信号作快时间移频调制，再将其变换到快时间频域，与噪声模板作频域乘积，最后将调制信号反变换到快时间域。在该阶段中，算法的主要计算量来自于快时间移频调制、快时间FFT变换、向量的点对点乘积以o(Nt) o(Nflog2Nf) o(Nf) o(Nflog2Nf) Nt Nf及快时间逆FFT变换，相应的计算复杂度分别为,  ,  和 ，其中 为干扰机对截获信号的采样点数， 为快时间FFT长度。在被雷达波束照射的合成孔径时间内，干扰机会不断截获到每个脉冲重复周期的SAR信号，干扰调制过程会重复进行，因此对实时性要求较高。鉴于在离线状态下已完成了对噪声模板的调制，因此剩余计算量控制在可接受的范围。4    仿真实验与结果分析4.1  干扰的有效性验证200 m
s 1在仿真实验中，SAR平台的飞行高度为5 km，飞行速度为   ，工作于正侧视模式；信号载频为10 GHz，脉冲时宽和带宽分别为5 µs和100 MHz，方位向波束宽度为0.03 rad，雷达下视角为45°，脉冲重复频率为800 Hz，成像区域的方位向和地距向范围分别为300 m和400 m。干扰机部署于场景中心，其固定转发时延为0.2 µs。2:3 dB图2所示为传统噪声卷积调制的干扰结果，其中雷达接收端的干信比为   ，噪声模板的快时间宽度采用式(16)所示的设置方法，压制长度分别设置为100 m和25 m。可以看出，干扰能量分布于整个方位向，形成了压制条带，但却在距离向始终滞后于干扰机(固定转发时延的滞后距离约42 m)，无法对其前方场景进行掩护。另一方面，图2(a)中的干扰压制范围较大，能量分布较为稀疏，对场景的掩护效果较差；但在图2(b)中，压制范围的缩小使干扰能量更加集中，对局部场景的掩护效果更好。2:3 dB 300m400m图3—图5所示为本文所提干扰方法的成像结果，雷达接收端的干信比依然设置为   。当干扰的距离向偏移量设置为–242 m，压制面积设置为 时，干扰实现了对全成像面的压  (a)  距离向压制长度100 m  (a) Suppression scope 100 m in range -150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200地距向 (m)(b)  距离向压制长度25 m(b) Suppression scope 25 m in range-150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200地距向 (m) 图 2 传统噪声卷积调制的干扰结果Fig. 2  The imaging results of the traditional noise convolution modulation jamming902 雷  达  学  报 第9卷制，如图3所示。然而，鉴于覆盖区域太大，能量分布稀疏，故干扰对场景的掩护效果较差。当同时减小距离向压制范围和噪声调制模板的慢时间频谱宽度时，如图4(a)所示，干扰压制面积将不断缩小，干扰能量愈发集中，对局部场景的掩护效果也会越好，如图4(b)所示。特别地，当压制面积仅为 (b) SAR 干扰结果(b) SAR jamming result-150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200地距向 (m)(a)  噪声模板的二维频域图(a) 2D spectrum of the noise template0100   50   0    -50 -100-400  -200    0    200   4000.20.4 幅度0.60.81.0慢时间频域 (Hz)快时间频域 (MHz) ∆Rshift = 242m 300m400m图 3 距离向偏移量 ，压制面积∆Rshift = 242m 300m400mFig. 3  Range offset  , suppression area  (b) SAR 干扰结果(b) SAR jamming result-150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200地距向 (m)(a)  噪声模板的二维频域图(a) 2D spectrum of the noise template0100   50   0    -50 -100-400  -200    0    200   4000.20.4 幅度0.60.81.0慢时间频域 (Hz)快时间频域 (MHz) ∆Rshift = 142m 100m200m图 4 距离向偏移量 ，压制面积∆Rshift = 142m 100m200mFig. 4  Range offset  , suppression area  (b) SAR 干扰结果(b) SAR jamming result-150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200地距向 (m)(a)  噪声模板的二维频域图(a) 2D spectrum of the noise template0100   50   0    -50 -100-400  -200    0    200   4000.20.4 幅度0.60.81.0慢时间频域 (Hz)快时间频域 (MHz) ∆Rshift = 167m 20m50m图 5 距离向偏移量 ，压制面积∆Rshift = 167m 20m50mFig. 5  Range offset  , suppression area 第5期 黄大通等：基于噪声卷积调制的SAR虚假信号生成新方法 90320m50m时，指定的局部场景则被彻底压制，如图5所示。通过以上实验结果对比可知，本文所提的干扰方法成功实现了在固定干扰机下对压制面距离向位置和方位向范围的灵活控制，有效弥补了传统噪声卷积调制干扰滞后于干扰机，方位向能量分散的缺陷，提高了干扰能量的利用率，与理论分析相一致。4.2  侦察误差对干扰效果的影响kr 30% vR0 50% kr vR0为反映干扰对侦察参数的敏感性和依赖性，这里以图5(b)为基准，考虑较大误差的恶劣情况：的相对侦察误差为 ,  ,  和 的相对侦察误差为 ，干扰结果如图6所示。相比于无误差的干扰结果，由式(20)可计算得 的误差将导致干扰超前约50 m，干扰的起始位置为4825 m，如图6(a)所示；由式(21)可计算得 的误差将导致压制面的方位向范围展宽到30 m，如图6(b)所示，而 和 的误差将导致压制面的方位向范围缩小到13 m，如图6(c)和图6(d)所示。由此可知，即使在侦察误差较大的情况下，本文所提的干扰方法依然能完成对局部场景的有效压制，仿真结果与理论分析相一致。4.3  干扰增益验证为进一步验证本文所提干扰方法相比于传统噪2:410 7s13 dB声卷积调制干扰的增益情况，表1记录了移频量为零时，噪声模板的不同快时间宽度和慢时间单边带宽对干扰所获增益的影响情况。可以看出，对于传统噪声卷积调制干扰，当不断减小噪声模板的快时间宽度时，其干扰所获增益不断增大，因此对局部场景的压制效果也将越好，与图2的实验结果相吻合。而对于本文所提的干扰方法，当不断减小噪声调制模板的快时间宽度和慢时间带宽时，其所获的增益也将不断增大，因此其对场景的压制效果也会越好，与图3—图5的实验结果相吻合。特别地，在相同的距离向压制范围(快时间宽度为   )条件下，由于减小了噪声模板的慢时间带宽，本文所提干扰方法获得的增益比传统噪声卷积调制干扰提高了约   ，干扰能量更加集中，有效提高了能量利用率，实验数值与理论数值几乎相等，验证了文中理论分析的正确性。5    结束语本文在传统噪声卷积调制干扰的基础上，提出了对SAR的干扰新方法：干扰方通过设置噪声调制模板的快时间宽度和慢时间带宽，直接控制干扰的压制面积；结合移频调制，实现对干扰距离向位置 -150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200地距向 (m)-150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200地距向 (m)(a)  ekr=0.3(c)  el=0.5-150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200-150 -100 -50 0 50 100 150方位向 (m)480048504900495050005050510051505200(b)  ev=0.5(d)  eR0=0.5地距向 (m) 地距向 (m) 图 6 侦察误差下的干扰结果Fig. 6  The jamming results under reconnaissance error904 雷  达  学  报 第9卷的控制，有效弥补了传统噪声卷积调制干扰存在的距离向位置滞后，方位向压制范围不可控的缺陷。文中首先建立了干扰的信号模型，推导出了其压缩结果的解析表达，并基于此对干扰的压制范围进行了详细分析；接着，文中从侦察误差影响、功率增益以及算法复杂度3个方面对所提干扰方法的性能情况作了理论论证。实验结果表明本文所提的干扰方法可有效实现对干扰的距离向位置和压制面大小的灵活控制，即使在较大侦察误差下仍能完成对局部场景的要地掩护，且相比于传统噪声卷积调制干扰，其在相同条件下所获的增益更大，对要地的掩护效果更好。参 考 文 献CUMMING I G and WONG F H. Digital Processing ofSynthetic Aperture Radar Data: Algorithms andImplementation[M]. Boston: Artech House, 2005.[1]纪朋徽, 代大海, 吴昊, 等. SAR成像电子对抗技术综述[J]. 无线电工程, 2019, 49(6): 508–513.JI Penghui, DAI Dahai, WU Hao, et al. Review of SARimaging electronic countermeasures[J]. Radio Engineering,2019, 49(6): 508–513.[2]SUN Hongbo, SHIMADA M, and XU Feng. Recentadvances in synthetic aperture radar remote sensing-systems, data processing, and applications[J]. IEEEGeoscience and Remote Sensing Letters, 2017, 14(11):2013–2016. doi: 10.1109/LGRS.2017.2747602.[3]CHANG Xin and DONG Chunxi. A barrage noise jammingmethod based on double jammers against three channelSAR GMTI[J]. IEEE Access, 2019, 7: 18755–18763. doi:10.1109/ACCESS.2019.2897043.[4]刘玉玲. SAR有源假目标精确位置欺骗干扰技术研究[D]. [硕士论文], 国防科学技术大学, 2012.LIU Yuling. Research on SAR precisely position deceptionjamming[D]. [Master dissertation], National University ofDefense Technology, 2012.[5]SCHLEHER D C. Electronic Warfare in the InformationAge[M]. Boston, MA: Artech House, 2000.[6]李欣, 王春阳, 原慧, 等. 雷达灵巧干扰样式建模及仿真研究[J].航天电子对抗, 2015, 31(5): 38–41, 59.LI Xin, WANG Chunyang, YUAN Hui, et al. The modelingand simulation of radar smart jamming types[J]. AerospaceElectronic Warfare, 2015, 31(5): 38–41, 59.[7]ABOUELFADL A, SAMIR A M, AHMED F M, et al.Performance analysis of LFM pulse compression radar undereffect of convolution noise jamming[C]. The 2016 33rdNational Radio Science Conference, Aswan, Egypt, 2016:282–289.[8]GONG Shixian, WEI Xizhang, LI Xiang, et al. Mathematicprinciple of active jamming against wideband LFM radar[J].Journal of Systems Engineering and Electronics, 2015,26(1): 50–60. doi: 10.1109/JSEE.2015.00008.[9]HAO Honglian, ZENG Dazhi, GE Ping, et al. Research onthe method of smart noise jamming on pulse radar[C]. The2015 5th International Conference on Instrumentation andMeasurement, Computer, Communication and Control,Qinhuangdao, China, 2015: 1339–1342.[10]邱杰. 灵巧噪声干扰本质含义探讨[J]. 海军航空工程学院学报,2011, 26(5): 481–484.QIU Jie. Research on essential signification of smart noisejamming[J].  Journal of Naval Aeronautical andAstronautical University, 2011, 26(5): 481–484.[11]曾云健, 陶建锋, 孙青. 灵巧噪声在SAR中的干扰应用研究[J].航天电子对抗, 2009, 25(2): 39–41.ZENG Yunjian, TAO Jianfeng, and SUN Qing. Applicationand research of smart noise in jamming SAR[J]. AerospaceElectronic Warfare, 2009, 25(2): 39–41.[12]吕波, 冯起, 张晓发, 等. 对SAR的随机脉冲卷积干扰研究[J].中国电子科学研究院学报, 2008, 3(3): 276–279.LV Bo, FENG Qi, ZHANG Xiaofa, et al. Study of randompulse convolution jamming to SAR[J]. Journal of CAEIT,2008, 3(3): 276–279.[13]表 1 干扰增益对比Tab. 1   The comparison of jamming gain干扰样式干扰参数 干扰增益(dB)噪声模板的快时间宽度(s) 噪声模板的慢时间单边带宽(Hz) 理论值 实验值传统噪声卷积调制干扰4:710 74007.5467 7.80622:410 740010.2765 10.74681:210 740013.0147 13.6788本文所提干扰18:910 62834.0968 3.94709:410 79411.2159 11.81952:410 71923.5414 23.2726第5期 黄大通等：基于噪声卷积调制的SAR虚假信号生成新方法 905降佳伟, 吴彦鸿, 王宏艳, 等. 一种对SAR-GMTI的密集假目标干扰技术[J]. 舰船电子对抗, 2016, 39(1): 9–14, 88.JIANG Jiawei, WU Yanhong, WANG Hongyan, et al. Adense false target jamming technique against SAR-GMTI[J].Shipboard Electronic Countermeasure, 2016, 39(1): 9–14, 88.[14]周传晟. 合成孔径雷达卷积干扰技术的研究[J]. 舰船电子对抗,2013, 36(2): 6–10.ZHOU Chuansheng. Study of convolution jammingtechnology of synthetic aperture radar[J]. ShipboardElectronic Countermeasure, 2013, 36(2): 6–10.[15]沈爱国, 姜秋喜. 脉冲卷积干扰技术在超宽带SAR中的应用分析[J]. 电光与控制, 2010, 17(2): 39–42.SHEN Aiguo and JIANG Qiuxi. Analysis of pulseconvolution jamming technique applied to counter UWBSAR[J]. Electronics Optics & Control, 2010, 17(2): 39–42.[16]黄洪旭, 黄知涛, 周一宇. 对合成孔径雷达的移频干扰研究[J].宇航学报, 2006, 27(3): 463–468.HUANG Hongxu, HUANG Zhitao, and ZHOU Yiyu. Astudy on the shift-frequency jamming to SAR[J]. Journal ofAstronautics, 2006, 27(3): 463–468.[17]黄洪旭, 黄知涛, 周一宇. 对合成孔径雷达的随机移频干扰[J].信号处理, 2007, 23(1): 41–45.HUANG Hongxu, HUANG Zhitao, and ZHOU Yiyu.Randomly-shift-frequency jamming style to SAR[J]. SignalProcessing, 2007, 23(1): 41–45.[18]黄洪旭, 黄知涛, 吴京, 等. 对合成孔径雷达的步进移频干扰[J].宇航学报, 2011, 32(4): 898–902.HUANG Hongxu, HUANG Zhitao, Wu Jing, et al. Stepped-shift-frequency jamming to SAR[J]. Journal of Astronautics,2011, 32(4): 898–902.[19]YANG Kaizhi, YE Wei, MA Fangfang, et al. A large-scenedeceptive jamming method for space-borne SAR based ontime-delay  and  frequency-shift  with  templatesegmentation[J]. Remote Sensing, 2020, 12(1): 53. doi: 10.3390/rs12010053.[20]和小冬. 弹载合成孔径雷达干扰技术研究[D]. [博士论文], 电子科技大学, 2015.HE Xiaodong. Research on the jamming technologies ofmissile-borne synthetic aperture radar[D]. [Ph. D.dissertation], University of Electronic Science andTechnology of China, 2015.[21]YE Wei, RUAN Hang, ZHANG Shuxian, et al. Study ofnoise jamming based on convolution modulation to SAR[C].2010 International Conference on Computer, Mechatronics,Control and Electronic Engineering, Changchun, China,[22]2010: 169–172.房明星, 王杰贵, 雷磊. SAR雷达二维噪声卷积调制干扰研究[J]. 现代防御技术, 2014, 42(2): 139–144, 160. doi:10.3969/j.issn.1009-086x.2014.02.025.FANG Mingxing, WANG Jiegui, and LEI Lei. Study on 2Dnoise convolution modulation jamming to SAR[J]. ModernDefense Technology, 2014, 42(2): 139–144, 160. doi:10.3969/j.issn.1009-086x.2014.02.025.[23]ZHANG Yanbin. Technology of smart noise jamming basedon multiplication modulation[C]. 2011 InternationalConference on Electric Information and ControlEngineering, Wuhan, China, 2011: 4557–4559.[24]朱宁龙, 刘海亮. 移频对卷积调制灵巧噪声干扰的影响研究[J].舰船电子对抗, 2016, 39(2): 17–20. doi: 10.16426/j.cnki.jcdzdk.2016.02.005.ZHU Ninglong and LIU Hailiang. Research into theinfluence of frequency shift on smart noise jamming withconvolution modulation[J]. Shipboard ElectronicCountermeasure, 2016, 39(2): 17–20. doi: 10.16426/j.cnki.jcdzdk.2016.02.005.[25]林欢, 甘德云, 李宏, 等. 基于移频卷积的合成孔径雷达干扰技术研究[J]. 军械工程学院学报, 2014, 26(4): 47–50.LIN Huan, GAN Deyun, LI Hong, et al. Interferencetechnology research of synthetic aperture radar based onfrequency-shift convolution[J]. Journal of OrdnanceEngineering College, 2014, 26(4): 47–50.[26]吴晓芳, 邢世其, 王雪松, 等. 对合成孔径雷达的脉间分段移频干扰[J]. 航天电子对抗, 2010, 26(1): 53–57.WU Xiaofang, XING Shiqi, WANG Xuesong, et al.Interpulse subsection-shift-frequency jamming to SAR[J].Aerospace Electronic Warfare, 2010, 26(1): 53–57.[27]栾琳. 灵巧噪声干扰的建模仿真研究[D]. [硕士论文], 西安电子科技大学, 2009.LUAN Lin. Modeling, simulation and study of smart noisejamming[D]. [Master dissertation], Xidian University, 2009.[28]LIU Yongcai, WANG Wei, DAI Shaoqi, et al. Transmitterpower comparison between incoherent noise and coherentdeceptive jamming against synthetic aperture radar[C]. 2017IEEE Microwaves, Radar and Remote Sensing Symposium,Kiev, Ukraine, 2017: 279–284.[29]刘永才. 基于卷积调制的SAR有源欺骗干扰技术[D]. [硕士论文], 国防科学技术大学, 2013.LIU Yongcai. Active deceptive jamming against SAR basedon convolutional modulation[D]. [Master dissertation],National University of Defense Technology, 2013.[30]906 雷  达  学  报 第9卷作 者 简 介 黄大通(1993–)，男，博士研究生，主要研究方向为合成孔径雷达信号处理与对抗。E-mail：huangdatong68@163.com邢世其(1984–)，男，副研究员，主要研究方向为极化雷达成像、雷达信号处理以及合成孔径雷达对抗。E-mail: xingshiqi_paper@163.com刘业民(1982–)，男，博士，工程师，主要研究方向为雷达干扰与反干扰、极化雷达信号处理。E-mail: liuyemin2@163.com李永祯(1977–)，男，研究员，博士生导师，主要研究方向为极化雷达与电子对抗。E-mail: e0061@sina.com肖顺平(1964–)，男，教授，博士生导师，主要研究方向为极化信息处理、电子仿真与雷达目标识别、雷达电子对抗。E-mail: xiaoshunping_nudt@163.com第5期 黄大通等：基于噪声卷积调制的SAR虚假信号生成新方法 907