CN117784040A - 一种时频域抗自卫式复合干扰方法 - Google Patents

一种时频域抗自卫式复合干扰方法 Download PDF

Info

Publication number
CN117784040A
CN117784040A CN202410202388.7A CN202410202388A CN117784040A CN 117784040 A CN117784040 A CN 117784040A CN 202410202388 A CN202410202388 A CN 202410202388A CN 117784040 A CN117784040 A CN 117784040A
Authority
CN
China
Prior art keywords
interference
data
time
frequency
target
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
CN202410202388.7A
Other languages
English (en)
Other versions
CN117784040B (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.)
Air Force Early Warning Academy
Original Assignee
Air Force Early Warning Academy
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 Air Force Early Warning Academy filed Critical Air Force Early Warning Academy
Priority to CN202410202388.7A priority Critical patent/CN117784040B/zh
Publication of CN117784040A publication Critical patent/CN117784040A/zh
Application granted granted Critical
Publication of CN117784040B publication Critical patent/CN117784040B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供一种时频域抗自卫式复合干扰方法,利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;最后得到目标的距离信息;由于发射信号的调频率是已知的,不仅可以避免了估计干扰调频率的复杂过程,从而可以采用先期估计滤波器,并进行存贮。

Description

一种时频域抗自卫式复合干扰方法
技术领域
本发明涉及各类发射线性调频信号技术领域,尤其涉及一种时频域抗自卫式复合干扰方法。
背景技术
由于近年来,由于电子对抗技术飞速发展,对抗雷达的各处干扰信号被开发出来,并得到广泛的运用。其中切片干扰和频谱弥散干扰是最常用的自卫式干扰,它们也是目前最有效的对抗线性调频信号的干扰。这类干扰充分利用了线性调频信号的时频特性,产生大量的虚假目标从而压制真实目标,对雷达抗干扰带来了很大的难度。当前,雷达最有效的对抗干扰的方法就是频率捷变,但由于雷达工作需要毫秒级的驻留时间,而切片转发式干扰可以在微秒量级跟上雷达的频率,从而导致雷达的这一对抗措施失效。
目前,对抗切片干扰最常用的是时域匿影方法,但该方法只适用于只存在切片转发式干扰的场景,且目标信号的能量会有严重损失,所以其使用条件受限。另外,对抗频谱弥散干扰则采用频域旋转的方法,再通过低通过滤波或带通滤波的方法来抑制干扰,这种方法需要首先估计出调制的干扰斜率。但在自卫式切片和频谱弥散混合干扰的场景下,由于切片干扰通常和噪声进行卷积处理,导致无法估计出干扰的调制斜率,同时由于频谱弥散干扰和切片干扰不同步,导致无法从时域找出切片干扰,这就导致现有的只能抗切片干扰或者抗频谱弥散干扰的方法无法适用于二种干扰均存在的混合干扰场景。
发明内容
本发明的目的在于针对上述现有技术的不足,提供了一种时频域抗自卫式复合干扰方法,通过构造时频域二维滤波器、通过滤波、搜索估计再重构的方式实现主瓣自卫式切片和频谱弥散混合干扰的抑制。
为实现上述目的,本发明采用了如下技术方案:
本发明提供了一种时频域抗自卫式复合干扰方法,包括以下步骤:
S1、利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;
S2、利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;
S3、最后得到目标的距离信息;
所述S1具体为:
S101、接收雷达的和通道的数据进行分析,从数据中提取有干扰的脉冲数据段:
(1);
其中, 为目标信号的回波幅度;为雷达的发射信号;为虚数;为数 据采样时间;为目标的多普勒频率;为自卫式干扰发射的切片干扰;为自 卫式干扰发射的频谱弥散干扰;为噪声;
雷达的发射信号的信号波形为:
(2);
雷达的发射信号能平分成份:
(3);
自卫式干扰发射的切片干扰的波形为:
(4);
自卫式干扰发射的频谱弥散干扰的波形为:
(5);
其中,为发射信号的起始频率;为线性调频信号的斜率;为发射脉冲的 周期;分别为二种干扰的幅度;为雷达发射信号的数据分块;为扰调制 噪声由干扰机按时序设置;为数据分块的采样时间;为频谱弥散信号的数据分块;为干扰机将干扰复制了份;为频谱弥散干扰的调频斜率;符号表示卷积运算;
S102、对雷达的发射信号和自卫式干扰发射的频谱弥散干扰进行短 时傅立叶变换:
(6);
对二维矩阵数据取绝对值,并进行归一化处理:
(7);
对二维矩阵数据进行归一化处理:
(8);
其中,均为二维矩阵,维数为为时 间维采样的长度,为频率维采样的长度,为已知信息;为门限;为短时 傅立叶变换;为所得到的滤波器;
所述S3具体为:对所述时域数据进行检测,最大值对应的距离门就是目标所在的距离门,将目标的距离信息输出。
进一步,所述S2具体为:
S201、根据已知目标的跟踪信息确定多普勒搜索范围和搜索步长,搜索范围为:
(9);
搜索步长为:
(10);
其中,多普勒估计值为为航迹的径向角,为目标的速度,为工作波长,为多普勒估计误差;其中,均为已知信息;为搜 索的次数;
对提取的数据进行次的时频域滤波:
(11);
其中,为点乘,为短时傅立叶变换的反变换,
S202、对组时频域滤波的数据取绝对值求和,得到个数据,通过 搜索找出其中的最大值,并得到最大值对应的多普勒频率,则重构的数据为:
(12);
S203、对重构的数据进行特异点剔除,即:
(13);
其中,为门限;
S204、对滤波后的数据进行脉冲压缩,得到时域数据为:
(14);
其中,为脉冲压缩。
进一步,所述多普勒估计误差能够按整数倍来设置。
进一步,所述最大值对应的多普勒频率能直接对最大值周围数据进行幅度 拟合得到。
本发明的有益效果为:(1)由于发射信号的调频率是已知的,不仅可以避免了估计干扰调频率的复杂过程,从而可以采用先期估计滤波器,并进行存贮;
(2)由于发射信号已知,其时频域滤波可以构造得很理想,则很容易抑制大部分的带外干扰,而不需要通过设计专门的滤波器,实现简单方便;
(3)通过带外的滤波之后只剩余少量的带内干扰,此时带内干扰在时域上占的区域比例少且幅度远大于信号,所以通过时域方法可以很容易滤除剩余的频谱弥散干扰;
(4)该方法可用于改造现有雷达的信号处理系统,不需要额外增加处理通道和设备,只需对数字化雷达的进行算法升级即可。
附图说明
图1 为本发明一种时频域抗自卫式复合干扰方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,下面结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
假设雷达发射线性调频信号,干扰机发射切片转发和频谱弥散混合式干扰,实施例中信号的归一化的调频率,干扰的归一化的调频率/>,/>Hz,,/>,/>Hz,/>米/秒,/>米,/>度,/>。下面结合实施例说明一下整个发明的详细步骤:
请参阅图1,一种时频域抗自卫式复合干扰方法,包括以下步骤:
S1、利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;
S2、利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;
S3、最后得到目标的距离信息。
所述S1具体为:
S101、接收雷达的和通道的数据进行分析,从数据中提取有干扰的脉冲数据段:
(1);
其中,为目标信号的回波幅度;为雷达的发射信号;为虚数;为数据 采样时间;为目标的多普勒频率;为自卫式干扰发射的切片干扰;为自卫 式干扰发射的频谱弥散干扰;为噪声;
雷达的发射信号的信号波形为:
(2);
雷达的发射信号能平分成份:
(3);
自卫式干扰发射的切片干扰的波形为:
(4);
自卫式干扰发射的频谱弥散干扰的波形为:
(5);
其中,为发射信号的起始频率;为线性调频信号的斜率;为发射脉冲的 周期;分别为二种干扰的幅度;为雷达发射信号的数据分块;为扰调制 噪声由干扰机按时序设置;为数据分块的采样时间;为频谱弥散信号的数据分块;为干扰机将干扰复制了份;为频谱弥散干扰的调频斜率;符号表示卷积运算
在一种具体的实施例中,,切片干扰/>,频谱弥散干扰/>,说明整个波形被切片干扰分成了50份进行调制,每一份卷积一个随机噪声,再构成切片干扰。频谱弥散干扰/>,分成了15份进行调制,每一份是相同的。将二种干扰波形相加则形成干扰的数据。
S102、对雷达的发射信号和自卫式干扰发射的频谱弥散干扰进行短 时傅立叶变换:
(6);
对二维矩阵数据取绝对值,并进行归一化处理:
(7);
对二维矩阵数据进行归一化处理:
(8);
其中,均为二维矩阵,维数为为时 间维采样的长度,为频率维采样的长度,为已知信息;为门限;为短时 傅立叶变换;为所得到的滤波器。
在一种具体的实施例中,假设,即大于0.5的被置 为1,其它情况均为0。另外,的维数均为
所述S2具体为:
S201、根据已知目标的跟踪信息确定多普勒搜索范围和搜索步长,
搜索范围为:
(9);
搜索步长为:
(10);
其中,多普勒估计值为为航迹的径向角,为目标的速度,为工作波长,为多普勒估计误差;其中,均为已知信息;为搜 索的次数;
对提取的数据进行次的时频域滤波:
(11);
其中,为点乘,为短时傅立叶变换的反变换,
在一种具体的实施例中,Hz,/>米/秒,/>米,/>度,则/>Hz,所以/>的搜索范围为/>,假设/>,则/>Hz,说明在280至320之间按1Hz步长进行搜索。
S202、对组时频域滤波的数据取绝对值求和,得到个数据,通过 搜索找出其中的最大值,并得到最大值对应的多普勒频率,则重构的数据为:
(12);
在一种具体的实施例中,共对41组数据进行取绝对值求和,并找出41个数据中最大值,由这个最大值确定多普勒频率,由于多普勒频率为310Hz,则第31个数据的幅度最大,对应的多普勒频率为310Hz,再利用这个值对数据进行重构,公式如下
S203、对重构的数据进行特异点剔除,即:
(13);
其中,为门限;
在一种具体的实施方式中,为/>所有数据的平均值的10倍,即大于平均值10倍的值被置零,其于值保留原值。
S204、对滤波后的数据进行脉冲压缩,得到时域数据为:
(14);
其中,为脉冲压缩。
所述S3具体为:对所述时域数据进行检测,最大值对应的距离门就是目标 所在的距离门,将目标的距离信息输出;其中,检测可以单元平均CFAR实现,也可以采用选 大、选小、排序等CFAR算法实现;实施例中采用单元平均CFAR方法。
所述多普勒估计误差能够按整数倍来设置。能够进一步增加多普勒的容差 范围;实施例中,是按1倍设置的,如容错范围需扩大至2倍,则可以设计为2;
所述最大值对应的多普勒频率能直接对最大值周围数据进行幅度拟合得 到,能够进一步增加估计的精度,实施例中,采用幅度最大的方式进行,没有采用拟合的方 式,如果需要提高精度,可以采用拟合的方式。
其中,S204中,脉冲压缩可以采用时域方法,也可以用频域方法。实施例中采用时域方法进行卷积处理得到。
以上所述实施例仅表达了本发明的实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求。

Claims (4)

1.一种时频域抗自卫式复合干扰方法,其特征在于,包括以下步骤:
S1、利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;
S2、利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;
S3、最后得到目标的距离信息;
所述S1具体为:
S101、接收雷达的和通道的数据进行分析,从数据中提取有干扰的脉冲数据段:
(1);
其中, 为目标信号的回波幅度;/>为雷达的发射信号;/>为虚数;/>为数据采样时间;/>为目标的多普勒频率;/>为自卫式干扰发射的切片干扰;/>为自卫式干扰发射的频谱弥散干扰;/>为噪声;
雷达的发射信号的信号波形为:
(2);
雷达的发射信号能平分成/>份:
(3);
自卫式干扰发射的切片干扰的波形为:
(4);
自卫式干扰发射的频谱弥散干扰的波形为:
(5);
其中,为发射信号的起始频率;/>为线性调频信号的斜率;/>为发射脉冲的周期;分别为二种干扰的幅度;/>为雷达发射信号的数据分块;/>为扰调制噪声由干扰机按时序设置;/>为数据分块的采样时间;/>为频谱弥散信号的数据分块;/>为干扰机将干扰复制了/>份;/>为频谱弥散干扰的调频斜率;符号/>表示卷积运算;
S102、对雷达的发射信号和自卫式干扰发射的频谱弥散干扰/>进行短时傅立叶变换:
(6);
对二维矩阵数据取绝对值,并进行归一化处理:
(7);
对二维矩阵数据进行归一化处理:
(8);
其中,、/>和/>均为二维矩阵,维数为/>,/>为时间维采样的长度,/>为频率维采样的长度,/>为已知信息;/>为门限;/>为短时傅立叶变换;/>为所得到的滤波器;
所述S3具体为:对所述时域数据进行检测,最大值对应的距离门就是目标所在的距离门,将目标的距离信息输出。
2.根据权利要求1所述的一种时频域抗自卫式复合干扰方法,其特征在于,所述S2具体为:
S201、根据已知目标的跟踪信息确定多普勒搜索范围和搜索步长,搜索范围为:
(9);
搜索步长为:
(10);
其中,多普勒估计值为;/>为航迹的径向角,/>为目标的速度,/>为工作波长,/>为多普勒估计误差;其中,/>、/>、/>和/>均为已知信息;/>为搜索的次数;
对提取的数据进行次的时频域滤波:
(11);
其中,为点乘,/>为短时傅立叶变换的反变换,,/>
S202、对组时频域滤波/>的数据取绝对值求和,得到/>个数据,通过搜索找出其中的最大值,并得到最大值对应的多普勒频率/>,则重构的数据为:
(12);
S203、对重构的数据进行特异点剔除,即:
(13);
其中,为门限;
S204、对滤波后的数据进行脉冲压缩,得到时域数据为:
(14);
其中,为脉冲压缩。
3.根据权利要求2所述的一种时频域抗自卫式复合干扰方法,其特征在于,所述多普勒估计误差能够按整数倍来设置。
4.根据权利要求2所述的一种时频域抗自卫式复合干扰方法,其特征在于,所述最大值对应的多普勒频率能直接对最大值周围数据进行幅度拟合得到。
CN202410202388.7A 2024-02-23 2024-02-23 一种时频域抗自卫式复合干扰方法 Active CN117784040B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410202388.7A CN117784040B (zh) 2024-02-23 2024-02-23 一种时频域抗自卫式复合干扰方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410202388.7A CN117784040B (zh) 2024-02-23 2024-02-23 一种时频域抗自卫式复合干扰方法

Publications (2)

Publication Number Publication Date
CN117784040A true CN117784040A (zh) 2024-03-29
CN117784040B CN117784040B (zh) 2024-05-10

Family

ID=90391219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410202388.7A Active CN117784040B (zh) 2024-02-23 2024-02-23 一种时频域抗自卫式复合干扰方法

Country Status (1)

Country Link
CN (1) CN117784040B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5790475A (en) * 1996-10-28 1998-08-04 Multispec Corporation Process and apparatus for improved interference suppression in echo-location and imaging systems
UA89945U (ru) * 2013-10-14 2014-05-12 Запорізький Національний Технічний Університет Способ помехозащиты когерентно-импульсной рлс от комбинированных помех
CN113176542A (zh) * 2021-03-24 2021-07-27 中国人民解放军空军预警学院 一种时域自适应抗主瓣自卫式干扰方法及系统
US20210231789A1 (en) * 2018-06-28 2021-07-29 Richwave Technology Corp. Doppler signal processing device and method thereof for interference suppression
CN113702921A (zh) * 2021-08-17 2021-11-26 西安电子科技大学 基于时频分析的干扰抑制方法
CN115453469A (zh) * 2022-09-21 2022-12-09 南京理工大学 针对fmcw毫米波雷达切片重构干扰的抑制方法及系统
CN116359854A (zh) * 2023-03-23 2023-06-30 哈尔滨工程大学 一种基于YOLOv5的防空警戒雷达复合干扰参数估计方法
CN117008064A (zh) * 2023-08-07 2023-11-07 青岛九维华盾科技研究院有限公司 一种基于改进型YOLOv5的雷达干扰自适应抑制方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5790475A (en) * 1996-10-28 1998-08-04 Multispec Corporation Process and apparatus for improved interference suppression in echo-location and imaging systems
UA89945U (ru) * 2013-10-14 2014-05-12 Запорізький Національний Технічний Університет Способ помехозащиты когерентно-импульсной рлс от комбинированных помех
US20210231789A1 (en) * 2018-06-28 2021-07-29 Richwave Technology Corp. Doppler signal processing device and method thereof for interference suppression
CN113176542A (zh) * 2021-03-24 2021-07-27 中国人民解放军空军预警学院 一种时域自适应抗主瓣自卫式干扰方法及系统
CN113702921A (zh) * 2021-08-17 2021-11-26 西安电子科技大学 基于时频分析的干扰抑制方法
CN115453469A (zh) * 2022-09-21 2022-12-09 南京理工大学 针对fmcw毫米波雷达切片重构干扰的抑制方法及系统
CN116359854A (zh) * 2023-03-23 2023-06-30 哈尔滨工程大学 一种基于YOLOv5的防空警戒雷达复合干扰参数估计方法
CN117008064A (zh) * 2023-08-07 2023-11-07 青岛九维华盾科技研究院有限公司 一种基于改进型YOLOv5的雷达干扰自适应抑制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
WANG XIAOGE 等: "Radar anti-false target jamming method based on phase modulation", SYSTEMS ENGINEERING AND ELECTRONICS, vol. 43, no. 9, 30 September 2021 (2021-09-30), pages 2476 - 2483 *
周必雷 等: "基于盲分离的空时联合处理抗复合干扰方法", 系统工程与电子技术, no. 11, 30 October 2018 (2018-10-30), pages 6 - 15 *
邓斌: "某型雷达性能参数测试仿真训练系统设计与实现", 2020中国系统仿真与虚拟现实技术高层论坛论文集, 31 December 2022 (2022-12-31), pages 98 - 102 *

Also Published As

Publication number Publication date
CN117784040B (zh) 2024-05-10

Similar Documents

Publication Publication Date Title
US5376939A (en) Dual-frequency, complementary-sequence pulse radar
CN109490851B (zh) 一种基于脉间伪随机码的抗干扰方法
CN106019243B (zh) 一种基于三次初相和frft的drfm干扰检测和对抗方法
CN112363121B (zh) 基于回波预处理和相参积累的c&i干扰抑制方法
CN113376601B (zh) 基于clean算法的捷变频雷达旁瓣抑制方法
CN110133632B (zh) 一种基于cwd时频分析的复合调制信号识别方法
CN114578296B (zh) 基于相位编码信号的间歇采样干扰抑制方法
Gulum et al. Parameter extraction of FMCW modulated radar signals using Wigner-Hough Transform
CN114966572A (zh) 基于脉内分段lfm波形的间歇采样转发干扰对抗方法
CN113176542B (zh) 一种时域自适应抗主瓣自卫式干扰方法及系统
CN109061626B (zh) 一种步进频相参处理检测低信杂比动目标的方法
CN113176541B (zh) 一种自适应抗频谱弥散干扰方法及系统
CN117784040B (zh) 一种时频域抗自卫式复合干扰方法
CN116400303A (zh) 一种脉内联合脉间时域编码抗isrj的波形设计方法
Li et al. A Radar Waveform Design of MCPC Method for Interrupted Sampling Repeater Jamming Suppression via Fractional Fourier Transform.
Li et al. A study on blanket noise jamming to LFM pulse compression radar
Zhang et al. Enhanced Signal Processing Method for Integrated Detection and Jamming System Considering the Complex Environment
CN112285656B (zh) 一种基于航迹特征的抗时敏目标欺骗干扰方法
Tian et al. An adaptive scheme for chopping and interleaving jamming suppression based on the mismatched filter
Rodriguez et al. A simple clutter suppression approach for OFDM-based passive radar exploiting reciprocal filter
Taboada et al. Intercept receiver signal processing techniques to detect low probability of intercept radar signals
CN118244210B (zh) 一种间歇采样转发干扰主动对抗方法、装置及系统
Fan et al. A New CCBFM/P4 RF Stealth Signal Waveform Design
Tian et al. An Improved Frequency Domain Interference Suppressor for DBT Processing of High Order BOC Signals
CN116774163A (zh) 一种基于复经验模态分解的抗密集假目标干扰方法

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