一种基于变分模态分解的间歇采样转发干扰的抑制方法技术

本发明专利技术涉及一种基于变分模态分解的间歇采样转发干扰的抑制方法，包括：根据发射信号得到目标回波信号；根据间歇采样直接转发干扰信号和目标回波信号得到雷达的接收信号；将雷达的接收信号与参考信号的共轭进行相乘处理，得到去斜后的信号；通过变分模态分解方法对去斜后的信号进行分解，经过多次迭代得到多个模态分量，将干扰与信号分离；基于分解的模态分量，使用阈值分割的方法将干扰剔除，再根据剔除干扰的模态分量重构出剔除干扰的去斜信号，然后实现脉冲压缩。本发明专利技术的方法直接将干扰与信号分离，对干扰的剔除比较干净，不会有太多脉压增益损失，栅瓣较低。栅瓣较低。栅瓣较低。

【技术实现步骤摘要】
一种基于变分模态分解的间歇采样转发干扰的抑制方法


[0001]本专利技术属于雷达
，涉及一种基于变分模态分解的间歇采样转发干扰的抑制方法。

技术介绍

[0002]数字射频存储器(DRFM，Digital Radio Frequency Memory)具有大容量的数字存储功能和高速率的输出转发功能，已成为目前产生相干干扰的主要实现方法和硬件基础。间歇采样转发干扰(Interrupted
‑
Sampling Repeater Jamming，ISRJ)利用数字射频存储器(DRFM)对雷达发射信号进行截获、延迟与转发，然后利用脉压雷达的匹配滤波特性，可以在信号处理时获得增益，产生高相干度的假目标串，传统的信号处理方法不能有效辨识与抑制间歇采样转发干扰，因此极大地影响了雷达对目标的探测以及对目标参数的估计。
[0003]本征模态函数(Intrinsic Mode Function，IMF)，可以通俗地将其理解为从原始信号中分解出的子信号，是调幅调频信号。VMD(变分模式分解，Variational Mode Decomp本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于变分模态分解的间歇采样转发干扰的抑制方法，其特征在于，所述抑制方法包括：步骤1、基于目标与雷达接收机的距离R和延迟时间，根据发射信号得到目标回波信号；步骤2、根据间歇采样直接转发干扰信号和所述目标回波信号得到雷达的接收信号；步骤3、将所述雷达的接收信号与参考信号的共轭进行相乘处理，得到去斜后的信号，其中，所述去斜后的信号包括目标回波去斜信号和干扰去斜信号；步骤4、通过变分模态分解方法对去斜后的信号进行分解，经过多次迭代得到多个模态分量，将干扰与信号分离；步骤5、基于分解的模态分量，使用阈值分割的方法将干扰剔除，再根据剔除干扰的模态分量重构出剔除干扰的去斜信号，然后实现脉冲压缩。2.根据权利要求1所述的基于变分模态分解的间歇采样转发干扰的抑制方法，其特征在于，所述发射信号表示为：其中，S(t)为发射信号，t为时间，T
p
为信号的脉宽，μ为信号的调频斜率；所述目标回波信号表示为：其中，S
r
(t)为目标回波信号，A
r
为目标回波幅度，τ为延迟时间，τ＝2R/c，c为光速。3.根据权利要求2所述的基于变分模态分解的间歇采样转发干扰的抑制方法，其特征在于，所述间歇采样直接转发干扰信号表示为：其中，S
j
(t)为间歇采样直接转发干扰信号，A
j
为干扰信号幅值，N为干扰的切片个数，0≤n≤N，T
j
为干扰的切片宽度；所述雷达的接收信号表示为：S
recv
(t)＝S
r
(t)+S
J
(t)其中，S
recv
(t)为雷达的接收信号。4.根据权利要求1所述的基于变分模态分解的间歇采样转发干扰的抑制方法，其特征在于，所述参考信号表示为：
其中，S
ref
(t)为参考信号，t为时间，T
ref
为参考信号脉宽，τ
ref
为参考信号延时，μ为信号的调频斜率；所述去斜后的信号表示为：S
if
(t)＝S
r_if
(t)+S
j_if
(t)其中，S
if
(t)为去斜后的信号，S
r_if
(t)为目标回波去斜信号，S
j_if
为干扰去斜信号；所述目标回波去斜信号表示为：其中，*为共轭，T
p
为信号的脉宽，A
r
为目标回波幅度，τ为延迟时间；所述干扰去斜信号表示为：其中，A
j
为干扰信号幅值，N为干扰的切片个数，0≤n≤N，T
j
为干扰的切片宽度。5.根据权利要求1所述的基于变分模态分解的间歇采样转发干扰的抑制方法，其特征在于，所述步骤4包括：步骤4.1、根据模态分量和所述去斜后的信号得到满足的约束条件，满足的约束条件为：其中，u
k
(t)为信号的第k个模态分量，ω
k
为第k个模态分量对应的中心频率，表示使得{
·
}最小时取得的{u
k
(t)}、{ω
k
}，表示二范数，s.t.表示使对应的式子得到满足，表示括号内的式子对时间t求偏导，S
if
(t)为去斜后的信号；步骤4.2、通过增广拉格朗日函数将模态分量和所述去斜后的信号满足的约束条件等效为无约束问题，等效的无约束问题表示为：其中，L表示增广拉格朗日函数，λ(...

【专利技术属性】
技术研发人员：吴耀君，方毅，全英汇，耿强，刘智星，邢孟道，
申请(专利权)人：北京华航无线电测量研究所，
类型：发明
国别省市：