CN110927693B - 一种匹配滤波结合稀疏反演的脉冲压缩方法 - Google Patents
一种匹配滤波结合稀疏反演的脉冲压缩方法 Download PDFInfo
- Publication number
- CN110927693B CN110927693B CN201911338589.5A CN201911338589A CN110927693B CN 110927693 B CN110927693 B CN 110927693B CN 201911338589 A CN201911338589 A CN 201911338589A CN 110927693 B CN110927693 B CN 110927693B
- Authority
- CN
- China
- Prior art keywords
- matched filtering
- target
- sparse inversion
- pulse compression
- echo signal
- 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
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种匹配滤波结合稀疏反演的脉冲压缩方法,属雷达信号处理技术领域。通过本发明包括四个步骤,首先建立匹配滤波雷达回波信号模型,然后在所述雷达回波信号模型基础上建立稀疏反演框架,实现在中等信噪比情况下有效抑制匹配滤波的主瓣展宽和副瓣抬升,利用稀疏反演压制在低信噪比情况下由匹配滤波的模糊调制造成的目标主瓣展宽和副瓣抬高现象,大大提高目标分辨能力;利用稀疏反演直接估计目标后向散射系数,实现在低信噪比情况下成功压制噪声对目标回波的干扰,有效提高目标检测的准确度。解决了现有单独采用匹配滤波法进行脉冲压缩处理,因匹配滤波的模糊调制和雷达接收机热噪声造成的目标主瓣展宽、副瓣抬升的问题。
Description
技术领域
本发明涉及一种匹配滤波结合稀疏反演的脉冲压缩方法,属雷达信号处理技术领域。
背景技术
脉冲压缩波形广泛应用于雷达系统,目前雷达系统普遍采用匹配滤波法进行脉冲压缩处理,该法以最大信噪比为准则实现最佳线性滤波和最大化积累信号能量。然而,依据模糊原理可知,匹配滤波法实际上是对真实目标的一种函数调制,以具有一定宽度的主瓣和高频副瓣的形式体现。在实际工程应用中,受限于元器件的稳定性,发射波形通常存在不同程度的幅相扰动而与设计波形产生差异,导致脉冲压缩处理信号的主瓣畸变和副瓣抬高;此外,发射波形与雷达接收机的热噪声往往无法做到完全正交,会进一步造成副瓣抬升,使得脉冲压缩处理结果无法达到理想的主副瓣比,严重降低对目标的分辨能力;且在低信噪比情况下无法实现最优检测。
发明内容
本发明的目的在于,针对上述现有技术的不足,提供一种匹配滤波结合稀疏反演的脉冲压缩方法,通过首先建立匹配滤波后雷达回波信号模型,再利用稀疏反演压制在低信噪比情况下由匹配滤波造成的目标主瓣展宽和副瓣抬高现象,提高目标信噪比和分辨能力;利用稀疏反演框架直接估计目标后向散射系数,实现在低信噪比情况下压制噪声对目标回波的干扰,有效提高目标检测的准确度。解决单独采用匹配率波法进行脉冲压缩处理,因匹配滤波的模糊调制及雷达接收机热噪声造成的目标主瓣展宽和副瓣抬升的问题。
本发明是通过如下的技术方案来实现上述目的的:
一种匹配滤波结合稀疏反演的脉冲压缩方法,包括匹配滤波器、稀疏反演框架,其特征在于:所述匹配滤波结合稀疏反演的脉冲压缩方法是通过包括如下步骤实现的:
步骤一、通过信号卷积模型和匹配滤波原理建立匹配滤波后雷达回波信号模型
式中:为经过匹配滤波后的输出回波信号,为算子,为考虑了球面扩散及电磁波传播损耗的目标后向散射复包络,为发射波形的复共轭,为卷积算子,为发射波形;当且仅当算子为狄拉克函数时,匹配滤波器满足,否则匹配滤波器输出为算子对真实目标后向散射系数的模糊调制,一般情况下算子的影响不可忽视;
步骤二、在步骤一基础上建立稀疏反演框架
1)根据目标后向散射复包络的估计问题可以描述为一个典型的反问题,定义如下目标函数:
3)为便于问题分析与数值实现,公式(2)可变形如下:
步骤三、对公式(3)进行正则化处理,引入对目标的稀疏性约束如下:
本发明与现有技术相比的有益效果在于:
该匹配滤波结合稀疏反演的脉冲压缩方法,通过包括四个步骤,实现在中等信噪比情况下有效抑制匹配滤波的主瓣展宽与副瓣抬升问题,大大提高目标信噪比和分辨能力;及利用稀疏反演框架直接估计目标后向散射系数,在低信噪比情况下有效压制噪声对目标回波的干扰,切实提高目标检测的准确度。很好地解决了现有单独采用匹配滤波法进行脉冲压缩处理,因匹配滤波的模糊调制和雷达接收机热噪声造成的目标主瓣展宽、副瓣抬升的问题。
附图说明
图1为一种匹配滤波结合稀疏反演的脉冲压缩方法的匹配滤波雷达发射波形示意图;
图2为一种匹配滤波结合稀疏反演的脉冲压缩方法的匹配滤波后雷达目标回波示意图;
图3为一种匹配滤波结合稀疏反演的脉冲压缩方法的信噪比SNR=-5dB时的匹配滤波效果示意图;
图4为一种匹配滤波结合稀疏反演的脉冲压缩方法的稀疏反演框架在信噪比SNR=-5dB时的脉冲压缩效果示意图;
图5为一种匹配滤波结合稀疏反演的脉冲压缩方法的稀疏反演框架在信噪比SNR=-15dB时的脉冲压缩与匹配滤波处理效果对比示意图;
图6为一种匹配滤波结合稀疏反演的脉冲压缩方法的稀疏反演框架在信噪比SNR=-18dB时的脉冲压缩与匹配滤波处理效果对比示意图;
图7为一种匹配滤波结合稀疏反演的脉冲压缩方法的稀疏反演框架在信噪比SNR=-20dB时的脉冲压缩与匹配滤波处理效果对比示意图。
具体实施方式
下面结合附图对该匹配滤波结合稀疏反演的脉冲压缩方法的实施方式作进一步详细说明:
(参见图1~4),通过图1~图4验证本发明匹配滤波结合稀疏反演的脉冲压缩方法对目标的分辨能力。图1所示的发射波形为线性调频信号,其脉宽为100 us,带宽4 MHz,设置两个目标间隔37.5 m,分别在100 km和100.0375 km处,匹配滤波处理选取hamming窗。图2所示为接收回波,加入了SNR = -5 dB的随机噪声。图3所示为匹配滤波处理结果,两个目标在脉压处理后无法识别,这是因为加窗造成的主瓣展宽效应使得脉压处理的分辨率低于理论分辨率37.5 m(目标间隔);此外,由于噪声与发射波形往往具有一定的相关性,匹配滤波处理后的旁瓣也无法达到理想的水平。图4为稀疏反演框架下的脉冲压缩技术的处理效果,由于稀疏反演框架下的最优化算法实施过程中不存在对目标的“模糊”调制,因此所述两个目标被清晰的分辨出来,且副瓣被压制到较低水平。本发明的稀疏反演框架下的最优化算法在matlab2016中的计算效率控制在10ms左右,达到了某些场景下处理的实时性要求。
(参见图5~7),图5~图7验证本发明的稀疏反演框架下脉冲压缩对噪声的抑制作用。图5、图6、图7分别给出了SNR为-15dB、-18dB和-20dB情况下,稀疏反演框架下脉冲压缩与匹配滤波处理的效果对比显示结果,显示结果经过归一化处理。如图5,可以看出,SNR=-15dB时,匹配滤波副瓣出现-5dB左右的噪声,整体噪底水平较高,易造成对目标的漏检,而稀疏反演框架下脉冲压缩处理的副瓣可控制在-40dB以下,体现出对噪声的强抑制能力。如图6,可以看出,SNR=-18dB时,匹配滤波副瓣进一步抬高,而稀疏反演框架下脉冲压缩处理仅局部出现-22dB左右的噪声,依然可以保证稳定检测。如图7,可以看出,当信噪比下降至-20dB时,稀疏反演框架下脉冲压缩处理的最大噪声水平依然比匹配滤波低6.8dB左右,体现出对噪声的强稳健性。
本发明一种匹配滤波结合稀疏反演的脉冲压缩方法,根据信号卷积模型和匹配滤波原理,首先通过匹配滤波器建立匹配滤波后的雷达回波信号模型,然后以该雷达回波信号模型为基础建立稀疏反演框架,利用稀疏反演框架估计目标后向散射系数,压制在低信噪比情况下由匹配滤波造成的目标主瓣展宽和高副瓣问题,最终实现超低副瓣脉冲压缩。下面对本发明包括的实现步骤进行详细阐述。
一、通过匹配滤波器建立匹配滤波后雷达回波信号模型
单个脉冲的回波信号模型可表达为如下卷积过程:
匹配滤波器是以最大信噪比为准则的最佳线性滤波器,然而考察上述雷达回波信号模型公式(1)可知,匹配滤波器的输出本质上是真实后向散射复包络与算子的卷积,当且仅当算子为狄拉克函数时,满足。对于一般情形,发射的脉冲波形为线性调频信号,算子为辛克函数,工程实现中考虑加窗的影响,匹配滤波的输出实际上是对真实目标的一种“模糊”调制,以平滑主瓣和高频副瓣的形式即辛克函数的平滑形式体现。
考虑实际工程应用场景,一方面受限于元器件的稳定性,发射波形往往存在不同程度的幅相扰动,易引起匹配滤波的输出信号出现主瓣畸变和副瓣抬高;另一方面,发射波形与雷达接收机的热噪声往往无法做到完全正交,一定程度上也会进一步抬升副瓣,无法达到理想的主副瓣比,在较低信噪比情况下无法实现最优检测。
二、在步骤一基础上建立稀疏反演框架
本发明申请人考虑到匹配滤波处理在脉冲压缩中的应用劣势主要在于精度不高,近年来,反演与最优化理论在雷达信号处理领域:雷达波形与发射方向图设计、自适应阵列信号处理、现代谱估计等得到广泛应用,体现出传统算法无法达到的精度优势。而雷达脉冲压缩处理的目标后向散射复包络的估计问题可以描述为一个典型的反问题,由此定义如下目标函数:
1)根据目标后向散射复包络的估计问题可以描述为一个典型的反问题,定义如下目标函数:
由于噪声的存在,直接求解所述变形公式(3)描述的反问题,通常存在数值不稳定性,因此需要对该反问题进行正则化处理。常见的正则化算子为:对解的平滑或阻尼约束,其可有效提升反问题求解的数值稳定性,但针对所述目标复包络的估计问题其物理意义不明确,且会引入解的平滑效应,影响参数估计精度。本发明申请人考虑到相对于整个回波信号的采样长度,目标的个数总是较少这一特性,引入对目标的稀疏性约束。
三、对所述变形公式(3)进行正则化处理,引入对目标的稀疏性约束如下:
目标函数的常用最优化数值求解算法包括正交匹配追踪、同伦算法、迭代收缩阈值算法等,本发明申请人经过多种算法仿真实验,最终选取迭代收缩阈值算法求解目标函数的最优化数值,该迭代收缩阈值算法对噪声具有良好的适应性,同时,迭代收敛速度较快。
本发明与现有技术相比的有益效果体现在两点:第一,中等信噪比情况下可有效抑制匹配滤波的主瓣展宽和副瓣抬升,从而提高目标分辨能力;第二,利用稀疏反演直接估计目标后向散射系数,在低信噪比情况下能压制噪声对目标回波的干扰,提高目标检测的准确度。
以上所述只是本发明的较佳实施例而已,上述举例说明不对本发明的实质内容作任何形式上的限制,所属技术领域的普通技术人员在阅读了本说明书后依据本发明的技术实质对以上具体实施方式所作的任何简单修改或变形,以及可能利用上述揭示的技术内容加以变更或修饰为等同变化的等效实施例,均仍属于本发明技术方案的范围内,而不背离本发明的实质和范围。
Claims (1)
1.一种匹配滤波结合稀疏反演的脉冲压缩方法,包括匹配滤波器、稀疏反演框架,其特征在于:所述匹配滤波结合稀疏反演的脉冲压缩方法是通过包括如下步骤实现的:
步骤一、通过信号卷积模型和匹配滤波原理建立匹配滤波后雷达回波信号模型如下:
式中:为经过匹配滤波后的输出回波信号,为算子,为考虑了球面扩散及电磁波传播损耗的目标后向散射复包络,为发射波形的复共轭,为卷积算子,为发射波形;当且仅当算子为狄拉克函数时,匹配滤波器满足,否则匹配滤波器输出为算子对真实目标后向散射系数的模糊调制,一般情况下算子的影响不可忽视;
步骤二、在步骤一基础上建立稀疏反演框架
1)根据目标后向散射复包络的估计问题可以描述为一个典型的反问题,定义如下目标函数:
3)为便于问题分析与数值实现,公式(2)可变形如下:
步骤三、对公式(3)进行正则化处理,引入对目标的稀疏性约束如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911338589.5A CN110927693B (zh) | 2019-12-23 | 2019-12-23 | 一种匹配滤波结合稀疏反演的脉冲压缩方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911338589.5A CN110927693B (zh) | 2019-12-23 | 2019-12-23 | 一种匹配滤波结合稀疏反演的脉冲压缩方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110927693A CN110927693A (zh) | 2020-03-27 |
CN110927693B true CN110927693B (zh) | 2021-07-27 |
Family
ID=69860687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911338589.5A Active CN110927693B (zh) | 2019-12-23 | 2019-12-23 | 一种匹配滤波结合稀疏反演的脉冲压缩方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110927693B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113835066B (zh) * | 2021-09-15 | 2023-12-26 | 中国人民解放军陆军工程大学 | 一种抗转发式干扰方法及装置 |
CN114779335B (zh) * | 2022-03-31 | 2023-05-02 | 吉林大学 | 基于各向异性全变分约束的弹性波直接包络反演方法 |
CN115453484B (zh) * | 2022-08-29 | 2023-08-25 | 电子科技大学 | 一种自适应散射中心估计的距离扩展目标检测方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2296345C2 (ru) * | 2004-12-30 | 2007-03-27 | Открытое акционерное общество "Корпорация "Фазотрон-Научно-исследовательский институт радиостроения" | Способ разрешения целей по дальности радиолокационной станцией и импульсная радиолокационная станция со сжатием импульсов и восстановлением сигналов |
CN103630897A (zh) * | 2012-08-28 | 2014-03-12 | 中国科学院电子学研究所 | 一种多通道合成孔径雷达成像的方法 |
CN103235295B (zh) * | 2013-04-02 | 2014-10-01 | 西安电子科技大学 | 基于压缩卡尔曼滤波的小场景雷达目标距离像估计方法 |
CN103869316A (zh) * | 2014-03-27 | 2014-06-18 | 西安电子科技大学 | 基于稀疏表征的前视阵列sar超分辨成像方法 |
EP3137856A4 (en) * | 2014-04-28 | 2017-12-27 | Windbidco Pty Ltd | Method and system for detecting aircraft induced wake turbulence |
CN106908768A (zh) * | 2017-02-24 | 2017-06-30 | 四川九洲空管科技有限责任公司 | 一种超低旁瓣的脉冲压缩方法 |
CN110441772B (zh) * | 2019-07-30 | 2021-03-30 | 北京环境特性研究所 | 大轨道分布范围下星载稀疏重航过sar三维成像方法 |
-
2019
- 2019-12-23 CN CN201911338589.5A patent/CN110927693B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110927693A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110927693B (zh) | 一种匹配滤波结合稀疏反演的脉冲压缩方法 | |
CN110515052B (zh) | 一种基于时间反演的超宽带频域非等间隔采样目标检测方法 | |
CN108120976A (zh) | 一种基于多普勒通道特性的地杂波谱泄露抑制方法 | |
CN106646406A (zh) | 基于改进小波阈值去噪的外弹道测速雷达功率谱检测方法 | |
CN104793194B (zh) | 基于改进的自适应多脉冲压缩的距离‑多普勒估计方法 | |
CN110208796B (zh) | 基于奇异值逆滤波的扫描雷达超分辨成像方法 | |
CN111796242B (zh) | 一种基于功率特征值提取的改进块稀疏贝叶斯抗干扰方法 | |
CN106772254A (zh) | 改进的基于数字自适应干扰对消的收发隔离方法 | |
CN110412552A (zh) | 一种干扰样本自适应选择的副瓣对消方法 | |
CN114660562A (zh) | 自适应滤波方法及滤波器 | |
CN107678029B (zh) | 一种基于随机参考平均互相关信息的后向投影成像方法 | |
CN112505665B (zh) | 适用于部分均匀混响环境的空时自适应检测方法及系统 | |
CN111398966A (zh) | 一种主动声呐目标检测方法 | |
CN116027280B (zh) | 一种低峰值旁瓣频率编码雷达波形设计方法 | |
CN111025253B (zh) | 用于雷达脉冲压缩处理的窗函数波形联合设计方法 | |
CN107976663A (zh) | 一种基于子空间投影的外辐射源雷达对目标压缩检测方法 | |
CN108152796A (zh) | 一种基于灰色卡尔曼滤波的主瓣移动干扰消除方法 | |
CN113625233B (zh) | 基于强化学习的雷达抗干扰智能决策方法 | |
CN108983227B (zh) | 一种基于白化滤波的极化mimo雷达检测方法 | |
CN110673118A (zh) | 一种主动声纳单频脉冲串波形设计及检测算法 | |
CN116165610A (zh) | 一种改进奇异值分解的海杂波抑制算法 | |
CN116243251A (zh) | 一种基于波形分集与子空间投影的主瓣欺骗干扰抑制方法 | |
CN111025258B (zh) | 一种针对雷达波形分集的联合失配滤波器及其设计方法 | |
CN116047488A (zh) | 一种针对压制式干扰源的主动声纳干扰抑制方法 | |
CN115542304A (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 |