雷达抗距离-速度联合欺骗干扰的自适应迭代滤波方法
【专利摘要】本发明提供一种雷达抗距离?速度联合欺骗干扰的自适应迭代滤波方法，分别对目标和干扰的距离维和脉冲维进行估计，然后通过迭代来提高算法的精度。在干扰相对于目标延迟已知个脉冲的前提下，将自适应滤波的思想运用到对干扰和目标的距离?多普勒平面的估计上，为了提高估计精度，同时引入了迭代步骤。本发明针对一种少有干扰形式：对抗距离?速度联合，不需要对波形进行设计，实施较为方便，并且在多目标多干扰的场景下仍能对干扰有效压制并能对目标正确检测。仿真结果表明，本发明采用的方法相对于传统的匹配滤波和多普勒处理有明显的性能优越性。
【专利说明】
雷达抗距离-速度联合欺骗干扰的自适应迭代滤波方法
技术领域
[0001] 本发明属于雷达抗干扰技术，特别设及自适应迭代滤波技术。
【背景技术】
[0002] 在电子战中，随着数字射频存储器DRFM技术在干扰机中的应用，使得电子反对抗 ECCM变得更加困难。装备了数字射频存储器的干扰机在截获了雷达发射信号后进行，复制， 存储，调制并转发。由于干扰信号复制了发射信号的模型，所W可W在雷达接收机中获得处 理增益，运样在雷达的显示屏上将出现多个目标，在不具备电子反对抗能力的雷达系统中， 运使得雷达无法辨别出真实目标，对传统雷达系统的生存构成了极大挑战。干扰机对截获 的雷达发射信号的调制可W改变时延，或者对多普勒频率变换，W分别产生距离假目标和 速度假目标，或者同时对两者进行调制，运样产生的即是距离-速度联合欺骗干扰。显然，后 者对雷达的威胁更大。为了确保雷达对真实目标的检测，W及后续的跟踪，提高雷达的抗距 离-速度联合欺骗干扰能力具有重要的理论价值和实际意义。
[0003] -种被广泛研究的抗干扰方法是对雷达发射信号的波形进行设计。例如设计雷达 发射的脉冲是互相正交的，假设干扰机需要几个脉冲重复周期来识别和捕获发射信号，那 么反射回来后的干扰信号经过当前脉冲重复周期的匹配滤波器后会衰减，因为设计波形为 脉冲间相互正交。但是，在动目标显示(MTI)和脉冲多普勒雷达的距离多普勒处理中，上述 方法会产生旁瓣变异，极大地限制了运一方法的有效性。文献[Jindong Zhang ,Daiyin Zhu,Gong Zhang,New Anti-velocity Deception Jamming Technique using Pulses with Adaptive Initial Phases. IEEE Transactions on Aerospace and Electronic Systems ,vol. 49,no. 2 ,pp. 1290-1300,2013]提出了一种通过设计自适应初始相位的编码 信号来抗速度欺骗干扰的方法，但是此方法仅考虑了单个目标，在多目标场景下，此方法将 丧失抗干扰能力。从公开发表的文献来看，针对距离-速度联合欺骗干扰抑制的方法还未有 研究。
【发明内容】
[0004] 本发明所要解决的技术问题是，提供一种通过自适应迭代滤波的方法来对抗距 离-速度联合欺骗干扰的方法。
[0005] 本发明为解决上述技术问题所采用的技术方案是，一种雷达抗距离-速度联合欺 骗干扰的自适应迭代滤波方法，包括W下步骤：
[0006] 1初始化步骤:设雷达对目标照射的驻留时间为M个脉冲，首先在快时间维对信号 进行匹配滤波，然后在慢时间维做傅里叶变换得到的距离-多普勒平面作为真实目标的距 离-多普勒平面初始值文|；*1和干扰信号的距离-多普勒平面初始值?皆>;
[0007] 2距离维估计迭代步骤：
[000引2-1 :构造真实目标的快时间变换矩阵GT,m与对角矩阵巧;,;,1) [d :
[0009]
[0010]
[0011] 其中U表示迭代次数，diag( ?)为求矩阵的对角线元素，ST,m=Sm为真实目标第m个 脉冲的回波信号，Sm=^m(I)，Sm(2)，…，Sm(N)]T是对信号模型的采样结果，N为采样点数， (.)T表示转置。，其中，C也；；1>表示真实目标第U-I次迭代距离像估计值，1 = 1，2,…，L，L表示 距离单元的个数，m= 1，…，M，M表示脉冲个数；
[0012] 构造干扰信号的快时间变换矩阵Gj,m与对角矩阵
[0013]
[0014]
[001引SJ,m表示干扰信号第m个脉冲的回波信号，其中，诚"1嗦示干扰信号第U-I次迭代距 离像估计值；
[0016] 2-2:得到真实目标在不同取值范围1里的距离像估计，/ = 1，:…，W-I， /二尼…，1-(媒-1)，/=1-的-2)，...，/.;得到干扰信号在不同取值范围1里的距离像估 计：I二V.'，'V-I，I二 N，…丄I 二 L-{N-2)，.-.，L'.
[0017] 步骤2-3:构造距离维上真实目标矩阵Al勺日干扰信号矩阵Af;
[0020] 3脉冲维估计迭代步骤：[0021] 3-1:估计真实目标第1个距离单元对应第k个多普勒单元的值乂
[001 引
[0019]
[0022]
[0023]
[0024] 估计干扰信号第I个距离单元对应第k个多普勒单元的值乂f ;
[0027]其中，悼')[/]和的'中]分别是!Alf和{碱聲的第巧U，F=[fi，f2，...，fM ]，fm=[exp(j23imfi)，exp( j23imf2)，…，exp( j化mfK)]T是KX 1的归一化多普勒相位矢量，Fk为F的第k行， I是单位矩阵，（.）H表示共辆转置；[00%] 3-3:令1 = 1，…，L得到真实目标的距离-多普勒平面估计：
[0025]
[0026]
[0029]
[0030] ;-多普勒平面估计：
[0031]
[0032] 其中
'
[0033] 步骤5:重复步骤2和步骤3，直到达到停止迭代条件。
[0034] 本发明分别对目标和干扰的距离维和脉冲维进行估计，然后通过迭代来提高算法 的精度。在干扰相对于目标延迟已知个脉冲的前提下，将自适应滤波的思想运用到对干扰 和目标的距离-多普勒平面的估计上，为了提高估计精度，同时引入了迭代步骤。
[0035] 本发明的有益效果是，提针对一种少有干扰形式:对抗距离-速度联合，不需要对 波形进行设计，实施较为方便，并且在多目标多干扰的场景下仍能对干扰有效压制并能对 目标正确检测。仿真结果表明，本发明采用的方法相对于传统的匹配滤波和多普勒处理有 明显的性能优越性。
【附图说明】
[0036] 图1为算法流程图。
[0037] 图2为传统匹配滤波抗距离-速度联合欺骗干扰估计结果。
[0038] 图3自适应迭代滤波法抗距离-速度联合欺骗干扰结果。
[0039] 3-1为自适应迭代滤波法对真实目标的距离-多普勒平面估计结果示意图。
[0040] 3-2为自适应迭代滤波法对干扰目标的距离-多普勒平面估计结果示意图。
【具体实施方式】
[0041] 本发明的具体实施步骤如图1所示：
[0042] 步骤1:脉冲多普勒雷达不断的发射脉冲信号，同时接收目标的回波信号，接收到 的回波信号不光包括目标回波，还包括干扰机产生的欺骗式干扰信号，并且干扰类型为距 离-速度联合欺骗式干扰，那么此时接着转到步骤2做下一步处理。
[0043] 步骤2:准备自适应迭代算法初始值。假设雷达对目标照射的驻留时间为M个脉冲， 首先在快时间维对信号进行匹配滤波，然后在慢时间维做傅里叶变换，运就是传统的脉冲 多普勒处理，此时可W得到的距离-多普勒品面作为下面迭代步骤的初始值
[0044] 其中，Xi;*哪則i分别表示真实目标和干扰的距离-多普勒平面初始值。
[0045] 步骤3:距离维估计步骤：
[0046] 步骤3-1:对脉冲维的估计做准备工作，构造真实目标的快时间变换矩阵Gt,根据 上一步的估计结果，构造对角矩阵11<;"4 W，其中：
[0047]
[004引
[0049] 对1 = 1，...，N-1，1=N，…，レ(N-1)，1 = レ(N-2)，...，L，在1位于不同的取值范围时 均做运样的处理。
[0050] 对于干扰信号的估计，重复上面一样的准备工作。
[0051] 步骤3-2:利用上一步已经得到的准备信息，分别得到在第U次迭代时，1在不同取 值范围里的距离像估计，即：
[0052] 4':,!,[/]，/ = 1，....jV-1，/ = W,…，1)，/ = L i.iV 2')，'..，么。同时也可^得到干 扰目标在不同取值范围里的距离像估计，即：[/]，/ =心..，W-1，/ = %…，王-(F-1)， /二I-(W-2)，...，I。对于真实目标的估计，具体如下：
[0化3]
[0化4]
[0化5]
[005引其中，gT,m,康示GT,m的第巧lj，ym[l] = [ym[l]，…，ym[l + (N-l)]]表示回波信号第m个 脉冲对应的从第1个距离单元到第1+(N-1)个距离单元的采样结果，I表示单位矩阵。
[0059]同样对于干扰目标，具体计算方式如下：
[0063] 其中，gj,m,康示Gj,m的第巧U，Sj,m=ST,m-康示干扰滞后信号i个脉冲重复周期。
[0064] 步骤3-3:利用步骤3-2的估计结果构造矩阵A护和。
[00化]其中；
[0068] 步骤4:脉冲维估计步骤：[0069] 步骤4-1:由步骤3的估计结果结合上一次的距离多普勒平面估计结果义!"
[0066]
[0067] 和来估计第1个距离单元对应第k个多普勒单元的值，即;同理也可W得
[0070] 到干扰目标的估计结果為")[a]，即：
[00对的'1 [/]和 P(;') [/]分别是{a巧『和仪f 的第1 列，F = [f 1，f2，...，fM]，fm= [exp(j如
[0071]
[0072]
[0073]
[0074] mf 1)，exp (j化mf2)，…，e邱（j23imfK) ]T是K X 1的目标多普勒相位矢量，Fk为F的第k行，I为单 位矩阵，（.）H表示转置共辆；
[0076] 步骤4-2:令k = 1，…，K，可W得到：
[0077]
[007引
[0079]
[0080] 步骤4-3:令1 = 1，…，L，可W得到真实目标与干扰的距离-多普勒平面估计：
[0081]
[0082]
[0083] 步骤5:重复步骤3和步骤4，直到达到停止迭代的要求，如迭代次数达到10次。
[0084] 仿真验证及分析 [00化]仿真场景；
[0086]假设一个CPI发射波形M+i = 33,干扰滞后发射信号i = l个PRT。发射信号的离散采 样点数N=32的随机相位波形。多普勒采样点数K=128,距离多普勒平面的估计范围是距离 单元个数L= 100。选取一个CPI的最后1 = 32个脉冲信号进行处理。中屯、频率f。= IG化W及 脉冲重复间隔PRT = Ims接收机引入的加性噪声功率为1。算法迭代次数为10次。
[0087]假设场景中目标为Swerling-O类型，真实目标和假目标的位置、速度和幅度如下 表：
[00则表1真实目标参数
[0089]
[0090]
[0091]
[0092]
[0093] 仿真分析：
[0094] 从图2可W看出，传统的匹配滤波和多普勒处理估计真实目标的距离-多普勒平面 结果中，不能有效识别真实目标。
[00M]从图3中利用自适应迭代滤波算法得到的距离-多普勒平面估计可W看出：
[0096] (1)自适应迭代滤波算法算法能够有效抑制距离-多普勒旁瓣和互干扰，估计真实 目标和干扰假目标的距离-多普勒平面。
[0097] (2)在图3-2中对干扰假目标的距离-多普勒估计中，距离维的旁瓣几乎没有扩展， 而多普勒维上的旁瓣扩展明显。运是由于距离维的分辨率与带宽有关，而多普勒维的分辨 率与脉冲的个数M有关，影响距离维和多普勒维的旁瓣扩展的参数不同，因此估计结果中距 离维和多普勒维的旁瓣扩展效果不同。
[0098] (3)图3-2中图中的干扰假目标估计比图3-1中图中真实目标估计的多普勒维旁瓣 扩展效果明显，因为干扰假目标的功率远大于真实目标，干扰假目标的估计结果中的旁瓣 的水平明显高于噪声基底，因此干扰假目标的旁瓣明显。而真实目标估计中的旁瓣和噪声 基底水平相近，因此旁瓣不明显。
【主权项】
1.雷达抗距离-速度联合欺骗干扰的自适应迭代滤波方法，其特征在于，包括以下步 骤： 1初始化步骤:设雷达对目标照射的驻留时间为Μ个脉冲，首先在快时间维对信号进行 匹配滤波，然后在慢时间维做傅里叶变换得到的距离-多普勒平面作为真实目标的距离-多 普勒平面初始值Xf和干扰信号的距离-多普勒平面初始值 2距离维估计迭代步骤： 2-1:构造真实目标的快时间变换矩阵GT,m与对角矩阵Π ^1其中u表示迭代次数，diag( ·)为求矩阵的对角线元素，ST,m=Sm为真实目标第m个脉冲 的回波信号，8"=[^(1)，^(2)，一，^(?]7是对信号模型的采样结果小为采样点数，（·， 表示转置。其中，C11表示真实目标第u-Ι次迭代距离像估计值，1 = 1，2，…，L，L表示距离单 元的个数，m=l，…，M，M表示脉冲个数； 构造干扰信号的快时间变换矩阵以》与对角矩阵!!；；，Μ :#?=8^表示干扰信号第m个脉冲的回波信号，i为干扰机捕获雷达信号需要的脉冲周 期数。其中，表示干扰信号第u-1次迭代距离像估计值； 2-2:得到真实目标在不同取值范围1里的距离像估计:α；5, |Τ|，1 = 1，…，N-1，1 = N，…， 得到干扰信号在不同取值范围丨里的距离像估计：c^m,1 = 1，···，Ν-1，1=Ν，···，?-(Ν-1)，1 = ?-(Ν-2)，···，Μ 步骤2-3:构造距离维上真实目标矩阵A#和干扰信号矩阵A；；1:3脉冲维估计迭代步骤： 3-1:估计真实目标第1个距离单元对应第k个多普勒单元的值;rp「/, A-1: k= 1，…，K估计干扰信号第1个距离单元对应第k个多普勒单元的值f/, Μ: k= 1，…，Κ其中，_ W和:_ Μ分别是丨摩f和！讀f的第1列，F= [fi，f2，…，fM]，fm= [eXp( j2JT mfi)，exp( j2Jimf2)，···，exp( j2JimfK)]T是KX 1的目标多普勒相位矢量，m=l，…，M，Fk为F的第 k行，I是单位矩阵，（·）H表示共辄转置； 3-3:令1 = 1，…，L得到真实目标的距离-多普勒平面估计：得到干扰信号的距离-多普勒平面估计：其中步骤5:重复步骤2和步骤3，直到达到停止迭代条件。
【文档编号】G01S7/36GK105954729SQ201610284089
【公开日】2016年9月21日
【申请日】2016年4月29日
【发明人】崔国龙, 张帅, 胡露, 卢金伟, 盛彪, 汪兵, 孔令讲, 杨晓波, 易伟
【申请人】电子科技大学