CN106199539A - 基于白化滤波器的地杂波抑制方法 - Google Patents
基于白化滤波器的地杂波抑制方法 Download PDFInfo
- Publication number
- CN106199539A CN106199539A CN201610701166.5A CN201610701166A CN106199539A CN 106199539 A CN106199539 A CN 106199539A CN 201610701166 A CN201610701166 A CN 201610701166A CN 106199539 A CN106199539 A CN 106199539A
- Authority
- CN
- China
- Prior art keywords
- clutter
- prewhitening filter
- filter
- range gate
- prewhitening
- 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.)
- Pending
Links
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
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
本发明公开了一种基于白化滤波器的地杂波抑制方法,包括:将含有杂波和噪声的雷达回波数据进行脉冲压缩处理,并按照距离门重排数据;将脉冲压缩后的第K个距离门设为参考单元,并将慢时间维的雷达回波数据建立为自回归过程,通过基于三阶累积量的Burg算法估计参考单元自回归模型参数;通过估计出的参考单元自回归模型参数构造白化滤波器;利用构造出的白化滤波器对当前检测单元的数据进行白化滤波处理;对滤波后的时域信号进行动目标检测处理。本发明对存在慢速目标情况下的杂波进行有效抑制,将原本在零频附近聚集的杂波有效抑制,分散其功率使得杂波在频域表现为近似白噪声的形式,功率分散后的杂波对于慢速目标检测的影响明显减小。
Description
技术领域
本发明涉及雷达信号处理中的杂波抑制技术,特别是一种基于白化滤波器的地杂波抑制方法。
背景技术
自20世纪30年代雷达诞生以来,雷达技术获得了高速的发展。近年来,各种雷达处理技术的理论与应用研究成为一大热门和关键性课题。在雷达进行目标检测时,目标周围存在着各种背景干扰,也称为杂波,即雷达所需探测的运动目标以外的一切无源反射物的反射回波,雷达目标回波信号中往往夹杂着强度更大的环境杂波,如果弱目标淹没在强杂波中,特别是当强杂波使接收机过载时,这些目标的发现将会变得非常困难。而即使目标并没有处于杂波背景内,想要在成片的杂波中很快地分辨出运动目标回波也是一项很难完成的工作。为了应付日益复杂的雷达工作环境,提高检测能力,降低虚警率,国内外很早就开始了在各种杂波背景下检测目标信号的研究工作,杂波抑制技术是其中的关键技术之一。
对地面侦察雷达而言,其主要功能是在复杂地杂波环境中探测敌方运动中的人员、车辆、坦克及低空直升机等活动目标,测定其位置、运动速度和运动方向,判明其性质。慢速目标检测中的杂波抑制一直是地面侦查雷达目标检测中一个难点,慢速目标频谱如果位于杂波的频谱附近,则目标的功率往往会低于杂波的功率,难以进行有效的信号检测。而且由于慢速目标的频率处于杂波谱展宽附近,传统的MTI、AMTI等杂波抑制技术由于地面侦查雷达的多普勒分辨率不高,在抑制杂波的同时也会抑制目标信号本身,难以起到效果。近年来,提出了许多新的杂波抑制技术,如基于小波分解和中值滤波,基于杂波图技术的杂波抑制方法,基于EMD分解和Hilbert变换的滤波方法等,但是对存在慢速目标情况下的杂波抑制效果并不是很理想。
发明专利201210211297.7公开了一种基于支持向量回归的水下混响抑制方法,主要针对水下混响,将水下混响建立为自回归过程,并对其进行白化滤波,实现混响干扰的消除。发明专利201110409933.2公开了一种自适应杂波抑制动目标信号处理技术及实现方法,主要针对地杂波干扰,将脉压后的信号进行一次对消处理,然后送入窄带FIR滤波器组,各滤波器组输出分别进行距离单元恒虚警处理,经过选大输出得到检测结果,但是并未专门对目标速度情况做出讨论。国内外有关地面侦查雷达存在慢速目标情况下的杂波抑制方法的专利尚未查到。
发明内容
本发明的目的在于提供一种基于白化滤波器的地杂波抑制方法。
实现本发明的技术方案为:一种基于白化滤波器的地杂波抑制方法,包括以下步骤:
步骤1,将含有杂波和噪声的雷达回波数据进行脉冲压缩处理,并按照距离门重排数据;
步骤2,将脉冲压缩后的第K个距离门设为参考单元,并将慢时间维的雷达回波数据建立为自回归过程,通过基于三阶累积量的Burg算法估计参考单元自回归模型参数;K为与当前处理单元相隔两个距离门的距离门数;
步骤3,通过参考单元自回归模型参数构造白化滤波器;
步骤4,利用构造出的白化滤波器对当前检测单元的数据进行白化滤波处理;
步骤5,对滤波后的时域信号进行动目标检测处理,完成地杂波抑制。
本发明与现有的杂波抑制技术相比,其显著优点为:
(1)本发明能够对存在慢速目标情况下的杂波进行抑制,将原本在零频附近聚集的杂波有效抑制,分散其功率使杂波在频域表现为近似白噪声的形式,功率分散后的杂波对于慢速目标检测的影响明显减小;
(2)现有的MTI、AMTI等技术在抑制杂波的同时会对目标信号也造成较大的损失,本发明对于目标信号造成的滤波损失较小;
(3)本发明中采用的Burg算法是一种参数化分析方法,利用前、后向预测误差平均功率最小准则来估计反射系数,再利用Levinson递推由低阶到高阶来求自回归模型的参数,避免了计算自相关矩阵和求解Yule-Walker方程,并且参数估计的精确度相比其他估计方法如自协方差法等在精度和运算量方面有一定的优势;
(4)针对地面侦查雷达分辨率不高的现实情况,本发明将杂波构建为瑞利杂波,并基于三阶累积量估计白化滤波器系数,相对于传统的白化滤波器滤波效果有所提升;
(5)在现有的时频结合杂波抑制技术中,其滤波过程往往需要较多人为的判断,本发明为自适应型滤波器,不会引入人为干扰。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是本发明基于白化滤波器的地杂波抑制方法流程图。
图2是常规MTD目标距离门多普勒切片图。
图3是本发明杂波抑制后的MTD多普勒切片图。
图4是本发明杂波抑制后MTD的目标检测效果距离-多普勒三围图。
具体实施方式
本发明针对战场侦察雷达目标检测中慢速目标多普勒频谱位于或靠近杂波频谱区域难以进行目标检测的情况提出了一种有效抑制地杂波的方法。
结合图1,本发明基于白化滤波器的地杂波抑制方法包括以下步骤:
步骤1,将含有杂波和噪声的雷达回波数据进行脉冲压缩处理,并按照距离门重排数据;
步骤2,将脉冲压缩后的第K个距离门设为参考单元,并将慢时间维的雷达回波数据建立为自回归过程,通过基于三阶累积量的Burg算法估计参考单元自回归模型参数;K为与当前检测单元相隔两个距离门的距离门数,即当前检测单元为第K+2个距离门;
步骤3,通过估计出的参考单元自回归模型参数构造白化滤波器;
步骤4,利用构造出的白化滤波器对当前检测单元的数据进行白化滤波处理;
步骤5,对滤波后的时域信号进行动目标检测处理,完成地杂波抑制。
进一步的,步骤2中通过基于三阶累积量的Burg算法及最速下降法估计参考距离单元自回归模型参数的具体过程为:
步骤2-1,初始化Burg算法所需的迭代参数:
其中Km为反射系数,y(n)为当前时刻的回波信号,n表示当前时刻,为第m阶前向误差功率,为第m阶后向误差功率;
步骤2-2,计算每次Burg算法迭代过程对应的基于三阶累积量的代价函数:
其中,p为白化滤波器的阶数,本具体实施方式中根据实际运算量和效果折中取阶数为5阶;
步骤2-3,通过最速下降法进行递推,定义迭代过程中间参数:
α为收敛因子,0<α<<1,i的取值为[-p,p];
步骤2-4,求解下一次迭代所需的最优反射系数:
其中μ为迭代步长,取值为[0,1],为当前代价函数关于当前反射系数Km的梯度;
步骤2-5,利用Burg算法计算下一次迭代前向误差功率、后向误差功率,利用格型滤波器结构进行递推:
步骤2-6,重复以上步骤2-2到步骤2-5过程完成整个迭代过程。
进一步的,步骤3中通过估计出的参考单元自回归模型参数构造白化滤波器的具体过程为:
步骤3-1,利用得到的反射系数Km并经过Levinson递推准则得出基于三阶累积量自回归模型的白化滤波器系数ap;
步骤3-2,利用得到的白化滤波器系数ap,构造白化滤波器,白化滤波器的频率特性为:
其中,为ap的共轭转置,w为频率。
下面结合具体实施例对本发明作进一步说明。
实施例
本实施例中仿真参数设置为:载频为10GHz,载波为线性调频信号,带宽10MHz,采样频率20MHz,脉冲时宽80μs,占空比为1:20,慢时间维积累次数为512次,目标沿径向做匀速直线运动,速度为0.4m/s,噪声为白噪声,信噪比为10dB,杂波为服从瑞利分布的地杂波,信杂比为0dB。
从图2中可以看出,在常规动目标检测处理的情况下,零频附近的杂波能量集中,在幅度上表现得跟目标能量近似,难以检测出目标,极其容易出现虚警。
从图3中可以看出,经过白化滤波预处理后再进行动目标检测,零频附近的杂波基本被抑制住,整个检测背景变成了类似“白噪声”的情况,目标清晰凸显出来,可以进行目标检测。
从图4中可以看出,经过白化滤波处理后,原本有大量杂波干扰的检测背景变成了类似白噪声的背景,原本在零频附近聚集的杂波不再对慢速目标造成强大的干扰,目标检测得以进行,说明了这种白化预处理是有效的。
Claims (3)
1.一种基于白化滤波器的地杂波抑制方法,其特征在于,包括以下步骤:
步骤1,将含有杂波和噪声的雷达回波数据进行脉冲压缩处理,并按照距离门重排数据;
步骤2,将脉冲压缩后的第K个距离门设为参考单元,并将慢时间维的雷达回波数据建立为自回归过程,通过基于三阶累积量的Burg算法估计参考单元自回归模型参数;K为与当前处理单元相隔两个距离门的距离门数;
步骤3,通过参考单元自回归模型参数构造白化滤波器;
步骤4,利用构造出的白化滤波器对当前检测单元的数据进行白化滤波处理;
步骤5,对滤波后的时域信号进行动目标检测处理,完成地杂波抑制。
2.根据权利要求1所述的基于白化滤波器的地杂波抑制方法,其特征在于,步骤2中通过基于三阶累积量的Burg算法估计参考距离单元自回归模型参数的具体过程为:
步骤2-1,初始化Burg算法所需的迭代参数:
其中Km为反射系数,y(n)为当前时刻的回波信号,n表示当前时刻,为第m阶前向误差功率,为第m阶后向误差功率;
步骤2-2,计算每次Burg算法迭代过程对应的基于三阶累积量的代价函数:
其中,p为白化滤波器的阶数;
步骤2-3,通过最速下降法进行递推,定义迭代过程中间参数:
α为收敛因子,0<α<<1,i的取值为[-p,p];
步骤2-4,求解下一次迭代所需的最优反射系数:
其中μ为迭代步长,取值为[0,1],为当前代价函数关于当前反射系数Km的梯度;
步骤2-5,利用Burg算法计算下一次迭代前向误差功率、后向误差功率,利用格型滤波器结构进行递推:
步骤2-6,重复以上步骤2-2到步骤2-5过程完成整个迭代过程。
3.根据权利要求1所述的基于白化滤波器的地杂波抑制方法,其特征在于,步骤3具体为:
步骤3-1,利用得到的反射系数Km并经过Levinson递推准则得出基于三阶累积量自回归模型的白化滤波器系数ap;
步骤3-2,利用得到的白化滤波器系数ap,构造白化滤波器,白化滤波器的频率特性为:
其中,为ap的共轭转置,w为频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610701166.5A CN106199539A (zh) | 2016-08-22 | 2016-08-22 | 基于白化滤波器的地杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610701166.5A CN106199539A (zh) | 2016-08-22 | 2016-08-22 | 基于白化滤波器的地杂波抑制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106199539A true CN106199539A (zh) | 2016-12-07 |
Family
ID=57524072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610701166.5A Pending CN106199539A (zh) | 2016-08-22 | 2016-08-22 | 基于白化滤波器的地杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106199539A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061626A (zh) * | 2018-07-19 | 2018-12-21 | 武汉滨湖电子有限责任公司 | 一种步进频相参处理检测低信杂比动目标的方法 |
CN110146851A (zh) * | 2019-05-17 | 2019-08-20 | 西安电子科技大学 | 一种基于数字统计特性分析改善雷达回波信噪比的方法 |
CN110927750A (zh) * | 2019-11-22 | 2020-03-27 | 中科院计算技术研究所南京移动通信与计算创新研究院 | 一种基于格型滤波Burg谱估计算法的低轨卫星多普勒频偏捕获方法 |
CN113759354A (zh) * | 2020-06-02 | 2021-12-07 | 中国科学院声学研究所 | 一种适用于侧扫声呐的自适应底混响抑制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004257761A (ja) * | 2003-02-24 | 2004-09-16 | Toshiba Corp | レーダ信号処理装置およびレーダ信号処理方法 |
CN102565763A (zh) * | 2011-12-12 | 2012-07-11 | 中国船舶重工集团公司第七二四研究所 | 一种自适应杂波抑制动目标信号处理技术及实现方法 |
CN103513240A (zh) * | 2012-06-20 | 2014-01-15 | 中国科学院声学研究所 | 一种基于支持向量回归的水下混响抑制方法 |
-
2016
- 2016-08-22 CN CN201610701166.5A patent/CN106199539A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004257761A (ja) * | 2003-02-24 | 2004-09-16 | Toshiba Corp | レーダ信号処理装置およびレーダ信号処理方法 |
CN102565763A (zh) * | 2011-12-12 | 2012-07-11 | 中国船舶重工集团公司第七二四研究所 | 一种自适应杂波抑制动目标信号处理技术及实现方法 |
CN103513240A (zh) * | 2012-06-20 | 2014-01-15 | 中国科学院声学研究所 | 一种基于支持向量回归的水下混响抑制方法 |
Non-Patent Citations (4)
Title |
---|
何平编著: "《功率谱估计基础》", 31 January 2016, 北京:气象出版社 * |
宁贝佳: "雷达杂波抑制技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
张媚等: "海面目标的自适应预白化检测前跟踪算法", 《计算机工程》 * |
赵知劲: "一种新的基于高阶统计量的自适应滤波算法", 《电子学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061626A (zh) * | 2018-07-19 | 2018-12-21 | 武汉滨湖电子有限责任公司 | 一种步进频相参处理检测低信杂比动目标的方法 |
CN109061626B (zh) * | 2018-07-19 | 2022-06-28 | 武汉滨湖电子有限责任公司 | 一种步进频相参处理检测低信杂比动目标的方法 |
CN110146851A (zh) * | 2019-05-17 | 2019-08-20 | 西安电子科技大学 | 一种基于数字统计特性分析改善雷达回波信噪比的方法 |
CN110146851B (zh) * | 2019-05-17 | 2022-12-23 | 西安电子科技大学 | 一种基于数字统计特性分析改善雷达回波信噪比的方法 |
CN110927750A (zh) * | 2019-11-22 | 2020-03-27 | 中科院计算技术研究所南京移动通信与计算创新研究院 | 一种基于格型滤波Burg谱估计算法的低轨卫星多普勒频偏捕获方法 |
CN113759354A (zh) * | 2020-06-02 | 2021-12-07 | 中国科学院声学研究所 | 一种适用于侧扫声呐的自适应底混响抑制方法 |
CN113759354B (zh) * | 2020-06-02 | 2024-02-09 | 中国科学院声学研究所 | 一种适用于侧扫声呐的自适应底混响抑制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guan et al. | Adaptive fractional Fourier transform-based detection algorithm for moving target in heavy sea clutter | |
Li et al. | Fast coherent integration for maneuvering target with high-order range migration via TRT-SKT-LVD | |
CN106468770B (zh) | K分布杂波加噪声下的近最优雷达目标检测方法 | |
CN102879766B (zh) | 一种高频地波雷达弱目标检测跟踪方法和装置 | |
CN106127110B (zh) | 一种基于uwb雷达和最优svm的人体细粒度运动识别方法 | |
CN106199539A (zh) | 基于白化滤波器的地杂波抑制方法 | |
Shui et al. | Range-spread target detection using consecutive HRRPs | |
CN106569193B (zh) | 基于前-后向收益参考粒子滤波的海面小目标检测方法 | |
WO2006080011A2 (en) | Using pulsed-wave ultrasonography for determining an aliasing-free radial velocity spectrum of matter moving in a region | |
CN107167802A (zh) | 一种基于超宽带雷达的呼吸信号检测算法 | |
CN106646406A (zh) | 基于改进小波阈值去噪的外弹道测速雷达功率谱检测方法 | |
CN105116388A (zh) | 基于鲁棒主成分分析的天波超视距雷达瞬态干扰抑制方法 | |
CN106772302A (zh) | 一种复合高斯背景下的知识辅助stap检测方法 | |
CN112255607B (zh) | 一种海杂波的抑制方法 | |
CN104614769A (zh) | 一种压制地震面波的聚束滤波方法 | |
CN103760540B (zh) | 基于重构信号和1-范数的动目标检测与参数估计方法 | |
CN106353743B (zh) | 匹配于等效形状参数的近最优雷达目标检测方法 | |
CN103169449B (zh) | 呼吸信号检测方法和装置 | |
CN106597445A (zh) | 一种基于自适应Chirp分解的SAR动目标检测方法 | |
CN103885044B (zh) | 一种基于clean算法的窄带雷达回波杂噪抑制方法 | |
CN116224277A (zh) | 一种基于时频谱能量分布与奇异值分解的海杂波抑制方法 | |
CN116719000A (zh) | 双通道并行检测机载相控阵雷达地面运动目标提取方法 | |
Juan et al. | A new Wavelet Prediction method for GPR clutter elimination Based on LSTM network | |
CN101849841A (zh) | 基于几何滤波器抑制超声彩色血流成像中杂波的方法 | |
Fossum | Exploration of micro-Doppler signatures associated with humans and dogs using UWB radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161207 |
|
RJ01 | Rejection of invention patent application after publication |