CN112904327A - 基于调频模糊函数的复合微动目标参数估计方法 - Google Patents
基于调频模糊函数的复合微动目标参数估计方法 Download PDFInfo
- Publication number
- CN112904327A CN112904327A CN202110071362.XA CN202110071362A CN112904327A CN 112904327 A CN112904327 A CN 112904327A CN 202110071362 A CN202110071362 A CN 202110071362A CN 112904327 A CN112904327 A CN 112904327A
- Authority
- CN
- China
- Prior art keywords
- frequency modulation
- fuzzy function
- motion
- target
- estimating
- 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
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达信号处理领域,涉及一种基于调频模糊函数的复合微动目标参数估计方法,包括以下步骤:首先,计算复合微动目标雷达信号的调频模糊函数;之后计算调频模糊函数的l4‑l2范数;在获得调频模糊函数的l4‑l2范数的基础上,估计调频率并生成调频模糊函数统计量;进一步估计目标的微动周期;在获得微动周期的基础上,估计多普勒中心;最后,估计目标的高阶平动参数。本发明的有益效果如下:第一,本发明提供的复合微动目标参数估计方法能够实现机动和微动复合运动条件下的微动周期估计,并且不需要进行平动补偿;第二,本发明提供的复合微动目标参数估计方法能够在估计微动周期的同时获得平动参数的估计。
Description
技术领域
本发明属于雷达信号处理领域,涉及雷达目标参数估计方法中的一种基于调频模糊函数的复合微动目标参数估计方法。
背景技术
雷达目标探测与识别是目前传感器信息处理领域中的一个重要研究方向,通过对雷达回波信号的处理可获得目标的运动轨迹、速度、形状、尺寸、电磁参数等信息。雷达观测目标通常为非合作目标,且随着隐身、伪装等目标特征控制技术和虚假目标成像等目标干扰技术的发展,使基于目标形状、尺寸等传统雷达特征的识别受到挑战,对目标特有运动信息的进一步挖掘,是应对该挑战的一个可行途径。
微动(Micro-Motion)是指目标或目标的组成部分除了主体平动之外的振动、转动等小幅往复性运动。目标微动来源于特定的受力作用下的目标特殊结构,微动状态具有独特性,反映了目标的几何结构信息与运动特性等精细特征,可作为目标识别的重要依据。目前,微动特征在雷达目标探测与识别的作用愈加凸显,被广泛用于空间、空中、地面、海面目标探测中。对于地/海面目标,微动特征可用于辨别人的行进状态(徒手行进、跑步等)、从防空阵地中检测带转动天线的雷达、分辨海上颠簸舰船等;对于空中目标,微动特征可用于旋翼飞机低空悬停检测、固定翼飞机和直升机的鉴别;对于空间目标,由于弹道导弹弹头和轻、重诱饵的质量分布具有明显差异,因此微动特征明显不同,可用作从诱饵群中选择弹头目标的重要依据。
目标运动产生雷达信号多普勒调制。对于复合微动目标,即同时包含平动和微动两种运动的目标,目标多普勒调制特性表现更加复杂,信号的微多普勒周期性遭到破坏,基于信号微多普勒周期性的估计方法失效。针对复合微动目标参数估计,对比文件“WenpengZhang,Yaowen Fu,Jiapeng Yin,Micro-Doppler Period Estimation based onConcentration Statistics of Ambiguity Function[J],IEEE Transactions onAerospace and Electronic Systems Vol.56,No.3 JUNE 2020”提出基于模糊函数聚集度的微动周期估计方法。该方法具有无需进行平动补偿,能够在较低信噪比条件下准确估计微动周期等优势。然而,该方法仅能够处理包含二阶多项式平动的复合微动目标,且无法实现其他目标运动参数的估计。
发明内容
本发明的目的是,针对同时包含微动和机动的复合微动目标参数估计问题,提出一种基于调频模糊函数的复合微动目标参数估计方法,首先对雷达信号处理中广泛使用的模糊函数进行改进,提出了调频模糊函数,相比模糊函数,调频模糊函数增加了调频率相乘项,能够处理更复杂的运动,同时,通过计算复合微动目标回波的调频模糊函数统计量,进一步估计微动周期,并实现目标平动参数的估计。
为达到上述目的,本发明是通过以下技术方案实现的:
一种基于调频模糊函数的复合微动目标参数估计方法,包括下列步骤:
S1计算复合微动目标雷达信号的调频模糊函数:
记复合微动目标的雷达信号为s(n),信号的采样时间间隔为ΔT,其中0≤n<N表示信号的离散时间,N为信号的长度,NΔT即为信号的持续时间。复合微动目标的雷达信号s(n)可建模为:
rT(n)=R0+a1(nΔT)+a2(nΔT)2+a3(nΔT)3
R0为目标到雷达的起始距离,a1、a2、a3分别为平动的一阶参数、二阶参数、三阶参数。通常将三阶多项式平动称为机动。该信号的调频模糊函数计算具体步骤如下:
其中,*为共轭运算符,m为调频瞬时自相关函数的延时,m,n满足条件:0≤m≤N-n+1,β为调频率,β的取值为其中Δβ为调频率采样间隔,::表示从第一个数开始以间隔Δβ进行采样直到最后一个数形成的数组。为了平衡精度和计算量,此处设置对于较大的延时m,满足条件0≤m≤N-n+1的n较少,将影响调频瞬时自相关函数的可靠性,因此一般设置0≤m<[rm×N],rm可设置为0.8,其中[]为取整运算。
其中Nf为离散傅里叶变换的点数,k=0,1,…,Nf-1表示离散频率;
S3估计调频率并生成调频模糊函数统计量,具体步骤如下:
S4.1计算调频模糊函数统计量Ts(m)的截尾平均值η:
S4.2寻找调频模糊函数统计量Ts(m)的峰值位置,即满足下述条件的延时m:
其中l表示延时偏移量,L为延时偏移量的范围,可取为5。用集合Θm={mi|1≤i≤P}记录所有的峰值位置,其中i为峰值位置的索引,P为峰值位置的数量;
S6估计目标的高阶平动参数:
与现有技术相比,本发明的有益效果如下:
第一,本发明提供的复合微动目标参数估计方法能够实现机动和微动复合运动条件下的微动周期估计,并且不需要进行平动补偿;
第二,本发明提供的复合微动目标参数估计方法能够在估计微动周期的同时获得平动参数的估计。
附图说明
图1复合微动目标参数估计流程图;
图2复合微动目标雷达信号的处理结果:(a)复合微动目标雷达信号时频图,(b)调频模糊函数切片,(c)模糊函数统计量,(d)调频模糊函数统计量。
具体实施方式
下面结合附图对本发明做进一步说明:
图1为本发明提出的复合微动目标参数估计流程图。首先,计算复合微动目标雷达信号的调频模糊函数;之后计算调频模糊函数的l4-l2范数;在获得调频模糊函数的l4-l2范数的基础上,估计调频率并生成调频模糊函数统计量;进一步估计目标的微动周期;在获得微动周期的基础上,估计多普勒中心;最后,估计目标的高阶平动参数。
图2为复合微动目标雷达信号的处理结果。在该实验中,目标的微动周期为2s,目标的平动为rT(t)=-t+0.4t2+0.02t3,通常也将三阶多项式的平动称为目标机动。图2(a)为复合微动目标雷达信号时频图,坐标横轴为时间,坐标纵轴为多普勒频率,从图中可看出,目标的频率包括两部分,一部分为正弦形式的微动引起的多普勒频率变化,另外一部分则为目标机动引起的趋势项。图2(b)为调频模糊函数的切片,即的幅度图像。图2(c)为根据对比文件“Wenpeng Zhang,Yaowen Fu,Jiapeng Yin,Micro-Doppler PeriodEstimation based on Concentration Statistics of Ambiguity Function[J],IEEETransactions on Aerospace and Electronic Systems Vol.56,No.3 JUNE 2020”计算得到的模糊函数统计量,在τ=2s附近,该统计量未出现明显的峰值,表明其无法估计微动周期。作为对比,图2(d)中的调频模糊函数统计量在在τ=2s处出现较为明显的峰值,据此能够估计得到微动周期为与真实的微动周期一致,表明本发明提出的方法能够处理机动条件下的微动目标信号,获取其周期估计。进一步,使用本发明能够得到高阶平动参数的估计,分别为和这两个参数和真实值0.4m/s2、0.02m/s3接近,表明本发明提出的方法能够同时获得复合微动目标的高阶平动参数估计。
Claims (7)
1.一种基于调频模糊函数的复合微动目标参数估计方法,其特征在于,该方法包括下列步骤:
S1计算复合微动目标雷达信号的调频模糊函数:
记复合微动目标的雷达信号为s(n),信号的采样时间间隔为ΔT,其中0≤n<N表示信号的离散时间,N为信号的长度,NΔT即为信号的持续时间;复合微动目标的雷达信号s(n)可建模为:
rT(n)=R0+a1(nΔT)+a2(nΔT)2+a3(nΔT)3
R0为目标到雷达的起始距离,a1、a2、a3分别为平动的一阶参数、二阶参数、三阶参数;该信号的调频模糊函数计算具体步骤如下:
其中,*为共轭运算符,m为调频瞬时自相关函数的延时,m,n满足条件:0≤m≤N-n+1,β为调频率;
其中Nf为离散傅里叶变换的点数,k=0,1,…,Nf-1表示离散频率;
S3估计调频率并生成调频模糊函数统计量,具体步骤如下:
S4.1计算调频模糊函数统计量Ts(m)的截尾平均值η:
S4.2寻找调频模糊函数统计量Ts(m)的峰值位置,即满足下述条件的延时m:
其中l表示延时偏移量,L为延时偏移量的范围;用集合Θm={mi|1≤i≤P}记录所有的峰值位置,其中i为峰值位置的索引,P为峰值位置的数量;
S6估计目标的高阶平动参数:
4.一种根据权利要求1所述基于调频模糊函数的复合微动目标参数估计方法,其特征在于:S1.1中,设置0≤m<[rm×N],其中[]为取整运算。
5.一种根据权利要求4所述基于调频模糊函数的复合微动目标参数估计方法,其特征在于:S1.1中,rm设置为0.8。
6.一种根据权利要求1所述基于调频模糊函数的复合微动目标参数估计方法,其特征在于:S2中,设置lp范数中的p=4。
7.一种根据权利要求1所述基于调频模糊函数的复合微动目标参数估计方法,其特征在于:S4.2中,延时偏移量的范围L取为5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110071362.XA CN112904327B (zh) | 2021-01-19 | 2021-01-19 | 基于调频模糊函数的复合微动目标参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110071362.XA CN112904327B (zh) | 2021-01-19 | 2021-01-19 | 基于调频模糊函数的复合微动目标参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112904327A true CN112904327A (zh) | 2021-06-04 |
CN112904327B CN112904327B (zh) | 2021-10-15 |
Family
ID=76116023
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110071362.XA Active CN112904327B (zh) | 2021-01-19 | 2021-01-19 | 基于调频模糊函数的复合微动目标参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112904327B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113609446A (zh) * | 2021-07-30 | 2021-11-05 | 北京果仁互动科技有限公司 | 一种运动参数确定方法和装置 |
CN116106857A (zh) * | 2023-04-13 | 2023-05-12 | 中国人民解放军国防科技大学 | 基于稀疏时间-频率-调频率表示的微动形式辨识方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698757A (zh) * | 2013-12-31 | 2014-04-02 | 中国人民解放军国防科学技术大学 | 低频段雷达目标微动特性估计方法 |
CN104007430A (zh) * | 2014-05-29 | 2014-08-27 | 西安电子科技大学 | 基于瞬时调频率估计的进动目标的微多普勒提取方法 |
CN107132536A (zh) * | 2017-04-10 | 2017-09-05 | 中国科学院国家空间科学中心 | 一种消除目标微动对雷达成像干扰的方法 |
CN107271997A (zh) * | 2017-06-05 | 2017-10-20 | 西北工业大学 | 机载多通道cssar地面运动目标运动参数估计方法 |
WO2017187242A1 (en) * | 2016-04-25 | 2017-11-02 | Uhnder, Inc. | On-demand multi-scan micro doppler for vehicle |
CN110221288A (zh) * | 2019-05-17 | 2019-09-10 | 常州第四无线电厂有限公司 | 一种穿墙雷达微动目标检测方法 |
-
2021
- 2021-01-19 CN CN202110071362.XA patent/CN112904327B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698757A (zh) * | 2013-12-31 | 2014-04-02 | 中国人民解放军国防科学技术大学 | 低频段雷达目标微动特性估计方法 |
CN104007430A (zh) * | 2014-05-29 | 2014-08-27 | 西安电子科技大学 | 基于瞬时调频率估计的进动目标的微多普勒提取方法 |
WO2017187242A1 (en) * | 2016-04-25 | 2017-11-02 | Uhnder, Inc. | On-demand multi-scan micro doppler for vehicle |
CN107132536A (zh) * | 2017-04-10 | 2017-09-05 | 中国科学院国家空间科学中心 | 一种消除目标微动对雷达成像干扰的方法 |
CN107271997A (zh) * | 2017-06-05 | 2017-10-20 | 西北工业大学 | 机载多通道cssar地面运动目标运动参数估计方法 |
CN110221288A (zh) * | 2019-05-17 | 2019-09-10 | 常州第四无线电厂有限公司 | 一种穿墙雷达微动目标检测方法 |
Non-Patent Citations (2)
Title |
---|
LAN DU ET AL.: "Micro-Doppler Feature Extraction Based on Time-Frequency Spectrogram for Ground Moving Targets Classification With Low-Resolution Radar", 《IEEE SENSORS JOURNAL》 * |
郭力仁等: "基于粒子滤波的高阶运动目标激光探测微动参数估计", 《光学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113609446A (zh) * | 2021-07-30 | 2021-11-05 | 北京果仁互动科技有限公司 | 一种运动参数确定方法和装置 |
CN116106857A (zh) * | 2023-04-13 | 2023-05-12 | 中国人民解放军国防科技大学 | 基于稀疏时间-频率-调频率表示的微动形式辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112904327B (zh) | 2021-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105425224B (zh) | 车载毫米波雷达系统多目标个数获取方法及装置 | |
CN112904327B (zh) | 基于调频模糊函数的复合微动目标参数估计方法 | |
CN105549001B (zh) | 车载毫米波雷达系统多目标检测方法 | |
CN106597437B (zh) | 高频机载调频连续波sar运动补偿及成像处理方法 | |
CN105549012B (zh) | 车载毫米波雷达系统多目标检测装置 | |
Liu et al. | Bistatic FMCW SAR signal model and imaging approach | |
CN107918123B (zh) | 太赫兹isar目标微动部件分离与特征参数估计方法 | |
Li et al. | Bistatic forward-looking SAR ground moving target detection and imaging | |
CN109324322A (zh) | 一种基于被动相控阵天线的测向与目标识别方法 | |
CN111273285B (zh) | 一种多旋翼无人机的微多普勒谱相关矩阵特征提取方法 | |
Park et al. | Performance analysis of the scenario-based construction method for real target ISAR recognition | |
Park et al. | Construction of ISAR training database for automatic target recognition | |
CN107942323B (zh) | 基于频域熵的进动目标时频曲线提取方法 | |
CN103064084A (zh) | 基于距离频域的解模糊方法 | |
CN110441771B (zh) | 基于方位时间重采样的高速机动飞行sar成像方法 | |
Li et al. | Joint clutter suppression and moving target indication in 2-D azimuth rotated time domain for single-channel bistatic SAR | |
CN111781589A (zh) | 提高微动锥体目标时频分辨率的时频分析方法 | |
Zhu et al. | Classification of ground vehicles based on micro-Doppler effect and singular value decomposition | |
CN113625245B (zh) | 一种弹道导弹中段目标进动参数估计方法 | |
CN106353751B (zh) | 一种高分辨率引信目标检测方法 | |
CN114565020A (zh) | 一种基于深度置信网络和扩展卡尔曼滤波的飞行器传感器信号融合方法 | |
CN110736988B (zh) | 双基地pfa运动目标参数估计和成像方法 | |
CN114488125A (zh) | 利用空间映射优化进动锥体目标参数的方法 | |
Feng et al. | Inverse synthetic aperture radar imaging of maneuvering targets based on joint time-frequency analysis | |
Liu et al. | Bistatic FMCW SAR raw signal simulator for extended scenes |
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 |