CN113203991A - 一种多基sar在多干扰机环境下的抗欺骗型干扰方法 - Google Patents

一种多基sar在多干扰机环境下的抗欺骗型干扰方法 Download PDF

Info

Publication number
CN113203991A
CN113203991A CN202110472443.0A CN202110472443A CN113203991A CN 113203991 A CN113203991 A CN 113203991A CN 202110472443 A CN202110472443 A CN 202110472443A CN 113203991 A CN113203991 A CN 113203991A
Authority
CN
China
Prior art keywords
jammer
interference
target
false
representing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110472443.0A
Other languages
English (en)
Other versions
CN113203991B (zh
Inventor
武俊杰
王雯璟
张潇文
裴季方
孙稚超
杨建宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202110472443.0A priority Critical patent/CN113203991B/zh
Publication of CN113203991A publication Critical patent/CN113203991A/zh
Application granted granted Critical
Publication of CN113203991B publication Critical patent/CN113203991B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/292—Extracting wanted echo-signals
    • G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9058—Bistatic or multistatic SAR
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/35—Details of non-pulse systems
    • G01S7/352—Receivers
    • G01S7/354—Extracting wanted echo-signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，包括以下步骤：S1、建立回波模型；S2、虚假目标的检测与判别：利用基于谱残差的视觉显著性检测方法与形态学方法相结合实现多基SAR图像中的目标检测，然后利用基于欧式距离准则的判别法判定目标是否虚假；S3、对虚假目标位置分布规律进行分析；S4、进行干扰机定位，并利用波束形成的方法对干扰信号进行抑制。本发明解决在多基SAR系统中存在多个转发不同时延干扰信号的干扰机时，对多个干扰机进行准确定位的问题，有利于后续通过干扰机位置信息对干扰信号进行有效抑制。

Description

一种多基SAR在多干扰机环境下的抗欺骗型干扰方法
技术领域
本发明属于雷达信号处理技术领域，具体涉及一种多基SAR在多干扰机环境下的抗欺骗型干扰方法。
背景技术
合成孔径雷达(SAR)是一种二维高分辨成像雷达，具有全天时、全天候的特点，可以大大提高雷达对目标的信息获取能力。现如今，SAR已经在灾害监测、资源勘察、地质测绘、军事侦察等领域得到了非常广泛的应用。为了获取目标更多的信息，多基SAR获得了越来越多的关注。
多基SAR由于具有灵活的构型，可以同时从多个角度观测感兴趣区域，获得目标不同角度的观测信息。在电子对抗领域，针对多基SAR系统干扰的研究已经成为近几年的研究热点。
SAR干扰技术通过多年的发展，按照干扰实施方式的不同，SAR干扰方式通常分为压制干扰和欺骗干扰。压制干扰是利用大功率的非相干或部分相干信号对SAR进行干扰，干扰SAR的信息获取。欺骗干扰主要是通过获取发射站的发射信号，然后分析发射信号的相关参数，最后通过参数调制，再发射干扰信号使SAR图像中出现虚假目标，该虚假目标很难被直接消除。由于，与压制干扰相比较，欺骗干扰对功率的需求更小，欺骗干扰技术则更容易实现。
随着转发式欺骗干扰技术的发展，对应的抗干扰研究也引起大家的关注。在不同SAR 体制情况下，对相关抗干扰技术进行了研究。分别有传统SAR、MIMO SAR、单通道SAR 以及多通道SAR等多个方面。然而，目前对多基SAR系统的抗干扰方法研究较少，目前已有的研究方法无法直接对多基SAR系统中的干扰机进行准确定位及干扰信号进行抑制。
发明内容
本发明的目的在于克服现有技术的不足，提供一种解决了在多基SAR系统中存在多个转发不同时延干扰信号的干扰机时，对多个干扰机进行准确定位的问题多基SAR在多干扰机环境下的抗欺骗型干扰方法。
本发明的目的是通过以下技术方案来实现的：一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，包括以下步骤：
S1、建立回波模型；
S2、虚假目标的检测与判别：利用基于谱残差的视觉显著性检测方法与形态学方法相结合实现多基SAR图像中的目标检测，然后利用基于欧式距离准则的判别法判定目标是否虚假；
S3、对虚假目标位置分布规律进行分析；
S4、进行干扰机定位，并利用波束形成的方法对干扰信号进行抑制。
进一步地，所述步骤S1具体实现方法为：多基SAR系统包括一个发射机、多个接收机和干扰机，发射机发射线性调频信号；
方位时刻为η时，第i个接收机接收到的真实环境中点目标回波表示如下：
Figure BDA0003045871370000021
其中，A0表示点目标雷达横截面面积，wr和wa分别表示距离向和方位向的包络函数，τ表示距离向时间，η表示方位向时间，ηc表示中心时刻；f0表示载波频率，Kr表示调频速率， c表示光速，R(i)(η)表示在η时刻发射机到点目标的距离RT(η)及第i个接收机到点目标的距离
Figure BDA0003045871370000022
的距离和；
第i个接收机的距离历史和R(i)(η)的表达式如下：
Figure BDA0003045871370000023
其中，
Figure BDA0003045871370000024
表示第i个接收机的位置，(xT(η),yT(η),zT(η))表示发射机的位置，(xq,yq,zq)表示成像环境中点目标的位置；
干扰机通过参数调制，进行时延处理后，再把干扰信号发射出去；由第m个干扰机进行第n个时延转发的干扰信号在η时刻在第i个接收机形成的干扰瞬时回波
Figure BDA0003045871370000025
表示如下：
Figure BDA0003045871370000031
其中，
Figure BDA0003045871370000032
表示第m个干扰机的第n个时延，
Figure BDA0003045871370000033
表示在η时刻第m个干扰机到发射机及第i个接收机的距离历史和，表示如下：
Figure BDA0003045871370000034
其中，
Figure BDA0003045871370000035
表示第m个干扰机的位置；
在多干扰机环境下，多基SAR的接收机会同时受到多个干扰机发射的干扰信号影响；所以，第i个接收机在η时刻收到的回波s(i)(τ,η)如下所示：
Figure BDA0003045871370000036
其中，M表示该干扰环境共有M个干扰机，N表示第M个干扰机发射N个不同时延的干扰信号；
利用后向投影算法将接收到的回波信号生成多基SAR图像。
进一步地，所述步骤S2中，基于谱残差的视觉显著性检测方法的实现步骤如下：
S21、给定SAR图像I(x,y)，其频谱表示为
Figure BDA0003045871370000037
其中
Figure BDA0003045871370000038
表示二维傅里叶变换；
对应的幅度谱和相位谱分别是：
Figure BDA0003045871370000039
Figure BDA00030458713700000310
其中，
Figure BDA00030458713700000311
表示输入频谱的幅度计算，
Figure BDA00030458713700000312
表示输入频谱的相位计算；
对SAR图像的幅度谱进行自然对数运算，得到对应的对数谱，如下表示：
Figure BDA0003045871370000041
S22、利用均值滤波器进行平滑滤波，获得SAR图像中的均值谱，表示为：
Figure BDA0003045871370000042
其中，h(fx,fy)是一个n×n的矩阵，h(fx,fy)为频谱均值滤波器，具体表示为：
Figure BDA0003045871370000043
S23、用对数谱减去均值谱，得到SAR图像的谱残差：
R(fx,fy)＝L(fx,fy)-B(fx,fy) (12)；
S24、经过二维逆傅里叶变换和高斯平滑滤波，将其还原到图像域，得到SAR图像中的新颖部分，即SAR图像的显著图：
Figure BDA0003045871370000044
S25、对SAR图像的显著图进行开运算，然后计算SAR图像的显著图中目标点的面积，与设定的面积阈值进行比较，若大于或等于面积阈值则将该目标作为感兴趣区域，否则无视该目标；
S26、利用基于欧式距离判别法实现对虚假目标的判定；首先，从这些图像中任选两幅不同的图像Mi和Mj；设Ci,k为第i幅图像中的第k个目标的重心，Cj,h表示第j幅图像中的第h个目标的重，γ为欧式距离门限；依次计算Ci,k和Cj,h的欧氏距离，若Ci,k可以在第 j幅图像中找到满足欧式距离小于γ的，说明Ci,k为真目标，反之为虚假目标。
进一步地，所述步骤S3具体实现方法为：假设干扰环境由M个发射多个不同时延的干扰机组成，第m个干扰机发射N个不同时延干扰信号，该干扰机在每幅SAR图像中会产生N个不同的虚假目标；根据距离欺骗干扰原理得到干扰机在η时刻的多普勒频率为：
Figure BDA0003045871370000045
其中，
Figure BDA0003045871370000051
和
Figure BDA0003045871370000052
分别表示发射机和第i个接收机的速度矢量，λ表示干扰机的发射波长，
Figure BDA0003045871370000053
和
Figure BDA0003045871370000054
分别表示发射机和第i个接收机的位置矢量，
Figure BDA0003045871370000055
表示第m个干扰机的位置矢量；
由该干扰机生成的虚假目标与该干扰机的瞬时多普勒相同，如下表示：
Figure BDA0003045871370000056
Figure BDA0003045871370000057
其中，
Figure BDA0003045871370000058
表示在η时刻，第i个接收机受到第m个干扰机第n个时延干扰信号产生的虚假目标的瞬时多普勒频率；
Figure BDA0003045871370000059
表示在第i个接收机所成的SAR图像中，第m个干扰机的第n个时延干扰信号产生的虚假目标位置矢量，n＝1,2,…,N；由此，N个不同时延干扰信号同时满足下列方程组：
Figure BDA00030458713700000510
推出：
Figure BDA00030458713700000511
同时，第i个接收机到由第m个干扰机的第n个时延干扰信号产生的虚假目标的距离历史和表示为：
Figure BDA00030458713700000512
其中，
Figure BDA00030458713700000513
表示第i幅SAR图像中该干扰机的第n个时延干扰信号生成的虚假目标位置；
得到
Figure BDA00030458713700000514
表示为：
Figure BDA00030458713700000515
其中，
Figure BDA0003045871370000061
表示第m个干扰机到发射机和第i个接收机的距离历史和；
通过上述分析发现第i个接收机成的SAR图像中，第m个干扰机第n个时延
Figure BDA0003045871370000062
转发的干扰信号产生的虚假目标位于
Figure BDA0003045871370000063
对应的等距离线与干扰机所在的等多普勒线的交点处；其中，
Figure BDA0003045871370000064
表示第m个干扰机到发射机和第i个接收机的距离历史和；
随着干扰时延的变化，虚假目标沿着干扰机所在的等多普勒线分布；所以，在多基SAR 图像中，将分布在同一条等多普勒线上的虚假目标视为同一干扰机生成。
进一步地，所述步骤S4中进行干扰机定位的具体实现方法为：
S41、虚假目标与干扰机对应关系组合：假设第m个干扰机发射N个不同时延的干扰信号，各接收机所成SAR图像中产生多个虚假目标；
组合规则为：在I个接收机接收到的I幅SAR图像中分别任意选择一个虚假目标，把这些虚假目标视为一个小组；已被选择过的虚假目标不能再被选择，继续选择虚假目标形成小组，直到所有虚假目标都被选完，每种组合包含N个小组，每个小组包含I个虚假目标，N个小组形成一种组合情况；
S42、定位方程组求解：虚假目标均满足公式(20)，将每个小组中每个SAR图像中的虚假目标参数带入公式(20)，得到有N个方程组成的方程组，通过求解方程组，得到N 个干扰机位置及N个时延大小；
S43、干扰机的位置判定：每种组合情况解得的N个干扰机位置分别表示为
Figure BDA0003045871370000066
其中
Figure BDA0003045871370000066
表示该组合情况中第n个方程组所解得的干扰机位置；当
Figure BDA0003045871370000067
时，则该组合情况为正确组合，即说明同一小组的虚假目标是由第m个干扰机的同一个时延转发干扰信号产生的，解得的干扰机位置为干扰机的正确位置；如果不满足，则继续进行下一个组合情况的判断，直到找到正确的组合为止。
本发明的有益效果是：本发明解决在多基SAR系统中存在多个转发不同时延干扰信号的干扰机时，对多个干扰机进行准确定位的问题，有利于后续通过干扰机位置信息对干扰信号进行有效抑制。
附图说明
图1为本发明的多基SAR系统结构示意图；
图2为本发明的多基SAR在多干扰机环境下的抗欺骗型干扰方法的流程图；
图3为本发明的目标场景示意图；
图4为本实施例含有虚假目标的成像结果，其中(a)、(b)、(c)和(d)分别对应接收站1、接收站2、接收站3和接收站4；
图5为消除虚假目标后的成像结果，其中(a)、(b)、(c)和(d)分别对应接收站1、接收站2、接收站3和接收站4。
具体实施方式
本发明的解决方案是，首先，多个接收站收到真实目标与虚假目标的回波，再利用后向投影算法(BP)生成多基SAR图像，然后，采用基于谱残差的视觉显著性检测方法对 SAR图像中感兴趣的目标进行有效检测后，使用欧式距离准则的判别法对虚假目标进行有效判定。最后，通过分析虚假目标的位置分布规律，将虚假目标与干扰机的对应关系进行排列组合，结合多基SAR构型信息实现多个干扰机的准确定位。在本发明中，多基SAR 系统至少需要4个接收机。
本发明主要采用仿真实验的方法进行验证，所有步骤、结论都在Matlab2016上验证正确。该具体实施方式是针对于一发四收的多基SAR系统，即一个发射站，四个接收站。干扰环境由两个干扰机组成，两个个干扰机发射多个不同时延的回波。具体结构如图1所示。
下面就具体实施方式和附图对本发明作进一步的详细描述。
如图2所示，本发明的一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，包括以下步骤：
S1、建立回波模型；具体实现方法为：多基SAR系统的结构如图1所示，其包括一个发射机、4个接收机和2个干扰机，发射机发射线性调频信号(LFM)；
多基SAR系统包括一个发射机、多个接收机和干扰机，发射机发射线性调频信号；
方位时刻为η时，第i个接收机接收到的真实环境中点目标回波表示如下：
Figure BDA0003045871370000071
其中，A0表示点目标雷达横截面面积，wr和wa分别表示距离向和方位向的包络函数，τ表示距离向时间，η表示方位向时间，ηc表示中心时刻；f0表示载波频率，Kr表示调频速率，c表示光速，R(i)(η)表示在η时刻发射机到点目标的距离RT(η)及第i个接收机到点目标的距离
Figure BDA0003045871370000081
的距离和；
第i个接收机的距离历史和R(i)(η)的表达式如下：
Figure BDA0003045871370000082
其中，
Figure BDA0003045871370000083
表示第i个接收机的位置，(xT(η),yT(η),zT(η))表示发射机的位置，(xq,yq,zq)表示成像环境中点目标的位置；
干扰机通过参数调制，进行时延处理后，再把干扰信号发射出去；由第m个干扰机进行第n个时延转发的干扰信号在η时刻在第i个接收机形成的干扰瞬时回波
Figure BDA0003045871370000084
表示如下：
Figure BDA0003045871370000085
其中，
Figure BDA0003045871370000086
表示第m个干扰机的第n个时延，
Figure BDA0003045871370000087
表示在η时刻第m个干扰机到发射机及第i个接收机的距离历史和，表示如下：
Figure BDA0003045871370000088
其中，
Figure BDA0003045871370000089
表示第m个干扰机的位置；
在多干扰机环境下，多基SAR的接收机会同时受到多个干扰机发射的干扰信号影响；所以，第i个接收机在η时刻收到的回波s(i)(τ,η)如下所示：
Figure BDA00030458713700000810
其中，M表示该干扰环境共有M个干扰机，N表示第M个干扰机发射N个不同时延的干扰信号；
利用后向投影算法将接收到的回波信号生成多基SAR图像，接收站同时受到多个干扰机的干扰，在所成的SAR图像中会出现多个不同的虚假目标，从而严重影响SAR图像的解释。
本实施例对经典海面场景分别进行成像，得到四幅双基SAR图像。系统仿真所需参数如表1所示，目标场景如图3所示。
表1
Figure BDA0003045871370000091
S2、虚假目标的检测与判别：利用基于谱残差的视觉显著性检测方法与形态学方法相结合实现多基SAR图像中的目标检测，然后利用基于欧式距离准则的判别法判定目标是否虚假；
基于谱残差的视觉显著性检测方法的实现步骤如下：
S21、给定SAR图像I(x,y)，其频谱表示为
Figure BDA0003045871370000092
其中
Figure BDA0003045871370000101
表示二维傅里叶变换；
对应的幅度谱和相位谱分别是：
Figure BDA0003045871370000102
Figure BDA0003045871370000103
其中，
Figure BDA0003045871370000104
表示输入频谱的幅度计算，
Figure BDA0003045871370000105
表示输入频谱的相位计算；
对SAR图像的幅度谱进行自然对数运算，得到对应的对数谱，如下表示：
L(fx,fy)＝In[A(fx,fy)] (9)；
S22、利用均值滤波器进行平滑滤波，获得SAR图像中的均值谱，表示为：
Figure BDA0003045871370000106
其中，h(fx,fy)是一个n×n的矩阵，h(fx,fy)为频谱均值滤波器，具体表示为：
Figure BDA0003045871370000107
S23、用对数谱减去均值谱，得到SAR图像的谱残差：
R(fx,fy)＝L(fx,fy)-B(fx,fy) (12)；
S24、经过二维逆傅里叶变换和高斯平滑滤波，将其还原到图像域，得到SAR图像中的新颖部分，即SAR图像的显著图：
Figure BDA0003045871370000108
S25、对SAR图像的显著图进行开运算，然后计算SAR图像的显著图中目标点的面积，与设定的面积阈值进行比较，若大于或等于面积阈值则将该目标作为感兴趣区域，否则无视该目标；
S26、利用基于欧式距离判别法实现对虚假目标的判定；首先，从这些图像中任选两幅不同的图像Mi和Mj；设Ci,k为第i幅图像中的第k个目标的重心，Cj,h表示第j幅图像中的第h个目标的重，γ为欧式距离门限；依次计算Ci,k和Cj,h的欧氏距离，若Ci,k可以在第j幅图像中找到满足欧式距离小于γ的，说明Ci,k为真目标，反之为虚假目标。
S3、对虚假目标位置分布规律进行分析；
具体实现方法为：假设干扰环境由M个发射多个不同时延的干扰机组成，第m个干扰机发射N个不同时延干扰信号，该干扰机在每幅SAR图像中会产生N个不同的虚假目标；根据距离欺骗干扰原理得到干扰机在η时刻的多普勒频率为：
Figure BDA0003045871370000111
其中，
Figure BDA0003045871370000112
和
Figure BDA0003045871370000113
分别表示发射机和第i个接收机的速度矢量，λ表示干扰机的发射波长，
Figure BDA0003045871370000114
和
Figure BDA0003045871370000115
分别表示发射机和第i个接收机的位置矢量，
Figure BDA0003045871370000116
表示第m个干扰机的位置矢量；
由该干扰机生成的虚假目标与该干扰机的瞬时多普勒相同，如下表示：
Figure BDA0003045871370000117
Figure BDA0003045871370000118
其中，
Figure BDA0003045871370000119
表示在η时刻，第i个接收机受到第m个干扰机第n个时延干扰信号产生的虚假目标的瞬时多普勒频率；
Figure BDA00030458713700001110
表示在第i个接收机所成的SAR图像中，第m个干扰机的第n个时延干扰信号产生的虚假目标位置矢量，n＝1,2,…,N；由此，N个不同时延干扰信号同时满足下列方程组：
Figure BDA00030458713700001111
推出：
Figure BDA00030458713700001112
同时，第i个接收机到由第m个干扰机的第n个时延干扰信号产生的虚假目标的距离历史和表示为：
Figure BDA0003045871370000121
其中，
Figure BDA0003045871370000122
表示第i幅SAR图像中该干扰机的第n个时延干扰信号生成的虚假目标位置；
得到
Figure BDA0003045871370000123
表示为：
Figure BDA0003045871370000124
其中，
Figure BDA0003045871370000125
表示第m个干扰机到发射机和第i个接收机的距离历史和；
通过上述分析发现第i个接收机成的SAR图像中，第m个干扰机第n个时延
Figure BDA0003045871370000126
转发的干扰信号产生的虚假目标位于
Figure BDA0003045871370000127
对应的等距离线与干扰机所在的等多普勒线的交点处；其中，
Figure BDA0003045871370000128
表示第m个干扰机到发射机和第i个接收机的距离历史和；
干扰时延不同，其对应的等距离线与干扰机所在等多普勒线的交点位置也就不同，即不同干扰时延在SAR图像中形成的虚假目标位置不同。根据距离欺骗干扰的原理可知，我们可以发现随着干扰时延的变化，虚假目标沿着干扰机所在的等多普勒线分布。所以，在多基SAR图像中，将分布在同一条等多普勒线上的虚假目标可视为同一干扰机生成。
S4、进行干扰机定位，并利用波束形成的方法对干扰信号进行抑制；具体实现方法为：
S41、虚假目标与干扰机对应关系组合：假设第m个干扰机发射N个不同时延的干扰信号，各接收机所成SAR图像中产生多个虚假目标；
组合规则为：在I个接收机接收到的I幅SAR图像中分别任意选择一个虚假目标，把这些虚假目标视为一个小组；已被选择过的虚假目标不能再被选择，继续选择虚假目标形成小组，直到所有虚假目标都被选完，每种组合包含N个小组，每个小组包含I个虚假目标，N个小组形成一种组合情况；
S42、定位方程组求解：虚假目标均满足公式(20)，将每个小组中每个SAR图像中的虚假目标参数带入公式(20)，得到有N个方程组成的方程组，通过求解方程组，得到N 个干扰机位置及N个时延大小；
S43、干扰机的位置判定：每种组合情况解得的N个干扰机位置分别表示为
Figure BDA0003045871370000131
其中
Figure BDA0003045871370000133
表示该组合情况中第n个方程组所解得的干扰机位置；当
Figure BDA0003045871370000132
时，则该组合情况为正确组合，即说明同一小组的虚假目标是由第m个干扰机的同一个时延转发干扰信号产生的，解得的干扰机位置为干扰机的正确位置；如果不满足，则继续进行下一个组合情况的判断，直到找到正确的组合为止。
图4为本实施例含有虚假目标的成像结果图，其中(a)、(b)、(c)和(d)分别对应接收站1、接收站2、接收站3和接收站4，图中，方框内为虚假目标。图5为消除虚假目标后的成像结果，其中(a)、(b)、(c)和(d)分别对应接收站1、接收站2、接收站3和接收站4。通过本发明具体实施方式可以看出，本发明可以实现针对多基SAR系统的抗欺骗性干扰，抗干扰效果明显。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (5)

1.一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，其特征在于，包括以下步骤：
S1、建立回波模型；
S2、虚假目标的检测与判别：利用基于谱残差的视觉显著性检测方法与形态学方法相结合实现多基SAR图像中的目标检测，然后利用基于欧式距离准则的判别法判定目标是否虚假；
S3、对虚假目标位置分布规律进行分析；
S4、进行干扰机定位，并利用波束形成的方法对干扰信号进行抑制。
2.根据权利要求1所述的一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，其特征在于，所述步骤S1具体实现方法为：多基SAR系统包括一个发射机、多个接收机和干扰机，发射机发射线性调频信号；
方位时刻为η时，第i个接收机接收到的真实环境中点目标回波表示如下：
Figure FDA0003045871360000011
其中，A0表示点目标雷达横截面面积，wr和wa分别表示距离向和方位向的包络函数，τ表示距离向时间，η表示方位向时间，ηc表示中心时刻；f0表示载波频率，Kr表示调频速率，c表示光速，R(i)(η)表示在η时刻发射机到点目标的距离RT(η)及第i个接收机到点目标的距离
Figure FDA0003045871360000012
的距离和；
第i个接收机的距离历史和R(i)(η)的表达式如下：
Figure FDA0003045871360000013
其中，
Figure FDA0003045871360000014
表示第i个接收机的位置，(xT(η),yT(η),zT(η))表示发射机的位置，(xq,yq,zq)表示成像环境中点目标的位置；
干扰机通过参数调制，进行时延处理后，再把干扰信号发射出去；由第m个干扰机进行第n个时延转发的干扰信号在η时刻在第i个接收机形成的干扰瞬时回波
Figure FDA0003045871360000015
表示如下：
Figure FDA0003045871360000021
其中，
Figure FDA0003045871360000022
表示第m个干扰机的第n个时延，
Figure FDA0003045871360000023
表示在η时刻第m个干扰机到发射机及第i个接收机的距离历史和，表示如下：
Figure FDA0003045871360000024
其中，
Figure FDA0003045871360000025
表示第m个干扰机的位置；
在多干扰机环境下，多基SAR的接收机会同时受到多个干扰机发射的干扰信号影响；所以，第i个接收机在η时刻收到的回波s(i)(τ,η)如下所示：
Figure FDA0003045871360000026
其中，M表示该干扰环境共有M个干扰机，N表示第M个干扰机发射N个不同时延的干扰信号；
利用后向投影算法将接收到的回波信号生成多基SAR图像。
3.根据权利要求1所述的一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，其特征在于，所述步骤S2中，基于谱残差的视觉显著性检测方法的实现步骤如下：
S21、给定SAR图像I(x,y)，其频谱表示为
Figure FDA0003045871360000027
其中
Figure FDA0003045871360000028
表示二维傅里叶变换；
对应的幅度谱和相位谱分别是：
Figure FDA0003045871360000029
Figure FDA0003045871360000031
其中，
Figure FDA0003045871360000032
表示输入频谱的幅度计算，
Figure FDA0003045871360000033
表示输入频谱的相位计算；
对SAR图像的幅度谱进行自然对数运算，得到对应的对数谱，如下表示：
L(fx,fy)＝In[A(fx,fy)] (9)；
S22、利用均值滤波器进行平滑滤波，获得SAR图像中的均值谱，表示为：
Figure FDA0003045871360000034
其中，h(fx,fy)是一个n×n的矩阵，h(fx,fy)为频谱均值滤波器，具体表示为：
Figure FDA0003045871360000035
S23、用对数谱减去均值谱，得到SAR图像的谱残差：
R(fx,fy)＝L(fx,fy)-B(fx,fy) (12)；
S24、经过二维逆傅里叶变换和高斯平滑滤波，将其还原到图像域，得到SAR图像中的新颖部分，即SAR图像的显著图：
Figure FDA0003045871360000036
S25、对SAR图像的显著图进行开运算，然后计算SAR图像的显著图中目标点的面积，与设定的面积阈值进行比较，若大于或等于面积阈值则将该目标作为感兴趣区域，否则无视该目标；
S26、利用基于欧式距离判别法实现对虚假目标的判定；首先，从这些图像中任选两幅不同的图像Mi和Mj；设Ci,k为第i幅图像中的第k个目标的重心，Cj,h表示第j幅图像中的第h个目标的几何中心，γ为欧式距离门限；依次计算Ci,k和Cj,h的欧氏距离，若Ci,k可以在第j幅图像中找到满足欧式距离小于γ的，说明Ci,k为真目标，反之为虚假目标。
4.根据权利要求1所述的一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，其特征在于，所述步骤S3具体实现方法为：假设干扰环境由M个发射多个不同时延的干扰机组成，第m个干扰机发射N个不同时延干扰信号，该干扰机在每幅SAR图像中会产生N 个不同的虚假目标；根据距离欺骗干扰原理得到干扰机在η时刻的多普勒频率为：
Figure FDA0003045871360000041
其中，
Figure FDA0003045871360000042
和
Figure FDA0003045871360000043
分别表示发射机和第i个接收机的速度矢量，λ表示干扰机的发射波长，
Figure FDA0003045871360000044
和
Figure FDA0003045871360000045
分别表示发射机和第i个接收机的位置矢量，
Figure FDA0003045871360000046
表示第m个干扰机的位置矢量；
由该干扰机生成的虚假目标与该干扰机的瞬时多普勒相同，如下表示：
Figure FDA0003045871360000047
Figure FDA0003045871360000048
其中，
Figure FDA0003045871360000049
表示在η时刻，第i个接收机受到第m个干扰机第n个时延干扰信号产生的虚假目标的瞬时多普勒频率；
Figure FDA00030458713600000410
表示在第i个接收机所成的SAR图像中，第m个干扰机的第n个时延干扰信号产生的虚假目标位置矢量，n＝1,2,…,N；由此，N个不同时延干扰信号同时满足下列方程组：
Figure FDA00030458713600000411
推出：
Figure FDA00030458713600000412
同时，第i个接收机到由第m个干扰机的第n个时延干扰信号产生的虚假目标的距离历史和表示为：
Figure FDA00030458713600000413
其中，
Figure FDA00030458713600000414
表示第i幅SAR图像中该干扰机的第n个时延干扰信号生成的虚假目标位置；
得到
Figure FDA0003045871360000051
表示为：
Figure FDA0003045871360000052
其中，
Figure FDA0003045871360000053
表示第m个干扰机到发射机和第i个接收机的距离历史和；
通过上述分析发现第i个接收机成的SAR图像中，第m个干扰机第n个时延
Figure FDA0003045871360000054
转发的干扰信号产生的虚假目标位于
Figure FDA0003045871360000055
对应的等距离线与干扰机所在的等多普勒线的交点处；其中，
Figure FDA0003045871360000056
Figure FDA0003045871360000057
表示第m个干扰机到发射机和第i个接收机的距离历史和；
随着干扰时延的变化，虚假目标沿着干扰机所在的等多普勒线分布；所以，在多基SAR图像中，将分布在同一条等多普勒线上的虚假目标视为同一干扰机生成。
5.根据权利要求4所述的一种多基SAR在多干扰机环境下的抗欺骗型干扰方法，其特征在于，所述步骤S4进行干扰机定位的具体实现方法为：
S41、虚假目标与干扰机对应关系组合：假设第m个干扰机发射N个不同时延的干扰信号，各接收机所成SAR图像中产生多个虚假目标；
组合规则为：在I个接收机接收到的I幅SAR图像中分别任意选择一个虚假目标，把这些虚假目标视为一个小组；已被选择过的虚假目标不能再被选择，继续选择虚假目标形成小组，直到所有虚假目标都被选完，每种组合包含N个小组，每个小组包含I个虚假目标，N个小组形成一种组合情况；
S42、定位方程组求解：虚假目标均满足公式(20)，将每个小组中每个SAR图像中的虚假目标参数带入公式(20)，得到有N个方程组成的方程组，通过求解方程组，得到N个干扰机位置及N个时延大小；
S43、干扰机的位置判定：每种组合情况解得的N个干扰机位置分别表示为
Figure FDA0003045871360000059
其中
Figure FDA0003045871360000059
表示该组合情况中第n个方程组所解得的干扰机位置；当
Figure FDA00030458713600000510
时，则该组合情况为正确组合，即说明同一小组的虚假目标是由第m个干扰机的同一个时延转发干扰信号产生的，解得的干扰机位置为干扰机的正确位置；如果不满足，则继续进行下一个组合情况的判断，直到找到正确的组合为止。
CN202110472443.0A 2021-04-29 2021-04-29 一种多基sar在多干扰机环境下的抗欺骗型干扰方法 Active CN113203991B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110472443.0A CN113203991B (zh) 2021-04-29 2021-04-29 一种多基sar在多干扰机环境下的抗欺骗型干扰方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110472443.0A CN113203991B (zh) 2021-04-29 2021-04-29 一种多基sar在多干扰机环境下的抗欺骗型干扰方法

Publications (2)

Publication Number Publication Date
CN113203991A true CN113203991A (zh) 2021-08-03
CN113203991B CN113203991B (zh) 2022-05-31

Family

ID=77027774

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110472443.0A Active CN113203991B (zh) 2021-04-29 2021-04-29 一种多基sar在多干扰机环境下的抗欺骗型干扰方法

Country Status (1)

Country Link
CN (1) CN113203991B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114895264A (zh) * 2022-06-23 2022-08-12 南京慧尔视智能科技有限公司 雷达虚假目标的抑制方法、装置、电子设备和存储介质

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITMI991154A1 (it) * 1999-05-25 2000-11-25 Milano Politecnico Procedimento per misure radar di spostamento di aere urbane e zone franose
US6650271B1 (en) * 1997-11-24 2003-11-18 Raytheon Company Signal receiver having adaptive interfering signal cancellation
EP2660623A2 (en) * 2012-09-03 2013-11-06 Institute of Electronics, Chinese Academy of Sciences Imaging method and device in SAB mobile bistatic SAR
CN103983974A (zh) * 2014-05-30 2014-08-13 电子科技大学 双站调频连续波合成孔径雷达成像方法
CN105005032A (zh) * 2015-07-21 2015-10-28 电子科技大学 一种基于级数反演成像算法的sar移频干扰方法
CN105354541A (zh) * 2015-10-23 2016-02-24 西安电子科技大学 基于视觉注意模型和恒虚警率的sar图像目标检测方法
CN105372637A (zh) * 2015-10-28 2016-03-02 西安电子科技大学 基于运动车辆目标微动特性的对雷达的窄带欺骗干扰方法
CN106291494A (zh) * 2016-07-21 2017-01-04 深圳大学 基于微分特征增强的sar欺骗干扰目标标识方法及系统
CN106772273A (zh) * 2016-11-11 2017-05-31 深圳大学 一种基于动态孔径的sar虚假目标干扰抑制方法及系统
CN106990396A (zh) * 2017-05-25 2017-07-28 电子科技大学 一种双基地合成孔径雷达定位误差计算方法
CN107247930A (zh) * 2017-05-26 2017-10-13 西安电子科技大学 基于cnn和选择性注意机制的sar图像目标检测方法
CN107274401A (zh) * 2017-06-22 2017-10-20 中国人民解放军海军航空工程学院 一种基于视觉注意机制的高分辨率sar图像舰船检测方法
CN107831475A (zh) * 2017-11-06 2018-03-23 电子科技大学 一种基于阴影特征的合成孔径雷达欺骗式干扰方法
CN108020836A (zh) * 2018-01-10 2018-05-11 电子科技大学 双基地合成孔径雷达运动目标定位方法
KR20190043478A (ko) * 2017-10-18 2019-04-26 서울대학교산학협력단 고해상도 수치 표고 모델 생성 시스템 및 방법
CN110058233A (zh) * 2019-04-28 2019-07-26 电子科技大学 一种多基地合成孔径雷达系统的抗欺骗性干扰方法
CN110109113A (zh) * 2019-05-07 2019-08-09 电子科技大学 一种基于级联对消的双基前视sar非平稳杂波抑制方法
CN110109117A (zh) * 2019-05-30 2019-08-09 电子科技大学 基于频控阵的星载合成孔径雷达卷积调制干扰方法
CN111650586A (zh) * 2020-04-03 2020-09-11 中国科学院电子学研究所苏州研究院 基于星载sar方位向nlcs成像的目标定位方法
CN112346087A (zh) * 2020-11-04 2021-02-09 上海交通大学 多峰检测和增益监测结合的gnss欺骗检测方法及系统

Patent Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6650271B1 (en) * 1997-11-24 2003-11-18 Raytheon Company Signal receiver having adaptive interfering signal cancellation
ITMI991154A1 (it) * 1999-05-25 2000-11-25 Milano Politecnico Procedimento per misure radar di spostamento di aere urbane e zone franose
EP2660623A2 (en) * 2012-09-03 2013-11-06 Institute of Electronics, Chinese Academy of Sciences Imaging method and device in SAB mobile bistatic SAR
CN103983974A (zh) * 2014-05-30 2014-08-13 电子科技大学 双站调频连续波合成孔径雷达成像方法
CN105005032A (zh) * 2015-07-21 2015-10-28 电子科技大学 一种基于级数反演成像算法的sar移频干扰方法
CN105354541A (zh) * 2015-10-23 2016-02-24 西安电子科技大学 基于视觉注意模型和恒虚警率的sar图像目标检测方法
CN105372637A (zh) * 2015-10-28 2016-03-02 西安电子科技大学 基于运动车辆目标微动特性的对雷达的窄带欺骗干扰方法
CN106291494A (zh) * 2016-07-21 2017-01-04 深圳大学 基于微分特征增强的sar欺骗干扰目标标识方法及系统
CN106772273A (zh) * 2016-11-11 2017-05-31 深圳大学 一种基于动态孔径的sar虚假目标干扰抑制方法及系统
CN106990396A (zh) * 2017-05-25 2017-07-28 电子科技大学 一种双基地合成孔径雷达定位误差计算方法
CN107247930A (zh) * 2017-05-26 2017-10-13 西安电子科技大学 基于cnn和选择性注意机制的sar图像目标检测方法
CN107274401A (zh) * 2017-06-22 2017-10-20 中国人民解放军海军航空工程学院 一种基于视觉注意机制的高分辨率sar图像舰船检测方法
KR20190043478A (ko) * 2017-10-18 2019-04-26 서울대학교산학협력단 고해상도 수치 표고 모델 생성 시스템 및 방법
CN107831475A (zh) * 2017-11-06 2018-03-23 电子科技大学 一种基于阴影特征的合成孔径雷达欺骗式干扰方法
CN108020836A (zh) * 2018-01-10 2018-05-11 电子科技大学 双基地合成孔径雷达运动目标定位方法
CN110058233A (zh) * 2019-04-28 2019-07-26 电子科技大学 一种多基地合成孔径雷达系统的抗欺骗性干扰方法
CN110109113A (zh) * 2019-05-07 2019-08-09 电子科技大学 一种基于级联对消的双基前视sar非平稳杂波抑制方法
CN110109117A (zh) * 2019-05-30 2019-08-09 电子科技大学 基于频控阵的星载合成孔径雷达卷积调制干扰方法
CN111650586A (zh) * 2020-04-03 2020-09-11 中国科学院电子学研究所苏州研究院 基于星载sar方位向nlcs成像的目标定位方法
CN112346087A (zh) * 2020-11-04 2021-02-09 上海交通大学 多峰检测和增益监测结合的gnss欺骗检测方法及系统

Non-Patent Citations (14)

* Cited by examiner, † Cited by third party
Title
BO ZHAO;FENG ZHOU;ZHENG BAO: "Deception Jamming for Squint SAR Based on Multiple Receivers", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》 *
BO ZHAO;FENG ZHOU;ZHENG BAO: "Deception Jamming for Squint SAR Based on Multiple Receivers", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》, vol. 8, 29 May 2014 (2014-05-29), pages 3988 - 3998, XP055312960, DOI: 10.1109/JSTARS.2014.2322612 *
QINGYANG SUN;TING SHU;KAI-BOR YU;WENXIAN YU: "A Novel Deceptive Jamming Method Against Two-Channel SAR-GMTI Based on Two Jammers", 《IEEE SENSORS JOURNAL》 *
QINGYANG SUN;TING SHU;KAI-BOR YU;WENXIAN YU: "A Novel Deceptive Jamming Method Against Two-Channel SAR-GMTI Based on Two Jammers", 《IEEE SENSORS JOURNAL》, vol. 19, 28 March 2019 (2019-03-28), pages 5600 - 5610, XP011731502, DOI: 10.1109/JSEN.2019.2908030 *
WENJING WANG;JUNJIE WU;JIFANG PEI;XINYU MAO;JIANYU YANG: "An Antideceptive Jamming Method for Multistatic", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》 *
WENJING WANG;JUNJIE WU;JIFANG PEI;XINYU MAO;JIANYU YANG: "An Antideceptive Jamming Method for Multistatic", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》, vol. 13, 31 December 2020 (2020-12-31), pages 2757 - 2768, XP011793528, DOI: 10.1109/JSTARS.2020.2997345 *
X. HE, J. ZHU, J. WANG, D. DU AND B. TANG: "False Target Deceptive Jamming for Countering Missile-Borne SAR", 《2014 IEEE 17TH INTERNATIONAL CONFERENCE ON COMPUTATIONAL SCIENCE AND ENGINEERING》 *
X. HE, J. ZHU, J. WANG, D. DU AND B. TANG: "False Target Deceptive Jamming for Countering Missile-Borne SAR", 《2014 IEEE 17TH INTERNATIONAL CONFERENCE ON COMPUTATIONAL SCIENCE AND ENGINEERING》, 19 January 2015 (2015-01-19), pages 974 - 1978 *
YUNPENG, ZHANG第: "Azimuth Intermittent Sampling Time-delay Repeater Jamming Method Against SAR-GMTI", 《JOURNAL OF DETECTION&CONTROL》 *
YUNPENG, ZHANG第: "Azimuth Intermittent Sampling Time-delay Repeater Jamming Method Against SAR-GMTI", 《JOURNAL OF DETECTION&CONTROL》, 12 December 2017 (2017-12-12), pages 84 - 90 *
裴季方: "多视角SAR目标识别方法研究", 《中国博士学位论文全文数据库》 *
裴季方: "多视角SAR目标识别方法研究", 《中国博士学位论文全文数据库》, 15 March 2019 (2019-03-15), pages 136 - 42 *
闫成章,刘畅: "基于显著性的SAR图像船舶目标检测方法", 《中国科学院大学学报》 *
闫成章,刘畅: "基于显著性的SAR图像船舶目标检测方法", 《中国科学院大学学报》, 15 May 2019 (2019-05-15), pages 401 - 409 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114895264A (zh) * 2022-06-23 2022-08-12 南京慧尔视智能科技有限公司 雷达虚假目标的抑制方法、装置、电子设备和存储介质

Also Published As

Publication number Publication date
CN113203991B (zh) 2022-05-31

Similar Documents

Publication Publication Date Title
CN107037410B (zh) 一种对雷达进行干扰的方法、装置及频控阵干扰机
CN110412559B (zh) 分布式无人机mimo雷达的非相参融合目标检测方法
CN103744068B (zh) 双通道调频连续波sar系统的动目标检测成像方法
Huang et al. Radar high speed small target detection based on keystone transform and linear canonical transform
CN106443671A (zh) 基于调频连续波的sar雷达动目标检测与成像方法
CN107728115A (zh) 一种雷达目标成像后基于svm的背景干扰抑制方法
Chabriel et al. Adaptive target detection techniques for OFDM-based passive radar exploiting spatial diversity
Daun et al. Tracking algorithms for multistatic sonar systems
CN104155632A (zh) 一种基于局部相关性的改进子空间海杂波抑制方法
CN103913725A (zh) 密集转发式干扰环境下的机载雷达地面动目标检测方法
CN104535972A (zh) 一种机载雷达相干转发式干扰抑制方法
CN104597434A (zh) 改进包络移位补偿和分数阶傅里叶变换的多帧相参tbd方法
CN114924236B (zh) 基于位置和速度信息的空地雷达协同抗欺骗式干扰方法
CN113203991B (zh) 一种多基sar在多干扰机环境下的抗欺骗型干扰方法
Lee et al. A preliminary study of machine-learning-based ranging with LTE channel impulse response in multipath environment
Wang et al. Antirange-deception jamming from multijammer for multistatic SAR
Chen et al. Improved target detection through DNN-based multi-channel interference mitigation in automotive radar
CN112835006A (zh) 一种基于帧间积累的跟踪雷达海上小目标检测方法及系统
Temiz et al. Improved target localization in multiwaveform multiband hybrid multistatic radar networks
CN113203992A (zh) 一种多基sar的抗欺骗性干扰方法
Wang et al. A novel anti-deceptive jamming method for multistatic sar
Bennaceur et al. Target localization in depth and range from passive sonar
Teng Fundamental aspects of netted radar performance
Makhoul et al. Evaluation of state-of-the-art GMTI techniques for future spaceborne SAR systems-Simulation validation
Nikaein et al. Multitarget detection in passive MIMO radar using block sparse recovery

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant