CN106405541A - 全相参连续波多普勒雷达及其测距测速方法 - Google Patents

全相参连续波多普勒雷达及其测距测速方法 Download PDF

Info

Publication number
CN106405541A
CN106405541A CN201610997913.4A CN201610997913A CN106405541A CN 106405541 A CN106405541 A CN 106405541A CN 201610997913 A CN201610997913 A CN 201610997913A CN 106405541 A CN106405541 A CN 106405541A
Authority
CN
China
Prior art keywords
target
frequency
distance
signal
doppler
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
CN201610997913.4A
Other languages
English (en)
Other versions
CN106405541B (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.)
Shandong Yinghe Road Freight Transport Co ltd
Original Assignee
Suzhou Tushi Electronic Technology 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 Suzhou Tushi Electronic Technology Co Ltd filed Critical Suzhou Tushi Electronic Technology Co Ltd
Priority to CN201610997913.4A priority Critical patent/CN106405541B/zh
Publication of CN106405541A publication Critical patent/CN106405541A/zh
Application granted granted Critical
Publication of CN106405541B publication Critical patent/CN106405541B/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/87Combinations of radar systems, e.g. primary radar and secondary radar

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+方位维MTD的二维联合检测方式,利用多个相参的高重频调制信号独立获得无模糊的径向速度信息,然后再将速度引入的多普勒信息在距离维频域信息中去除,消除距离‑速度耦合的影响,从而进一步得到准确的距离信息。本发明同时具备相参雷达测速和连续波雷达无距离盲区测距的特点,满足无距离盲区、同时测距测速、实时性、小型化、低功耗的要求。

Description

全相参连续波多普勒雷达及其测距测速方法
技术领域
本发明涉及连续波测量雷达领域,具体地说,是一种全相参连续波多普勒雷达及其测距测速方法。
背景技术
传统的可同时测距测速的全相参雷达采用脉冲多普勒体制,通过发射相干脉冲串测量目标反射信号与发射脉冲的时间差确定目标的距离,脉冲间采用FFT分析回波信号的多普勒频谱获得目标的速度信息,从而具有距离与速度的两维测量能力。但是该体制存在距离盲区,并且分辨率要求较高的情况下,脉冲体制需要同时兼顾窄脉冲和大带宽的要求,从而加大了信号产生的难度,并且后端信号处理需要高采样率,这大大增加了系统的复杂度。
传统的连续波雷达不存在距离盲区,一方面可以采用相参体制,通过对雷达接收回波信号进行频谱分析,可以完成测速功能但是不具备测距能力;另一方面可以采用非相参体制,通过发射线性射调频信号,利用接收频差信号的特性与目标距离映射关系提取距离信息,具有极高距离分辨率的同时大幅降低信号处理带宽以及不存在距离盲区的优势,但是由于采用的是非相参体制,在测距的同时没有多普勒信息而无法测速,并且其速度带来的多普勒信息耦合到信号中,将导致距离测量误差较大。
目前常用的三角调频连续波可以同时测距和测速,通过正负扫频完成距离-速度解耦,但是需要经过频域配对。在复杂的高分辨多目标环境背景下,频域精确配对难以实现,容易导致目标信息测量错误。并且三角波调制法在相同距离分辨率下会使中频带宽增大一倍,从而增加了中频采样率和信号处理运算量,难以满足系统实时性要求。
因此需要一种将全相参脉冲多普勒雷达和连续波雷达相结合,可以无距离盲区、同时测距测速并且可以有效解决距离-速度耦合的全相参连续波多普勒雷达。
发明内容
1、本发明的目的。
本发明将全相参脉冲多普勒雷达和连续波雷达相结合,提出了一种全相参连续波多普勒雷达及其测距测速方法,具备同时测距测速功能,并且可以消除距离-速度耦合的影响,满足无距离盲区测量、实时性、小型化、低功耗的要求。
2、本发明所采用的技术方案。
本发明提出的全相参连续波多普勒雷达测距测速方法,按照如下步骤进行:
步骤1、根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱;
步骤2、根据距离向频谱计算目标方位向的频谱,获取多普勒频率和目标径向速度,即利用多个相参的高重频调制信号独立获得无模糊的径向速度信息;
步骤3、在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离,即将速度引入的多普勒信息在距离维频域信息中去除,消除距离-速度耦合的影响,从而进一步得到距离信息。
优选的,所述的步骤1中根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱:
1.1目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号;
1.2对差频信号进行距离FFT变换,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合。
优选的,所述的步骤2中根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率和目标径向速度具体为:对距离向频谱进行方位向FFT变换,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标,其频谱峰值受到多普勒频率调制。
优选的,所述的步骤1中具体为:
设第K个调频周期内发射信号为S(k,t)可以表示为:
S(k,t)=A0cos{2π[f0(t-kTr)+0.5kr(t-kTr)2]+jφ0} (1)
其中,t为信号时间,Kr为调频率,Tr为调频周期,f0为载波频率,φ0为初始相位;
则对于t=0时刻目标距离为R0,径向速度为v的回波信号可以表示为:
S(k,t)=A1cos{2π[f0(t-kTr-τ(t))+0.5kr(t-kTr-τ(t))2]+jφ1} (2)
其中,点目标回波延时τ(t)=2(R0+vt)/c;
目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号,则第k个CPI周期的差拍信号表示为:
其中,v为目标径向速度,φ为固定相位;
对差频信号进行距离向FFT后的频谱可以表示为:
可以得到,S(k,fr)峰值点处的中心频率fc=fb+fd=2krR0/c+2vf0/c,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合。
优选的,所述的步骤2中:
根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率fd和目标径向速度:
假设方位维MTD处理所需调频周期数为M,对距离向频谱进行方位向FFT后的频谱可以表示为:
从上式可以看出,对于多周期回波差拍信号频谱峰值的变化,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标来说,回波差拍信号频谱峰值只有exp{j4πkvf0Tr/c}项是变化的,即频谱峰值是受到多普勒频率fd=2vf0/c调制的。
优选的,所述的步骤3中:
在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离和径向速度:
设经过MTD处理后,在两维频谱幅度最高的谱线对应的单元号为(m,n),第m个速度通道对应的目标多普勒频率和第n个距离通道对应的中心频率分别为:
fd=(m-M/2-1)fr/M (6)
fc=(n-1)fs/N (7)
其中,N为距离维FFT点数,fs为距离向采样率;
则该通道对应的速度可以表示为vt=cfd/2f0,从而可以获得每一速度通道的实际距离对应的频率为fb=fc-fd,根据频率和距离的对应关系,进而可以得到目标的实际距离为Rt=cfb/2kr
一种全相参连续波多普勒测距测速雷达,包括:
用于根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱的装置;
用于根据距离向频谱计算目标方位向的频谱,获取多普勒频率和目标径向速度的装置,即利用多个相参的高重频调制信号独立获得无模糊的径向速度信息;
用于在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离的装置,即将速度引入的多普勒信息在距离维频域信息中去除,消除距离-速度耦合的影响,从而进一步得到距离信息。
优选的,所述的计算目标差拍回波信号的距离向频谱的装置包括:
用于目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号的装置
用于对差频信号进行距离FFT变换,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合的装置。
优选的,所述的根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率和目标径向速度的装置具体为:对距离向频谱进行方位向FFT变换,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标,其频谱峰值受到多普勒频率调制。
优选的,用于根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱的装置具体为:
S(k,t)=A0cos{2π[f0(t-kTr)+0.5kr(t-kTr)2]+jφ0} (1)
其中,t为信号时间,Kr为调频率,Tr为调频周期,f0为载波频率,φ0为初始相位;
则对于t=0时刻目标距离为R0,径向速度为v的回波信号可以表示为:
S(k,t)=A1cos{2π[f0(t-kTr-τ(t))+0.5kr(t-kTr-τ(t))2]+jφ1} (2)
其中,点目标回波延时τ(t)=2(R0+vt)/c;
目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号,则第k个CPI周期的差拍信号表示为:
其中,v为目标径向速度,φ为固定相位;
对差频信号进行距离向FFT后的频谱可以表示为:
可以得到,S(k,fr)峰值点处的中心频率fc=fb+fd=2krR0/c+2vf0/c,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合。
用于根据距离向频谱计算目标方位向的频谱,获取多普勒频率和目标径向速度,即利用多个相参的高重频调制信号独立获得无模糊的径向速度信息的装置具体为:
设第K个调频周期内发射信号为S(k,t)可以表示为:
根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率fd和目标径向速度:
假设方位维MTD处理所需调频周期数为M,对距离向频谱进行方位向FFT后的频谱可以表示为:
从上式可以看出,对于多周期回波差拍信号频谱峰值的变化,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标来说,回波差拍信号频谱峰值只有exp{j4πkvf0Tr/c}项是变化的,即频谱峰值是受到多普勒频率fd=2vf0/c调制的。
用于在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离,即将速度引入的多普勒信息在距离维频域信息中去除,消除距离-速度耦合的影响,从而进一步得到距离信息的装置具体为:
在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离和径向速度:
设经过MTD处理后,在两维频谱幅度最高的谱线对应的单元号为(m,n),第m个速度通道对应的目标多普勒频率和第n个距离通道对应的中心频率分别为:
fd=(m-M/2-1)fr/M (6)
fc=(n-1)fs/N (7)
其中,N为距离维FFT点数,fs为距离向采样率;
则该通道对应的速度可以表示为vt=cfd/2f0,从而可以获得每一速度通道的实际距离对应的频率为fb=fc-fd,根据频率和距离的对应关系,进而可以得到目标的实际距离为Rt=cfb/2kr
3、本发明的有益效果。
(1)本发明采用距离维FFT和方位维MTD的二维联合检测方式,相比于传统的测距测速雷达,同时具备相参雷达测速和连续波雷达无距离盲区测距的特点;
(2)本发明充分考虑了由于距离-速度耦合导致的平台与目标间相对距离速度的调制特性,能够综合全相参脉冲体制和连续波体制的特点,在方位维采用多调频周期的相参锯齿波进行MTD处理独立获取无模糊速度信息,并消除差频信号中速度引入的多普勒频率在距离频率上的影响,从而进一步得到准确的距离信息。
附图说明
图1为本发明的脉冲信号和一次调频信号示意图。
图2为本发明的全相参准连续波发射信号、接收回波和差频信号关系示意图。
图3为距离维FFT和方位维MTD处理后对应的二维数据结构示意图。
图4为本发明方法计算的二维频谱示意图。
具体实施方式
实施例
本发明提供了一种全相参连续波多普勒雷达及其测距测速方法,下面结合附图并举实施例,对本发明进行详细描述。
本发明采用全相参准锯齿连续波发射波形,综合了脉冲体制和连续波体制的特点,每一次锯齿的发射、接收和数据采集受脉冲重复信号控制,获取多个调频周期的全相参准连续波信号,其信号形式如图1和图2所示。在该技术体制下,首先获得回波信号与发射信号混频得到的差频复信号,并根据距离频率-多普勒频率与目标距离-速度的调制特性,采用距离维FFT+方位维MTD的二维联合检测方式,如图3所示,利用多个相参的高重频调制信号独立获得无模糊的径向速度信息,然后再将速度引入的多普勒信息在距离维频域信息中去除,消除距离-速度耦合的影响,从而进一步得到准确的距离信息。具体包括:
步骤一:根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱:
设第K个调频周期内发射信号为S(k,t)可以表示为:
S(k,t)=A0cos{2π[f0(t-kTr)+0.5kr(t-kTr)2]+jφ0} (1)
其中,t为信号时间,Kr为调频率,Tr为调频周期,f0为载波频率,φ0为初始相位。
则对于t=0时刻目标距离为R0,径向速度为v的回波信号可以表示为:
S(k,t)=A1cos{2π[f0(t-kTr-τ(t))+0.5kr(t-kTr-τ(t))2]+jφ1} (2)
其中,点目标回波延时τ(t)=2(R0+vt)/c。
目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号,则第k个CPI周期的差拍信号表示为:
其中,v为目标径向速度,φ为固定相位。
对差频信号进行距离向FFT后的频谱可以表示为:
可以得到,S(k,fr)峰值点处的中心频率fc=fb+fd=2krR0/c+2vf0/c,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合。
步骤二:根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率fd和目标径向速度:
假设方位维MTD处理所需调频周期数为M,对距离向频谱进行方位向FFT后的频谱可以表示为:
从上式可以看出,对于多周期回波差拍信号频谱峰值的变化,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标来说,回波差拍信号频谱峰值只有exp{j4πkvf0Tr/c}项是变化的,即频谱峰值是受到多普勒频率fd=2vf0/c调制的。
步骤三:在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离和径向速度:
设经过MTD处理后,在两维频谱幅度最高的谱线对应的单元号为(m,n),第m个速度通道对应的目标多普勒频率和第n个距离通道对应的中心频率分别为:
fd=(m-M/2-1)fr/M (6)
fc=(n-1)fs/N (7)
其中,N为距离维FFT点数,fs为距离向采样率。
则该通道对应的速度可以表示为vt=cfd/2f0。从而可以获得每一速度通道的实际距离对应的频率为fb=fc-fd,根据频率和距离的对应关系,进而可以得到目标的实际距离为Rt=cfb/2kr
至此目标距离和速度计算完毕。
实施例2
下面给出一个实施例,进一步说明本发明技术方案的实施方式。
假设雷达工作频率f0为32.4GHz,发射信号调制带宽B为120MHz,发射重复周期Tr为3us,采样率fs为70MHz,处理的调频周期数为64,平台与目标距离设定为R0=10m,相对速度为v=700m/s。
根据本发明所述的技术方案,对多个调频周期接收的差频信号进行数字采样,然后希尔伯特滤波获得目标的差频复信号,对该信号进行64次128点FFT,64点MTD处理,得到目标二维频谱分布,频谱分布仿真结果如图4所示。
按以上处理步骤,在二维频域检测得到目标所在的多普勒和距离通道为(62,11),目标多普勒频率为fd=151.04kHz,fb=2.58MHz,得到目标的速度测量值为vt=699.27m/s,距离测量值为Rt=9.69m,由于存在测量噪声的影响,测量值与实际目标设定的距离和速度值基本吻合,满足同时测距测速的精度要求。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种全相参连续波多普勒雷达测距测速方法,其特征在于按照如下步骤进行:
步骤1、根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱;
步骤2、根据距离向频谱计算目标方位向的频谱,获取多普勒频率和目标径向速度,即利用多个相参的高重频调制信号独立获得无模糊的径向速度信息;
步骤3、在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离,即将速度引入的多普勒信息在距离维频域信息中去除,消除距离-速度耦合的影响,从而进一步得到距离信息。
2.根据权利要求1所述的全相参连续波多普勒雷达测距测速方法,其特征在于所述的步骤1中根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱:
1.1目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号;
1.2对差频信号进行距离FFT变换,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合。
3.根据权利要求1所述的全相参连续波多普勒雷达测距测速方法,其特征在于所述的步骤2中根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率和目标径向速度具体为:对距离向频谱进行方位向FFT变换,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标,其频谱峰值受到多普勒频率调制。
4.根据权利要求1或2任一所述的全相参连续波多普勒雷达测距测速方法,其特征在于:所述的步骤1中具体为:
设第K个调频周期内发射信号为S(k,t)可以表示为:
S(k,t)=A0cos{2π[f0(t-kTr)+0.5kr(t-kTr)2]+jφ0} (1)
其中,t为信号时间,Kr为调频率,Tr为调频周期,f0为载波频率,φ0为初始相位;
则对于t=0时刻目标距离为R0,径向速度为v的回波信号可以表示为:
S(k,t)=A1cos{2π[f0(t-kTr-τ(t))+0.5kr(t-kTr-τ(t))2]+jφ1} (2)
其中,点目标回波延时τ(t)=2(R0+vt)/c;
目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号,则第k个CPI周期的差拍信号表示为:
S ( k , t ) = A exp { j 2 π [ ( 2 R 0 c k r + 2 v c f 0 ) t + k 2 v c f 0 T r ] + j φ } - - - ( 3 )
其中,v为目标径向速度,φ为固定相位;
对差频信号进行距离向FFT后的频谱可以表示为:
S ( k , f r ) = A 1 sin c [ π ( f r - 2 R 0 c k r - 2 v c f 0 ) ] exp { j 2 π k 2 v c f 0 T r + j φ } - - - ( 4 )
可以得到,S(k,fr)峰值点处的中心频率fc=fb+fd=2krR0/c+2vf0/c,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合。
5.根据权利要求1所述的全相参连续波多普勒雷达测距测速方法,其特征在于所述的步骤2中:
根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率fd和目标径向速度:
假设方位维MTD处理所需调频周期数为M,对距离向频谱进行方位向FFT后的频谱可以表示为:
S ( f d , f r ) = A 1 sin c [ π ( f r - 2 R 0 c k r - 2 v c f 0 ) ] sin c [ π ( f d - 2 v c f 0 ) MT r ] exp { j φ } - - - ( 5 )
从上式可以看出,对于多周期回波差拍信号频谱峰值的变化,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标来说,回波差拍信号频谱峰值只有exp{j4πkvf0Tr/c}项是变化的,即频谱峰值是受到多普勒频率fd=2vf0/c调制的。
6.根据权利要求1所述的全相参连续波多普勒雷达测距测速方法,其特征在于所述的步骤3中:
在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离和径向速度:
设经过MTD处理后,在两维频谱幅度最高的谱线对应的单元号为(m,n),第m个速度通道对应的目标多普勒频率和第n个距离通道对应的中心频率分别为:
fd=(m-M/2-1)fr/M (6)
fc=(n-1)fs/N (7)
其中,N为距离维FFT点数,fs为距离向采样率;
则该通道对应的速度可以表示为vt=cfd/2f0,从而可以获得每一速度通道的实际距离对应的频率为fb=fc-fd,根据频率和距离的对应关系,进而可以得到目标的实际距离为Rt=cfb/2kr
7.一种全相参连续波多普勒测距测速雷达,其特征在于:
用于根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱的装置;
用于根据距离向频谱计算目标方位向的频谱,获取多普勒频率和目标径向速度的装置,即利用多个相参的高重频调制信号独立获得无模糊的径向速度信息;
用于在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离的装置,即将速度引入的多普勒信息在距离维频域信息中去除,消除距离-速度耦合的影响,从而进一步得到距离信息。
8.根据权利要求7所述的全相参连续波多普勒测距测速雷达,其特征在于所述的计算目标差拍回波信号的距离向频谱的装置包括:
用于目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号的装置
用于对差频信号进行距离FFT变换,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合的装置。
9.根据权利要求7所述的全相参连续波多普勒测距测速雷达,其特征在于所述的根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率和目标径向速度的装置具体为:对距离向频谱进行方位向FFT变换,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标,其频谱峰值受到多普勒频率调制。
10.根据权利要求7或8任一所述的全相参连续波多普勒测距测速雷达,其特征在于:
用于根据全相参准锯齿调频连续波发射波形,计算目标差拍回波信号的距离向频谱的装置具体为:
S(k,t)=A0cos{2π[f0(t-kTr)+0.5kr(t-kTr)2]+jφ0} (1)
其中,t为信号时间,Kr为调频率,Tr为调频周期,f0为载波频率,φ0为初始相位;
则对于t=0时刻目标距离为R0,径向速度为v的回波信号可以表示为:
S(k,t)=A1cos{2π[f0(t-kTr-τ(t))+0.5kr(t-kTr-τ(t))2]+jφ1} (2)
其中,点目标回波延时τ(t)=2(R0+vt)/c;
目标回波信号经过与发射信号混频获得差频信号,通过数字希尔伯特正交变换和低通滤波后,得到两路差频信号形成的复信号,则第k个CPI周期的差拍信号表示为:
S ( k , t ) = A exp { j 2 π [ ( 2 R 0 c k r + 2 v c f 0 ) t + k 2 v c f 0 T r ] + j φ } - - - ( 3 )
其中,v为目标径向速度,φ为固定相位;
对差频信号进行距离向FFT后的频谱可以表示为:
S ( k , f r ) = A 1 sin c [ π ( f r - 2 R 0 c k r - 2 v c f 0 ) ] exp { j 2 π k 2 v c f 0 T r + j φ } - - - ( 4 )
可以得到,S(k,fr)峰值点处的中心频率fc=fb+fd=2krR0/c+2vf0/c,包含两部分:第一项是由目标初始距离所产生的频移,第二项由运动目标的多普勒效应产生的距离-速度耦合;
用于根据距离向频谱计算目标方位向的频谱,获取多普勒频率和目标径向速度,即利用多个相参的高重频调制信号独立获得无模糊的径向速度信息的装置具体为:
设第K个调频周期内发射信号为S(k,t)可以表示为:
根据距离向频谱表达式,计算目标方位向的频谱,获取多普勒频率fd和目标径向速度:
假设方位维MTD处理所需调频周期数为M,对距离向频谱进行方位向FFT后的频谱可以表示为:
S ( f d , f r ) = A 1 sin c [ π ( f r - 2 R 0 c k r - 2 v c f 0 ) ] sin c [ π ( f d - 2 v c f 0 ) MT r ] exp { j φ } - - - ( 5 )
从上式可以看出,对于多周期回波差拍信号频谱峰值的变化,由于目标速度相对于光速而言很小,多周期时间也很短,则周期间由于目标距离变化产生的雷达接收信号强弱变化可以忽略不计,对于一个确知信号和目标来说,回波差拍信号频谱峰值只有exp{j4πkvf0Tr/c}项是变化的,即频谱峰值是受到多普勒频率fd=2vf0/c调制的;
用于在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离,即将速度引入的多普勒信息在距离维频域信息中去除,消除距离-速度耦合的影响,从而进一步得到距离信息的装置具体为:
在两维频域进行目标检测,获取目标所在的距离多普勒单元,计算目标距离和径向速度:
设经过MTD处理后,在两维频谱幅度最高的谱线对应的单元号为(m,n),第m个速度通道对应的目标多普勒频率和第n个距离通道对应的中心频率分别为:
fd=(m-M/2-1)fr/M (6)
fc=(n-1)fs/N (7)
其中,N为距离维FFT点数,fs为距离向采样率;
则该通道对应的速度可以表示为vt=cfd/2f0,从而可以获得每一速度通道的实际距离对应的频率为fb=fc-fd,根据频率和距离的对应关系,进而可以得到目标的实际距离为Rt=cfb/2kr
CN201610997913.4A 2016-11-14 2016-11-14 全相参连续波多普勒雷达及其测距测速方法 Active CN106405541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610997913.4A CN106405541B (zh) 2016-11-14 2016-11-14 全相参连续波多普勒雷达及其测距测速方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610997913.4A CN106405541B (zh) 2016-11-14 2016-11-14 全相参连续波多普勒雷达及其测距测速方法

Publications (2)

Publication Number Publication Date
CN106405541A true CN106405541A (zh) 2017-02-15
CN106405541B CN106405541B (zh) 2021-06-01

Family

ID=59230210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610997913.4A Active CN106405541B (zh) 2016-11-14 2016-11-14 全相参连续波多普勒雷达及其测距测速方法

Country Status (1)

Country Link
CN (1) CN106405541B (zh)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107861117A (zh) * 2017-10-17 2018-03-30 电子科技大学 一种适用于连续波周界监视雷达的多目标参数测量方法
CN108614267A (zh) * 2018-07-09 2018-10-02 南京普陆康电子科技有限公司 一种基于77GHz毫米波技术的防撞雷达设计方法
CN108693523A (zh) * 2017-12-06 2018-10-23 深圳大学 基于锯齿波线性调频雷达多目标测距测速的方法及系统
CN108761415A (zh) * 2018-05-25 2018-11-06 中国人民解放军国防科技大学 基于特定切换选通时序的雷达运动目标相位恢复的方法
CN109658715A (zh) * 2019-01-31 2019-04-19 厦门精益远达智能科技有限公司 多车道的车流量统计方法、装置、设备和存储介质
CN109703770A (zh) * 2018-12-12 2019-05-03 国耀量子雷达科技有限公司 基于测风激光雷达和cfd数据库的舰载机助降方法
CN110361726A (zh) * 2019-06-21 2019-10-22 广东工业大学 一种毫米波雷达测速方法
CN110412558A (zh) * 2019-07-03 2019-11-05 南京理工大学 基于tdm mimo的解车载fmcw雷达速度模糊方法
CN110488263A (zh) * 2018-05-14 2019-11-22 杭州海康威视数字技术股份有限公司 一种雷达设备的测量方法及雷达设备
CN110632587A (zh) * 2019-10-10 2019-12-31 南京凌云科技发展有限公司 一种基于快速fmcw雷达的弱运动物体监测方法
CN110873877A (zh) * 2019-04-25 2020-03-10 北京航空航天大学 目标运动轨迹的确定方法及装置
CN110927717A (zh) * 2019-12-11 2020-03-27 深圳大学 调频连续波雷达的成像方法、装置及成像系统
CN111505618A (zh) * 2020-04-14 2020-08-07 东南大学 一种适用于车载毫米波雷达领域的基于频率估计的去耦合校正测距方法
CN111630406A (zh) * 2017-12-07 2020-09-04 法雷奥开关和传感器有限责任公司 使用机动车辆雷达系统的不同传输器的相位编码fmcw斜坡进行改善目标分离的方法
CN112285684A (zh) * 2020-12-29 2021-01-29 南昌大学 相对运动速度跟踪方法、装置、可读存储介质及电子设备
CN113189577A (zh) * 2021-04-13 2021-07-30 电子科技大学 一种基于快速斜坡模式lfmcw车载雷达测距测速方法
CN113534125A (zh) * 2021-06-04 2021-10-22 惠州市德赛西威汽车电子股份有限公司 一种估算目标模糊速度的方法
WO2024061206A1 (zh) * 2022-09-19 2024-03-28 电子科技大学长三角研究院(衢州) 一种线性调频脉冲雷达的目标探测方法
CN113534125B (zh) * 2021-06-04 2024-06-07 惠州市德赛西威汽车电子股份有限公司 一种估算目标模糊速度的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4916452A (en) * 1985-11-16 1990-04-10 Aeg Process for range measurement with a pulse radar of high pulse repetition frequency
CN102121989A (zh) * 2010-01-08 2011-07-13 清华大学 一种同时测量高速运动目标的速度和距离的方法
CN102183762A (zh) * 2011-03-15 2011-09-14 北京航空航天大学 一种压缩感知合成孔径雷达数据获取与成像方法
CN102565763A (zh) * 2011-12-12 2012-07-11 中国船舶重工集团公司第七二四研究所 一种自适应杂波抑制动目标信号处理技术及实现方法
CN102565793A (zh) * 2011-12-30 2012-07-11 北京华航无线电测量研究所 一种单天线阵列全方位扫描的毫米波成像系统
CN104991248A (zh) * 2015-06-25 2015-10-21 苏州途视电子科技有限公司 脉冲全相参宽带通用化实验雷达

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4916452A (en) * 1985-11-16 1990-04-10 Aeg Process for range measurement with a pulse radar of high pulse repetition frequency
CN102121989A (zh) * 2010-01-08 2011-07-13 清华大学 一种同时测量高速运动目标的速度和距离的方法
CN102183762A (zh) * 2011-03-15 2011-09-14 北京航空航天大学 一种压缩感知合成孔径雷达数据获取与成像方法
CN102565763A (zh) * 2011-12-12 2012-07-11 中国船舶重工集团公司第七二四研究所 一种自适应杂波抑制动目标信号处理技术及实现方法
CN102565793A (zh) * 2011-12-30 2012-07-11 北京华航无线电测量研究所 一种单天线阵列全方位扫描的毫米波成像系统
CN104991248A (zh) * 2015-06-25 2015-10-21 苏州途视电子科技有限公司 脉冲全相参宽带通用化实验雷达

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘勇: "基于DRFM的PD雷达干扰技术及其实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107861117B (zh) * 2017-10-17 2020-11-17 电子科技大学 一种适用于连续波周界监视雷达的多目标参数测量方法
CN107861117A (zh) * 2017-10-17 2018-03-30 电子科技大学 一种适用于连续波周界监视雷达的多目标参数测量方法
CN108693523A (zh) * 2017-12-06 2018-10-23 深圳大学 基于锯齿波线性调频雷达多目标测距测速的方法及系统
CN111630406A (zh) * 2017-12-07 2020-09-04 法雷奥开关和传感器有限责任公司 使用机动车辆雷达系统的不同传输器的相位编码fmcw斜坡进行改善目标分离的方法
CN111630406B (zh) * 2017-12-07 2023-10-13 法雷奥开关和传感器有限责任公司 确定物体信息项的方法和雷达系统、车辆驾驶员辅助系统
CN110488263B (zh) * 2018-05-14 2021-07-20 杭州海康威视数字技术股份有限公司 一种雷达设备的测量方法及雷达设备
CN110488263A (zh) * 2018-05-14 2019-11-22 杭州海康威视数字技术股份有限公司 一种雷达设备的测量方法及雷达设备
CN108761415A (zh) * 2018-05-25 2018-11-06 中国人民解放军国防科技大学 基于特定切换选通时序的雷达运动目标相位恢复的方法
CN108761415B (zh) * 2018-05-25 2020-09-29 中国人民解放军国防科技大学 基于特定切换选通时序的雷达运动目标相位恢复的方法
CN108614267A (zh) * 2018-07-09 2018-10-02 南京普陆康电子科技有限公司 一种基于77GHz毫米波技术的防撞雷达设计方法
CN109703770A (zh) * 2018-12-12 2019-05-03 国耀量子雷达科技有限公司 基于测风激光雷达和cfd数据库的舰载机助降方法
CN109658715A (zh) * 2019-01-31 2019-04-19 厦门精益远达智能科技有限公司 多车道的车流量统计方法、装置、设备和存储介质
US11333749B2 (en) 2019-04-25 2022-05-17 Beihang University Method and device for determining motion trajectory of target
CN110873877A (zh) * 2019-04-25 2020-03-10 北京航空航天大学 目标运动轨迹的确定方法及装置
CN110361726B (zh) * 2019-06-21 2022-12-16 广东工业大学 一种毫米波雷达测速方法
CN110361726A (zh) * 2019-06-21 2019-10-22 广东工业大学 一种毫米波雷达测速方法
CN110412558A (zh) * 2019-07-03 2019-11-05 南京理工大学 基于tdm mimo的解车载fmcw雷达速度模糊方法
CN110412558B (zh) * 2019-07-03 2022-05-17 南京理工大学 基于tdm mimo的解车载fmcw雷达速度模糊方法
CN110632587A (zh) * 2019-10-10 2019-12-31 南京凌云科技发展有限公司 一种基于快速fmcw雷达的弱运动物体监测方法
CN110927717B (zh) * 2019-12-11 2022-01-11 深圳大学 调频连续波雷达的成像方法、装置及成像系统
CN110927717A (zh) * 2019-12-11 2020-03-27 深圳大学 调频连续波雷达的成像方法、装置及成像系统
CN111505618A (zh) * 2020-04-14 2020-08-07 东南大学 一种适用于车载毫米波雷达领域的基于频率估计的去耦合校正测距方法
CN112285684B (zh) * 2020-12-29 2021-03-16 南昌大学 相对运动速度跟踪方法、装置、可读存储介质及电子设备
CN112285684A (zh) * 2020-12-29 2021-01-29 南昌大学 相对运动速度跟踪方法、装置、可读存储介质及电子设备
CN113189577A (zh) * 2021-04-13 2021-07-30 电子科技大学 一种基于快速斜坡模式lfmcw车载雷达测距测速方法
CN113534125A (zh) * 2021-06-04 2021-10-22 惠州市德赛西威汽车电子股份有限公司 一种估算目标模糊速度的方法
CN113534125B (zh) * 2021-06-04 2024-06-07 惠州市德赛西威汽车电子股份有限公司 一种估算目标模糊速度的方法
WO2024061206A1 (zh) * 2022-09-19 2024-03-28 电子科技大学长三角研究院(衢州) 一种线性调频脉冲雷达的目标探测方法

Also Published As

Publication number Publication date
CN106405541B (zh) 2021-06-01

Similar Documents

Publication Publication Date Title
CN106405541A (zh) 全相参连续波多普勒雷达及其测距测速方法
RU2628566C1 (ru) Способ работы радиолокационной станции с повышенными допплеровскими характеристиками
US8305257B2 (en) Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents
CN105093227B (zh) 一种交通流量测量装置及车辆运行信息获得方法
CN107861117A (zh) 一种适用于连续波周界监视雷达的多目标参数测量方法
CN107966688A (zh) 基于相位干涉技术的宽带雷达目标速度解模糊方法
CN103760545A (zh) 合成孔径雷达中子带串扰对成像性能影响的分析方法
CN104614713A (zh) 一种适合于艇载雷达系统的雷达回波信号模拟器
CN102914773B (zh) 一种多航过圆周sar三维成像方法
US20220187158A1 (en) Method and apparatus of full-field vibration measurement via microwave sensing
CN107843892A (zh) 一种基于最小二乘法的高速目标多普勒测速方法
CN109738894A (zh) 一种雷达的成像方法
DE112017007143T5 (de) Radarvorrichtung
CN110850400A (zh) 基于干涉仪测向的lfmcw雷达多目标检测方法
CN105467369A (zh) 一种目标回波仿真方法和装置
CN103869298A (zh) 一种分布式mimo天波超视距雷达海杂波仿真方法
CN105676212A (zh) 一种近程测距雷达系统和基于该系统的目标测量方法
Sediono Method of measuring Doppler shift of moving targets using FMCW maritime radar
CN105116408A (zh) 一种舰船isar图像结构特征提取方法
CN109001729A (zh) 太赫兹成像中调频连续波线性度实时校准方法及其系统
CN104181521B (zh) 一种发射多频载波的高距离分辨率雷达
CN106054154A (zh) 一种基于icpf的机动目标的步进频率信号高分辨率成像方法
WO2011021262A1 (ja) レーダ装置
CN103412302A (zh) 基于先验知识的多载频miso雷达目标定位方法
CN108594216A (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
TR01 Transfer of patent right

Effective date of registration: 20240206

Address after: 264000 meters west of Housun Jia Bao Village, Gaogezhuang Town, Laiyang City, Yantai City, Shandong Province

Patentee after: Shandong Yinghe Road Freight Transport Co.,Ltd.

Country or region after: China

Address before: Room 603, building 1, kechuangyuan, 11 Sihai Road, Changshu Economic and Technological Development Zone, Suzhou City, Jiangsu Province

Patentee before: SUZHOU TUSHI ELECTRONIC TECHNOLOGY Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right