CN108415010B - 一种基于梯形lfmcw调制的雷达多目标检测方法 - Google Patents

一种基于梯形lfmcw调制的雷达多目标检测方法 Download PDF

Info

Publication number
CN108415010B
CN108415010B CN201810404328.8A CN201810404328A CN108415010B CN 108415010 B CN108415010 B CN 108415010B CN 201810404328 A CN201810404328 A CN 201810404328A CN 108415010 B CN108415010 B CN 108415010B
Authority
CN
China
Prior art keywords
frequency
spectrum peak
frequency band
peak
spectrum
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
CN201810404328.8A
Other languages
English (en)
Other versions
CN108415010A (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.)
Xidian University
Xian Cetc Xidian University Radar Technology Collaborative Innovation Research Institute Co Ltd
Original Assignee
Xidian University
Xian Cetc Xidian University Radar Technology Collaborative Innovation Research Institute Co Ltd
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 Xidian University, Xian Cetc Xidian University Radar Technology Collaborative Innovation Research Institute Co Ltd filed Critical Xidian University
Priority to CN201810404328.8A priority Critical patent/CN108415010B/zh
Publication of CN108415010A publication Critical patent/CN108415010A/zh
Application granted granted Critical
Publication of CN108415010B publication Critical patent/CN108415010B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/583Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
    • G01S13/584Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements

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

本发明属于雷达目标检测技术领域,公开了一种基于梯形LFMCW调制的雷达多目标检测方法,包括:雷达目标发射梯形LFMCW信号,接收目标回波信号,产生数字差频信号;计算上扫频段谱峰、恒频段谱峰、下扫频段谱峰,将上扫频段和下扫频段谱峰两两组合得到两谱峰组合,并计算对应的频率和与幅度差,按幅度差升序排列得到序列S;按频率优先原则,标记满足第一频率匹配条件的三谱峰组合;若存在未成功匹配的恒频段谱峰,按频率优先原则,标记满足第二频率匹配条件的三谱峰组合;若P不等于L,按幅度信息优先原则,标记满足第二频率匹配条件的L‑P对三谱峰组合;计算得到目标的距离估计值和速度估计值。本发明能够减少虚假目标,提高检测可靠性。

Description

一种基于梯形LFMCW调制的雷达多目标检测方法
技术领域
本发明涉及雷达目标检测技术领域,尤其涉及一种基于梯形线性调频连续波(Linear Frequency Modulation Continuous Wave,LFMCW)调制的雷达多目标检测方法,可适用于复杂环境条件下的多目标雷达检测。
背景技术
调频连续波雷达,是指发射频率受特定信号调制的连续波雷达,它是一种发射线性调频连续波的雷达。调频连续波雷达通过比较回波信号频率与发射信号的频率之差的方法来获得目标的距离信息。其他测距测速雷达相比,调频连续波雷达具有分辨率高、无测距盲区、低截获率、低成本等一系列优点,被广泛应用于汽车防撞、电子警察、无接触测量、高度计等领域。
目前,调频连续波雷达的目标检测方法主要有基于频移键控连续波调制(Frequency-shift Keying Continuous Wave,FSKCW)的目标检测方法、基于对称三角波调制的目标检测方法、基于锯齿波调制的目标检测方法及基于梯形波调制的目标检测方法。其中,基于频移键控连续波调制的目标检测方法通过检测目标多普勒频移以及中频信号相位差得到运动目标参数。对于静止目标,由于其多普勒频移为零,因此使用频移键控连续波无法获得静止目标的目标信息;而基于对称三角波调制的目标检测方法使用对称三角波能够在一定程度上克服该问题,但是在多目标条件下存在上下扫频段频谱配对难的问题,因此该方法难以分辨多目标;基于锯齿波调制的目标检测方法通过对多个调频周期差频信号做二维快速傅里叶变换处理可以较好地解决速度距离耦合以及多目标分辨问题,然而通常由于硬件调频斜率限制以及测距精度要求,调频周期过长,从而导致速度模糊;基于梯形波调制的目标检测传统方法通过对上、下扫频段和恒频段的差拍信号频谱进行定量分析,能够解决上述目标检测方法存在的静止目标检测、上下扫频配对、速度模糊等一系列问题,但其实际应用中,例如交通流量统计、汽车防撞等多目标复杂环境下,易产生虚假目标。
发明内容
有鉴于此,本发明提供一种基于梯形LFMCW调制的雷达多目标检测方法,以至少解决在多目标复杂环境下现有技术易产生虚假目标的问题,能够减少虚假目标,提高多目标检测的可靠性。
为达到上述目的,本发明采用如下技术方案:
提供一种基于梯形LFMCW调制的雷达多目标检测方法,包括以下步骤:
步骤1,雷达向其检测范围内的目标发射梯形LFMCW信号,并接收目标回波信号;将所述目标回波信号与所述梯形LFMCW信号进行混频,得到混频信号;对所述混频信号经低通滤波器滤波后,再经模拟/数字A/D转换,得到数字差频信号;其中,所述数字差频信号包括上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down
步骤2,计算得到所述数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱,进而对所述数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱分别进行恒虚警检测,得到所述数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱中包含的谱峰,分别记为上扫频段谱峰、恒频段谱峰、下扫频段谱峰;
其中,第i个上扫频段谱峰记为pi up,i=1,2,…,M,pi up对应的频率为fi up、幅度为Ai up,M表示上扫频段差频Sb,up的频谱中包含的谱峰个数;第k个恒频谱峰记为pk const,k=1,2,…,P,pk const对应的频率为fk const、幅度为Ak const,P表示恒频段差频Sb,const的频谱中包含的谱峰个数;第j个下扫频谱峰记为pi down,j=1,2,…,N,pj down对应的频率为fj down、幅度为Aj down,N表示下扫频段差频Sb,down的频谱中包含的谱峰个数;
步骤3,将每个上扫频段谱峰分别与各上扫频段谱峰进行组合,得到M×N对两谱峰组合(pi up,pj down);计算每对两谱峰组合对应的频率和与归一化幅度差,并按照幅度差从小到大的顺序对所述M×N对两谱峰组合进行排序得到序列S;
步骤4,对于第k个恒频段谱峰,将其对应的频率与序列S中两谱峰组合对应的频率和进行对比,确定序列S中是否存在与所述第k个恒频段谱峰满足第一频率匹配条件的第一两谱峰组合;若存在,则将所述第一两谱峰组合与所述第k个恒频段谱峰进行匹配,得到对应的三谱峰组合,并将所述第一两谱峰组合从序列S中删除;若不存在,则确定本阶段第k个恒频段谱峰未成功匹配到两谱峰组合;其中,所述第一频率匹配条件为:
Figure BDA0001646485010000031
fres表示预设的频谱分辨率,k依次取1、2、3…P;
判断P个恒频段谱峰中是否存在未成功匹配到两谱峰组合的恒频段谱峰:若是,则转至步骤5;否则,转至步骤6;
步骤5,对于步骤4中未成功匹配到两谱峰组合的各恒频段谱峰,确定经步骤4更新后的序列S中与该恒频段谱峰满足第二频率匹配条件的第二两谱峰组合,将所述第二两谱峰组合与该恒频段谱峰进行匹配,得到对应的三谱峰组合,并将所述第二两谱峰组合从更新后的序列S中删除,转至步骤6;
其中,所述第二频率匹配条件为:
Figure BDA0001646485010000032
步骤6,判断P是否等于L,L表示雷达最多可分辨的目标数目,L=max(M,N,P);若P不等于L,则转至步骤7;若P等于L,则转至步骤8;
步骤7,对经步骤5更新后的序列S中的各对两谱峰组合,将其对应的频率和与各恒频段谱峰对应的频率分别进行对,确定P个恒频段谱峰中与该对两谱峰组合满足第二频率匹配条件的恒频段谱峰,将所述恒频段谱峰与该对两谱峰组合进行匹配,得到对应的三谱峰组合,转至步骤8;
步骤8,对于步骤4至步骤7得到的每对三谱峰组合,根据其中的上扫频段谱峰、恒频段谱峰以及下扫频段谱峰对应的频率,计算得到目标的距离估计值和速度估计值。
基于本发明实施例提供的基于梯形LFMCW调制的雷达多目标检测方法,在对谱峰进行配对时,利用差拍信号的频率和幅度信息,根据频率为主,幅度为辅的匹配原则,能够属于同一目标的谱峰进行正确配对,因此能够降低虚假目标产生的概率,有利于获取准确的目标速度和距离信息,具有较高的可靠性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种基于梯形LFMCW调制的雷达多目标检测方法流程图;
图2为本发明实施例中雷达发射的梯形LFMCW信号、雷达接收的目标回波信号以及对应的差频信号的时间-频率关系图;
图3为梯形LFMCW解算原理示意图;
图4为传统梯形LFMCW解算方法出现虚假目标的情况示意图,图中白色椭圆表示虚假目标,黑色椭圆表示真实目标;
图5为对仿真数据采用传统LFMCW算法得到的20dB信噪比下解算虚假目标个数图;
图6为对仿真数据采用本发明方法得到的20dB信噪比下解算虚假目标个数图;
图7为对仿真数据采用传统LFMCW算法得到信噪比与解算虚假目标个数关系图;
图8为对仿真数据采用本发明方法得到的得到信噪比与解算虚假目标个数关系图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1所示为本发明实施例提供的一种基于梯形LFMCW调制的雷达多目标检测方法的流程示意图;
如图1所示,本发明实施例提供的基于梯形LFMCW调制的雷达多目标检测方法,具体可以包括以下步骤:
步骤1,雷达向其检测范围内的目标发射梯形LFMCW信号,并接收目标回波信号;将目标回波信号与梯形LFMCW信号进行混频,得到混频信号;对混频信号经低通滤波器滤波后,再经模拟/数字A/D转换,得到数字差频信号。
其中,数字差频信号包括上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down
示意性的,雷达向其检测范围内的目标发射的梯形LFMCW信号、雷达接收的目标回波信号以及对应的差频信号的时间-频率关系图具体可参见图2。图2中上半部分中实线梯形信号为发射信号,虚线梯形信号为接收信号,图2中下半部分所示为差频信号。观察图2,梯形LFMCW信号包括上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down,调频带宽为B,每段信号持续时常分别为Tup、Tconst和Tdown,且Tup=Tdown=Tconst;类似的,目标回波信号对应包含上扫频段Sr,up、恒频段Sr,const和下扫频段Sr,down,目标距离引起的时间延迟为
Figure BDA0001646485010000051
目标运动速度v引起的多普勒频移为
Figure BDA0001646485010000052
其中C为光速,λ为雷达发射的梯形LFMCW信号的载波波长;差频信号在上扫频段、恒频段和下扫频段的频偏分别为fup、fconst、fdown
步骤2,计算得到数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱,进而对数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱分别进行恒虚警检测,得到数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱中包含的谱峰,分别记为上扫频段谱峰、恒频段谱峰、下扫频段谱峰。
其中,第i个上扫频段谱峰记为pi upi=1,2,…,M,pi up对应的频率为fi up、幅度为Ai up,M表示上扫频段差频Sb,up的频谱中包含的谱峰个数;第k个恒频谱峰记为pk const,k=1,2,…,P,pk const对应的频率为fk const、幅度为Ak const,P表示恒频段差频Sb,const的频谱中包含的谱峰个数;第j个下扫频谱峰记为pj down,j=1,2,…,N,pj down对应的频率为fj down、幅度为Aj down,N表示下扫频段差频Sb,down的频谱中包含的谱峰个数。
步骤3,将每个上扫频段谱峰分别与各上扫频段谱峰进行组合,得到M×N对两谱峰组合(pi up,pj down);计算每对两谱峰组合对应的频率和与归一化幅度差,并按照幅度差从小到大的顺序对M×N对两谱峰组合进行排序得到序列S。
具体的,步骤3中,计算每对两谱峰组合对应的频率和与归一化幅度差,具体包括:
根据第一预设公式:
Figure BDA0001646485010000061
计算得到两谱峰组合(pi up,pj down)对应的频率和;以及,
根据第二预设公式:
Figure BDA0001646485010000062
计算得到两谱峰组合(pi up,pj down)对应的归一化幅度差。
步骤4,对于第k个恒频段谱峰,按照频率优先原则,将其对应的频率与序列S中两谱峰组合对应的频率和进行对比,确定序列S中是否存在与第k个恒频段谱峰满足第一频率匹配条件的第一两谱峰组合;若存在,则将第一两谱峰组合与第k个恒频段谱峰进行匹配,得到对应的三谱峰组合,并将第一两谱峰组合从序列S中删除;若不存在,则确定本阶段第k个恒频段谱峰未成功匹配到两谱峰组合;判断P个恒频段谱峰中是否存在未成功匹配到两谱峰组合的恒频段谱峰:若是,则转至步骤5;否则,转至步骤6。
其中,第一频率匹配条件为:
Figure BDA0001646485010000063
fres表示预设的频谱分辨率,k依次取1、2、3…P。
步骤5,对于步骤4中未成功匹配到两谱峰组合的各恒频段谱峰,按照频率优先原则,确定经步骤4更新后的序列S中与该恒频段谱峰满足第二频率匹配条件的第二两谱峰组合,将第二两谱峰组合与该恒频段谱峰进行匹配,得到对应的三谱峰组合,并将第二两谱峰组合从更新后的序列S中删除,转至步骤6。
其中,第二频率匹配条件为:
Figure BDA0001646485010000064
步骤6,判断P是否等于L,L表示雷达最多可分辨的目标数目,L=max(M,N,P);若P不等于L,则转至步骤7;若P等于L,则转至步骤8。
步骤7,对经步骤4、步骤5更新后的序列S中的前L-P对两谱峰组合,按照幅度优先原则,将其对应的频率和与各恒频段谱峰对应的频率分别进行对比,确定P个恒频段谱峰中与该对两谱峰组合满足第二频率匹配条件的恒频段谱峰,将恒频段谱峰与该对两谱峰组合进行匹配,得到对应的L-P对三谱峰组合,转至步骤8。
需要说明的是,从步骤4至步骤7,将产生共L对三谱峰组合。具体来说,假设步骤4中P个恒频段谱峰均成功匹配到两谱峰组合,也即产生了P对三谱峰组合,则转至步骤6:若P等于L,三谱峰组合的数量即是L个;若P不等于L,则转至步骤7,步骤7中又产生L-P对三谱峰组合,一共是L对三谱峰组合;假设步骤4中P个恒频段谱峰中有Y个恒频段谱峰成功匹配到两谱峰组合,则仍有P-Y个恒频段谱峰未成功匹配到两谱峰组合,对于这P-Y个恒频段谱峰,经步骤5将产生P-Y对三谱峰组合,加起来仍然是P对三谱峰组合,再转至步骤6,与前述情况相同:若P等于L,三谱峰组合的数量即是L个;若P不等于L,则转至步骤7,步骤7中又产生L-P对三谱峰组合,一共是L对三谱峰组合。
步骤8,对于步骤4至步骤7得到的每对三谱峰组合,根据其中的上扫频段谱峰、恒频段谱峰以及下扫频段谱峰对应的频率,计算得到目标的距离估计值和速度估计值。
如前所述,步骤4至步骤7将得到的L对三谱峰组合,因此步骤8中是对每对三谱峰组合,根据其中的上扫频段谱峰、恒频段谱峰以及下扫频段谱峰对应的频率,计算得到对应的一个目标的距离估计值和速度估计值,最终将得到L个目标的距离估计值和速度估计值,L是雷达最多可分辨的目标数目。
步骤8中,根据上扫频段谱峰、恒频段谱峰以及下扫频段谱峰对应的频率,计算得到目标的距离估计值和速度估计值,具体包括:
根据恒频段对应的频率,计算得到目标的距离估计值
Figure BDA0001646485010000071
以及,
根据上扫频段谱峰和下扫频段谱峰对应的频率,计算得到目标的速度估计值
Figure BDA0001646485010000072
其中,Rl表示第l个目标的距离估计值,vl表示第l个目标的速度估计值,l=1,2,…L,fl,up、fl,down、fl,const分别为第l对三谱峰组合中上扫频谱峰、下扫频谱峰、恒频段谱峰对应的频率,C为光速,λ为雷达发射的LFMCW信号的载波波长。
至此,本发明实施例提供的基于梯形LFMCW调制的雷达多目标检测方法即结束。
为了清楚地说明本发明实施例提供的基于梯形LFMCW调制的雷达多目标检测方法的有益效果,以下先对传统LFMCW算法在多目标复杂场景下产生虚假目标的原因进行分析如下:
根据梯形信号目标解算理论,差频信号在上扫频段、恒频段和下扫频段的频偏分别可以表示为fup=fd-uτ、fdown=fd+uτ、fconst=fd,其中u表示调谐频率。反映在距离-速度坐标轴上,fup、fconst、fdown即为三条直线,而这三条直线的交点对应的速度和距离即为目标的信息,如图3所示。当在多目标复杂环境下,如图4所示,由于环境中有多个目标,对应的fup、fconst、fdown即为多条直线,而传统的LFMCW算法在对这些频偏进行配对时,会将不同目标对应的fup、fconst、fdown进行错误配对,这样就会产生虚假目标。
而基于本发明实施例提供的基于梯形LFMCW调制的雷达多目标检测方法,在对谱峰进行配对时,利用差拍信号的频率和幅度信息,根据频率为主,幅度为辅的匹配原则,能够属于同一目标的谱峰进行正确配对,因此能够降低虚假目标产生的概率,有利于获取准确的目标速度和距离信息,具有较高的可靠性。
以下通过仿真实验对本发明上述效果作进一步验证说明:
仿真实验1:本仿真实验中每次随机产生5个目标,在信噪比为20dB的条件下重复进行1000次蒙特卡洛实验,统计每次产生的虚假目标个数。
实验结果参见图5-图6,其中,图5为对仿真数据采用传统LFMCW算法得到的20dB信噪比下解算虚假目标个数图,图6为对仿真数据采用本发明方法得到的20dB信噪比下解算虚假目标个数图;图5和图6中,横轴表示信噪比,单位为dB,纵轴为解算虚假目标个数。
比较图5和图6,可以发现,经过本发明方法处理之后,虚假目标产生的次数及数量均减少。
仿真实验2:本仿真实验中每次随机产生5个目标,在不同的信噪比下重复进行1000次蒙特卡洛实验,统计总虚假目标个数。
实验结果参见图7-图8,其中图7为对仿真数据采用传统梯形LFMCW算法得到的信噪比与解算虚假目标个数关系图,图8为对仿真数据采用本发明方法得到的信噪比与解算虚假目标个数关系图。
比较图7和图8,可以发现,在不同信噪比下,本发明方法降低虚假目标出现概率的效果基本保持稳定,相比现有技术能够减少约50%的虚假目标。
综上,本发明实施例提供的基于梯形LFMCW调制的雷达多目标检测方法,根据LFMCW波形特点,通过分析同一目标在不同信号调制段的回波差拍频谱特点,结合频率和幅度信息确定信号频谱进行匹配原则,从而显著提高了频谱匹配的准确度,算法稳定有效,可解决多目标复杂环境下传统梯形信号处理方法易产生虚假目标的问题,能够提高多目标检测的可靠性。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (3)

1.一种基于梯形LFMCW调制的雷达多目标检测方法,其特征在于,包括以下步骤:
步骤1,雷达向其检测范围内的目标发射梯形LFMCW信号,并接收目标回波信号;将所述目标回波信号与所述梯形LFMCW信号进行混频,得到混频信号;对所述混频信号经低通滤波器滤波后,再经模拟/数字A/D转换,得到数字差频信号;其中,所述数字差频信号包括上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down
步骤2,计算得到所述数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱,进而对所述数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱分别进行恒虚警检测,得到所述数字差频信号的上扫频段差频Sb,up、恒频段差频Sb,const以及下扫频段差频Sb,down的频谱中包含的谱峰,分别记为上扫频段谱峰、恒频段谱峰、下扫频段谱峰;
其中,第i个上扫频段谱峰记为pi up,i=1,2,…,M,pi up对应的频率为fi up、幅度为Ai up,M表示上扫频段差频Sb,up的频谱中包含的谱峰个数;第k个恒频谱峰记为pk const,k=1,2,…,P,pk const对应的频率为fk const、幅度为Ak const,P表示恒频段差频Sb,const的频谱中包含的谱峰个数;第j个下扫频谱峰记为pj down,j=1,2,…,N,pj down对应的频率为fj down、幅度为Aj down,N表示下扫频段差频Sb,down的频谱中包含的谱峰个数;
步骤3,将每个上扫频段谱峰分别与各下扫频段谱峰进行组合,得到M×N对两谱峰组合(pi up,pj down);计算每对两谱峰组合对应的频率和与归一化幅度差,并按照幅度差从小到大的顺序对所述M×N对两谱峰组合进行排序得到序列S;
步骤4,对于第k个恒频段谱峰,按照频率优先原则,将其对应的频率与序列S中两谱峰组合对应的频率和进行对比,确定序列S中是否存在与所述第k个恒频段谱峰满足第一频率匹配条件的第一两谱峰组合;若存在,则将所述第一两谱峰组合与所述第k个恒频段谱峰进行匹配,得到对应的三谱峰组合,并将所述第一两谱峰组合从序列S中删除;若不存在,则确定本阶段第k个恒频段谱峰未成功匹配到两谱峰组合;其中,所述第一频率匹配条件为:
Figure FDA0003517921100000021
Figure FDA0003517921100000024
表示两谱峰组合(pi up,pj down)对应的频率和,fres表示预设的频谱分辨率,k依次取1、2、3…P;
判断P个恒频段谱峰中是否存在未成功匹配到两谱峰组合的恒频段谱峰:若是,则转至步骤5;否则,转至步骤6;
步骤5,对于步骤4中未成功匹配到两谱峰组合的各恒频段谱峰,按照频率优先原则,确定经步骤4更新后的序列S中与该恒频段谱峰满足第二频率匹配条件的第二两谱峰组合,将所述第二两谱峰组合与该恒频段谱峰进行匹配,得到对应的三谱峰组合,并将所述第二两谱峰组合从更新后的序列S中删除,转至步骤6;
其中,所述第二频率匹配条件为:
Figure FDA0003517921100000022
步骤6,判断P是否等于L,L表示雷达最多可分辨的目标数目,L=max(M,N,P);若P不等于L,则转至步骤7;若P等于L,则转至步骤8;
步骤7,对经步骤4、步骤5更新后的序列S中的前L-P对两谱峰组合,按照幅度优先原则,将对应的频率和与各恒频段谱峰对应的频率分别进行对比,确定P个恒频段谱峰中与该对两谱峰组合满足第二频率匹配条件的恒频段谱峰,将所述恒频段谱峰与该对两谱峰组合进行匹配,得到对应的L-P对三谱峰组合,转至步骤8;
步骤8,对于步骤4至步骤7得到的每对三谱峰组合,根据其中的上扫频段谱峰、恒频段谱峰以及下扫频段谱峰对应的频率,计算得到目标的距离估计值和速度估计值。
2.根据权利要求1所述的方法,其特征在于,步骤3中,所述计算每对两谱峰组合对应的频率和与归一化幅度差,包括:
根据第一预设公式:
Figure FDA0003517921100000025
计算得到两谱峰组合(pi up,pj down)对应的频率和;
根据第二预设公式:
Figure FDA0003517921100000023
计算得到两谱峰组合(pi up,pj down)对应的归一化幅度差;其中,Amax表示全部上扫频段谱峰和下扫频谱峰对应的归一化幅度中的最大值。
3.根据权利要求1所述的方法,其特征在于,步骤8中,根据上扫频段谱峰、恒频段谱峰以及下扫频段谱峰对应的频率,计算得到目标的距离估计值和速度估计值,包括:
根据恒频段对应的频率,计算得到目标的距离估计值
Figure FDA0003517921100000031
根据上扫频段谱峰和下扫频段谱峰对应的频率,计算得到目标的速度估计值
Figure FDA0003517921100000032
其中,Rl表示第l个目标的距离估计值,vl表示第l个目标的速度估计值,l=1,2,…L,fl,up、fl,down、fl,const分别为第l对三谱峰组合中上扫频谱峰、下扫频谱峰、恒频段谱峰对应的频率,B为调频带宽,T为每段信号持续时常,C为光速,λ为雷达发射的LFMCW信号的载波波长。
CN201810404328.8A 2018-04-28 2018-04-28 一种基于梯形lfmcw调制的雷达多目标检测方法 Active CN108415010B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810404328.8A CN108415010B (zh) 2018-04-28 2018-04-28 一种基于梯形lfmcw调制的雷达多目标检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810404328.8A CN108415010B (zh) 2018-04-28 2018-04-28 一种基于梯形lfmcw调制的雷达多目标检测方法

Publications (2)

Publication Number Publication Date
CN108415010A CN108415010A (zh) 2018-08-17
CN108415010B true CN108415010B (zh) 2022-04-19

Family

ID=63137394

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810404328.8A Active CN108415010B (zh) 2018-04-28 2018-04-28 一种基于梯形lfmcw调制的雷达多目标检测方法

Country Status (1)

Country Link
CN (1) CN108415010B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109471095B (zh) * 2018-11-06 2023-02-14 哈尔滨工程大学 一种基于快速迭代插值的fmcw雷达距离估计方法
CN109521417B (zh) * 2018-12-07 2023-01-03 哈尔滨工程大学 基于fmcw雷达波形的多目标检测计算方法
CN110988881B (zh) * 2019-10-31 2021-12-31 惠州市德赛西威智能交通技术研究院有限公司 一种车载毫米波雷达的目标检测方法及其车载雷达系统
CN110850400B (zh) * 2019-11-13 2022-10-14 南京理工大学 基于干涉仪测向的lfmcw雷达多目标检测方法
CN111103588B (zh) * 2019-12-31 2022-12-13 哈尔滨工程大学 一种利用信号能量的三角波多目标识别方法
CN111983627A (zh) * 2020-08-31 2020-11-24 宁波飞芯电子科技有限公司 一种探测装置及探测方法
CN112083405B (zh) * 2020-09-11 2023-11-14 广东工业大学 一种基于混合波形的目标检测方法及相关装置
CN112099038B (zh) * 2020-09-17 2023-11-10 上海波汇科技有限公司 一种基于fmcw激光雷达的多物体识别方法及装置
CN114690171A (zh) * 2020-12-31 2022-07-01 上海禾赛科技有限公司 利用调频连续波进行探测的方法和激光雷达
CN112882006B (zh) * 2021-01-20 2024-02-13 南京理工大学 基于复合调频连续波的毫米波近程目标探测方法
CN112859061B (zh) * 2021-03-12 2021-08-24 兰州理工大学 一种基于调频连续波雷达的多目标检测方法
CN115685099A (zh) * 2021-07-30 2023-02-03 北京万集科技股份有限公司 一种运动方向测量方法及激光雷达系统
CN114185046B (zh) * 2021-12-08 2024-08-16 西安电子科技大学 基于发射波形参数可调的视频sar成像系统及方法
CN114217301B (zh) * 2021-12-10 2022-08-05 中大智能科技股份有限公司 一种高精度边坡监测雷达目标检测及测距方法
CN114428255B (zh) * 2022-04-06 2022-07-08 北京摩尔芯光科技有限公司 用于激光雷达的距离速度确定方法和装置

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6320531B1 (en) * 1999-12-09 2001-11-20 Denso Corporation FM-CW radar system for measuring distance to and relative speed of target
CN103257346A (zh) * 2013-05-15 2013-08-21 桂林电子科技大学 一种汽车防撞雷达多目标探测方法与系统
CN105182341A (zh) * 2015-09-29 2015-12-23 西安知几天线技术有限公司 基于lfm三角波与恒频波组合波形的车辆防撞雷达多目标频率匹配方法
CN105549012A (zh) * 2015-12-02 2016-05-04 大连楼兰科技股份有限公司 车载毫米波雷达系统多目标检测装置
CN105629235A (zh) * 2015-12-29 2016-06-01 大连楼兰科技股份有限公司 多目标检测的组合波形汽车变道辅助系统的信号处理装置
CN106338727A (zh) * 2016-09-05 2017-01-18 南京理工大学 一种车载辅助驾驶雷达的目标检测方法
CN106842182A (zh) * 2016-08-31 2017-06-13 时艳玲 基于对称三角lfmcw雷达的多目标测速测距方法
CN107132534A (zh) * 2017-06-27 2017-09-05 西安电子科技大学 一种高速雷达目标频域检测的优化方法
CN107255814A (zh) * 2017-07-31 2017-10-17 西安电子科技大学 一种基于lfmsk波形的雷达目标检测方法
CN107783128A (zh) * 2016-08-25 2018-03-09 大连楼兰科技股份有限公司 基于毫米波雷达的固定翼无人机多目标防撞系统
CN107861117A (zh) * 2017-10-17 2018-03-30 电子科技大学 一种适用于连续波周界监视雷达的多目标参数测量方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6320531B1 (en) * 1999-12-09 2001-11-20 Denso Corporation FM-CW radar system for measuring distance to and relative speed of target
CN103257346A (zh) * 2013-05-15 2013-08-21 桂林电子科技大学 一种汽车防撞雷达多目标探测方法与系统
CN105182341A (zh) * 2015-09-29 2015-12-23 西安知几天线技术有限公司 基于lfm三角波与恒频波组合波形的车辆防撞雷达多目标频率匹配方法
CN105549012A (zh) * 2015-12-02 2016-05-04 大连楼兰科技股份有限公司 车载毫米波雷达系统多目标检测装置
CN105629235A (zh) * 2015-12-29 2016-06-01 大连楼兰科技股份有限公司 多目标检测的组合波形汽车变道辅助系统的信号处理装置
CN107783128A (zh) * 2016-08-25 2018-03-09 大连楼兰科技股份有限公司 基于毫米波雷达的固定翼无人机多目标防撞系统
CN106842182A (zh) * 2016-08-31 2017-06-13 时艳玲 基于对称三角lfmcw雷达的多目标测速测距方法
CN106338727A (zh) * 2016-09-05 2017-01-18 南京理工大学 一种车载辅助驾驶雷达的目标检测方法
CN107132534A (zh) * 2017-06-27 2017-09-05 西安电子科技大学 一种高速雷达目标频域检测的优化方法
CN107255814A (zh) * 2017-07-31 2017-10-17 西安电子科技大学 一种基于lfmsk波形的雷达目标检测方法
CN107861117A (zh) * 2017-10-17 2018-03-30 电子科技大学 一种适用于连续波周界监视雷达的多目标参数测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LFMCW radar multi-target acceleration and velocity estimation method;Rongquan Zhang et al.;《Proceedings 7th International Conference on Signal Processing》;20050627;全文 *
基于改进梯形波的连续波雷达多目标差频匹配算法*;田茂泰 等;《贵州师范大学学报( 自然科学版)》;20120930;第30卷(第5期);全文 *

Also Published As

Publication number Publication date
CN108415010A (zh) 2018-08-17

Similar Documents

Publication Publication Date Title
CN108415010B (zh) 一种基于梯形lfmcw调制的雷达多目标检测方法
CN111352102B (zh) 一种基于调频连续波雷达的多目标个数检测方法及装置
CN107861117B (zh) 一种适用于连续波周界监视雷达的多目标参数测量方法
CN106338727B (zh) 一种车载辅助驾驶雷达的目标检测方法
CN101738606B (zh) 基于广义多普勒滤波器组的雷达目标相参积累检测方法
CN109375202B (zh) 一种基于车载毫米波雷达的车辆测距测速方法
CN107255814B (zh) 一种基于lfmsk波形的雷达目标检测方法
CN112834992B (zh) 一种脉冲多普勒雷达的信号处理方法、装置及存储介质
CN107037432A (zh) 用于车辆的雷达装置及其确定目标的方法
CN112882006B (zh) 基于复合调频连续波的毫米波近程目标探测方法
CN110850400A (zh) 基于干涉仪测向的lfmcw雷达多目标检测方法
CN113009439A (zh) 一种基于调频连续波毫米波雷达的多目标探测方法
CN105487060A (zh) 一种双通道四斜率调制的多目标提取方法
CN110320509A (zh) 一种毫米波雷达目标检测方法及系统
CN111505631A (zh) 基于lfmcw雷达的心率估计算法
Kuptsov et al. Features of multi-target detection algorithm for automotive FMCW radar
Liu et al. A PELT-KCN algorithm for FMCW radar interference suppression based on signal reconstruction
CN110673130A (zh) 一种基于航迹关联的运动目标航迹跟踪方法
WO2022000333A1 (zh) 一种雷达探测方法及相关装置
CN114002654A (zh) 一种毫米波雷达解速度模糊方法
CN109085568B (zh) 一种基于二次混频的调频连续波多目标检测方法
CN110609263B (zh) 一种同时计算脉冲激光雷达目标回波时延和频偏的方法
CN108120975A (zh) 基于梯形连续波的雷达测速测距方法
CN110441749A (zh) 一种频率步进雷达目标运动参数估计方法
Willetts et al. Optimal time-frequency distribution selection for LPI radar pulse classification

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