CN113406592B - 一种高频地波雷达弱目标积累检测方法、计算设备 - Google Patents
一种高频地波雷达弱目标积累检测方法、计算设备 Download PDFInfo
- Publication number
- CN113406592B CN113406592B CN202110663010.3A CN202110663010A CN113406592B CN 113406592 B CN113406592 B CN 113406592B CN 202110663010 A CN202110663010 A CN 202110663010A CN 113406592 B CN113406592 B CN 113406592B
- Authority
- CN
- China
- Prior art keywords
- wave radar
- ground wave
- frequency ground
- target
- detection
- 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
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- 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
- G01S7/414—Discriminating targets with respect to background clutter
-
- 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
- G01S7/418—Theoretical aspects
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
本发明提供了一种高频地波雷达弱目标积累检测方法、计算设备,其中方法包括:获取高频地波雷达的回波信号,进行脉冲压缩得到高频地波雷达脉压信号;依据高频地波雷达参数设置径向距离、径向速度和真实速度三维联合搜索的搜索参数;根据搜索参数进行高频地波雷达的目标轨迹信号抽取;沿高频地波雷达慢时间维构建相位补偿函数;沿着高频地波雷达脉冲序列,补偿抽取的目标轨迹信号的在脉冲间的相位波动,将补偿后的目标轨迹信号的回波包络同相相加,完成相参积累;基于相参积累的结果进行三维联合恒虚警率检测,判断是否检测到目标;对检测出的目标进行参数估计,并输出。本发明能够提升高频地波雷达对高速、高机动弱目标的检测性能。
Description
技术领域
本发明涉及目标检测技术领域,特别涉及一种高频地波雷达弱目标积累检测方法、计算设备及计算机可读存储介质。
背景技术
高频地波雷达(High Frequency Surface Wave Radar)利用长电磁波在地球表面(包括海面)的绕射效应使电磁波沿地球曲面传播。基于高频段垂直极化电磁波沿海面超水平线传播损耗低的特点,高频地波雷达可以检测超视距的船、飞机、导弹等目标,作用距离大致为300~400km,因此高频地波雷达可以填补微波视距雷达与天波超视距雷达所探测不到的盲区。此外,由于传播稳定、工作波长较长等特点,高频地波雷达能有效抑制现有隐身技术的隐身作用,具有抗隐身的能力,能发现从海平面直到高空的各种目标。
虽然高频地波雷达可对其作用范围内实施全天候连续监视,及时发现作用范围内的海面船只和空中运动目标,测定目标的距离、方位和径向速度,并进行目标的航迹处理,给出目标的基本属性特征;但是,仅根据海面船只和空中运动目标的距离、方位和径向速度,无法正确地对目标属性进行划分。当目标运动方向与雷达视线方向夹角很大时,目标的真实运动速度与径向速度存在较大差距,例如,当超音速战斗机以接近垂直雷达视线的方向飞行时,雷达测得的其径向速度将趋近于零,这给正确划分目标属性造成了巨大的困难,因此,有必要通过技术手段获得目标的真实速度,而目前常规雷达解算方法都无法正确获得高频地波雷达监测目标真实的运动速度。
此外,随着各类高速、高机动性目标不断出现,例如超高音速无人机、高空高音速巡航导弹、弹道导弹与空间高机动飞行器等,特别是高速运动的弱目标,高频地波雷达能接收到的回波信号极其微弱,信噪比(Signal to Noise Ratio,SNR)极低,现有的高频地波雷达脉冲积累法遭遇了巨大挑战,难以获得准确的测量结果。
发明内容
本发明提供了一种高频地波雷达弱目标积累检测方法、计算设备及计算机可读存储介质,能够对高频地波雷达探测结果进行距离徙动校正与多普勒相位补偿,同时能够测量目标真实速度,实现提高雷达回波的积累增益,抑制噪声与杂波干扰,从而有效提升对高速弱目标的检测性能。
第一方面,本发明实施例提供了一种高频地波雷达弱目标积累检测方法,包括:
步骤S1、获取高频地波雷达的回波信号,进行脉冲压缩,得到高频地波雷达脉压信号;
步骤S2、根据高频地波雷达的参数,设置径向距离、径向速度和真实速度三维联合搜索的搜索参数;
步骤S3、对所述高频地波雷达脉压信号进行径向距离、径向速度和真实速度三维联合搜索;
步骤S4、根据搜索参数进行高频地波雷达的目标轨迹信号抽取;
步骤S5、沿高频地波雷达慢时间维构建相位补偿函数;
步骤S6、沿着高频地波雷达脉冲序列,补偿抽取的目标轨迹信号在脉冲间的相位波动,将补偿后目标轨迹信号的回波包络同相相加,完成相参积累;
步骤S7、设置检测参数,基于相参积累的结果进行三维恒虚警率检测,判断是否检测到目标,是则执行下一步骤;
步骤S8、对检测出的目标进行参数估计,根据目标的径向速度与真实速度划分目标属性,并输出目标检测结果与目标属性。
优选地,所述步骤S7包括:
步骤S7-1、初始化检测参数,设定虚警概率Pfa,分别在径向距离域、径向速度域和真实速度域三个检测域中,根据完成相参积累的目标峰值在检测域所占单元数,设定参考单元与保护单元的数目;
步骤S7-2、分别在径向距离域、径向速度域和真实速度域三个检测域中,根据参考单元数目及虚警概率Pfa计算对应的CFAR门限因子,表达式为:
其中,T表示在一个检测域中对应的CFAR门限因子,Nc表示在检测域中参考单元的数目;
步骤S7-3、分别在径向距离域、径向速度域和真实速度域三个检测域中,根据检测单元及多个参考单元,计算噪声基底,噪声基底等于一个检测域中每个参考单元所对应的相参积累回波幅度求和后除以参考单元的数目;
步骤S7-4、分别在径向距离域、径向速度域和真实速度域三个检测域中,计算检测单元对应的检测门限,表达式为:
S=T*N
其中,S为一个检测域中对应的检测门限,T为检测域中对应的CFAR门限因子,N为检测域中对应的噪声基底;
步骤S7-5、判断检测单元在径向速度域、径向距离域和真实速度域三个检测域中的相参积累回波幅度是否均过对应的检测门限,若是,且该检测单元是峰值单元,则判定为一个目标,若否,则判定为非目标;
步骤S7-6、判断是否遍历完所有检测单元,是则执行下一步骤,否则检测单元步进,并返回步骤S7-2。
优选地,所述步骤S8,对检测出的目标进行参数估计时,通过三点二次插值法得到各检测域中峰值目标参数估值大小;
当目标对应的峰值出现在第J个单元,利用峰值所在的第J个单元的相参积累幅度、峰值相邻的第J-1个单元的相参积累幅度及峰值相邻的第J+1个单元的相参积累幅度进行三点二次插值,估算出最大值所在单元对应的参数值,作为一个检测域中峰值目标参数估值大小,表达式为:
其中,xJ表示第J个单元的参数值,f(J)表示第J个单元的相参积累幅度。
优选地,所述步骤S1中,得到高频地波雷达脉压信号的表达式为:
其中,tm=mTr(m=0,1,...,M-1)表示慢时间,M为高频地波雷达脉冲积累数,Tr为高频地波雷达脉冲周期,τ表示快时间,K为目标数,Ai1为第i个目标高频地波雷达脉压信号幅度,B表示高频地波雷达发射信号带宽,c为光速,λ=c/fc表示高频地波雷达波长,fc为高频地波雷达载波频率,ri(tm)为第i个目标的径向距离方程。
优选地,所述步骤S2中,设置径向距离、径向速度和真实速度三维联合搜索的搜索参数,包括:
步骤S2-1、根据高频地波雷达的威力范围0~Rmax,设定径向距离搜索范围rmin~rmax:
rmin=0,rmax=Rmax
根据高频地波雷达探测目标类型最大速度Vmax,设定真实速度搜索范围vmin~vmax:
vmin=0,vmax=Vmax
步骤S2-2、根据高频地波雷达距离波门大小设置径向距离搜索间隔:
其中fs为高频地波雷达采样率大小,c为光速;
根据高频地波雷达多普勒分辨率设置径向速度搜索间隔:
根据高频地波雷达多普勒分辨率设置速度搜索间隔:
步骤S2-3、确定目标的距离、径向速度和真实速度的搜索数目,分别表示为:
其中round(·)表示取整操作,取整方法是四舍五入;
步骤S2-4、得到离散化的搜索距离、搜索径向速度和搜索真实速度,分别表示为:
rn=nΔr,n=0,1,…,Nr-1
vq=qΔv,q=0,1,…,Nv-1
优选地,所述步骤S3中,对所述高频地波雷达脉压信号进行径向距离、径向速度和真实速度三维联合搜索,包括:
得到径向距离、径向速度和真实速度三维联合搜索的定义式为:
代入高频地波雷达脉压信号sc(tm,τ)表达式,有:
优选地,所述步骤S4中,根据搜索参数进行高频地波雷达的目标轨迹信号抽取,包括:
将每个脉冲周期回波包络对齐到初始径向距离所在的距离单元,完成高频地波雷达的目标轨迹信号抽取,表达式为:
s(tm)=sc(tm,k(tm))
其中,s(tm)表示抽取的tm时刻的目标轨迹信号。
优选地,所述步骤S5中,沿高频地波雷达慢时间维构建相位补偿函数,表达式为:
所述步骤S6中,将补偿后目标轨迹信号的回波包络同相相加,完成相参积累,表达式为:
第二方面,本发明实施例还提供了一种计算设备,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器执行所述计算机程序时,实现本说明书任一实施例所述的高频地波雷达弱目标积累检测方法。
第三方面,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,当所述计算机程序在计算机中执行时,令计算机执行本说明书任一实施例所述的高频地波雷达弱目标积累检测方法。
本发明实施例提供了一种高频地波雷达弱目标积累检测方法、计算设备及计算机可读存储介质,本发明针对高频地波雷达,根据雷达系统的参数,设置径向距离、径向速度和真实速度三维联合搜索的搜索范围与搜索间隔,对高频地波雷达脉压信号进行径向距离、径向速度和真实速度三维联合搜索,基于高频地波雷达探测结果进行距离徙动校正与多普勒相位补偿,合理利用包含径向速度、径向距离、真实速度三维信息的匀速运动模型,可以在完成高频地波雷达目标能量积累的同时,具备高信噪比增益的特性,并实现对目标真实速度的估计,以便对各类目标正确划分目标属性,同时通过增加目标真实速度域检测信息,进一步提高了高速弱目标的检测概率,提升高频地波雷达对微弱高速、高机动目标的检测性能。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一实施例提供的一种高频地波雷达弱目标积累检测方法流程图;
图2(a)是本发明一实施例提供的一种高频地波雷达弱目标积累检测方法所得的RD图;图2(b)为MTD方法所得的RD图;
图3示出了本发明一实施例提供的一种高频地波雷达弱目标积累检测方法与MTD方法在径向速度域检测结果;
图4示出了本发明一实施例提供的一种高频地波雷达弱目标积累检测方法与MTD方法在径向距离域检测结果;
图5示出了本发明一实施例提供的一种高频地波雷达弱目标积累检测方法与MTD方法在真实速度域检测结果;
图6示出了本发明一实施例提供的一种高频地波雷达弱目标积累检测方法所绘制的目标真实速度谱;
图7示出了本发明一实施例提供的一种高频地波雷达弱目标积累检测方法与现有径向匀速运动相参积累方法性能对比。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
如前所述,目标的高速特性导致传统相参积累过程中将产生距离徙动(RangeMigration,RM)现象,目标的高机动特性也将导致多普勒频率徙动(Doppler FrequencyMigration,DFM)问题,使得传统的高频地波雷达相参积累方法在高速弱目标探测中面临着严峻的挑战。本发明借鉴广义Radon傅里叶变换(Generalized Radon Fourier Transform,GRFT)方法多维联合搜索的思想,对高频地波雷达脉压后的回波数据,按照径向与非径向统一的匀速模型,进行径向速度、径向距离以及真实速度三维联合搜索,通过抽取目标轨迹信号、构造相应的相位补偿函数,对抽取出的目标轨迹信号进行相位补偿,将补偿后的目标轨迹信号同相相加,完成对目标能量的相参积累,本发明提出一种高频地波雷达弱目标积累检测方法,该方法基于GRFT方法且具备测速功能,记为CVGRFT方法。
下面描述以上构思的具体实现方式。
请参考图1,本发明实施例提供了一种高频地波雷达弱目标积累检测方法,该方法包括:
步骤S1、获取高频地波雷达的回波信号,进行脉冲压缩,得到高频地波雷达脉压信号;
步骤S2、根据高频地波雷达的参数,设置径向距离、径向速度和真实速度三维联合搜索的搜索参数;
步骤S3、对高频地波雷达脉压信号进行径向距离、径向速度和真实速度三维联合搜索;
步骤S4、根据搜索参数进行高频地波雷达的目标轨迹信号抽取;
步骤S5、沿高频地波雷达慢时间维构建相位补偿函数,用于补偿抽取的目标轨迹信号在脉冲间的相位波动;
步骤S6、沿着高频地波雷达脉冲序列,补偿抽取的目标轨迹信号在脉冲间的相位波动,将补偿后的目标轨迹信号的回波包络同相相加,完成相参积累;
步骤S7、设置检测参数,基于相参积累的结果进行三维恒虚警率检测,判断是否检测到目标,若检测到目标,则执行下一步骤,即继续执行步骤S8,否则检测单元步进,直至遍历所有检测单元;
步骤S8、对检测出的目标进行参数估计,根据目标的径向速度与真实速度划分目标属性,并输出目标检测结果与目标属性。
本发明实施例中,在完成高频地波雷达回波信号相参积累之后,进行径向距离、径向速度与真实速度三维搜索空间的联合CFAR(恒虚警率)检测与目标参数估计,通过增加真实速度谱信息来进一步提高真实目标的检测概率,具有较大的性能优势,可以有效地提高高频地波雷达高速运动目标的检测概率,减小漏警概率。
下面描述图1所示的各个步骤的执行方式。
可选地,针对步骤S1中,得到高频地波雷达脉压信号sc(tm,τ)的表达式为:
其中,tm=mTr(m=0,1,...,M-1)表示慢时间,M为高频地波雷达脉冲积累数,Tr为高频地波雷达脉冲周期,τ表示快时间,K为目标数,Ai1为第i个目标高频地波雷达脉压信号幅度,B表示高频地波雷达发射信号带宽,c为光速,λ=c/fc表示高频地波雷达波长,fc为高频地波雷达载波频率,ri(tm)为第i个目标的径向距离方程。
可选地,针对步骤S2中,设置径向距离、径向速度和真实速度三维联合搜索的搜索参数,包括:
步骤S2-1、依据高频地波雷达参数初始化搜索参数,确定目标的径向距离、径向速度和真实速度的搜索范围。
根据高频地波雷达的威力范围0~Rmax,设定径向距离搜索范围rmin~rmax:
rmin=0,rmax=Rmax
其中λ为高频地波雷达波长,Tr为高频地波雷达脉冲周期,M为高频地波雷达脉冲积累数;
根据高频地波雷达探测目标类型(包括海面舰船和海面飞行目标)最大速度Vmax,设定真实速度搜索范围vmin~vmax:
vmin=0,vmax=Vmax
步骤S2-2、根据高频地波雷达参数确定径向距离、径向速度和真实速度的搜索间隔。
根据高频地波雷达距离波门大小设置径向距离的搜索间隔,表示为:
其中,fs为高频地波雷达采样率,c为光速;
根据高频地波雷达多普勒分辨率设置径向速度搜索间隔,表示为:
真实速度的搜索间隔与径向速度的搜索间隔保持一致,即,根据高频地波雷达多普勒分辨率设置速度搜索间隔,表示为:
步骤S2-3、确定目标的距离、径向速度和真实速度的搜索数目,分别表示为:
其中round(·)表示取整操作,取整方法是四舍五入,M为高频地波雷达脉冲积累数;
步骤S2-4、得到离散化的搜索距离、搜索径向速度和搜索真实速度,分别表示为:
rn=nΔr,n=0,1,…,Nr-1
vq=qΔv,q=0,1,…,Nv-1
可选地,针对步骤S3中,对高频地波雷达脉压信号进行径向距离、径向速度和真实速度三维联合搜索,包括:
得到径向距离、径向速度和真实速度三维联合搜索的定义式为:
代入高频地波雷达脉压信号sc(tm,τ)表达式,有:
由上式可知,只有在每个脉冲周期tm时刻,r(tm)与ri(tm)恒等时,有积累峰值;同时,相位项中ri(tm)-r(tm)利用平方差公式展开有:
可得到,r(tm)与ri(tm)恒等时有积累峰值,等价于搜索检测域目标参数相符时有积累峰值。
可选地,针对步骤S4中,根据搜索参数进行高频地波雷达的目标轨迹信号抽取,包括:
round(·)表示取整操作,取整方法是四舍五入;
将每个脉冲周期回波包络对齐到初始径向距离所在的距离单元,完成高频地波雷达的目标轨迹信号抽取,表达式为:
s(tm)=sc(tm,k(tm))
其中,s(tm)表示抽取的tm时刻的目标轨迹信号,也即抽取的高频地波雷达目标tm时刻的目标轨迹。通过把回波包络对齐在初始径向距离单元,也即避免了距离徙动现象。
可选地,针对步骤S5中,沿高频地波雷达慢时间维构建相位补偿函数,表达式为:
其中,φ(tm)表示tm时刻对应的相位补偿函数;
步骤S6中,将补偿后目标轨迹信号的回波包络同相相加,完成相参积累,表达式为:
可选地,步骤S7中,设置检测参数,基于相参积累的结果进行恒虚警率检测,判断是否检测到目标,是则执行下一步骤,包括:
步骤S7-1、初始化检测参数,设定虚警概率Pfa,分别在径向距离域、径向速度域和真实速度域三个检测域中,根据完成相参积累的目标峰值在检测域所占单元数,设定参考单元与保护单元的数目;保护单元为目标单元在检测域相邻的左、右单元,参考单元在保护单元的外侧,即远离检测单元的两侧,一般选择2的整数倍,如8、16、32;
步骤S7-2、分别在径向距离域、径向速度域和真实速度域三个检测域中,根据参考单元数目及虚警概率Pfa计算对应的CFAR门限因子,表达式为:
其中,T表示在一个检测域中对应的CFAR门限因子,Nc表示在检测域中参考单元的数目;
步骤S7-3、分别在径向距离域、径向速度域和真实速度域三个检测域中,根据检测单元及多个参考单元,计算噪声基底,噪声基底等于一个检测域中每个参考单元所对应的相参积累回波幅度求和后除以参考单元的数目,即,分别在径向速度域、径向距离域与真实速度域以其参考单元的平均作为该域的噪声基底大小,表达式为:
步骤S7-4、分别在径向距离域、径向速度域和真实速度域三个检测域中,计算检测单元对应的检测门限,表达式为:
S=T*N
其中,S为一个检测域中对应的检测门限,T为检测域中对应的CFAR门限因子,N为检测域中对应的噪声基底;
步骤S7-5、判断检测单元在径向速度域、径向距离域和真实速度域三个检测域中的相参积累回波幅度是否均过对应的检测门限,若是,且该检测单元是峰值单元,则判定为一个目标,若否,则判定为非目标;
步骤S7-6、判断是否遍历完所有检测单元,是则执行下一步骤,否则检测单元步进,并返回步骤S7-2。
可选地,步骤S8中,对检测出的目标进行参数估计时,通过三点二次插值法得到各检测域中峰值目标参数估值大小。目标相参积累形成的尖峰通常会占据若干单元,因此若在积累记结果中遇到一个峰值,则其可能为目标,这就是文中所说的峰值目标;而目标峰值是在这若干个单元中,幅度最大的单元的幅度信息,即目标峰值。
在一个检测域中通过三点二次插值法进行参数估计,包括:
当目标对应的峰值出现在第J个单元,利用峰值所在的第J个单元的相参积累幅度、峰值相邻的第J-1个单元的相参积累幅度及峰值相邻的第J+1个单元的相参积累幅度进行三点二次插值,估算出最大值所在单元对应的参数值,作为一个检测域中峰值目标参数估值大小,表达式为:
采用上述方式在分别在径向距离域、径向速度域和真实速度域三个检测域计算峰值目标参数估值大小,即可得到估计的径向距离、估计的径向速度和估计的真实速度。
可选地,步骤S8中,对检测出的目标进行参数估计,包括:
对于被检测为目标的三维单元(i,j,k),可以计算得到目标径向距离、径向速度与真实速度的初步估计值为:
得到初步估计值后,可通过三点二次插值法进行修正。
本发明还对提供的高频地波雷达弱目标积累检测方法性能进行了验证。如图2(a)至图5所示,对5120个高频地波雷达实测脉冲数据,分别进行MTD(动目标检测)方法、本发明提供的高频地波雷达弱目标积累检测方法(简称CVGRFT方法)处理,得到RD图分别如图2(b)、2(a)所示,对积累的结果进行基于CVGRFT三个检测域的CFAR检测,在径向距离域、径向速度域以及真实速度域三个检测域中分别取定参考单元、保护单元以及CFAR门限因子,检测结果如图3至图5所示,从检测结果可以看到,MTD方法共检测出10个目标,而CVGRFT方法积累下共检测出16个目标,其中被MTD方法积累后检测出来的10个目标均能够在本发明方法积累情况下被检测出来,同时有6个目标在MTD积累后信噪比未超过检测门限,而在本发明方法情况下能够被正确检测出来,由此可见,本发明方法对于高频地波雷达系统而言,在具备目标真实速度测量功能的同时,能够有效地提高高速弱目标的相参积累性能,进一步提高高速弱目标的检测概率。图6示出了高频地波雷达弱目标检测方法所绘制的目标真实速度谱,对应测速功能。测速结果大小为262.9m/s,径向速度峰值为-129.2m/s,目标远离雷达运动。在真实速度与径向距离切片下,对比了径向速度谱,图7示出了高频地波雷达弱目标积累检测方法与现有径向匀速运动相参积累方法(即MTD方法、基于Keystone变换的相参积累方法,简称KT方法、基于Radon傅里叶变换的相参积累方法,简称RFT方法)性能对比,在同样积累时长条件下,CVGRFT方法积累峰值113.7dB,RFT方法峰值均为109.5dB,KT积累峰值为109.2dB,MTD方法积累峰值109.3dB,可知本发明提供的CVGRFT方法积累尖峰最为尖锐,杂波基底有所减低,信杂比提升。
本发明实施例还提供了一种计算设备,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器执行所述计算机程序时,实现本发明任一实施例中的一种高频地波雷达弱目标积累检测方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序在被处理器执行时,使所述处理器执行本发明任一实施例中的一种高频地波雷达弱目标积累检测方法。
具体地,可以提供配有存储介质的系统或者装置,在该存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码,且使该系统或者装置的计算机(或CPU或MPU)读出并执行存储在存储介质中的程序代码。
在这种情况下,从存储介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能,因此程序代码和存储程序代码的存储介质构成了本发明的一部分。
用于提供程序代码的存储介质实施例包括软盘、硬盘、磁光盘、光盘(如CD-ROM、CD-R、CD-RW、DVD-ROM、DVD-RAM、DVD-RW、DVD+RW)、磁带、非易失性存储卡和ROM。可选择地,可以由通信网络从服务器计算机上下载程序代码。
此外,应该清楚的是,不仅可以通过执行计算机所读出的程序代码,而且可以通过基于程序代码的指令使计算机上操作的操作系统等来完成部分或者全部的实际操作,从而实现上述实施例中任意一项实施例的功能。
此外,可以理解的是,将由存储介质读出的程序代码写到插入计算机内的扩展板中所设置的存储器中或者写到与计算机相连接的扩展模块中设置的存储器中,随后基于程序代码的指令使安装在扩展板或者扩展模块上的CPU等来执行部分和全部实际操作,从而实现上述实施例中任一实施例的功能。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种高频地波雷达弱目标积累检测方法,其特征在于,包括:
步骤S1、获取高频地波雷达的回波信号,进行脉冲压缩,得到高频地波雷达脉压信号;
步骤S2、根据高频地波雷达的参数,设置径向距离、径向速度和真实速度三维联合搜索的搜索参数;
步骤S3、对所述高频地波雷达脉压信号进行径向距离、径向速度和真实速度三维联合搜索;
步骤S4、根据搜索参数进行高频地波雷达的目标轨迹信号抽取;
步骤S5、沿高频地波雷达慢时间维构建相位补偿函数;
步骤S6、沿着高频地波雷达脉冲序列,补偿抽取的目标轨迹信号在脉冲间的相位波动,将补偿后目标轨迹信号的回波包络同相相加,完成相参积累;
步骤S7、设置检测参数,基于相参积累的结果进行三维恒虚警率检测,判断是否检测到目标,是则执行下一步骤;
步骤S8、对检测出的目标进行参数估计,根据目标的径向速度与真实速度划分目标属性,并输出目标检测结果与目标属性;
所述步骤S7包括:
步骤S7-1、初始化检测参数,设定虚警概率Pfa,分别在径向距离域、径向速度域和真实速度域三个检测域中,根据完成相参积累的目标峰值在检测域所占单元数,设定参考单元与保护单元的数目;
步骤S7-2、分别在径向距离域、径向速度域和真实速度域三个检测域中,根据参考单元数目及虚警概率Pfa计算对应的CFAR门限因子,表达式为:
其中,T表示在一个检测域中对应的CFAR门限因子,Nc表示在检测域中参考单元的数目;
步骤S7-3、分别在径向距离域、径向速度域和真实速度域三个检测域中,根据检测单元及多个参考单元,计算噪声基底,噪声基底等于一个检测域中每个参考单元所对应的相参积累回波幅度求和后除以参考单元的数目;
步骤S7-4、分别在径向距离域、径向速度域和真实速度域三个检测域中,计算检测单元对应的检测门限,表达式为:
S=T*N
其中,S为一个检测域中对应的检测门限,T为检测域中对应的CFAR门限因子,N为检测域中对应的噪声基底;
步骤S7-5、判断检测单元在径向速度域、径向距离域和真实速度域三个检测域中的相参积累回波幅度是否均过对应的检测门限,若是,且该检测单元是峰值单元,则判定为一个目标,若否,则判定为非目标;
步骤S7-6、判断是否遍历完所有检测单元,是则执行下一步骤,否则检测单元步进,并返回步骤S7-2。
4.根据权利要求3所述的高频地波雷达弱目标积累检测方法,其特征在于,所述步骤S2中,设置径向距离、径向速度和真实速度三维联合搜索的搜索参数,包括:
步骤S2-1、根据高频地波雷达的威力范围0~Rmax,设定径向距离搜索范围rmin~rmax:
rmin=0,rmax=Rmax
根据高频地波雷达探测目标类型最大速度Vmax,设定真实速度搜索范围vmin~vmax:
vmin=0,vmax=Vmax
步骤S2-2、根据高频地波雷达距离波门大小设置径向距离搜索间隔:
其中fs为高频地波雷达采样率大小,c为光速;
根据高频地波雷达多普勒分辨率设置径向速度搜索间隔:
根据高频地波雷达多普勒分辨率设置速度搜索间隔:
步骤S2-3、确定目标的距离、径向速度和真实速度的搜索数目,分别表示为:
Nr=M
其中round(·)表示取整操作,取整方法是四舍五入;
步骤S2-4、得到离散化的搜索距离、搜索径向速度和搜索真实速度,分别表示为:
rn=nΔr,n=0,1,…,Nr-1
vq=qΔv,q=0,1,…,Nv-1
8.一种计算设备,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器执行所述计算机程序时,实现如权利要求1-7中任一项所述的高频地波雷达弱目标积累检测方法。
9.一种计算机可读存储介质,其上存储有计算机程序,当所述计算机程序在计算机中执行时,令计算机执行权利要求1-7中任一项所述的高频地波雷达弱目标积累检测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110663010.3A CN113406592B (zh) | 2021-06-15 | 2021-06-15 | 一种高频地波雷达弱目标积累检测方法、计算设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110663010.3A CN113406592B (zh) | 2021-06-15 | 2021-06-15 | 一种高频地波雷达弱目标积累检测方法、计算设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113406592A CN113406592A (zh) | 2021-09-17 |
CN113406592B true CN113406592B (zh) | 2022-06-10 |
Family
ID=77684059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110663010.3A Active CN113406592B (zh) | 2021-06-15 | 2021-06-15 | 一种高频地波雷达弱目标积累检测方法、计算设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113406592B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115356694B (zh) * | 2022-08-26 | 2023-08-22 | 哈尔滨工业大学(威海) | 高频地波雷达抗冲击干扰方法、系统、计算机设备及应用 |
CN117268496B (zh) * | 2023-11-21 | 2024-01-30 | 无锡航征科技有限公司 | 基于雷达波的水位测量方法、装置及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005111655A2 (en) * | 2004-05-03 | 2005-11-24 | Raytheon Canada Limited | System and method for concurrent operation of multiple radar or active sonar systems on a common frequency |
WO2008093092A2 (en) * | 2007-01-31 | 2008-08-07 | Cambridge Consultants Limited | Adaptive radar |
CN101738606A (zh) * | 2008-11-21 | 2010-06-16 | 清华大学 | 基于广义多普勒滤波器组的雷达目标相参积累检测方法 |
CN103344949A (zh) * | 2013-06-18 | 2013-10-09 | 中国人民解放军海军航空工程学院 | 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法 |
CN105652256A (zh) * | 2015-12-28 | 2016-06-08 | 哈尔滨工业大学 | 一种基于极化信息的高频地波雷达tbd方法 |
CN107561508A (zh) * | 2017-08-24 | 2018-01-09 | 电子科技大学 | 一种用于匀加速运动目标的相参积累检测方法 |
CN109799488A (zh) * | 2019-01-10 | 2019-05-24 | 中国人民解放军海军航空大学 | 一种非参数搜索的雷达机动目标长时间相参积累方法 |
CN111505599A (zh) * | 2020-05-07 | 2020-08-07 | 哈尔滨工业大学 | 一种基于精确距离演化模型的相参积累检测方法 |
CN112255609A (zh) * | 2020-10-14 | 2021-01-22 | 哈尔滨工业大学 | 一种恒加速机动目标相参积累检测方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140121123A1 (en) * | 2010-10-29 | 2014-05-01 | Kevin Caili Wang | Methods for diversifying antibodies, antibodies derived therefrom and uses thereof |
US9638793B2 (en) * | 2013-01-15 | 2017-05-02 | Raytheon Canada Limited | Virtual antenna extension for sampled aperture arrays |
CN103176178B (zh) * | 2013-02-04 | 2014-11-19 | 中国人民解放军海军航空工程学院 | 雷达动目标Radon-分数阶傅里叶变换长时间相参积累检测方法 |
CN103323829B (zh) * | 2013-06-04 | 2015-09-09 | 中国人民解放军海军航空工程学院 | 基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法 |
WO2015173891A1 (ja) * | 2014-05-13 | 2015-11-19 | 三菱電機株式会社 | レーダ装置 |
CN106249216B (zh) * | 2016-08-08 | 2018-08-21 | 哈尔滨工业大学 | 基于高频地波雷达的静止目标双路径回波信息实现电离层高度估计的方法 |
CN106970371B (zh) * | 2017-04-28 | 2019-05-14 | 电子科技大学 | 一种基于Keystone和匹配滤波的目标检测方法 |
CN111123214B (zh) * | 2019-12-18 | 2023-09-01 | 南京理工大学 | 多项式旋转-多项式傅里叶变换的高速高机动目标检测方法 |
CN112924959B (zh) * | 2021-02-02 | 2022-05-31 | 哈尔滨工业大学 | 基于带宽合成的分布式高频地波雷达目标定位和测速方法 |
-
2021
- 2021-06-15 CN CN202110663010.3A patent/CN113406592B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005111655A2 (en) * | 2004-05-03 | 2005-11-24 | Raytheon Canada Limited | System and method for concurrent operation of multiple radar or active sonar systems on a common frequency |
WO2008093092A2 (en) * | 2007-01-31 | 2008-08-07 | Cambridge Consultants Limited | Adaptive radar |
CN101738606A (zh) * | 2008-11-21 | 2010-06-16 | 清华大学 | 基于广义多普勒滤波器组的雷达目标相参积累检测方法 |
CN103344949A (zh) * | 2013-06-18 | 2013-10-09 | 中国人民解放军海军航空工程学院 | 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法 |
CN105652256A (zh) * | 2015-12-28 | 2016-06-08 | 哈尔滨工业大学 | 一种基于极化信息的高频地波雷达tbd方法 |
CN107561508A (zh) * | 2017-08-24 | 2018-01-09 | 电子科技大学 | 一种用于匀加速运动目标的相参积累检测方法 |
CN109799488A (zh) * | 2019-01-10 | 2019-05-24 | 中国人民解放军海军航空大学 | 一种非参数搜索的雷达机动目标长时间相参积累方法 |
CN111505599A (zh) * | 2020-05-07 | 2020-08-07 | 哈尔滨工业大学 | 一种基于精确距离演化模型的相参积累检测方法 |
CN112255609A (zh) * | 2020-10-14 | 2021-01-22 | 哈尔滨工业大学 | 一种恒加速机动目标相参积累检测方法 |
Non-Patent Citations (6)
Title |
---|
adon-Fourier Transform for Radar Target Detection, I: Generalized Doppler Filter Bank;Jia Xu; Ji Yu; Ying-Ning Peng; Xiang-Gen Xia;《IEEE 》;20110415;第47卷(第2期);第1186-1202页 * |
Anna Dzvonkovskaya ; Klaus-Werner Gurgel ; Hermann Rohling ; Thomas Sc.Low power High Frequency Surface Wave Radar application for ship detection and tracking.《2008 International Conference on Radar》.2008,第627-632页. * |
MIMO雷达微弱目标长时积累技术综述;陈小龙,黄勇,关键,宋伟健,薛永华;《信号处理》;20201125;第36卷(第12期);第1947-1964页 * |
Position accuracy for 3-D curved surface in bistatic high frequency surface wave radar;Zhao Kong-rui; Zhou Gong-jian;Quan Tai-fan;《IEEE International Geoscience and Remote Sensing Symposium (IGARSS)》;20121112;第34卷(第7期);第1344-1349页 * |
Weighted distance iteration matrix CFAR detector in Sea clutter;Ye, L; Jian, FS;Chen, QS;《IET RADAR SONAR AND NAVIGATION》;20200902;第14卷(第10期);第1631-1639页 * |
高速目标的长时间相参积累算法研究;杨宇超;《中国优秀硕士学位论文全文数据库信息科技辑 》;20210215(第2期);第I136-1413页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113406592A (zh) | 2021-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110412559B (zh) | 分布式无人机mimo雷达的非相参融合目标检测方法 | |
CN108761419B (zh) | 基于组合空时主通道自适应处理的低空风切变风速估计方法 | |
CN113406592B (zh) | 一种高频地波雷达弱目标积累检测方法、计算设备 | |
CN107976660B (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN103529437B (zh) | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 | |
CN107450055B (zh) | 基于离散线性调频傅立叶变换的高速机动目标检测方法 | |
CN107607943B (zh) | 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法 | |
CN107607938B (zh) | 一种频率分集阵列雷达低空目标检测方法 | |
CN104155632A (zh) | 一种基于局部相关性的改进子空间海杂波抑制方法 | |
CN107607914B (zh) | 一种弹载pd体制雷达超低空目标与多径回波建模方法 | |
CN110146873B (zh) | 分布式非相参雷达的目标位置和速度的估计方法 | |
Gingras et al. | Electromagnetic matched-field processing: Basic concepts and tropospheric simulations | |
CN110673130A (zh) | 一种基于航迹关联的运动目标航迹跟踪方法 | |
CN113885006B (zh) | 一种机动转弯目标相参积累检测方法 | |
CN113311398B (zh) | 强杂波复杂背景高机动弱小目标的跟踪方法 | |
CN108490425B (zh) | 一种双基地mimo雷达的测角方法 | |
Fang et al. | E 2 DTF: An End-to-End Detection and Tracking Framework for Multiple Micro-UAVs With FMCW-MIMO Radar | |
CN113189554A (zh) | 雷达实测回波数据的处理方法、电子设备及存储介质 | |
Chen et al. | Radar coherent detection for maneuvering target based on product‐scaled integrated cubic phase function | |
Stinco et al. | Non‐cooperative target recognition in multistatic radar systems | |
CN113985376B (zh) | 一种雷达综合显控激励系统 | |
CN114355324A (zh) | 一种航迹生成方法 | |
CN113687321B (zh) | 雷达目标探测距离评估方法和装置 | |
CN112630744B (zh) | 一种多相参积累方法融合的海上小目标检测方法及系统 | |
CN111679270B (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 |