一种用于合成孔径雷达大场景欺骗干扰的逆距离徙动方法技术

一种用于合成孔径雷达大场景欺骗干扰的逆距离徙动方法，采取的技术方案如下：1、根据干扰方拟在雷达图像中产生的虚假电磁特性，构建干扰机欺骗模板：2、计算补余聚焦滤波器HHdiff(k,kx)在二维频域(k,kx)平面上的等间隔采样值。3、计算一致聚焦滤波器HHbulk(k,kx)在二维频域(k,kx)平面上的等间隔采样值，并与步骤二得到的补余聚焦滤波器HHdiff(k,kx)在二维频域(k,kx)平面上的等间隔采样值相乘。4、将步骤三得到的结果进行方位向逆傅里叶变换，并与干扰机消隐滤波器Hhelim(k,x)在距离频域(k,x)平面上的等间隔采样值相乘。本发明专利技术克服现有技术的不足，提供一种SAR大场景欺骗干扰方法，精确、高效地计算干扰机系统频率响应，使得雷达具有高分辨力、大斜视角和长合成孔径时，能够兼顾计算量、聚焦效果两个方面的要求。

【技术实现步骤摘要】
【专利说明】
本专利技术属于电子对抗和信号处理领域，具体涉及对合成孔径雷达的对抗技术，特 别是一种合成孔径雷达大场景欺骗干扰方法。 【
技术介绍
】 SAR(Synthetic Aperture Radar，合成孔径雷达)具有全天时、全天候、远距离、高 分辨成像能力，使信息安全的防护受到前所未有的挑战，对SAR的电子对抗技术因此受到广 泛的研究和注意。作为SAR电子对抗领域的热点难点问题之一，大场景欺骗干扰的目的是通 过在雷达图像中形成虚假的电磁散射特性而不引起雷达操作员的注意，从而迷惑雷达的目 标识别过程、破坏雷达的信息认知环节。 转发式干扰机实现SAR欺骗干扰的一般流程包括侦收雷达信号、干扰机调制、转发 干扰信号三个步骤。文献1:王盛利，于立，倪晋鳞，张光义合成孔径雷达的有源欺骗干扰 方法研究，"电子学报，2003(12): 1900-1902.中公开的技术将干扰机调制建模为一个线性 系统，并由线性系统的频率响应来表征。由于大场景欺骗干扰中，干扰机系统频率响应的计 算量与干扰实时性要求的矛盾突出，多项技术设计了优化的方法计算干扰机系统频率响 应，以降低计算量、提高干扰实时性。文献2:F.Zhou，B.Zhao，M.L.Tao，X.R.Bai，B.Chen，and G.C.Sun，"A large scene deceptive jamming method for space-borne SAR，''IEEE Trans.Geosci .Remote Sens.，vol .51，no.8，pp.4486-4495,Aug.2013.中公开的技术分解 干扰机频率响应，通过离线计算缓解在线计算的运算量;文献3: Y. C. Liu，W. Wang，X. Y. Pan， D.H.Dai,and D.J.Feng，"A frequency-domain three-stage algorithm for active deception jamming against synthetic aperture radar，''IET Radar，Sonar & Navig., vol.8，no.6，pp.639-646，Jul. 2014.中公开的技术通过二维频域的计算避免迭代的积分计 算。然而，为了换取计算量的优化，现有方法不同程度地牺牲干扰信号的聚焦效果。目前，当 雷达具有较高的分辨力、明显的斜视角、或较大的合成孔径时，现有方法产生的大场景虚假 目标将不能够良好的聚焦。 【
技术实现思路
】 本专利技术的目的在于克服现有技术的上述不足，提供一种SAR大场景欺骗干扰方法， 精确、高效地计算干扰机系统频率响应，使得雷达具有高分辨力、大斜视角和长合成孔径 时，仍能够兼顾计算量、聚焦效果两个方面的要求。 为了达到上述目的，本专利技术采取的技术方案如下： 步骤一:根据干扰方拟在雷达图像中产生的虚假电磁特性，构建干扰机欺骗模板： 以欺骗模板坐标系的原点0dt为中心，将二维空域(u，v)平面等间隔地采样。根据干 扰方拟在雷达图像中产生的虚假电磁特性，在每个采样位置上设定该位置所对应虚假散射 中心的雷达后向散射系数。将所有采样位置的众多虚假散射中心通过线性叠加组合在一 起，即可构建出任意虚假电磁特性，得到干扰机欺骗模板〇(u，V): 其中，u，v分别代表虚假散射中心在欺骗模板坐标系中的距离向位置和方位向位 置，m，n为整数，分别代表虚假散射中心在距离向和方位向的序号，δ( ·，·）为二维狄拉克 (Dirac)函数，Au，Δν分别代表欺骗模板坐标系在距离向和方位向的采样间隔，Σ为求和 符号。 步骤二:计算补余聚焦滤波器HHdiff(k，kx)在二维频域(k，k x)平面上的等间隔采样 值。补余聚焦滤波器HHdlff(k，kx)的数学表达式为：其中，（k，kx)代表二维频域坐标，k，kx分别代表雷达发射信号平面波频率和方位频 率，rect( ·)为归一化的矩形窗，ko为二倍的雷达载频除以电磁波传播速度，Br为二倍的雷 达带宽除以电磁波传播速度，《 a( ·)代表雷达天线方位向上的方向图，arctan( ·)为反正 切函数，9sq为雷达斜视角，Π代表双重积分运算，exp( ·)为指数函数，j为虚数单位。步骤二有两种实现方式，第一种实现方式的具体步骤为： 第一，对欺骗模板进行二维快速傅里叶变换，得到欺骗模板二维频谱HHsigma(kr， k x)在二维频域(kr，kx)平面上的等间隔采样值。其中 HHsigma(kr,kx) = JJ〇(u,v) · exp(-j23Tkru-j23Tkxv)dudv (3) 其中，kr代表距离频率。 第二，幅度加权，得到补余滤波器变换函数在二维频域(kr，kx)平 面上的等间隔采样值。其中， 为此，应在二维频域(kr，kx)平面上选取由式(5)所围成的扇形区域其中，0bw为雷达天线方位向波束宽度。截取HHsigma(kr，k x)在该扇形区域范围内的 数据，并根据式(4)中的rect( ·)项和coa( ·)项进行幅度加权。根据式(5)可知，扇形区域 的范围由雷达信号载频和带宽，雷达天线方位向波束宽度9bw和雷达斜视角0sq共四个雷达 参数确定。第三，通过对////,"", (/?.Λ )沿距离频率kr轴进行插值，得到补余聚焦滤波器 (k，kx)在二维频域(k，kx)平面上的等间隔采样值。插值运算中，应采用基带插值的核函数(如Sine函数），并根据斯托尔特(Stolt)变选取插值位置。具体而言，为获得补余聚焦滤波器HHdlff(k，k x)在二维频域 平面(k，kx)某一点处的采样值，应通过插值获得在二维频域平面(kr，kx)上位处的值。 步骤二的第二种实现方式的具体步骤为： 第一，对欺骗模板沿方位向v进行一维快速傅里叶变换。此时，欺骗模板将被变换 到方位频域，并且沿着方位频率kx轴的方向被等间隔地采样。 第二，沿距离向u进行Chirp-Z变换。设补余聚焦滤波器HHdlff(k，kx)在二维频域平 面(k，kx)上沿着距离频率k轴方向的等间隔采样点的初始位置Sks，采样间隔为Ak，那么 Chirp-Z变换运算中，距离频率kr的起始位置k rs和采样间隔Akr应根据式(7)所示的近似 Stolt变换选取，具体为 其中，D(kx)为距离徙动因子第三，幅度加权。在二维频域(kr，kx)平面上选取由式(5)所围成的圆环区域，截取 上一步所得结果在该圆环范围内的数据，并根据式(4)中的rect( ·)项和coa( ·)项进行幅 度加权。 步骤三:计算一致聚焦滤波器HHbulk(k，kx)在二维频域(k，kx)平面上的等间隔采样 值，并与步骤二得到的补余聚焦滤波器HHdlff(k，kx)在二维频域(k，kx)平面上的等间隔采样 值相乘。 其中，一致聚焦滤波器HHbulk(k，kx)由如下数学表达式给出步骤四：将步骤三得到的结果进行方位向逆傅里叶变换，并与干扰机消隐滤波器 Hhelim(k，x)在距离频域(k，x)平面上的等间隔采样值相乘。其中，干扰机消隐滤波器Hhelim (k，x)由如下数学表达式给出其中，x代表雷达的瞬时方位向位置，rj代表雷达与干扰机之间的最近距离，'代表雷达本文档来自技高网...

【技术保护点】
一种用于合成孔径雷达大场景欺骗干扰的逆距离徙动方法，其特征在于，包括如下步骤：步骤一：根据干扰方拟在雷达图像中产生的虚假电磁特性，构建干扰机欺骗模板：以欺骗模板坐标系的原点Odt为中心，将二维空域(u,v)平面等间隔地采样；根据干扰方拟在雷达图像中产生的虚假电磁特性，在每个采样位置上设定该位置所对应虚假散射中心的雷达后向散射系数；将所有采样位置的众多虚假散射中心通过线性叠加组合在一起，构建出任意虚假电磁特性，得到干扰机欺骗模板σ(u,v)：σ(u,v)=ΣmΣnσ[m,n]·δ(u-mΔu,v-nΔv)---(1)]]>其中，u,v分别代表虚假散射中心在欺骗模板坐标系中的距离向位置和方位向位置，m,n为整数，分别代表虚假散射中心在距离向和方位向的序号，δ(·,·)为二维狄拉克Dirac函数，Δu,Δv分别代表欺骗模板坐标系在距离向和方位向的采样间隔，Σ为求和符号；步骤二：计算补余聚焦滤波器HHdiff(k,kx)在二维频域(k,kx)平面上的等间隔采样值；补余聚焦滤波器HHdiff(k,kx)的数学表达式为：HHdiff(k,kx)=rect(k-k0Br)·ωa(arctan(kxk2-kx2)-θsq)·∫∫σ(u,v)·exp(-j2πk2-kx2u-j2πkxv)dudv---(2)]]>其中，(k,kx)代表二维频域坐标，k,kx分别代表雷达发射信号平面波频率和方位频率，rect(·)为归一化的矩形窗，k0为二倍的雷达载频除以电磁波传播速度，Br为二倍的雷达带宽除以电磁波传播速度，ωa(·)代表雷达天线方位向上的方向图，arctan(·)为反正切函数，θsq为雷达斜视角，∫∫代表双重积分运算，exp(·)为指数函数，j为虚数单位；步骤三：计算一致聚焦滤波器HHbulk(k,kx)在二维频域(k,kx)平面上的等间隔采样值，并与步骤二得到的补余聚焦滤波器HHdiff(k,kx)在二维频域(k,kx)平面上的等间隔采样值相乘；其中，一致聚焦滤波器HHbulk(k,kx)由如下数学表达式给出HHbulk(k,kx)=rJ+pck0·exp(-j2πk2-kx2(rJ+pc)-j2πkxqc)---(3)]]>步骤四：将步骤三得到的结果进行方位向逆傅里叶变换，并与干扰机消隐滤波器Hhelim(k,x)在距离频域(k,x)平面上的等间隔采样值相乘；其中，干扰机消隐滤波器Hhelim(k,x)由如下数学表达式给出Hhelim(k,x)=1ωa(-arctan(x/rJ)-θsq)·exp(j2πkrJ2+x2)---(4)]]>其中，x代表雷达的瞬时方位向位置，rJ代表雷达与干扰机之间的最近距离，代表雷达与干扰机之间的瞬时距离。...

【技术特征摘要】

【专利技术属性】
技术研发人员：王伟，刘永才，冯德军，潘小义，傅其祥，顾赵宇，饶彬，李程，陆珊珊，张文明，
申请(专利权)人：中国人民解放军国防科学技术大学，
类型：发明
国别省市：湖南;43