CN113702921B - 基于时频分析的干扰抑制方法 - Google Patents

基于时频分析的干扰抑制方法 Download PDF

Info

Publication number
CN113702921B
CN113702921B CN202110944340.XA CN202110944340A CN113702921B CN 113702921 B CN113702921 B CN 113702921B CN 202110944340 A CN202110944340 A CN 202110944340A CN 113702921 B CN113702921 B CN 113702921B
Authority
CN
China
Prior art keywords
time
frequency
signal
interference
frequency domain
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.)
Active
Application number
CN202110944340.XA
Other languages
English (en)
Other versions
CN113702921A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202110944340.XA priority Critical patent/CN113702921B/zh
Publication of CN113702921A publication Critical patent/CN113702921A/zh
Application granted granted Critical
Publication of CN113702921B publication Critical patent/CN113702921B/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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

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

Abstract

本发明公开了一种基于时频分析的干扰抑制方法，主要解决在低干信比条件下，干扰类型未知或多种干扰组合产生的干扰抑制问题。方案包括：1)雷达接收机接收目标与多种干扰形成的混合回波信号；2)对混合回波信号进行解线调；3)对解线调之后的信号进行短时傅里叶变换；4)设置时频筛选门限值，获取时频域筛选矩阵；5)根据获取到的时频域筛选矩阵求解时频域滤波器；6)利用时频滤波器对变换后解线调信号进行时频域滤波；7)对滤波后信号进行恢复与去解线调处理，并对处理结果进行脉冲压缩得到最终的干扰抑制结果。本发明能够对抗组合干扰，并使干扰有效对抗率大幅度提升，可用于雷达频谱弥散与切片重复转发组合干扰等多种干扰的抑制。

Description

基于时频分析的干扰抑制方法
技术领域
本发明属于信号处理技术领域，进一步涉及干扰抑制技术，具体为一种基于时频分析的干扰抑制方法，可用于在低干信比背景下对未知干扰进行抑制。
背景技术
随着干扰装备性能的不断提升，使用方式愈加灵活多样，致使干扰场景日趋复杂，对雷达探测目标参数构成了更大的挑战。电子干扰的存在极大的制约了雷达作战效能的发挥，如何提高雷达适应复杂战场环境的能力已成为当前国内外该领域研究热点问题。
虽然国内外采用脉冲压缩、频率捷变等技术，使得雷达的抗干扰性能得到了显著提高，但在对抗基于数字射频存储DRFM(Digital Radio Frequency Memory)技术的有源干扰时，效果仍不明显。
传统的干扰抑制方法在对抗基于DRFM技术的有源干扰时，通常分为两类，一是利用DRFM调制带来的量化特性差异去分辨目标与干扰信号，二是在某个域或多个域联合对目标与干扰进行区分；随着干扰机器件的发展进步，量化位数逐渐增大，量化特性差异几乎可以忽略此类抗干扰技术的抗干扰性能，无法有效的对抗新型干扰。近些年来，在多域进行联合对抗干扰的方法，如利用短时傅里叶变换、小波变换、极化域等联合多种域进行干扰目标区分的大量研究成果涌现。
频谱弥散SMSP(Smeared Spectrum)干扰是一种基于DRFM技术的针对线性调频信号的有源干扰。通过将一个完整的信号切片形成多个片段，进而对切片信号的调频率进行调制，使每段信号能够拓展到雷达信号的整个带宽。频谱弥散SMSP干扰在经过脉冲压缩处理后可以获得信号增益，干扰机可以通过改变干扰信号能量与调频率的调制，对雷达回波信号可以同时实现压制与欺骗干扰。由于SMSP等新型干扰样式与雷达信号在多域耦合，对其进行窄带滤波与匹配滤波等方法不可避免的会损失目标信息，在耦合过强时甚至会失效，所以研究一种的对抗新型干扰的抑制技术是当前领域亟待解决的问题。
张亮等人在“联合时频重排和双正交傅里叶变换对抗频谱弥散干扰[J].信号处理,2020,36(4):495-501.DOI:10.16798/j.issn.1003-0530.2020.04.003.”中，提出了一种利用联合时频重排和双正交傅里叶变换的SMSP干扰抑制算法，对回波重排矩阵进行Radon变换，估计干扰子信号个数，再将回波中的SMS P干扰调制成一个LFM信号，对调制后回波进行BFT峰值滤波处理，该方法虽然通过减去重构的干扰信号实现了干扰抑制，但对于多种干扰同时存在的组合干扰情况并不适用。
发明内容
本发明目的在于针对上述现有技术的不足，提出一种基于时频分析的干扰抑制方法，用于解决在干扰类型未知时或多种干扰组合下，在低干信比条件下的干扰抑制问题。可适用于雷达频谱弥散SMSP干扰与切片重复转发组合干扰等多种干扰在低干信比背景下的抑制。
实现本发明方案的思路是：首先根据发射信号的先验信息，对干扰与目标混合回波信号进行解线调处理，通过对解线调处理后的信号进行短时傅里叶变换STFT得到其时频域信息，然后通过计算出合适的时频域能量阈值求解出解线调后的回波信号的时频域信息的时频域筛选矩阵，将时频域筛选矩阵中大于时频域能量阈值的时频能量点设为1；利用求解的时频域筛选矩阵计算出时频滤波器，利用时频滤波器对解线调后的回波信号进行时频滤波；最后对滤波后的信号恢复解线调，再对其进行脉冲压缩处理，最终得到干扰抑制后的信号。
本发明实现上述目的具体步骤如下：
(1)雷达接收机接收由频谱弥散SMSP干扰、切片重复转发干扰与真实目标回波构成的混合回波信号y(n)：
y(n)＝x(n)+j1(n)+j2(n)，
其中，x(n)表示真实目标回波信号，j1(n)表示频谱弥散干扰信号，j2(n)表示切片重复转发干扰信号；
(2)对混合回波信号y(n)进行解线调，得到解线调之后的信号r(n)：
Figure BDA0003216247080000031
其中，μ表示雷达发射线性调频信号的调频率，exp()表示e的指数函数，j表示复数，t表示雷达混合回波信号的时长，Ts表示采样周期；
(3)对解线调之后的信号r(n)进行短时傅里叶变换STFT，得到变换后解线调信号STFT|r(n)(k,n)：
Figure BDA0003216247080000032
其中，k∈{0,1,…,L-1}表示窗内信号离散频率，n∈{0,1,…,N-1}表示离散时间，N表示接收到回波信号的采样点数，L表示矩形窗的窗长；
Figure BDA0003216247080000033
表示解线调信号经过短时傅里叶变换得到的第k个离散频率、第n个离散时间对应的时频能量点；
(4)对变换后解线调信号STFT|r(n)(k,n)的时频域进行筛选，步骤如下：
(4.1)设置时频点筛选门限值为p；
(4.2)构建与变换后解线调信号STFT|r(n)(k,n)维度相同的矩阵Ψ：
Figure BDA0003216247080000041
其中，αab为矩阵Ψ中第a行、第b列的时频筛选点，a∈{1,2,…,L}，b∈{1,2,…,N}；
(4.3)根据时频点筛选门限值p对变换后解线调信号STFT|r(n)(k,n)的时频域进行筛选，即对其时频能量点进行如下判断：
若
Figure BDA0003216247080000042
则令矩阵Ψ中的参数αab值为0；
若
Figure BDA0003216247080000043
则令矩阵Ψ中的参数αab值为1；
(4.4)得到由0和1组成的时频筛选矩阵Ψ；
(5)构造时频滤波器：
(5.1)对时频域筛选矩阵Ψ中的频率维度求和，表示如下：
Figure BDA0003216247080000044
其中，L表示频率维度；
(5.2)根据下式得到频滤波器的中心频率
Figure BDA0003216247080000045
Figure BDA0003216247080000046
(5.3)构建与时频域筛选矩阵Ψ维度相同的矩阵Ψ'：
Figure BDA0003216247080000051
(5.4)令Ψ'中
Figure BDA0003216247080000052
其余元素均为0；得到由0和1组成的时频滤波器Ψ'；
(6)利用时频滤波器Ψ'对变换后解线调信号STFT|r(n)(k,n)进行时频域滤波，得到滤波后时频域混合回波信号y1(k,n)；
(7)对滤波后时频域混合回波信号y1(k,n)进行恢复与去解线调处理，得到干扰抑制后的信号
Figure BDA0003216247080000053
Figure BDA0003216247080000054
其中，ξn表示第n个离散时间点干扰抑制后的信号；
(8)对干扰抑制后的信号
Figure BDA0003216247080000055
进行脉冲压缩处理，得到最终的干扰抑制结果。
本发明与现有技术相比具有以下优点：
第一、通常情况下干扰信号的能量强度往往远大于目标信号的能量强度，但在低干信比的时候，若利用此特征对干扰进行剔除的话，存在将目标信号也剔除出去的可能性，因此本发明选择利用目标信号持续时间比干扰持续时间长的特性，在时频域对信号进行处理，从而可在低干信比时对干扰实现有效抑制；
第二、由于本发明充分考虑了低干信比背景下干扰与目标信号难以剔除的实际情况，面对切片干扰与频谱弥散干扰的组合干扰形式，本发明无需获取干扰方的先验信息，便可使干扰有效对抗率大幅度提升。
附图说明
图1为本发明方法的实现流程图；
图2为真实目标回波与干扰的时域仿真图；
图3为对真实目标回波与干扰分别进行短时傅里叶变换后得到的时频仿真图；
图4为对混合回波信号进行时频点筛选后的时频仿真图；
图5为受干扰前后的回波信号分别进行脉冲压缩的仿真结果对比图；
图6为使用本发明方法对混合回波信号进行干扰抑制与未进行干扰抑制的脉冲压缩仿真对比图。
具体实施方式
下面结合附图对本发明的实现步骤及效果作进一步详细描述。
实施例一：
参照图1，本发明提出的一种基于时频分析的干扰抑制方法，具体实现步骤如下：
步骤1，雷达接收机接收目标与干扰形成的混合回波信号，该回波信号由由频谱弥散SMSP干扰、切片重复转发干扰与真实目标回波构成；混合回波信号y(n)可表示如下：
y(n)＝x(n)+j1(n)+j2(n)，
其中，x(n)表示真实目标回波信号，j1(n)表示频谱弥散SMSP干扰信号，j2(n)表示切片重复转发干扰信号；
步骤2，对混合回波信号y(n)进行解线调，得到解线调之后的信号r(n)：
Figure BDA0003216247080000061
其中，μ表示雷达发射线性调频信号的调频率，该调频率μ对于雷达来说是已知的；exp()表示e的指数函数，j表示复数，t表示雷达混合回波信号的时长，Ts表示采样周期；
步骤3，对解线调之后的信号r(n)进行短时傅里叶变换STFT，按如下公式实现：
Figure BDA0003216247080000071
其中，n∈{0,1,…,N-1}表示离散时间，k∈{0,1,…,L-1}表示窗内信号离散频率，l∈{0,1,…,L-1}表示滑窗内信号的离散时间，N表示接收到回波信号的采样点数，L表示矩形窗的窗长。
对r(n)进行短时傅里叶变换后，得到如下：
Figure BDA0003216247080000072
其中，
Figure BDA0003216247080000073
表示r(n)在离散时间为n时的离散傅里叶变换；进一步得到变换后解线调信号STFT|r(n)(k,n)：
Figure BDA0003216247080000074
其中，
Figure BDA0003216247080000075
表示解线调信号经过短时傅里叶变换得到的第k个离散频率、第n个离散时间对应的时频能量点；
步骤4，对变换后解线调信号STFT|r(n)(k,n)的时频域进行筛选，步骤如下：
(4.1)设置时频点筛选门限值为p；
本发明实施例中，所述时频点筛选门限值p为：
p＝2*mean(mean(STFT|r(n)(k,n)))，
其中，mean(·)表示对·取平均；当然此处也可根据实际情况选用其它方式确定该门限；在低干信比的电磁背景下，干扰信号的时频能量相比真实目标回波信号的时频能量值不会大很多，对所有时频点求均值可以获得一个高于噪声且小于干扰和目标信号的时频能量均值。
(4.2)构建与变换后解线调信号STFT|r(n)(k,n)维度相同的矩阵Ψ：
Figure BDA0003216247080000081
其中，αab为矩阵Ψ中第a行、第b列的时频筛选点，a∈{1,2,…,L}，b∈{1,2,…,N}；
(4.3)根据时频点筛选门限值p对变换后解线调信号STFT|r(n)(k,n)的时频域进行筛选，即对其时频能量点进行如下判断：
若
Figure BDA0003216247080000082
则令矩阵Ψ中的参数αab值为1；
若
Figure BDA0003216247080000083
则令矩阵Ψ中的参数αab值为0；
(4.4)若变换后解线调信号STFT|r(n)(k,n)时频点的能量大于这一均值则令所对应时频筛选矩阵Ψ中的时频点能量为1，反之则为0，得到由0和1组成的时频筛选矩阵Ψ；
步骤5，构造时频滤波器：
(5.1)对时频域筛选矩阵Ψ中的频率维度求和，表示如下：
Figure BDA0003216247080000091
其中，L表示频率维度；
(5.2)根据下式得到频滤波器的中心频率
Figure BDA0003216247080000092
Figure BDA0003216247080000093
(5.3)构建与时频域筛选矩阵Ψ维度相同的矩阵Ψ'：
Figure BDA0003216247080000094
(5.4)令Ψ'中
Figure BDA0003216247080000095
其余元素均为0；得到由0和1组成的时频滤波器Ψ'；
本发明针对未知干扰类型进行抑制，由于在低干信比的电磁背景下，对所有时频点求均值获得的时频能量均值，有可能在目标信号能量均值附近。若直接对混合信号进行干扰剔除的话，会将大量的目标信号时频点剔除。可能会导致抑制效果不佳甚至抑制失败。所以本发明提出的技术方案根据目标持续时间比干扰时间长来进行判断，即根据目标信号在时域的时间长度比干扰信号在时域的时间长度长这一特性，在时频域通过求每一频率维度所在的时域的长度，便可以得到真实目标所在的频域中心维度；通过设置合理的滤波器从而对干扰进行抑制。
步骤6，利用时频滤波器Ψ'对变换后解线调信号STFT|r(n)(k,n)进行时频域滤波，得到滤波后时频域混合回波信号y1(k,n)；具体如下：
Figure BDA0003216247080000101
其中，εkn表示第k个离散频率、第n个离散时间对应的滤波后时频域混合回波信号时频点。
步骤7，对筛选后时频域混合回波信号y1(k,n)进行恢复与去解线调处理，得到干扰抑制后的信号
Figure BDA0003216247080000102
Figure BDA0003216247080000103
其中，ξn表示第n个离散时间点干扰抑制后的信号，其按如下公式得到：
Figure BDA0003216247080000104
步骤8，对干扰抑制后的信号
Figure BDA0003216247080000105
进行脉冲压缩处理，得到最终的干扰抑制结果。
实施例二：参照图1，本发明提出的一种基于时频分析的干扰抑制方法，其实现步骤同实施例一，另外对于步骤5作进一步描述，具体如下：
在低干信比背景下，干扰信号时频点的能量值比目标时频点的能量值并不会大很多，在进行选取阈值之后，可以取得一个大于噪声时频能量而小于干扰回波与目标回波时频能量的阈值，令变换后解线调信号STFT|r(n)(k,n)大于阈值所对应时频筛选矩阵Ψ中的时频点能量为1，由于目标回波的持续时间在实际中是大于干扰回波的持续时间的，干扰与目标回波在解线调之后为单频信号，可以根据时频筛选矩阵中持续时间最长的频率维度来确定目标回波所在的频域维度，通过设计时频滤波器滤除干扰回波信号，从而可以有效进行干扰抑制。
对时频域筛选矩阵Ψ中的频率维度求和，表示如下：
Figure BDA0003216247080000111
其中，L表示频率维度；
根据下式得到频滤波器的中心频率
Figure BDA0003216247080000112
Figure BDA0003216247080000113
由上述两步可以求解出解线调后目标回波所存在的频域维度信息，通过设置合适的时频滤波器宽度，便可以对除目标信号以外的信号进行剔除，从而抑制点干扰回波信号。
本发明的效果可通过以下仿真实验进一步说明：
1)仿真实验环境和内容
实验环境：MATLAB R2015b,Intel(R)Core(TM)i7-10700 CPU 2.90GHz，Windows10。
仿真实验内容：实验使用了matlab仿真雷达线性调频信号与两个不同调制方式的频谱弥散SMSP干扰混合回波信号。雷达发射带宽设置为40MHz,时宽为8μs。雷达接收机采样频率为80MHz。干扰机发射SMSP干扰子脉冲数设定为n1＝4,切片重复转发干扰转发数量n2＝7；干信比设置为10dB。信号的短时傅里叶变换采用的长度M＝31的Hamming窗。时频点筛选条件设为在每个时频点能量小于p时为0，反之为1。实验中对干扰混合信号的时域仿真结果如图2所示；对混合信号进行短时傅里叶变换而得到的时频仿真结果如图3所示；在没有进行干扰抑制与进行了干扰抑制方面，实验中对混合回波信号进行时频点筛选之后的时频点筛选矩阵时频仿真结果如图4所示。实验中关于干扰抑制前后脉冲压缩信号的对比进行了仿真，结果如图5、图6所示。
2)实验结果：
结果如图5、图6，可以看到本发明提出的干扰抑制方法在低干信比的背景下，有着良好的抗干扰性能。本发明利用目标信号时间比干扰信号时间长，干扰信号强度比目标信号强度大的两个特性，可以对于频谱弥散干扰、切片干扰、距离拖引等干扰可以实现有效的抑制。
实验结果表明，本发明在低干信比的背景下，能够有效的抑制掉多种干扰组成的混合干扰，如频谱弥散与切片重复转发干扰，经过脉冲压缩之后能够准确的测量目标的参数信息；上述仿真分析证明了本发明所提方法的正确性与有效性。
本发明未详细说明部分属于本领域技术人员公知常识。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，显然对于本领域的专业人员来说，在了解了本发明内容和原理后，都可能在不背离本发明原理、结构的情况下，进行形式和细节上的各种修正和改变，但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。

Claims (5)

1.一种基于时频分析的干扰抑制方法，其特征在于，包括如下步骤：
(1)雷达接收机接收由频谱弥散SMSP干扰、切片重复转发干扰与真实目标回波构成的混合回波信号y(n)：
y(n)＝x(n)+j1(n)+j2(n)，
其中，x(n)表示真实目标回波信号，j1(n)表示频谱弥散干扰信号，j2(n)表示切片重复转发干扰信号；
(2)对混合回波信号y(n)进行解线调，得到解线调之后的信号r(n)：
r(n)＝y(n)·ref(n)
Figure FDA0003216247070000011
其中，μ表示雷达发射线性调频信号的调频率，exp()表示e的指数函数，j表示复数，t表示雷达混合回波信号的时长，Ts表示采样周期；
(3)对解线调之后的信号r(n)进行短时傅里叶变换STFT，得到变换后解线调信号STFT|r(n)(k,n)：
Figure FDA0003216247070000012
其中，k∈{0,1,…,L-1}表示窗内信号离散频率，n∈{0,1,…,N-1}表示离散时间，N表示接收到回波信号的采样点数，L表示矩形窗的窗长；
Figure FDA0003216247070000013
表示解线调信号经过短时傅里叶变换得到的第k个离散频率、第n个离散时间对应的时频能量点；
(4)对变换后解线调信号STFT|r(n)(k,n)的时频域进行筛选，步骤如下：
(4.1)设置时频点筛选门限值为p；
(4.2)构建与变换后解线调信号STFT|r(n)(k,n)维度相同的矩阵Ψ：
Figure FDA0003216247070000021
其中，αab为矩阵Ψ中第a行、第b列的时频筛选点，a∈{1,2,…,L}，b∈{1,2,…,N}；
(4.3)根据时频点筛选门限值p对变换后解线调信号STFT|r(n)(k,n)的时频域进行筛选，即对其时频能量点进行如下判断：
若
Figure FDA0003216247070000025
则令矩阵Ψ中的参数αab值为0；
若
Figure FDA0003216247070000026
则令矩阵Ψ中的参数αab值为1；
(4.4)得到由0和1组成的时频筛选矩阵Ψ；
(5)构造时频滤波器：
(5.1)对时频域筛选矩阵Ψ中的频率维度求和，表示如下：
Figure FDA0003216247070000022
其中，L表示频率维度；
(5.2)根据下式得到频滤波器的中心频率
Figure FDA0003216247070000023
Figure FDA0003216247070000024
(5.3)构建与时频域筛选矩阵Ψ维度相同的矩阵Ψ'：
Figure FDA0003216247070000031
Figure FDA0003216247070000032
其余元素均为0；得到由0和1组成的时频滤波器Ψ'；
(6)利用时频滤波器Ψ'对变换后解线调信号STFT|r(n)(k,n)进行时频域滤波，得到滤波后时频域混合回波信号y1(k,n)；
(7)对滤波后时频域混合回波信号y1(k,n)进行恢复与去解线调处理，得到干扰抑制后的信号
Figure FDA0003216247070000035
Figure FDA0003216247070000033
其中，ξn表示第n个离散时间点干扰抑制后的信号；
(8)对干扰抑制后的信号
Figure FDA0003216247070000036
进行脉冲压缩处理，得到最终的干扰抑制结果。
2.根据权利要求1所述的方法，其特征在于：步骤(3)中所述对解线调之后的信号r(n)进行短时傅里叶变换STFT，按如下公式实现：
Figure FDA0003216247070000034
其中，l∈{0,1,…,L-1}表示滑窗内信号的离散时间。
3.根据权利要求1所述的方法，其特征在于：步骤(4.1)中所述时频点筛选门限值p，根据下式进行设置：
p＝2*mean(mean(STFT|r(n)(k,n)))，
其中，mean(·)表示对·取平均。
4.根据权利要求1所述的方法，其特征在于：步骤(6)中所述滤波后时频域混合回波信号y1(k,n)，具体如下：
Figure FDA0003216247070000041
其中，εkn表示第k个离散频率、第n个离散时间对应的滤波后时频域混合回波信号时频点。
5.根据权利要求1所述的方法，其特征在于：步骤(7)中第n个离散时间点干扰抑制后的信号，按如下公式得到：
Figure FDA0003216247070000042
CN202110944340.XA 2021-08-17 2021-08-17 基于时频分析的干扰抑制方法 Active CN113702921B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110944340.XA CN113702921B (zh) 2021-08-17 2021-08-17 基于时频分析的干扰抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110944340.XA CN113702921B (zh) 2021-08-17 2021-08-17 基于时频分析的干扰抑制方法

Publications (2)

Publication Number Publication Date
CN113702921A CN113702921A (zh) 2021-11-26
CN113702921B true CN113702921B (zh) 2023-06-30

Family

ID=78653078

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110944340.XA Active CN113702921B (zh) 2021-08-17 2021-08-17 基于时频分析的干扰抑制方法

Country Status (1)

Country Link
CN (1) CN113702921B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117784040B (zh) * 2024-02-23 2024-05-10 中国人民解放军空军预警学院 一种时频域抗自卫式复合干扰方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014048193A1 (zh) * 2012-09-28 2014-04-03 北京理工大学 一种用于舰艇编队情况下同型雷达同频干扰抑制方法
CN105891789A (zh) * 2016-03-30 2016-08-24 西安电子科技大学 联合时频分布和压缩感知对抗雷达频谱弥散干扰的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014048193A1 (zh) * 2012-09-28 2014-04-03 北京理工大学 一种用于舰艇编队情况下同型雷达同频干扰抑制方法
CN105891789A (zh) * 2016-03-30 2016-08-24 西安电子科技大学 联合时频分布和压缩感知对抗雷达频谱弥散干扰的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
联合时频重排和双正交傅里叶变换对抗频谱弥散干扰;张亮;王国宏;李思文;;信号处理(第04期);全文 *

Also Published As

Publication number Publication date
CN113702921A (zh) 2021-11-26

Similar Documents

Publication Publication Date Title
CN110519003B (zh) 一种基于信号特征差异的无人机上下行通信链路识别方法及装置
CN104657948B (zh) 一种用于海洋探测的激光水下成像图像去噪和增强方法
CN101631102A (zh) 一种跳频系统干扰样式识别技术
CN108880604B (zh) 一种扩频通信系统中多分量调频干扰抑制方法和装置
CN113702921B (zh) 基于时频分析的干扰抑制方法
CN106680790B (zh) 一种伪码调相连续波雷达干扰检测方法
CN106877903A (zh) 一种基于混合跳频的二维信息隐匿抗干扰通信方法
CN111323794A (zh) 一种基于周期重采样的周期调频干扰消除方法
CN112800863B (zh) 用于通信信号调制模式识别的时频特征提取方法及系统
CN106154236B (zh) 一种基于cemd抑制高频地波雷达射频干扰的方法
Zhao et al. The RFI suppression method based on STFT applied to SAR
Song et al. Estimation and mitigation of time-variant RFI in low-frequency ultra-wideband radar
CN113176542A (zh) 一种时域自适应抗主瓣自卫式干扰方法及系统
Abdoush et al. Time-frequency interference rejection based on the S-transform for GNSS applications
CN113702922B (zh) 基于多维协同的干扰抑制方法
Silva Lorraine et al. An efficient interference mitigation approach for NavIC receivers using improved variational mode decomposition and wavelet packet decomposition
Ruan et al. Wide band noise interference suppression for SAR with dechirping and eigensubspace filtering
Das et al. Interference cancellation schemes in UWB systems used in wireless personal area network based on wavelet based pulse spectral shaping and transmitted reference UWB using AWGN channel model
Ozdemir et al. Narrowband interference resilient receiver design for unknown UWB signal detection
Meng et al. An intelligent anti-jamming frequency hopping system
CN104777497A (zh) 单天线北斗卫星信号小波分解抗干扰算法
Aldirmaz et al. Broadband interference excision in spread spectrum communication systems based on short-time Fourier transformation
Khezzar et al. New Thresholding Technique in DCT Domain for Interference Mitigation in GNSS Receivers.
Jiang et al. Time-frequency-sparsity-based interference mitigation using order statistics for automotive radar
Niemelä et al. IEEE 802.15. 4a UWB receivers' performance in different body area network channels

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