雷达抗距离-速度联合欺骗干扰的自适应迭代滤波方法技术

本发明专利技术提供一种雷达抗距离‑速度联合欺骗干扰的自适应迭代滤波方法，分别对目标和干扰的距离维和脉冲维进行估计，然后通过迭代来提高算法的精度。在干扰相对于目标延迟已知个脉冲的前提下，将自适应滤波的思想运用到对干扰和目标的距离‑多普勒平面的估计上，为了提高估计精度，同时引入了迭代步骤。本发明专利技术针对一种少有干扰形式：对抗距离‑速度联合，不需要对波形进行设计，实施较为方便，并且在多目标多干扰的场景下仍能对干扰有效压制并能对目标正确检测。仿真结果表明，本发明专利技术采用的方法相对于传统的匹配滤波和多普勒处理有明显的性能优越性。

【技术实现步骤摘要】

本专利技术属于雷达抗干扰技术，特别涉及自适应迭代滤波技术。
技术介绍
在电子战中，随着数字射频存储器DRFM技术在干扰机中的应用，使得电子反对抗ECCM变得更加困难。装备了数字射频存储器的干扰机在截获了雷达发射信号后进行，复制，存储，调制并转发。由于干扰信号复制了发射信号的模型，所以可以在雷达接收机中获得处理增益，这样在雷达的显示屏上将出现多个目标，在不具备电子反对抗能力的雷达系统中，这使得雷达无法辨别出真实目标，对传统雷达系统的生存构成了极大挑战。干扰机对截获的雷达发射信号的调制可以改变时延，或者对多普勒频率变换，以分别产生距离假目标和速度假目标，或者同时对两者进行调制，这样产生的即是距离-速度联合欺骗干扰。显然，后者对雷达的威胁更大。为了确保雷达对真实目标的检测，以及后续的跟踪，提高雷达的抗距离-速度联合欺骗干扰能力具有重要的理论价值和实际意义。一种被广泛研究的抗干扰方法是对雷达发射信号的波形进行设计。例如设计雷达发射的脉冲是互相正交的，假设干扰机需要几个脉冲重复周期来识别和捕获发射信号，那么反射回来后的干扰信号经过当前脉冲重复周期的匹配滤波器后会衰减，因为设计波形为脉冲间相互正交。但是，在动目标显示(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]提出了一种通过设计自适应初始相位的编码信号来抗速度欺骗干扰的方法，但是此方法仅考虑了单个目标，在多目标场景下，此方法将丧失抗干扰能力。从公开发表的文献来看，针对距离-速度联合欺骗干扰抑制的方法还未有研究。
技术实现思路
本专利技术所要解决的技术问题是，提供一种通过自适应迭代滤波的方法来对抗距离-速度联合欺骗干扰的方法。本专利技术为解决上述技术问题所采用的技术方案是，一种雷达抗距离-速度联合欺骗干扰的自适应迭代滤波方法，包括以下步骤：1初始化步骤：设雷达对目标照射的驻留时间为M个脉冲，首先在快时间维对信号进行匹配滤波，然后在慢时间维做傅里叶变换得到的距离-多普勒平面作为真实目标的距离-多普 勒平面初始值和干扰信号的距离-多普勒平面初始值2距离维估计迭代步骤：2-1：构造真实目标的快时间变换矩阵GT,m与对角矩阵其中u表示迭代次数，diag(·)为求矩阵的对角线元素，sT,m＝sm为真实目标第m个脉冲的回波信号，sm＝[sm(1),sm(2),…,sm(N)]T是对信号模型的采样结果，N为采样点数，(.)T表示转置。，其中，表示真实目标第u-1次迭代距离像估计值，l＝1,2,…,L，L表示距离单元的个数，m＝1,…,M，M表示脉冲个数；构造干扰信号的快时间变换矩阵GJ,m与对角矩阵sJ,m表示干扰信号第m个脉冲的回波信号，其中，表示干扰信号第u-1次迭代距离像估计值；2-2：得到真实目标在不同取值范围l里的距离像估计：得到干扰信号在不同取值范围l里的距离像估计： 步骤2-3：构造距离维上真实目标矩阵和干扰信号矩阵3脉冲维估计迭代步骤：3-1：估计真实目标第l个距离单元对应第k个多普勒单元的值估计干扰信号第l个距离单元对应第k个多普勒单元的值其中，和分别是和的第l列，F＝[f1,f2,…,fM]，fm＝[exp(j2πmf1),exp(j2πmf2),…,exp(j2πmfK)]T是K×1的归一化多普勒相位矢量，Fk为F的第k行，I是单位矩阵，(.)H表示共轭转置；3-3：令l＝1,…,L得到真实目标的距离-多普勒平面估计：得到干扰信号的距离-多普勒平面估计：其中，步骤5：重复步骤2和步骤3，直到达到停止迭代条件。本专利技术分别对目标和干扰的距离维和脉冲维进行估计，然后通过迭代来提高算法的精度。 在干扰相对于目标延迟已知个脉冲的前提下，将自适应滤波的思想运用到对干扰和目标的距离-多普勒平面的估计上，为了提高估计精度，同时引入了迭代步骤。本专利技术的有益效果是，提针对一种少有干扰形式：对抗距离-速度联合，不需要对波形进行设计，实施较为方便，并且在多目标多干扰的场景下仍能对干扰有效压制并能对目标正确检测。仿真结果表明，本专利技术采用的方法相对于传统的匹配滤波和多普勒处理有明显的性能优越性。附图说明图1为算法流程图。图2为传统匹配滤波抗距离-速度联合欺骗干扰估计结果。图3自适应迭代滤波法抗距离-速度联合欺骗干扰结果。3-1为自适应迭代滤波法对真实目标的距离-多普勒平面估计结果示意图。3-2为自适应迭代滤波法对干扰目标的距离-多普勒平面估计结果示意图。具体实施方式本专利技术的具体实施步骤如图1所示：步骤1：脉冲多普勒雷达不断的发射脉冲信号，同时接收目标的回波信号，接收到的回波信号不光包括目标回波，还包括干扰机产生的欺骗式干扰信号，并且干扰类型为距离-速度联合欺骗式干扰，那么此时接着转到步骤2做下一步处理。步骤2：准备自适应迭代算法初始值。假设雷达对目标照射的驻留时间为M个脉冲，首先在快时间维对信号进行匹配滤波，然后在慢时间维做傅里叶变换，这就是传统的脉冲多普勒处理，此时可以得到的距离-多普勒品面作为下面迭代步骤的初始值和其中，和分别表示真实目标和干扰的距离-多普勒平面初始值。步骤3：距离维估计步骤：步骤3-1：对脉冲维的估计做准备工作，构造真实目标的快时间变换矩阵GT,m；根据上一步的估计结果，构造对角矩阵其中：对l＝1,…,N-1，l＝N,…,L-(N-1)，l＝L-(N-2),…,L，在l位于不同的取值范围时均做这样的处理。对于干扰信号的估计，重复上面一样的准备工作。步骤3-2：利用上一步已经得到的准备信息，分别得到在第u次迭代时，l在不同取值范围里的距离像估计，即：同时也可以得到干扰目标在不同取值范围里的距离像估计，即：对于真实目标的估计，具体如下：其中，gT,m,l表示GT,m的第l列，ym[l]＝[ym[l],…,ym[l+(N-1)]]表示回波信号第m个脉冲对应的从第l个距离单元到第l+(N-1)个距离单元的采样结果，I表示单位矩阵。同样对于干扰目标，具体计算方式如下：其中，gJ,m,l表示GJ,m的第l列，sJ,m＝sT,m-i表示干扰滞后信号i个脉冲重复周期。步骤3-3：利用步骤3-2的估计结果构造矩阵和其中：步骤4：脉冲维估计步骤：步骤4-1：由步骤3的估计结果结合上一次的距离多普勒平面估计结果和 来估计第l个距离单元对应第k个多普勒单元的值，即同理也可以得到干扰目标的估计结果即：其中：和分别是和的第l列，F＝[f1,f2,…,fM]，fm＝[exp(j2πmf1),exp(j2πmf2),…,exp(j2πmfK)]T是K×1的目标多普勒相位矢量,Fk为F的第k行，I为单位矩阵，(.)H表示转置共轭；步骤4-2：令k＝1,…,K，可以得到：对干扰亦然，即：步骤4-3：令l＝1,…,L，可以得到真实目本文档来自技高网...

【技术保护点】
雷达抗距离‑速度联合欺骗干扰的自适应迭代滤波方法，其特征在于，包括以下步骤：1初始化步骤：设雷达对目标照射的驻留时间为M个脉冲，首先在快时间维对信号进行匹配滤波，然后在慢时间维做傅里叶变换得到的距离‑多普勒平面作为真实目标的距离‑多普勒平面初始值和干扰信号的距离‑多普勒平面初始值2距离维估计迭代步骤：2‑1：构造真实目标的快时间变换矩阵GT,m与对角矩阵ΠT,m(u-1)[l]=diag(|αT,m(u-1)[l-(N-1)]|2,...|αT,m(u-1)[l+(N-1)]|2)]]>其中u表示迭代次数，diag(·)为求矩阵的对角线元素，sT,m＝sm为真实目标第m个脉冲的回波信号，sm＝[sm(1),sm(2),…,sm(N)]T是对信号模型的采样结果，N为采样点数，(·)T表示转置。其中，表示真实目标第u‑1次迭代距离像估计值，l＝1,2,…,L，L表示距离单元的个数，m＝1,…,M，M表示脉冲个数；构造干扰信号的快时间变换矩阵GJ,m与对角矩阵ΠT,m(u-1)[l]=diag(|αJ,m(u-1)[l-(N-1)]|2,...|αJ,m(u-1)[l+(N-1)]|2)]]>sJ,m＝sm‑i表示干扰信号第m个脉冲的回波信号，i为干扰机捕获雷达信号需要的脉冲周期数。其中，表示干扰信号第u‑1次迭代距离像估计值；2‑2：得到真实目标在不同取值范围l里的距离像估计：l＝1,…,N‑1，l＝N,…,L‑(N‑1)，l＝L‑(N‑2),…,L；得到干扰信号在不同取值范围l里的距离像估计：l＝1,…,N‑1，l＝N,…,L‑(N‑1)，l＝L‑(N‑2),…,L；步骤2‑3：构造距离维上真实目标矩阵和干扰信号矩阵AT(u)=[αT,1(u),...,αT,M(u)]=αT,1(u)[1]...αT,M(u)[1]······αT,1(u)[L]...αT,M(u)[L]]]>AJ(u)=[αJ,1(u),...,αJ,M(u)]=αJ,1(u)[1]...αJ,M(u)[1]······αJ,1(u)[L]...αJ,M(u)[L]]]>3脉冲维估计迭代步骤：3‑1：估计真实目标第l个距离单元对应第k个多普勒单元的值XT(u)[l,k]=FkH{FDT(u-1)[l]FH+I/smHsm}-1βT(u)[l]FkH{FDT(u-1)[l]FH+I/smHsm}-1Fk]]>k＝1,…,KDT(u-1)[l]=diag(|XT(u-1)[l,1]|2,...,|XT(u-1)[l,K]|2)]]>估计干扰信号第l个距离单元对应第k个多普勒单元的值XJ(u)[l,k]=FkH{FDJ(u-1)[l]FH+I/smHsm}-1βJ(u)[l]FkH{FDJ(u-1)[l]FH+I/smHsm}-1Fk]]>k＝1,…,KDJ(u-1)[l]=diag(|XJ(u-1)[l,1]|2,...,|XJ(u-1)[l,K]|2)]]>其中，和分别是和的第l列，F＝[f1,f2,…,fM]，fm＝[exp(j2πmf1),exp(j2πmf2),…,exp(j2πmfK)]T是K×1的目标多普勒相位矢量，m＝1,…,M,Fk为F的第k行，I是单位矩阵，(·)H表示共轭转置；3‑3：令l＝1,…,L得到真实目标的距离‑多普勒平面估计：XT(u)={xT(u)[1]}T···{xT(u)[L]}T]]...

【技术特征摘要】
1.雷达抗距离-速度联合欺骗干扰的自适应迭代滤波方法，其特征在于，包括以下步骤：1初始化步骤：设雷达对目标照射的驻留时间为M个脉冲，首先在快时间维对信号进行匹配滤波，然后在慢时间维做傅里叶变换得到的距离-多普勒平面作为真实目标的距离-多普勒平面初始值和干扰信号的距离-多普勒平面初始值2距离维估计迭代步骤：2-1：构造真实目标的快时间变换矩阵GT,m与对角矩阵 Π T , m ( u - 1 ) [ l ] = d i a g ( | α T , m ( u - 1 ) [ l - ( N - 1 ) ] | 2 , ... | α T , m ( u - 1 ) [ l + ( N - 1 ) ] | 2 ) ]]>其中u表示迭代次数，diag(·)为求矩阵的对角线元素，sT,m＝sm为真实目标第m个脉冲的回波信号，sm＝[sm(1),sm(2),…,sm(N)]T是对信号模型的采样结果，N为采样点数，(·)T表示转置。其中，表示真实目标第u-1次迭代距离像估计值，l＝1,2,…,L，L表示距离单元的个数，m＝1,…,M，M表示脉冲个数；构造干扰信号的快时间变换矩阵GJ,m与对角矩阵 Π T , m ( u - 1 ) [ l ] = d i a g ( | α J , m ( u - 1 ) [ l - ( N - 1 ) ] | 2 , ... | α J , m ( u - 1 ) [ l + ( N - 1 ) ] | 2 ) ]]>sJ,m＝sm-i表示干扰信号第m个脉冲的回波信号，i为干扰机捕获雷达信号需要的脉冲周期数。其中，表示干扰信号第u-1次迭代距离像估计值；2-2：得到真实目标在不同取值范围l里的距离像估计：l＝1,…,N-1，l＝N,…,L-(N-1)，l＝L-(N-2),…,L；得到干扰信号在不同取值范围l里的距离像估计：l＝1,…,N-1，l＝N,…,L-(N-1)，l＝L-(N-2),…,L；步骤2-3：构造距离维上真实目标矩阵和干扰信号矩阵 A T ( u ) = [ α T , 1 ( u ) , ... , α T , M ( u ) ] = α T , 1 ( u ) [ 1 ] ... α T , M ( u ) [ 1 ] · · · · · · α T ...

【专利技术属性】
技术研发人员：崔国龙，张帅，胡露，卢金伟，盛彪，汪兵，孔令讲，杨晓波，易伟，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51