CN105549002A - 一种基于组合波形的调频连续波雷达测量方法 - Google Patents

一种基于组合波形的调频连续波雷达测量方法 Download PDF

Info

Publication number
CN105549002A
CN105549002A CN201610073333.6A CN201610073333A CN105549002A CN 105549002 A CN105549002 A CN 105549002A CN 201610073333 A CN201610073333 A CN 201610073333A CN 105549002 A CN105549002 A CN 105549002A
Authority
CN
China
Prior art keywords
frequency
signal
target
speed
signal section
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
Application number
CN201610073333.6A
Other languages
English (en)
Other versions
CN105549002B (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.)
Xiamen University
Original Assignee
Xiamen University
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 Xiamen University filed Critical Xiamen University
Priority to CN201610073333.6A priority Critical patent/CN105549002B/zh
Publication of CN105549002A publication Critical patent/CN105549002A/zh
Application granted granted Critical
Publication of CN105549002B publication Critical patent/CN105549002B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

本发明一种基于组合波形的调频连续波雷达测量方法,通过发射调制信号为单频信号段、两个不同周期三角波信号段和若干个相同周期的锯齿波信号段组成的调频连续波信号,将回波信号与发射端本振信号去斜得到去斜回波信号,对去斜回波信号进行傅里叶变换,提取得到目标对应的频率值,该三角波信号段可以完成对目标的距离速度解耦得到目标的距离速度值,利用单频信号段得到的速度可以去除三角波信号段的大部分的虚假目标,从而实现真实目标的快速匹配提取,然后结合三角波信号段得到的无模糊粗速度和若干个锯齿波信号段得到的有模糊精速度,计算出相应的速度模糊周期数,进而求得目标的无模糊精速度,实现在低FFT点数下完成对目标速度的高精度测量。

Description

一种基于组合波形的调频连续波雷达测量方法
技术领域
本发明属于调频连续波雷达领域,涉及一种基于组合波形的调频连续波雷达测量方法。
背景技术
调频连续波测距测速的基本原理是发射信号频率按照某种调制信号变化,发射信号在遇到目标后产生滞后的回波信号,回波信号与发射端本振信号进行去斜处理得到去斜回波信号,该去斜回波信号和目标的距离速度存在对应关系,利用该去斜回波信号即可对目标实现测速测距。
目前调频连续波的调制方式主要有线性调频和非线性调频两大类。对于线性调频连续波,由于存在着严重的距离-速度耦合问题,通常采用三角波调制,利用正、负斜率调制区域下运动目标关于真实距离镜像的特点,完成对目标距离和速度的解耦运算。但是在多目标场景中,存在目标频率正确配对困难的问题。目前,解决该问题的常用方法是变周期三角波匹配法。对于实际工程中,测距测速的精度与频率估计的精度直接相关,特别是测速精度,受频率估计的精度影响很大。频率估计的精度受限于FFT的点数,FFT的点数越多,所得到的目标频率估计精度越高。但是在实际工程中,由于信号处理单元计算能力有限,难以承载大量的FFT运算;另外,如果直接采用变周期三角波匹配法进行目标关联匹配来提取真实目标,当目标个数较多的时候,计算量较大,耗时多。
发明内容
本发明的目的在于提供一种基于组合波形的调频连续波雷达测量方法,利用单频信号段可以极大地提高匹配速度,去除大部分虚假目标,并结合三角波信号段速度粗估计和锯齿波信号段速度精估计,可以提高在较低FFT点数条件下和多目标场景下的速度估计精度。
本发明一种基于组合波形的调频连续波雷达测量方法,包括如下步骤:
步骤1、雷达发射频率由组合波形调制的调频连续波信号:
该组合波形调制的调频连续波信号是周期性的,一个周期T内的组合波形主要由三段信号组成,包括一时间为Ta的单频信号段、一时间为Tb的三角波和时间为Tc的三角波组成的三角波信号段,以及若干个时间为TS的锯齿波组成的锯齿波信号段,其中,锯齿波信号的时间TS应为Tb或Tc的一半;
步骤2、发射信号在碰到目标后产生滞后的目标回波信号,利用目标回波信号与发射端的本振信号进行去斜处理,得到去斜回波信号:
该去斜回波信号与发射信号是一一对应的,一个周期T内的去斜回波信号包括:一时间为Ta的单频信号段去斜回波信号sa(t)、一时间为Tb三角波信号段上调频去斜回波信号sbu(t)和下调频去斜回波信号sbd(t)、一时间为Tc三角波信号段上调频去斜回波信号scu(t)和下调频去斜回波信号scd(t)和若干个时间为TS的锯齿波信号段去斜回波信号分别为ssd1(t)、ssd2(t)、ssd3(t)、ssd4(t);
步骤3、对单频信号段去斜回波信号sa(t)与三角波信号段去斜回波信号sbu(t)、sbd(t)、scu(t)和scd(t),分别进行傅里叶变换,得到单频信号段与三角波信号段去斜回波频域信号,然后分别对所得的各段去斜回波频域信号进行目标检测(如CFAR等),提取出目标所在谱峰对应的频率,得到各信号段目标频率集如下:时间Ta的单频信号段目标频率集Ca(f);时间Tb的三角波信号段目标频率集Cbu(f)、Cbd(f);时间Tc的三角波信号段目标频率集Ccu(f)、Ccd(f);
步骤4、对步骤3所得的时间Ta的单频信号段目标频率集Ca(f),利用单频回波信号的频率只包含多普勒频率的原理,根据频率与速度的换算公式,可以得到时间Ta的单频信号段所对应的单频信号段目标速度集Ca(v);对时间Tb的三角波信号段目标谱峰频率集Cbu(f)、Cbd(f),分别逐个取出两集合中的频率元素进行两两搭配组成频率对,利用三角波信号上调频段与下调频段可以联合完成对目标距离速度解耦的原理,根据频率对与目标距离、速度的换算公式,得到时间Tb的三角波信号段目标距离速度集Cb(r,v);同理,对时间Tc的三角波信号段执行同样的操作,得到时间Tc的三角波信号段目标距离速度集Cc(r,v);
步骤5、分别取出Ca(v)和Cb(r,v)的元素,以预置准则一对这两个元素中的速度进行关联判断,提取出Cb(r,v)中符合准则要求的所有元素,即为真实目标距离速度集U1(r,v);同理,对Ca(v)和Cc(r,v)以预置准则一进行关联提取,得到真实目标距离速度集U2(r,v);再对U1(r,v)与U2(r,v)以预置准则二进行关联判断,得到真实目标距离速度集U(r,v);
步骤6、根据步骤5关联得到的真实目标的距离速度集U(r,v),取出U(r,v)的元素,选取与锯齿波信号段斜率一致的三角波信号段去斜回波频域信号并计算该元素在该信号中所对应谱峰的距离单元PTcd
对步骤2所得的四个时间TS的锯齿波信号段去斜回波信号ssd1(t)、ssd2(t)、ssd3(t)、ssd4(t)分别做傅里叶变换,分别得到对应的锯齿波信号段去斜回波频域信号Ssd1(f)、Ssd2(f)、Ssd3(f)、Ssd4(f);
根据距离单元PTcd提取锯齿波信号段去斜回波频域信号Ssd1(f)、Ssd2(f)、Ssd3(f)、Ssd4(f)对应的距离单元的频谱值,然后对这些值做傅里叶变换,提取出谱峰位置,计算出该真实目标所对应的有模糊精速度va-pre
步骤7、利用步骤5所得到的真实目标距离速度集U(r,v)中的无模糊粗速度v和TS所对应的无模糊速度范围Vunam,可以求得该目标的速度模糊周期数N,利用求得N和Vunam对va-pre进行解模糊,从而得到该目标对应的无模糊精速度vpre
步骤8、对其他真实目标距离速度对重复进行步骤6和步骤7的处理,可以求得所有目标的无模糊精速度。
所述的步骤4具体包括:
步骤41、设共有个N目标,则时间Ta的单频信号段目标频率集Ca(f)含有N个真实目标所对应的N个频率,用fi代表Ca(f)里的频点编号,则Ca(f)如下:Ca(f)=[f1f2…fN],然后利用目标速度与频率的换算公式:其中λ为波长,求得单频信号段目标速度集Ca(v)如下:Ca(v)=[v1v2…vN];
步骤42、时间Tb的三角波信号段目标谱峰频率集Cbu(f)与Cbd(f)各自均含有N个目标所对应的频点,用fm代表Cbu(f)里的频点编号,fn代表Cbd(f)里的频点编号,利用目标距离和速度与频率的换算公式:求得目标距离与目标速度,其中T为对应的周期,B为带宽,c为光速,λ为波长;
步骤43、将Cbu(f)与Cbd(f)中的fm与fn依次配对代入公式计算,得到N2组(r,v),即时间Tb的三角波信号段目标距离速度集Cb(r,v):
步骤44、将Ccu(f)与Ccd(f)中的fm与fn依次配对代入公式计算得到N2组(r,v),即时间Tc的三角波信号段目标距离速度集Cc(r,v):
所述的步骤5具体包括:
步骤51、取出步骤4得的Ca(v)的元素Ca(vi),将其跟Cb(r,v)中的元素Cb(rmn,vmn)逐个以预置准则一进行关联判断计算,预置准则一如下:
计算速度误差Δv Δ v = ( v i - v m n ) 2 ,
其中,vi∈Ca(v),vmn∈Cb(r,v),i,m,n=1,2,…N;
步骤52、判断是否存在Δv≤ΔV-TH,若存在,则认为这两个元素是重合的,完成对同一目标的关联,元素的速度和距离即为真实目标的速度和距离;匹配后得到三角波周期Tb所对应真实目标距离速度集U1(r,v);同理,对时间Tc的三角波信号段目标距离速度集Cc(r,v)重复以上步骤,得到三角波周期Tc所对应真实目标距离速度集U2(r,v);
步骤53、分别逐个取出U1(r,v)、U2(r,v)中的元素U1(rmn,vmn)、U2(rm'n',vm'n')以预置准则二进行关联判断计算,预置准则二如下:
计算距离速度误差Δrv Δ r v = ( r m n - r m ′ n ′ ) 2 + ( v m n - v m ′ n ′ ) 2 ,
其中,rmn,vmn∈U1(r,v),rm'n',vm'n'∈U2(r,v),m,n,m',n'=1,2,…N;
步骤54、判断是否存在Δrv≤ΔRV-TH,若存在,则认为这两个元素是重合的,完成对同一目标的关联,元素的速度和距离即为真实目标的速度和距离,得到真实目标距离速度集U(r,v)。
本发明通过发射调制信号为单频信号段、两个不同周期三角波信号段和若干个相同周期的锯齿波信号段组成的调频连续波信号,回波信号与发射端本振信号进行去斜处理得到去斜回波信号,对去斜回波信号进行傅里叶变换,提取得到目标对应的频率估计值,该三角波信号段可以完成对目标的距离速度解耦得到目标的距离速度值,利用单频信号段得到的速度可以去除三角波信号段的大部分的虚假目标,从而实现真实目标的快速匹配提取,然后结合三角波信号段得到的无模糊粗速度和若干个锯齿波信号段得到的有模糊精速度,计算出相应的速度模糊周期数,进而求得目标的无模糊精速度。
本发明通过对调频连续波调制信号的改善设计,利用单频信号段得到的速度来去除大部分假目标,结合三角波信号段的速度的粗估计与锯齿波信号段的速度精估计,从而极大地提高了匹配速度,实现在低FFT点数下完成对目标速度的高精度测量。
附图说明
图1为已有的线性调频方式(以三角波为例)的发射信号调制信号;
图2为本发明组合波形调制信号的波形示意图;
图3为本发明的工作流程示意图;
图4为本发明的仿真结果示意。
以下结合附图和具体实施例对本发明作进一步详述。
具体实施方式
如图3所示,本发明一种基于组合波形的调频连续波雷达测量方法,包括如下步骤:
步骤1、雷达发射频率由组合波形调制的调频连续波信号:
该组合波形调制的调频连续波信号是周期性的,一个周期T内的组合波形主要由三段信号组成,包括一时间为Ta的单频信号段、一时间为Tb的三角波和时间为Tc的三角波组成的三角波信号段,以及若干个时间为TS的锯齿波组成的锯齿波信号段,如图2所示,其中,其中,锯齿波信号的时间TS应为Tb或Tc的一半;
步骤2、发射信号在碰到目标后产生滞后的目标回波信号,利用目标回波信号与发射端的本振信号进行去斜处理,得到去斜回波信号:
所述的去斜回波信号与发射信号是一一对应的,一个周期T内的去斜回波信号包括:一时间为Ta的单频信号段去斜回波信号sa(t)、一时间为Tb三角波信号段上调频去斜回波信号sbu(t)和下调频去斜回波信号sbd(t)、一时间为Tc三角波信号段上调频去斜回波信号scu(t)和下调频去斜回波信号scd(t)和若干个时间为TS的锯齿波信号段去斜回波信号分别为ssd1(t)、ssd2(t)、ssd3(t)、ssd4(t);
步骤3、对单频信号段去斜回波信号sa(t)与三角波信号段去斜回波信号sbu(t)、sbd(t)、scu(t)和scd(t),分别进行傅里叶变换,得到单频信号段与三角波信号段去斜回波频域信号,然后分别对所得的各段去斜回波频域信号进行目标检测,提取出目标所在谱峰对应的频率,得到各信号段目标频率集如下:时间Ta的单频信号段目标频率集Ca(f);时间Tb的三角波信号段目标频率集Cbu(f)、Cbd(f);时间Tc的三角波信号段目标频率集Ccu(f)、Ccd(f);
步骤4、对步骤3所得的时间Ta的单频信号段目标频率集Ca(f),利用单频回波信号的频率只包含多普勒频率的原理,根据频率与速度的换算公式,得到时间Ta的单频信号段所对应的单频信号段目标速度集Ca(v);对时间Tb的三角波信号段目标谱峰频率集Cbu(f)、Cbd(f),分别逐个取出两集合中的频率元素进行两两搭配组成频率对,利用三角波信号上调频段与下调频段可以联合完成对目标距离速度解耦的原理,根据频率对与目标距离、速度的换算公式,可以得到时间Tb的三角波信号段目标距离速度集Cb(r,v);同理,对时间Tc的三角波信号段执行同样的操作,可以得到时间Tc的三角波信号段目标距离速度集Cc(r,v),具体为:
步骤41、设共有个N目标,则时间Ta的单频信号段目标频率集Ca(f)含有N个真实目标所对应的N个频率,用fi代表Ca(f)里的频点编号,则Ca(f)如下:Ca(f)=[f1f2…fN]。然后利用目标速度与频率的换算公式:其中λ为波长,可求得单频信号段目标速度集Ca(v)如下:Ca(v)=[v1v2…vN];
步骤42、时间Tb的三角波信号段目标谱峰频率集Cbu(f)与Cbd(f)各自均含有N个目标所对应的频点,用fm代表Cbu(f)里的频点编号,fn代表Cbd(f)里的频点编号,利用目标距离和速度与频率的换算公式:求得目标距离与目标速度,其中T为对应的周期,B为带宽,c为光速,λ为波长;
步骤43、将Cbu(f)与Cbd(f)中的fm与fn依次配对代入公式计算,得到N2组(r,v),即时间Tb的三角波信号段目标距离速度集Cb(r,v):
步骤44、将Ccu(f)与Ccd(f)中的fm与fn依次配对代入公式计算得到N2组(r,v),即时间Tc的三角波信号段目标距离速度集Cc(r,v):
步骤5、分别取出Ca(v)和Cb(r,v)的元素,以预置准则一对这两个元素中的速度进行关联判断,提取出Cb(r,v)中符合准则要求的所有元素,即为真实目标距离速度集U1(r,v);同理,对Ca(v)和Cc(r,v)以预置准则一进行关联提取,得到真实目标距离速度集U2(r,v);再对U1(r,v)与U2(r,v)以预置准则二进行关联判断,得到真实目标距离速度集U(r,v),具体为:
步骤51、取出步骤4得的Ca(v)的元素Ca(vi),将其跟Cb(r,v)中的元素Cb(rmn,vmn)逐个以预置准则一进行关联判断计算,预置准则一如下:
计算速度误差Δv Δ v = ( v i - v m n ) 2 ,
其中,vi∈Ca(v),vmn∈Cb(r,v),i,m,n=1,2,…N;
步骤52、判断是否存在Δv≤ΔV-TH,若存在,则认为这两个元素是重合的,完成对同一目标的关联,元素的速度和距离即为真实目标的速度和距离;匹配后得到三角波周期Tb所对应真实目标距离速度集U1(r,v);同理,对时间Tc的三角波信号段目标距离速度集Cc(r,v)重复以上步骤,得到三角波周期Tc所对应真实目标距离速度集U2(r,v);
步骤53、分别逐个取出U1(r,v)、U2(r,v)中的元素U1(rmn,vmn)、U2(rm'n',vm'n')以预置准则二进行关联判断计算,预置准则二如下:
计算距离速度误差Δrv Δ r v = ( r m n - r m ′ n ′ ) 2 + ( v m n - v m ′ n ′ ) 2 ,
其中,rmn,vmn∈U1(r,v),rm'n',vm'n'∈U2(r,v),m,n,m',n'=1,2,…N;
步骤54、判断是否存在Δrv≤ΔRV-TH,若存在,则认为这两个元素是重合的,完成对同一目标的关联,元素的速度和距离即为真实目标的速度和距离,得到真实目标距离速度集U(r,v);
步骤6、根据步骤5关联得到的真实目标的距离速度集U(r,v),取出U(r,v)的元素,选取与锯齿波信号段斜率一致的三角波信号段去斜回波频域信号并计算该元素在该信号中所对应谱峰的距离单元PTcd
对步骤2所得的四个时间TS的锯齿波信号段去斜回波信号ssd1(t)、ssd2(t)、ssd3(t)、ssd4(t)分别做傅里叶变换,分别得到对应的锯齿波信号段去斜回波频域信号Ssd1(f)、Ssd2(f)、Ssd3(f)、Ssd4(f);
根据距离单元PTcd提取锯齿波信号段去斜回波频域信号Ssd1(f)、Ssd2(f)、Ssd3(f)、Ssd4(f)对应的距离单元的频谱值,然后对这些值做傅里叶变换,提取出谱峰位置,计算出该真实目标所对应的有模糊精速度va-pre
步骤7、利用步骤5所得到的真实目标的距离速度集U(r,v)中的无模糊粗速度v和TS所对应的无模糊速度范围Vunam,可以求得该目标的速度模糊周期数N,利用求得N和Vunam对va-pre进行解模糊,从而得到该目标对应的无模糊精速度vpre
步骤8、对其他真实目标距离速度对重复进行步骤6和步骤7的处理,可以求得所有目标的无模糊精速度。
本发明主要是利用单频信号段来去除大部分虚假目标,并通过三角波信号段与锯齿波信号段的速度粗、精估计方法相结合,来提高目标匹配速度和速度估计精度。为了便于说明,这里设雷达载频fc=24GHz,B=100MHz,三个调制周期分别为TTa=2ms,TTb=4ms,TTc=2ms,TS=2ms,采样频率fs=0.15MHz,FFT点数NFFT=512,雷达静止,假设目标个数N=1,其中目标1:初始距离R1=50m,初始速度V1=10m/s,加速度a1=0.1m/s2;仿真100个周期。仿真结果如图4所示。
以上所述仅为本发明的较佳实施用例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换以及改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于组合波形的调频连续波雷达测量方法,其特征在于包括如下步骤:
步骤1、雷达发射频率由组合波形调制的调频连续波信号:
该组合波形调制的调频连续波信号是周期性的,一个周期T内的组合波形主要由三段信号组成,包括一时间为Ta的单频信号段、一时间为Tb的三角波和时间为Tc的三角波组成的三角波信号段,以及若干个时间为TS的锯齿波组成的锯齿波信号段,其中,锯齿波信号的时间TS应为Tb或Tc的一半;
步骤2、发射信号在碰到目标后产生滞后的目标回波信号,利用目标回波信号与发射端的本振信号进行去斜处理,得到去斜回波信号:
该去斜回波信号与发射信号是一一对应的,一个周期T内的去斜回波信号包括:一时间为Ta的单频信号段去斜回波信号sa(t)、一时间为Tb三角波信号段上调频去斜回波信号sbu(t)和下调频去斜回波信号sbd(t)、一时间为Tc三角波信号段上调频去斜回波信号scu(t)和下调频去斜回波信号scd(t)和若干个时间为TS的锯齿波信号段去斜回波信号分别为ssd1(t)、ssd2(t)、ssd3(t)、ssd4(t);
步骤3、对单频信号段去斜回波信号sa(t)与三角波信号段去斜回波信号sbu(t)、sbd(t)、scu(t)和scd(t),分别进行傅里叶变换,得到单频信号段与三角波信号段去斜回波频域信号,然后分别对所得的各段去斜回波频域信号进行目标检测,提取出目标所在谱峰对应的频率,得到各信号段目标频率集如下:时间Ta的单频信号段目标频率集Ca(f);时间Tb的三角波信号段目标频率集Cbu(f)、Cbd(f);时间Tc的三角波信号段目标频率集Ccu(f)、Ccd(f);
步骤4、对步骤3所得的时间Ta的单频信号段目标频率集Ca(f),利用单频回波信号的频率只包含多普勒频率的原理,根据频率与速度的换算公式,得到时间Ta的单频信号段所对应的单频信号段目标速度集Ca(v);对时间Tb的三角波信号段目标谱峰频率集Cbu(f)、Cbd(f),分别逐个取出两集合中的频率元素进行两两搭配组成频率对,利用三角波信号上调频段与下调频段可以联合完成对目标距离速度解耦的原理,根据频率对与目标距离、速度的换算公式,得到时间Tb的三角波信号段目标距离速度集Cb(r,v);同理,对时间Tc的三角波信号段执行同样的操作,得到时间Tc的三角波信号段目标距离速度集Cc(r,v);
步骤5、分别取出Ca(v)和Cb(r,v)的元素,以预置准则一对这两个元素中的速度进行关联判断,提取出Cb(r,v)中符合准则要求的所有元素,即为真实目标距离速度集U1(r,v);同理,对Ca(v)和Cc(r,v)以预置准则一进行关联提取,得到真实目标距离速度集U2(r,v);再对U1(r,v)与U2(r,v)以预置准则二进行关联判断,得到真实目标距离速度集U(r,v);
步骤6、根据步骤5关联得到的真实目标的距离速度集U(r,v),取出U(r,v)的元素,选取与锯齿波信号段斜率一致的三角波信号段去斜回波频域信号并计算该元素在该信号中所对应谱峰的距离单元PTcd
对步骤2所得的四个时间TS的锯齿波信号段去斜回波信号ssd1(t)、ssd2(t)、ssd3(t)、ssd4(t)分别做傅里叶变换,分别得到对应的锯齿波信号段去斜回波频域信号Ssd1(f)、Ssd2(f)、Ssd3(f)、Ssd4(f);
根据距离单元PTcd提取锯齿波信号段去斜回波频域信号Ssd1(f)、Ssd2(f)、Ssd3(f)、Ssd4(f)对应的距离单元的频谱值,然后对这些值做傅里叶变换,提取出谱峰位置,计算出该真实目标所对应的有模糊精速度va-pre
步骤7、利用步骤5所得到的真实目标距离速度集U(r,v)中的无模糊粗速度v和单个锯齿波的时间TS所对应的无模糊速度范围Vunam,求得该目标的速度模糊周期数N,利用求得N和Vunam对va-pre进行解模糊,从而得到该目标对应的无模糊精速度vpre
步骤8、对其他真实目标距离速度对重复进行步骤6和步骤7的处理,可以求得所有目标的无模糊精速度。
2.根据权利要求1所述的一种基于组合波形的调频连续波雷达测量方法,其特征在于所述的步骤4具体包括:
步骤41、设共有个N目标,则时间Ta的单频信号段目标频率集Ca(f)含有N个真实目标所对应的N个频率,用fi代表Ca(f)里的频点编号,则Ca(f)如下:Ca(f)=[f1f2…fN],然后利用目标速度与频率的换算公式:其中λ为波长,求得单频信号段目标速度集Ca(v)如下:Ca(v)=[v1v2…vN];
步骤42、时间Tb的三角波信号段目标谱峰频率集Cbu(f)与Cbd(f)各自均含有N个目标所对应的频点,用fm代表Cbu(f)里的频点编号,fn代表Cbd(f)里的频点编号,利用目标距离和速度与频率的换算公式:求得目标距离与目标速度,其中T为对应的周期,B为带宽,c为光速,λ为波长;
步骤43、将Cbu(f)与Cbd(f)中的fm与fn依次配对代入公式计算,得到N2组(r,v),即时间Tb的三角波信号段目标距离速度集Cb(r,v):
步骤44、将Ccu(f)与Ccd(f)中的fm与fn依次配对代入公式计算得到N2组(r,v),即时间Tc的三角波信号段目标距离速度集Cc(r,v):
3.根据权利要求1所述的一种基于组合波形的调频连续波雷达测量方法,其特征在于所述的步骤5具体包括:
步骤51、取出步骤4得的Ca(v)的元素Ca(vi),将其跟Cb(r,v)中的元素Cb(rmn,vmn)逐个以预置准则一进行关联判断计算,预置准则一如下:
计算速度误差Δv
其中,vi∈Ca(v),vmn∈Cb(r,v),i,m,n=1,2,…N;
步骤52、判断是否存在Δv≤ΔV-TH,若存在,则认为这两个元素是重合的,完成对同一目标的关联,元素的速度和距离即为真实目标的速度和距离;匹配后得到三角波周期Tb所对应真实目标距离速度集U1(r,v);同理,对时间Tc的三角波信号段目标距离速度集Cc(r,v)重复以上步骤,得到三角波周期Tc所对应真实目标距离速度集U2(r,v);
步骤53、分别逐个取出U1(r,v)、U2(r,v)中的元素U1(rmn,vmn)、U2(rm'n',vm'n')以预置准则二进行关联判断计算,预置准则二如下:
计算距离速度误差Δrv Δ rv = ( r mn - r m ′ n ′ ) 2 + ( v mn - v m ′ n ′ ) 2 ,
其中,rmn,vmn∈U1(r,v),rm'n',vm'n'∈U2(r,v),m,n,m',n'=1,2,…N;
步骤54、判断是否存在Δrv≤ΔRV-TH,若存在,则认为这两个元素是重合的,完成对同一目标的关联,元素的速度和距离即为真实目标的速度和距离,得到真实目标距离速度集U(r,v)。
CN201610073333.6A 2016-02-02 2016-02-02 一种基于组合波形的调频连续波雷达测量方法 Active CN105549002B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610073333.6A CN105549002B (zh) 2016-02-02 2016-02-02 一种基于组合波形的调频连续波雷达测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610073333.6A CN105549002B (zh) 2016-02-02 2016-02-02 一种基于组合波形的调频连续波雷达测量方法

Publications (2)

Publication Number Publication Date
CN105549002A true CN105549002A (zh) 2016-05-04
CN105549002B CN105549002B (zh) 2018-04-20

Family

ID=55828311

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610073333.6A Active CN105549002B (zh) 2016-02-02 2016-02-02 一种基于组合波形的调频连续波雷达测量方法

Country Status (1)

Country Link
CN (1) CN105549002B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106443660A (zh) * 2016-09-09 2017-02-22 电子科技大学 一种线性调频连续波雷达多运动目标匹配方法
CN109039497A (zh) * 2017-06-09 2018-12-18 中兴通讯股份有限公司 单频信号组合分配方法、数据发射方法、接收方法及设备
CN110109046A (zh) * 2019-04-22 2019-08-09 天津大学 一种水声信标的示位方法及装置
CN110865350A (zh) * 2019-11-21 2020-03-06 清华大学 基于双波段三角lfmw雷达的多目标测距测速方法及系统
CN111220975A (zh) * 2020-01-15 2020-06-02 深圳大学 多目标检测方法、装置、设备及存储介质
CN111562407A (zh) * 2020-04-26 2020-08-21 武汉拓宝科技股份有限公司 一种非接触式行驶车辆加速度测量方法
CN111722212A (zh) * 2020-06-29 2020-09-29 江苏集萃深度感知技术研究所有限公司 基于锯齿波及单频率信号的雷达测速方法
WO2021170133A1 (zh) * 2020-02-28 2021-09-02 加特兰微电子科技(上海)有限公司 提升目标探测精度的方法、集成电路、无线电器件及电子设备
CN117991282A (zh) * 2024-04-03 2024-05-07 厦门大学 FMCW-Lidar距离和速度的估计方法、目标检测装置和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103630888A (zh) * 2013-02-27 2014-03-12 中国科学院电子学研究所 基于对称三角lfmcw雷达的高精度实时微波测速测距装置
CN104155649A (zh) * 2014-09-02 2014-11-19 中国电子科技集团公司第五十四研究所 一种基于三周期调频连续波相参雷达距离速度去耦合方法
CN105182341A (zh) * 2015-09-29 2015-12-23 西安知几天线技术有限公司 基于lfm三角波与恒频波组合波形的车辆防撞雷达多目标频率匹配方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103630888A (zh) * 2013-02-27 2014-03-12 中国科学院电子学研究所 基于对称三角lfmcw雷达的高精度实时微波测速测距装置
CN104155649A (zh) * 2014-09-02 2014-11-19 中国电子科技集团公司第五十四研究所 一种基于三周期调频连续波相参雷达距离速度去耦合方法
CN105182341A (zh) * 2015-09-29 2015-12-23 西安知几天线技术有限公司 基于lfm三角波与恒频波组合波形的车辆防撞雷达多目标频率匹配方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
丁顺宝 等: "一种LFMCW雷达多目标距离-速度配对新方法", 《现代雷达》 *
孙昱 等: "汽车防撞雷达系统设计", 《军事交通学院学报》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106443660A (zh) * 2016-09-09 2017-02-22 电子科技大学 一种线性调频连续波雷达多运动目标匹配方法
CN106443660B (zh) * 2016-09-09 2019-03-01 电子科技大学 一种线性调频连续波雷达多运动目标匹配方法
CN109039497A (zh) * 2017-06-09 2018-12-18 中兴通讯股份有限公司 单频信号组合分配方法、数据发射方法、接收方法及设备
CN109039497B (zh) * 2017-06-09 2021-08-17 中兴通讯股份有限公司 单频信号组合分配方法、数据发射方法、接收方法及设备
CN110109046A (zh) * 2019-04-22 2019-08-09 天津大学 一种水声信标的示位方法及装置
CN110865350A (zh) * 2019-11-21 2020-03-06 清华大学 基于双波段三角lfmw雷达的多目标测距测速方法及系统
CN110865350B (zh) * 2019-11-21 2022-04-01 清华大学 基于双波段三角lfmw雷达的多目标测距测速方法及系统
CN111220975A (zh) * 2020-01-15 2020-06-02 深圳大学 多目标检测方法、装置、设备及存储介质
WO2021170133A1 (zh) * 2020-02-28 2021-09-02 加特兰微电子科技(上海)有限公司 提升目标探测精度的方法、集成电路、无线电器件及电子设备
CN111562407A (zh) * 2020-04-26 2020-08-21 武汉拓宝科技股份有限公司 一种非接触式行驶车辆加速度测量方法
CN111722212A (zh) * 2020-06-29 2020-09-29 江苏集萃深度感知技术研究所有限公司 基于锯齿波及单频率信号的雷达测速方法
CN117991282A (zh) * 2024-04-03 2024-05-07 厦门大学 FMCW-Lidar距离和速度的估计方法、目标检测装置和存储介质

Also Published As

Publication number Publication date
CN105549002B (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN105549002A (zh) 一种基于组合波形的调频连续波雷达测量方法
CN102628937B (zh) 基于广义keystone变换和非相参积累的雷达检测方法
CN101089653B (zh) 近程调频连续波fmcw雷达抗干扰方法
CN1932551B (zh) 雷达
CN102162848B (zh) 基于伪随机编码与线性调频连续波的汽车主动防撞雷达
CN105738889B (zh) 一种调频连续波测速测距方法
CN105093227B (zh) 一种交通流量测量装置及车辆运行信息获得方法
CN107688178A (zh) 一种基于77GHz毫米波雷达的锯齿波测距测速方法
CN105487060B (zh) 一种双通道四斜率调制的多目标提取方法
CN105425224A (zh) 车载毫米波雷达系统多目标个数获取方法及装置
CN103529444A (zh) 一种车载毫米波雷达动目标识别器及识别方法
CN101980046A (zh) 调频步进雷达复合测速运动补偿方法
CN106842166A (zh) 一种适用于线性调频连续波雷达系统的解速度模糊方法
CN107843892A (zh) 一种基于最小二乘法的高速目标多普勒测速方法
CN105301591A (zh) 一种道路交通监控雷达及其实现方法
CN105093215A (zh) 基于多普勒信息的雷达对低空慢速小目标的跟踪方法
CN102798866A (zh) 激光雷达系统及采用该系统的正弦调幅-脉冲相位编码调制的复合测距和测速方法
CN106170713B (zh) 车载用雷达装置
CN105549001A (zh) 车载毫米波雷达系统多目标检测方法
CN108693524A (zh) 基于线性调频连续波雷达的多运动目标匹配方法及其系统
CN102650689A (zh) 一种步进频率脉冲雷达位移测量方法
CN104991247B (zh) 一种低截获测速方法及雷达装置
CN105549012A (zh) 车载毫米波雷达系统多目标检测装置
CN102288942A (zh) 一种毫米波雷达信号波形的设计方法
CN203012135U (zh) 一种调频连续波雷达系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant