基于滑动多周期FFT的调频引信抗扫频干扰方法_周文.pdf

第 44 卷第 6 期2 0 2 3 年 6 月兵工学报ACTA AMAMENTAIIVol 44 No 6Jun2023DOI:10 12382/bgxb 2022 0119基于滑动多周期 FFT 的调频引信抗扫频干扰方法周文，郝新红，董二娃，陈彦君(北京理工大学 机电动态控制重点实验室，北京 100081)摘要:针对调频引信快速傅里叶变换(FFT)谐波包络提取方法对抗扫频干扰效果差的问题，提出一种滑动多周期 FFT 处理方法。在分析扫频干扰作用下调频引信的失效机理的基础上，提出通过多周期 FFT 进行回波能量的相干积累，实现对干扰的有效抑制，并在相邻处理窗口间滑动更新处理数据，进一步提取预定炸高对应的差频谐波包络进行门限判决。仿真及实测结果表明:滑动多周期 FFT 处理方法能在扫频干扰下有效提取目标谐波包络特征，相较于传统单周期 FFT 处理方法，显著改善峰值旁瓣比，提高了调频引信抗扫频干扰的能力。关键词:调频引信;扫频干扰;相干积累;谐波包络;峰值旁瓣比中图分类号:TJ43+4.1文献标志码:A文章编号:1000-1093(2023)06-1744-10收稿日期:2022-03-02基金项目:国家自然科学基金项目(61871414);国防科技基础加强计划项目(2019-JCJQ-ZD-324)Anti Frequency Sweeping Jamming Method for FMFuze Using Sliding Multi cycle FFTZHOU Wen，HAO Xinhong，DONG Erwa，CHEN Yanjun(Science and Technology on Electromechanical Dynamic Control Laboratory，Beijing Institute ofTechnology，Beijing 100081，China)Abstract:This study proposes a sliding multi-cycle fast Fourier transform(FFT)processing method toaddress the problem of poor performance of the FFT harmonic envelope extraction method used in FM fuzeunder frequency sweeping jamming conditionsThe proposed method aims to achieve effectiveinterference suppression by coherently accumulating echo energy through multi-period FFT processingThe processing data is updated by sliding between adjacent processing windows，and the differencefrequency harmonic envelope corresponding to the predetermined blast height is further extracted forthreshold judgment The simulated and measured results show that compared with traditional single-cycleFFT processing method，the new method can effectively extract the target harmonic envelopecharacteristics under frequency sweeping jamming，significantly improve the peak-to-side lobe ratio，andenhance the anti-frequency sweeping jamming performance of FM fuzeKeywords:FM fuze;frequency sweeping jamming;coherent accumulation;harmonic envelope;peak-to-side lobe ratio0引言引信作为一种弹药终端毁伤控制系统，直接决定了武器系统能否发挥出最佳毁伤效能，在现代战争及地区防御中占据重要地位1。作为引信的一个重要组成类型，调频引信由于其具有定距精度高、第 6 期基于滑动多周期 FFT 的调频引信抗扫频干扰方法结构简单等特点，已在常规弹药中得到广泛应用2。而在日趋复杂的战场电磁环境中，诸如扫频干扰3 4 等人为电子干扰手段对调频引信等无线电体制引信造成了巨大威胁5 7。在无线电引信抗干扰方面，近些年已有学者展开了相关研究8 14。文献 8 提出基于熵特征与支持向量机的信号识别方法，通过提取检波信号的香农熵和奇异谱熵，结合支持向量机分类器，达到了较高的目标信号识别正确率;文献 10 利用调频信号在分数阶域的聚集性，提出一种基于分数阶傅里叶变换的调频引信定距方法。但上述方法所需计算量较大，对计算资源要求较高。文献 13提出基于时序及相关检测的抗干扰方法，依据弹目接近过程中谐波峰值从高到低依序出现的规律，进行两路谐波时序逻辑判决，但在低信噪比条件下表现不够理想。针对调频引信抗扫频干扰问题，本文提出一种滑动多周期快速傅里叶变换(FFT)谐波包络提取的处理方法。该方法通过多周期 FFT 相干积累抑制扫频干扰信号，并在相邻处理窗口间滑动更新一个调制周期长度的数据，降低差频信号不规则区带来的影响，能在扫频干扰下有效提取目标差频信号的谐波包络，明显优于调频引信单周期 FFT 谐波包络提取方法。1扫频干扰下调频引信失效机理1.1调频引信模型调频引信的基本工作过程为:调制信号控制压控振荡器产生频率周期线性变化的射频信号，一路通过发射天线辐射出去，另一路送入本地混频器。接收天线接收到的回波信号通过低噪声放大器后与混频器中的本地参考信号进行混频处理，经模数转换器(ADC)采样后通过 FFT 提取差频的不同谐波包络，根据差频频率 f 的大小与距离的对应关系，完成预定炸高的设置。常见调频引信的主要工作过程如图 1 所示。以锯齿波线性调频引信为例，设发射信号幅值为 At，初始频率为 f0，调制周期为 Tm，调制频率为fm，调制带宽为 B，调频斜率为，发射信号 xt(t)可表示为xt(t)=Atn=1rectt (n 1)TmTmej 2f0t+t (n 1)Tm2(1)图 1调频引信工作原理Fig 1Working mechanism of the FM fuze式中:n 为当前时刻对应的周期数。设真实目标回波经过的路径时延为，信号幅值为 Ar，则回波信号 xr(t)可表示为xr(t)=Arn=1rectt (n 1)TmTmej 2f0(t )+t (n 1)Tm2(2)图 2单周期 FFT 处理结果Fig 2Single-cycle FFT processing results忽略非规则区影响，经过混频器处理，单个调制周期内的差频信号 xf(t)可表示为xf(t)=AtAr2ej2(f0+t 12)2(3)对单个周期内的差频信号进行 FFT 处理，结果如图 2 所示，路径时延 在单个周期内可视为常数，所得各次谐波系数的模可近似表示为5471兵工学报第 44 卷|XTm(f)|=AtAr(21 T)msinc(f)(Tm)(4)进一步可得距离分辨率 为=c2B(5)根据式(4)，差频频率随路径时延 线性变化，在不同距离处，离对应差频频率最近的谐波系数幅值包络最大;对于单次谐波系数，其距离维包络表现为 sinc 函数，sinc 函数零点为相关峰的峰值点，位于差频频率所对应弹目距离处。由式(5)可知，相关峰宽度与距离分辨率 相关，其大小由发射信号调制带宽决定。1.2扫频干扰模型干扰机扫频带宽能覆盖引信工作频带，扫频干扰信号的频率在扫频带宽内按一定规律步进，并在每个步进频点上驻留一定时间，当驻留频点进入引信工作带宽内，将对引信产生干扰效果15 18。如图 3 所示，干扰机采取扫频干扰方式，f 为信号频率，t 为当前时间，扫频初始频率为 fj1，相邻频点间步进频率为 fj，总扫频点数为 G，每个频点的驻留时间为 tj，扫频周期为 Tj，则 1 个扫频周期内的干扰信号频率 fj可表示为fj=fj1+(g 1)fj，g=1，2，G(6)进一步可得干扰信号 xj(t):xj(t)=Ajej(2fjgt)Pg(t)(7)Pg(t)=1，(g 1)tj tgtj0，其他(8)图 3扫频干扰模型Fig 3Frequency sweep jamming model1.3扫频干扰下失效机理分析一般来说，扫频干扰的驻留时间 tj均大于引信信号调制周期，扫频干扰信号进入引信接收通道，并于本地参考信号混频后输出信号在 1 个调制周期内可表示为xfj(t)=AtAj2ej(2(f0 fj)t+(t (n 1)Tm)2)(9)fj=f0 fj+(t (n 1)Tm)(10)设低通滤波器的截止频率为 fLP，经过低通滤波器后，能进入信号处理系统的干扰信号频带范围为 fLPfjfLP(11)如图 4 所示，当干扰信号频率进入引信信号带宽内，干扰信号与本地参考信号混频后，所得差频信号可看作一段线性调频信号，有效长度由低通滤波器的截止频率 fLP决定。图 4扫频干扰下失效机理Fig 4Failure mechanism under frequencysweeping jamming设有效干扰段的起始时刻为 tjL，结束时刻为tj，有效干扰时间 T=tj tjL，进行变量代换 t=t1+T2，则有效干扰段的频谱(见图 5)可表示为U(f)=AtAj21TT2T2ej2(t1)2 e j2ft1dt1=AtAj212槡Te jf2/c(v1)+c(v2)+j s(v1)+s(v2)(12)式中:c(v)=v0(cosx2)2dx，s(v)=v0(sinx2)2dx(13)为菲涅尔积分公式，积分限为v1=2槡(T2f)，v2=2槡(T2+f)(14)由图 5 可以看出，扫频干扰信号差频频谱覆盖整个引信信号处理带宽，由于干扰信号能量一般高于真实目标回波信号，当扫频干扰信号进入后，将对真实目标回波差频信号的频谱产生压制式遮盖效6471第 6 期基于滑动多周期 FFT 的调频引信抗扫频干扰方法图 5干扰信号差频频谱Fig 5Spectrum of beat signals under frequencysweeping jamming果，影响调频引信基于提取谐波包络特征的炸高判定效果，从而引起早炸或瞎火。2滑动多周期 FFT 处理方法2.1多周期 FFT 处理方法利用发射信号在各个调制周期内具有强相干性，而干扰信号与本地参考信号非相干的特点，提出多周期 FFT 处理方法。考虑扫频干扰信号的影响，位于相关器两端的本地参考信号 Sl(t)与混合接收信号 Sr(t)可分别表示为Sl(t)=xt(t)(15)Sr(t)=xr(t)+xj(t)(16)对 Q 个周期长度的相关器输出信号做傅里叶变换，各次谐波的系数由 XQTm(f)表示:XQTm(f)=1QTmQTm0Sl(t)Sr(t)e j2ftdt=1QTmQTm0 xt(t)xr(t)e j2ftdt+QTm0 xt(t)xj(t)e j2ftdt(17)进行变量替换 t=t2+(q 1)Tm，根据积分的可加性，式(17)可改写为XQTm(f)=1QTmQq=1Tm0 xt(t2+(q 1)Tm)xr(t2+(q 1)Tm)e j2f(t2+(q 1)Tm)dt2+Tm0 xt(t2+(q 1)Tm)xj(t2+(q 1)Tm)e j2f(t2+(q 1)Tm)dt2(18)设式(18)中第 1 项为 Xtr(f)，由于发射信号与真实回波信号在各个调制周期内具有强相干性，Xtr(f)取模进一步简化为|Xtr(f)|=AtAr(21 T)msin c(f)(Tm)sin fQTmQsin(fTm)|XTm(f)|(19)设式(18)中第 2 项为 Xtj(f)，由于干扰信号与本地参考信号非相干，不同周期内的同一频率分量间存在相位差，因此在进行多周期积累时存在失配损失，Xtj(f)取模进一步简化为|Xtj(f)|1QTmQq=1Uq(f)|U(f)|(20)当且仅当 Q=1 时，等号成立，Uq(f)表示第 q 个周期的差频信号频谱。由以上分析可知，经过多周期 FFT 处理，真实目标回波差频信号分量的主要特征没有变化，固定次谐波的距离维包络依旧表现为 sin c 包络，但在频率轴方向上频率分辨率提高到单周期 FFT 处理结果的 Q 倍;干扰对应的差频分量在不同周期间存在相位失配，因此会抑制干扰效果。2.2滑动数据更新方法考虑到在高速相对运动过程中，一个距离分辨单元内的信号调制周期数有限，若每次处理取多个周期数据，下一次处理数据全部更新，这会对谐波包络提取造成降采样的效果，将不利于后续判决处理。为了解决这个问题，借鉴短时傅里叶变换时频分析方法的思想，提出滑动多周期 FFT 谐波包络提取方法。主要过程如图 6 所示，由于每个周期中的差频信号均存在一段非规则区，为保证每个数据处理窗口中非规则区带来的影响水平相同，因此将滑动更新数据的时间长度设置为调制周期的整数倍。为了得到更为细致的谐波包络，滑动更新数据的时间长度应尽可能短，因此，相邻处理窗口中原始数据滑动更新 1 个调制周期长度的数据较为合理。当 Q=1时，滑动多周期 FFT 方法退化为传统单周期 FFT方法。设整个交会过程中采集的目标回波中频输出信号长度为 Tall，进行 Q 个调制周期长度的滑动多周期 FFT 处理后(见图 7)，数据处理窗口总数 i 由式(21)决定:i=Tall Tm(Q 1)Tm(21)7471兵工学报第 44 卷图 6滑动多周期 FFT 处理流程Fig 6Sliding multi-cycle FFT processing flow式中:?表示向下取整。整个数据处理过程所得谐波频率矩阵 X(f)可表示为X(f)=X1(f)X2(f)X3(f)Xi(f)(22)图 7滑动多周期 FFT 处理结果Fig 7Sliding multi-cycle FFT processing results在 Q 取值较小时，可将数据处理窗口内的时延 视为常数，此时 Xi(f)可由式(19)近似表示。当Q 值较大或弹目相对速度较大时，式(19)的准确性将变差。由图 7(a)可以看出，真实目标回波中频输出信号经过 8 个周期长度滑动多周期 FFT 处理后，所得时频矩阵的分布特征与式(19)描述一致，距离分辨单元在时频矩阵上线性分布，呈现出一个明显的时频脊。随着相对距离减小，对应差频频率随之减小，相邻分辨单元对应的谐波频率相差 8 次谐波;由图 7(b)可以看出，在未施加扫频干扰的条件下，8 个周期长度滑动多周期 FFT 处理所得距离维谐波包络依然为 sin c 函数的形状。3仿真及实验验证3.1仿真验证为了验证上述分析，增加扫频干扰进行仿真验证，仿真参数设置如表 1 所示。表 1仿真参数Table 1Simulation parameters参数数值调制频率/kHz100调制带宽/MHz200弹目距离/m15 0JS/dB15扫频带宽/MHz400注:JS 为干信比。对仿真所得差频数据分别进行单周期 FFT 及本文所提滑动多周期 FFT 方法处理，分别绘制仅目标信号作用下、仅干扰信号作用下及目标与干扰信号同时作用下的交会过程时频矩阵和距离维包络，所得结果如图 8 和图 9 所示。由图 8(a)及图 8(b)可以看出，仅目标信号作用下，差频信号由单周期 FFT 处理后，所得时频脊较为明显，距离维相关峰高于背景噪声;由图 9(a)及图 9(b)可以看出，进行 8 个周期长度滑动多周期FFT 处理后，背景噪声得到了有效抑制，距离维相关峰的 sin c 包络特征更加显著。由图 8(c)及图 8(d)可以看出，仅干扰信号作用下，差频信号由单周期 FFT 处理后，全局时频分布中包含大量异常峰值，峰值大小高于目标信号相关峰水平;由图 9(c)及图 9(d)可以看出，进行 8 个周期长度滑动多周期 FFT 处理后，干扰信号带来的异常峰值已低于目标信号相关峰水平。由图 8(e)可以看出，在目标信号及干扰信号的同时作用下，由单周期 FFT 处理后所得时频矩阵，其各次谐波中均混入大量由干扰信号带来的8471第 6 期基于滑动多周期 FFT 的调频引信抗扫频干扰方法图 8仿真数据单周期 FFT 处理结果Fig 8Single-cycle FFT processing results of simulation data虚假峰值，真实目标回波所对应的时频脊被遮盖;由图 8(f)可以看出，9 m 对应谐波包络中真实目标相关峰掩盖在干扰信号中，对门限判决产生不利影响。由图 9(e)可以看出，进行 8 个周期长度滑动多周期 FFT 处理后，时频矩阵中由干扰信号带来的虚假峰值被大幅压低，可以清楚观察到真实目标回波所对应的时频脊;由图 9(f)可以看出，通过滑动多周期 FFT 有效压制干扰信号后，距离维谐波包络中再次出现明显的相关峰。峰值旁瓣比(PSL)是衡量无线电信号抗干扰能力的一项重要指标19 20，其数值越小，表示抗干扰性能越好。为进一步进行抗干扰效果指标度量，采用 PSL 对在不同 JS 条件下不同周期长度滑动多周期 FFT 处理结果的抗干扰性能(见图 10)进行表征:9471兵工学报第 44 卷图 9仿真数据滑动多周期 FFT 处理结果Fig 9Sliding multi-cycle FFT processing results of simulated dataPPSL(=10lgPsideP)peak(23)式中:PPSL为 PSL;Pside为最大旁瓣功率;Ppeak为峰值响应功率。由图 10 可以看出:当 JS 较低时(如JS=9 dB)，在一定范围内，随着周期数 Q 的增加，在不同 JS 条件下的 PSL 均有明显下降;当 Q 值较大时，PSL 下降幅度有限;当 JS 较高时(如JS=21 dB)，随着周期数 Q 的增加，PSL 的改善幅度较为有限，当 Q 取值高于 8 时，PSL 没有得到明显改善。从前文分析可知，当 Q 取值过大，Q 个周期内的回波数据相干性变差，会影响滑动多周期 FFT 方法的处理效果，根据仿真结果，建议周期数 Q 的取值范围为 2 16。Q 的增大也会带来计算量的增加，对实时性会造成一定影响，如表 2 所示，给出了不同处理周期数 Q 对 应 的 计 算 时 间(基 于 美 国 Xilinx 公 司Spartan-6 系列 FPGA 的测试结果)。一般来说，引0571第 6 期基于滑动多周期 FFT 的调频引信抗扫频干扰方法图 10滑动多周期 FFT 处理后 PSL 结果Fig 10PSL results of sliding multi-cycleFFT processing表 2不同周期数所需计算时间Table 2Calculation time required for different cycles参数Q12481632数据长度/B2565121 0242 0484 0968 192计算时间/ms0.060.130.280.571.122.18信信号处理的延迟应控制在毫秒量级，因此，综合考虑最优抗干扰性能与实时性的平衡，本文选用 Q 值为 8。3.2实验验证为了进一步验证基于公式推导及仿真分析的正确性，现基于某引信样机进行实测数据采集及分析。实验场景如图 11 所示。主要实验过程为:设置目标散射截面积(CS)为 1 m2的金属平板，运动平台从12 m 运动到8 m位置处，交会速度 1 m/s，采集交会过程引信差频输出信号。图 11实验场景示意图Fig 11Schematic diagram of the test setup将采集的差频数据分别进行单周期 FFT 及本文所提滑动多周期 FFT 方法处理，提取并绘制交会过程时频矩阵及 9 m 炸高对应谐波的距离维包络，所得结果如图 12 和图 13 所示。图 12实测数据单周期 FFT 处理结果Fig 12Single-cycle FFT processing results ofmeasured data由图 12 和图 13 可以看出:由单周期 FFT 处理后所提取的时频矩阵无法看出主要频率分量随距离的线性变化关系，相关峰不明显，低于由扫频干扰带来的干扰信号谱线高度。采用本文所提滑动多周期FFT 处理方法，提取的时频分布图可以明显看出时频线性对应关系，并且 9 m 炸高谐波包络的相关峰特征明显，有利于下一级判决处理。4结论为提高调频引信对抗扫频干扰的能力，本文在分析扫频干扰下调频引信失效机理的基础上，提出一种滑动多周期 FFT 处理方法，通过仿真及实测结果验证了该方法的有效性。得出的主要结论如下:1)扫频干扰信号能覆盖调频引信工作带宽，对真实目标回波差频信号的频谱产生压制效果，对目标判决造成不利影响。2)通过在数据处理窗口进行多周期 FFT 相参积累，可以实现对扫频干扰的有效抑制，PSL 得到1571兵工学报第 44 卷图 13实测数据滑动多周期 FFT 处理结果Fig 13Sliding multi-cycle FFT processingresults of measured data显著改善。3)在相邻处理窗口中滑动更新部分数据，可以获得更为细致的距离维包络，有利于下一步目标判决处理。4)本文所提方法的算法复杂度低，能够满足引信的实时性要求。参考文献(eferences)1崔占忠，宋世和，徐立新 近炸引信原理M 北京:北京理工大学出版社，2009:68 75CUI Z Z，SONG S H，XUN L X Principle of proximity fuze M Beijing:Beijing Institute of Technology Press，2009:68 75(inChinese)2赵惠昌 无线电引信设计原理与方法M 北京:国防工业工业出版社，2012:35 46ZHAO H C Principles and methods of radio fuze design M National Defense Industry Press，2012:35 46(in Chinese)3CHOI D Y，KIM W K，KIM J H，et al Performance of analogand digital modulation schemes under sweep jammingCProceedings of the 2016Eighth International Conference onUbiquitous and Future Networks Vienna，Austria:IEEE，2016 4CHOI Y S，LEE S Y，CHOI H H，et al The improved spatialnuller with frequency swept jammerCProceedings of 2014IEEE/IONPosition，LocationNavigationSymposiumMonterey，CA，US:IEEE，2014 5金文，武海东，沙作金，等 无线电引信的抗干扰性能评估方法研究J 制导与引信，2020，41(2):11 16JIN W，WU H D，SHA Z J，et al esearch on anti-interferenceperformance evaluation method of radio fuze J Guidance Fuze，2020，41(2):11 16(in Chinese)6熊久良，武占成，孙永卫，等 无线电引信互扰效应规律研究 J 电波科学学报，2015，30(5):903 909XIONG J L，WU Z C，SUN Y W，et alStudy on mutualinterference effect of radio fuzeJ Journal of adio Science，2015，30(5):903 909(in Chinese)7张星宇，齐杏林，李黎明 国外无线电引信发展历程及未来研究方向综述J 飞航导弹，2021(12):152 157ZHANG X Y，QI X L，LI L M eview on the developmenthistory and future research direction of foreign radio fuzeJAeronautical Missile，2021(12):152 157(in Chinese)8黄莹，郝新红，孔志杰，等，基于熵特征的调频引信目标与干扰信号识别J 兵工学报，2017，38(2):254 260HUANG Y，HAO X H，KONG Z J，et al Target and jammingsignal recognition of FM fuze based on entropy featureJ ActaArmamentarius，2017，38(2):254 260(in Chinese)9DAI J，HAO X H，LI Z，et al Adaptive target and jammingrecognition for the pulse doppler radar fuze based on a time-frequency joint feature and an online-updated naive bayesianclassifier with minimal riskJ Defence Technology，2021，18(3):457 466 10岳凯，郝新红，栗苹，等 基于分数阶傅里叶变换的线性调频引信定距方法J 兵工学报，2015，36(5):801 808YUE K，HAO X H，LI P，et al A range method for chirp fuzebased on fractional Fourier transformJ Acta Armamentarii，2015，36(5):801 808(in Chinese)11DAI J，HAO X H，LIU Q，et al epeater jamming suppressionmethod for pulse Doppler fuze based on identity recognition andchaotic encryptionJ Defence Technology，2020，17(3):1002 1012 12黄海燕 连续波无线电引信信号的识别D 成都:电子科技大学，2014HUANG H Y ecognition of continuous wave radio fuze signal D Chengdu:University of Electronic Science and Technologyof China，2014(in Chinese)13孔志杰，郝新红，栗苹，等 基于时序及相关检测的调频引信抗扫频干扰方法J 兵工学报，2017，38(8):1483 1490KONG Z J，HAO X H，LI P，et al Anti-sweep interferencemethod of FM fuze based on timing and correlation detectionJ2571第 6 期基于滑动多周期 FFT 的调频引信抗扫频干扰方法Acta Armamentarii，2017，38(8):1483 1489(in Chinese)14李志强 连续波多普勒无线电引信目标信号识别方法研究D 北京:北京理工大学，2014LI Z Qesearch on target signal identification method ofcontinuous Wave Doppler radio fuze DBeijing:BeijingInstitute of Technology，2014(in Chinese)15李泽，闫晓鹏，栗苹，等 扫频式干扰对调频多普勒引信的干扰机理研究J 兵工学报，2017，38(9):55 61LI Z，YAN X P，LI P，et al Jamming mechanism of frequencysweep jamming to FM doppler fuzeJ Acta Armamentarii，2017，38(9):55 61(in Chinese)16夏彩杰，王爱华，安建平 扫频干扰的自适应时变抑制研究J 计算机工程，2008，34(4):37 39XIA C J，WANG A H，AN J P esearch on adaptive time-varying suppression of frequency sweeping interferenceJ Computer Engineering，2008，34(4):37 39(in Chinese)17ALLOND S，STAK W，ALI M，et alInterference inAutomotiveadarSystems:Characteristics，mitigationtechniques，and current and future researchJ IEEE SignalProcessing Magazine，2019，36(5):45 59 18陈齐乐，郝新红，闫晓鹏，等 基于谐波系数幅值平均的复合调制引信抗扫频式干扰方法J 北京航空航天大学学报，2020(7):1317 1324CHEN Q L，HAO X H，YAN X P，et al Anti sweep jammingmethod of hybrid modulation fuze based on harmonic coefficientamplitude averaging JJournal of Beijing University ofAeronautics and Astronautics，2020(7):1317 1324(inChinese)19SALEH M，OMA S M，GIVEL E，et al A modified steppedfrequency phase coding radar waveform designed for the frequencydomain algorithmJ Digital Signal Processing，2019，88:116129 20JI J，ZHENG S L，ZHANG X M Pre-distortion compensationfor optical-based broadband LFM signal generation systemJOptics Communications，2019，435:277 282作者简介:周文(1997)，男，博士研究生，研究方向为无线电引信信号处理和信息对抗理论。E-mail:zwdavidbit163 com郝新红(1974)，通信作者，女，教授，博士生导师，研究方向为目标探测理论、引信信号处理和信息对抗理论。E-mail:haoxinhong bit edu cn3571