CN109799488B - 一种非参数搜索的雷达机动目标长时间相参积累方法 - Google Patents

一种非参数搜索的雷达机动目标长时间相参积累方法 Download PDF

Info

Publication number
CN109799488B
CN109799488B CN201910024363.1A CN201910024363A CN109799488B CN 109799488 B CN109799488 B CN 109799488B CN 201910024363 A CN201910024363 A CN 201910024363A CN 109799488 B CN109799488 B CN 109799488B
Authority
CN
China
Prior art keywords
time
distance
radar
long
fourier transform
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
CN201910024363.1A
Other languages
English (en)
Other versions
CN109799488A (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.)
Naval Aeronautical University
China Academy of Civil Aviation Science and Technology
Original Assignee
Naval Aeronautical University
China Academy of Civil Aviation Science and Technology
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 Naval Aeronautical University, China Academy of Civil Aviation Science and Technology filed Critical Naval Aeronautical University
Priority to CN201910024363.1A priority Critical patent/CN109799488B/zh
Publication of CN109799488A publication Critical patent/CN109799488A/zh
Application granted granted Critical
Publication of CN109799488B publication Critical patent/CN109799488B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种非参数搜索的雷达机动目标长时间相参积累方法,属于雷达信号处理技术领域。首先,沿距离向对脉压后的雷达数据进行傅里叶变换,得到距离频率‑脉间慢时间二维数据;然后,进行非均匀采样降阶运算和变标尺度变换;最后,进行距离频率逆傅里叶变换,非均匀采样慢时间傅里叶变换,则机动目标将在二维平面形成峰值,实现长时间相参积累。本发明仅通过一次非均匀采样降阶运算,将高阶相位信号降阶至一阶相位信号,相比传统逐次降阶法降低了交叉项的影响,提高了参数估计精度;无需参数搜索匹配计算,降低了运算量,适合工程应用;能够同时补偿距离和高阶多普勒走动,实现长时间相参积累,提高雷达强杂波背景下的机动目标检测能力。

Description

一种非参数搜索的雷达机动目标长时间相参积累方法
技术领域
本发明属于雷达信号处理技术领域,更具体地,本发明涉及一种非参数搜索的雷达机动目标长时间相参积累方法,可用于雷达对机动目标检测和估计。
背景技术
具有低可观测性的运动目标在距离和多普勒分辨单元中往往具有较低的信噪比或信杂比(Signal-to-Noise/Clutter Ratio,SNR/SCR),降低了雷达的检测性能。在雷达信号处理中,通常可以延长积累时间以增加目标的能量,达到提高信号SNR/SCR的目的。目前,对机动目标的长时间积累主要面临以下两个方面的问题:一方面由于雷达距离分辨力的不断提高和目标的高速运动,目标回波包络在不同脉冲周期之间走动和弯曲,产生距离徙动效应,使目标能量在距离向分散,例如由目标匀速运动产生的跨距离单元走动,称为一阶距离走动(First-order Range Migration,FRM);由目标加速或减速运动产生的跨距离单元走动,称为二阶距离走动(Second-order Range Migration,SRM);由目标高阶运动,如加速度变化或微动,产生的跨距离单元走动称为三阶距离走动(Third-order RangeMigration,TRM),本专利主要讨论二阶或二阶以上的目标运动。另一方面,目标的加速、减速、高阶运动以及转动等会引起回波相位变化,使雷达回波信号具有时变特性并表现为高阶多项式相位形式,目标的多普勒频率将跨越多个多普勒单元,产生多普勒徙动效应,使得目标能量在频域分散,降低了相参积累增益。因此,对机动目标的长时间相参积累需补偿距离和多普勒走动。目前,雷达机动目标长时间积累检测的主要思路有:
1)基于运动轨迹搜索的长时间非相参积累方法,最为典型的是检测前跟踪技术(Track-Before-Detection,TBD),通过Hough变换(Hough Transform,HT)等轨迹搜索方法,把可能是同一运动轨迹的回波能量进行幅值积累,也称长时间非相参积累,对系统没有严格的相参要求,该方法在工程实现上比较简单,但其信号积累效率和SNR改善均明显低于相参积累方法,不适于复杂环境下微弱动目标的检测;
2)分段长时间积累方法,即非相参和相参积累结合的混合积累方法,首先将距离分段,在段内完成脉冲间的相参积累,然后通过距离包络对齐或Keystone变换(KeystoneTransform,KT)法实现回波的非相参积累,但该方法的积累增益仍非常有限,无法应对强杂波背景和机动目标的复杂运动形式。
3)分步长时间积累方法,即先采用二阶KT和包络插值补偿距离走动,然后,在此基础上采用变换将多普勒走动进行补偿,如相位匹配法、分数阶傅里叶变换(FractionalFourier Transform,FRFT)、Chirplet变换等,但后续多普勒走动补偿的效果受距离走动补偿结果影响,容易造成目标多普勒能量扩散,运动参数估计精度较差。
4)基于参数搜索的长时间相参积累(Long Time Coherent Integration,LTCI)方法,如Radon傅里叶变换(Radon FT,RFT)、Radon-FRFT(RFRFT)等,该方法根据预先设定的目标运动参数(初始距离、速度和加速度)搜索范围,提取位于距离-慢时间二维平面中的目标观测值,然后在相应的变换域选择合适的变换参数对该观测值进行匹配和积累,实现对动目标能量的长时间相参积累,该方法的主要问题在于需要多维参数搜索,导致运算量较大。
5)基于相关及相位差分降阶的相参积累方法,由于机动目标脉压后回波相位可近似为多项式函数,目标机动性越强多项式相位越高,因此,先采用相关或相位差分方法将高阶多项式降阶为低阶多项式,再采用FRFT或FT进行积累,从而降低运算量,但该方法多为逐次降阶,每一次降阶均会产生交叉项,从而影响了积累效果和参数估计精度。
因此,亟需设计一种无需参数搜索,适合机动目标复杂运动,能够同时补偿距离和多普勒走动的长时间相参积累快速实现方法,用于雷达机动目标的检测和参数估计。
发明内容
本发明的目的在于改善雷达机动目标积累增益,提高雷达检测和参数估计性能,实现长时间相参,并且提高运算效率,提出一种非参数搜索的雷达机动目标长时间相参积累方法。其中要解决的技术问题包括:
(1)机动目标的回波SCR低,且具有时变特性,随着积累时间增加,机动目标回波将产生距离走动和多普勒走动效应,导致目标能量在距离维和多普勒维能量发散,检测性能下降;
(2)已有的分段或分步积累方法,积累增益有限,无法应对强杂波背景和机动目标的复杂运动形式;
(3)基于参数搜索的长时间相参积累方法,需要多维参数搜索,以确定最佳参数匹配动目标回波信号,算法复杂度高;
(4)逐次降阶机动目标积累方法,每一次降阶均会产生交叉项,影响积累效果和参数估计精度。
本发明所述的一种非参数搜索的雷达机动目标长时间相参积累方法,其特征在于包括以下技术措施:
步骤一、沿距离向对脉压后的雷达数据进行傅里叶变换,得到距离频率-脉间慢时间二维数据;
步骤二、对脉间慢时间进行非均匀采样降阶运算,并进行变标尺度变换;
步骤三、距离频率维逆傅里叶变换,变标尺度变换后的时间变量维傅里叶变换,若存在机动目标,将在二维平面形成峰值,完成长时间相参积累。
对比现有技术,本技术方案所述的一种非参数搜索的雷达机动目标长时间相参积累方法,有益效果在于:
(1)该方法能够有效处理具有线性调频和二次调频特性的机动目标信号(匀加速和变加速目标),同时补偿距离和高阶多普勒走动,实现长时间相参积累;
(2)该方法仅通过一次非均匀采样降阶运算,即可实现高阶相位信号降至低阶相位信号,
降低了交叉项的影响;
(3)该方法无需运动参数搜索匹配计算,极大降低了运算量,适合工程应用。
附图说明
附图是非参数搜索的雷达机动目标长时间相参积累方法的实施流程图。
具体实施方式
以下结合说明书附图对本发明作进一步详细描述。参照说明书附图,本发明的处理流程分以下步骤:
1)沿距离向对脉压后的雷达数据进行傅里叶变换,得到距离频率-脉间慢时间二维数据
设雷达发射线性调频信号,接收的基带回波信号表示为
Figure BDA0001941928210000031
式中,t为脉内快时间,tm为脉间慢时间,sr(t,tm)为基带回波信号,Ar是回波幅度,K为发射的线性调频信号调频斜率,Tp为脉冲长度,Rs(tm)为雷达与目标的径向距离,c0代表光速,2Rs(tm)/c为回波延迟,λ为信号波长,对式(1)进行距离脉冲压缩,得到脉内积累后的雷达回波数据sPC(t,tm),
Figure BDA0001941928210000032
式中,APC为信号幅度,B为发射信号带宽,假设目标朝向雷达运动,且仅考虑径向速度分量,则目标的距离走动为慢时间的多项式函数,三阶多项式可描述大部分的机动目标的运动,
Figure BDA0001941928210000033
式中,r0表示初始距离,v0、as、gs为矢量,分别代表目标运动初速度、加速度和急动度。
由式(2)可知,第一项sinc函数表示距离走动,在长时间观测条件下,目标包络的峰值位置会随慢时间变化而偏移,当偏移量大于雷达距离单元时,将产生距离徙动效应,目标能量将部分泄露到相邻的距离单元中;第二个指数项表示多普勒信息,机动目标的高阶调频分量将会引起回波多普勒展宽,当多普勒频率跨越多个多普勒单元时,便会产生多普勒徙动效应。
沿距离向对脉压后的雷达数据进行傅里叶变换,得到距离频率-脉间慢时间二维数据,
Figure BDA0001941928210000034
式中,A'PC为信号幅度,fc为发射信号载频,f为距离频率。
2)对脉间慢时间进行非均匀采样降阶运算
非均匀采样降阶运算定义为
S(f,t'm)=SPC(f,tp+t'm)SPC(f,tp-t'm) (4)
式中,t'm为非均匀采样后的慢时间,tp为非均匀采样时间区间的中心,
Figure BDA0001941928210000041
c为尺度因子,控制非均匀采样的密度,τm为新的时间变量。则通过式(4)的运算,式(3)可转变为
Figure BDA0001941928210000042
由式(5)可知,机动目标回波相位已由关于tm的三阶多项式,降阶为关于新的时间变量τm的一次项。
3)变标尺度变换
由于式(5)中距离频率f和τm具有线性耦合关系,仍存在距离走动效应,为消除两者之间的耦合,对时间坐标轴进行尺度变换,称为变标尺度变换,
Figure BDA0001941928210000043
式中,q为尺度变换因子,为保证参数无模糊估计,通常qfc=1,τ'm为变标尺度变换后的时间变量,因此,式(5)可改写为
Figure BDA0001941928210000044
由式(6)可知,经过慢时间非均匀采样降阶运算以及变标尺度变换后,得到了关于(f,τ'm)的二维指数函数,不存在耦合,且为一次项。
4)长时间相参积累
对式(6)分别进行时间变量τ'm傅里叶变换,距离频率f逆傅里叶变换,实现机动目标能量的二维积累
Figure BDA0001941928210000045
式中,
Figure BDA0001941928210000046
表示长时间相参积累后的二维数据,
Figure BDA0001941928210000047
表示对τ'm进行傅里叶变换,IFT()|f表示对f进行逆傅里叶变换,ALTCI表示长时间相参积累幅值,Tn为积累时长。
由式(7)可知,若回波信号存在机动目标,将在
Figure BDA0001941928210000048
二维平面形成峰值,从而实现了机动目标的距离和多普勒走动补偿,完成长时间相参积累。

Claims (4)

1.一种非参数搜索的雷达机动目标长时间相参积累方法,其特征在于,包括以下步骤:
步骤一、沿距离向对脉压后的雷达数据进行傅里叶变换,得到距离频率-脉间慢时间二维数据;
步骤二、对脉间慢时间进行非均匀采样降阶运算,并进行变标尺度变换;
所述非均匀采样降阶运算为:
非均匀采样降阶运算定义为
S(f,t'm)=SPC(f,tp+t'm)SPC(f,tp-t'm)
式中,SPC表示距离频率-脉间慢时间二维数据,f为距离频率,t'm为非均匀采样后的慢时间,tp为非均匀采样时间区间的中心,
Figure FDA0003601728130000011
c为尺度因子,控制非均匀采样的密度,τm为新的时间变量;
所述变标尺度变换为:
为消除距离频率f和τm之间的耦合关系,对时间坐标轴进行尺度变换,
Figure FDA0003601728130000012
式中,q为尺度变换因子,fc为发射信号载频,τ'm为变标尺度变换后的时间变量;
步骤三、距离频率维逆傅里叶变换,变标尺度变换后的时间变量维傅里叶变换,若存在机动目标,将在二维平面形成峰值,完成长时间相参积累。
2.根据权利要求1所述的一种非参数搜索的雷达机动目标长时间相参积累方法,其特征在于,步骤一所述的距离频率-脉间慢时间二维数据计算方法为:
对雷达回波信号进行距离脉冲压缩,沿距离向对脉压后的雷达数据进行傅里叶变换,得到距离频率-脉间慢时间二维数据,
Figure FDA0003601728130000013
式中,A'PC为信号幅度,tm为脉间慢时间,B为发射信号带宽,c0代表光速,Rs(tm)为雷达与目标的径向距离,假设目标朝向雷达运动,且仅考虑径向速度分量,则目标的距离走动为慢时间的多项式函数,
Figure FDA0003601728130000014
式中,r0表示初始距离,v0、as、gs为矢量,分别代表目标运动初速度、加速度和急动度。
3.根据权利要求2所述的一种非参数搜索的雷达机动目标长时间相参积累方法,其特征在于,步骤一的输出经过非均匀采样降阶运算后变为
Figure FDA0003601728130000021
4.根据权利要求3所述的一种非参数搜索的雷达机动目标长时间相参积累方法,其特征在于,步骤三所述的机动目标将在二维平面形成峰值,完成长时间相参积累方法为:
对步骤二的输出分别进行时间变量τ'm傅里叶变换,距离频率f逆傅里叶变换,实现机动目标能量的二维积累
Figure FDA0003601728130000022
式中,
Figure FDA0003601728130000023
表示长时间相参积累后的二维数据,
Figure FDA0003601728130000024
表示对τ'm进行傅里叶变换,IFT()|f表示对f进行逆傅里叶变换,ALTCI表示长时间相参积累幅值,Tn为积累时长,若回波信号存在机动目标,将在
Figure FDA0003601728130000025
二维平面形成峰值。
CN201910024363.1A 2019-01-10 2019-01-10 一种非参数搜索的雷达机动目标长时间相参积累方法 Active CN109799488B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910024363.1A CN109799488B (zh) 2019-01-10 2019-01-10 一种非参数搜索的雷达机动目标长时间相参积累方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910024363.1A CN109799488B (zh) 2019-01-10 2019-01-10 一种非参数搜索的雷达机动目标长时间相参积累方法

Publications (2)

Publication Number Publication Date
CN109799488A CN109799488A (zh) 2019-05-24
CN109799488B true CN109799488B (zh) 2022-07-15

Family

ID=66557109

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910024363.1A Active CN109799488B (zh) 2019-01-10 2019-01-10 一种非参数搜索的雷达机动目标长时间相参积累方法

Country Status (1)

Country Link
CN (1) CN109799488B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112180342A (zh) * 2020-09-29 2021-01-05 中国船舶重工集团公司第七二四研究所 一种海上机动目标长时积累观测参数补偿方法
CN112327285B (zh) * 2020-10-30 2023-10-03 北京理工大学 基于混合广义瑞登傅里叶变换的雷达目标检测方法
CN112346022B (zh) * 2020-11-03 2022-10-25 中国人民解放军海军航空大学 基于二维精简分数阶傅里叶变换的密集假目标干扰抑制方法
CN112363120B (zh) * 2020-11-03 2022-10-25 中国人民解放军海军航空大学 基于二维分数阶傅里叶变换的移频干扰鉴别方法
CN112748415B (zh) * 2020-12-28 2022-08-12 电子科技大学 一种基于广义rft的长时间能量积累方法
CN113406592B (zh) * 2021-06-15 2022-06-10 哈尔滨工业大学 一种高频地波雷达弱目标积累检测方法、计算设备
CN114415122B (zh) * 2022-01-27 2023-03-28 电子科技大学 一种基于频域分段处理的高速目标积累检测方法
CN117784077B (zh) * 2024-02-27 2024-05-24 南京奥联智驾科技有限公司 一种基于频率积累的弱小目标检测方法、终端及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102012503A (zh) * 2010-07-29 2011-04-13 中国人民解放军海军航空工程学院 基于秩和局部估计的选大和选小非参量恒虚警检测器
CN102778674A (zh) * 2012-05-25 2012-11-14 安徽理工大学 非均匀采样的Chirp脉冲时延估计方法
CN103364768A (zh) * 2012-03-31 2013-10-23 中国科学院电子学研究所 压缩感知雷达重构方法
WO2014012828A1 (de) * 2012-07-19 2014-01-23 Deutsches Zentrum für Luft- und Raumfahrt e.V. Methode zur prozessierung von hochauflösenden weltraumgestützt erhaltenen spotlight-sar rohdaten
CN108549067A (zh) * 2018-07-27 2018-09-18 电子科技大学 一种应用于三阶机动目标的相参积累检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2003213762A1 (en) * 2002-03-13 2003-09-29 Raytheon Canada Limited An adaptive system and method for radar detection

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102012503A (zh) * 2010-07-29 2011-04-13 中国人民解放军海军航空工程学院 基于秩和局部估计的选大和选小非参量恒虚警检测器
CN103364768A (zh) * 2012-03-31 2013-10-23 中国科学院电子学研究所 压缩感知雷达重构方法
CN102778674A (zh) * 2012-05-25 2012-11-14 安徽理工大学 非均匀采样的Chirp脉冲时延估计方法
WO2014012828A1 (de) * 2012-07-19 2014-01-23 Deutsches Zentrum für Luft- und Raumfahrt e.V. Methode zur prozessierung von hochauflösenden weltraumgestützt erhaltenen spotlight-sar rohdaten
CN108549067A (zh) * 2018-07-27 2018-09-18 电子科技大学 一种应用于三阶机动目标的相参积累检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"Spatial Compressive Sensing for MIMO Radar";Marco Rossi et al.;《IEEE Transactions on Signal Processing》;20150115;第62卷(第2期);419-430 *
"海杂波FRFT域分形特征在动目标检测中的应用";陈小龙 等;《海军航空工程学院学报》;20100630;第25卷(第6期);611-616 *
"雷达机动目标长时间积累信号处理算法研究";田静;《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》;20150715(第7期);51-64 *
IgorDjurovićet al.."Parameter estimation of non-uniform sampled polynomial-phase signals Using the HOCPF-WD".《Signal Processing》.2015,第106卷253-258. *

Also Published As

Publication number Publication date
CN109799488A (zh) 2019-05-24

Similar Documents

Publication Publication Date Title
CN109799488B (zh) 一种非参数搜索的雷达机动目标长时间相参积累方法
CN109581318B (zh) 基于时间反转非均匀采样的雷达高机动目标相参积累检测方法
CN110109078B (zh) 基于相参频率捷变雷达的速度解模糊方法
CN109541568B (zh) 一种雷达机动目标跨距离和多普勒单元快速相参积累检测方法
CN109633622B (zh) 基于时间反转二阶Keystone变换的雷达机动目标跨单元积累方法
CN102628937B (zh) 基于广义keystone变换和非相参积累的雷达检测方法
CN109521417B (zh) 基于fmcw雷达波形的多目标检测计算方法
CN107907878B (zh) 高精度获取fmcw雷达距离测量值的方法
CN107132534B (zh) 一种高速雷达目标频域检测的优化方法
CN108089171B (zh) 一种针对无人机目标的雷达快速检测方法
CN108919249B (zh) 一种基于二维局部插值的雷达目标距离联合估计方法
CN110824439B (zh) 一种雷达目标快速长时间相参积累方法
CN109613507B (zh) 一种针对高阶机动目标雷达回波的检测方法
CN113093120B (zh) 基于capon算法的PRI捷变雷达目标参数估计方法
CN111551922A (zh) 三维空间双/多基地雷达高速目标检测方法
CN110850384A (zh) 一种基于扫频数据产生宽带去斜回波的方法
CN113805166A (zh) 一种雷达物位计的目标跟踪测距方法及系统
CN108776332B (zh) 利用高频地波雷达检测海上机动目标的方法
CN113253279A (zh) 一种正反线性调频信号测距测速方法
CN113900088A (zh) 一种针对匀加速机动目标的长时间相参积累方法及系统
CN110082748B (zh) 一种无源雷达目标检测方法及装置
CN110609264B (zh) 一种针对脉冲激光雷达的目标回波多普勒频率估计方法
CN116736297B (zh) 一种异构多帧联合相参积累方法
Song et al. LRS: A Positive and Negative Linear Frequency Modulation Method for Ranging and Speed Measurement
Pei et al. Long-time coherent integration method combining pulse compression and Radon fractional Fourier transform

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