CN112965067A - 一种适用于fmcw汽车雷达目标速度扩展方法 - Google Patents
一种适用于fmcw汽车雷达目标速度扩展方法 Download PDFInfo
- Publication number
- CN112965067A CN112965067A CN201911278548.1A CN201911278548A CN112965067A CN 112965067 A CN112965067 A CN 112965067A CN 201911278548 A CN201911278548 A CN 201911278548A CN 112965067 A CN112965067 A CN 112965067A
- Authority
- CN
- China
- Prior art keywords
- speed
- signal
- frequency
- distance
- target
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000008878 coupling Effects 0.000 claims abstract description 30
- 238000010168 coupling process Methods 0.000 claims abstract description 30
- 238000005859 coupling reaction Methods 0.000 claims abstract description 30
- 238000001514 detection method Methods 0.000 claims abstract description 25
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 abstract description 3
- 230000005540 biological transmission Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- CLOMYZFHNHFSIQ-UHFFFAOYSA-N clonixin Chemical compound CC1=C(Cl)C=CC=C1NC1=NC=CC=C1C(O)=O CLOMYZFHNHFSIQ-UHFFFAOYSA-N 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
- G01S13/931—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种适用于FMCW汽车雷达目标速度扩展方法,该方法包括:确定发射信号波形;发射信号波第一部分为最小周期为Tc的N个连续锯齿波;第二部分为一个持续时间为Tcheck的锯齿波;将发射信号与接收信号进行混频、滤波再进行A/D转换;对第一部分时间信号进行至少两次FFT变换,求解出目标对应的距离以及模糊速度;对第二部分时间信号进行FFT变换求解出目距离和速度的第一耦合频率;使用模糊速度和最高检测速度计算出所有可能的实际速度,将实际速度以及距离代入耦合方程计算出对应的第二耦合频率;当判定第一耦合频率与第二耦合频率之间的差值小于预设门限值,将第二耦合频率对应的实际速度选定为运行速度。该方法能够有效处理目标速度模糊问题,同时兼具适用面广、运算简单、处理速度快的特点。
Description
技术领域
本发明涉及FMCW雷达信号处理技术领域,尤其涉及使用24GHz或77GHz频段采用FMCW汽车雷达目标速度扩展方法。
背景技术
近年来,毫米波雷达在高级驾驶员辅助系统(ADAS)中的使用呈现了爆发式的增长。毫米波雷达通过发射并处理接收到的电磁波,可以精确地测量出目标的距离,速度和角度,受天气,光线等环境因素的影响小,非常适合汽车雷达应用。雷达的探测距离,速度和角度的最大值,最小值和分辨率是其最重要的性能指标。这些指标受波形制式,算法,中射频电路性能,DSP处理能力和内存大小等因素综合影响。一般来说汽车雷达需要达到100到200Km/h的最高检测速度,现有的FMCW雷达受中射频电路性能的限制,很难在满足其它性能指标的前提下支持这么高的检测速度。如何在现有的硬件平台上,使用合适的波形设计或者信号处理方法解决速度模糊问题是当前汽车雷达信号处理的关键技术之一。
目前存在基于中国剩余定理,基于多普勒相偏补偿假设和基于目标跟踪假设等三种适合FMCW雷达的速度扩展方法,可通过相应的波形设计或数字信号处理方法将雷达的最高检测速度提高至车厂要求。上述三种方法各有特点,其中基于中国剩余定理的速度扩展算法适用面广,速度估计的准确度高,但是需要结合Fast-slow chirp的特殊波形,数据处理运算量较大。基于Doppler相偏补偿假设的速度扩展算法非常适合TDM-MIMO的波形配置,不需要改变雷达的波形,但是当前只能支持最大两倍的速度扩展,并且仅适用于MIMO雷达制式。基于目标跟踪假设的速度扩展算法充分利用了目标跟踪的历史信息,对目标检测层的要求低,但是加大了目标跟踪的计算量,对目标速度突变的响应会受跟踪周期的影响。如何提供一种速度扩展方法既能够有效处理目标速度模糊问题,同时兼具适用面广、运算简单、更新速度和处理速度快的特点是当前亟待解决的一个问题。
发明内容
有鉴于此,本发明提供一种适用于FMCW汽车雷达目标速度扩展方法用于解决现有技术的不足。
为达到上述目的,本发明采用如下技术方案:
本发明实施例提供了一种适用于FMCW汽车雷达目标速度扩展方法,其特征在于,包括:
根据FMCW雷达芯片最高中频频率、最大检测距离、距离分辨率以及速度分辨率确定实现目标速度扩展的发射信号波形;所述发射信号波形由两部分组成,第一部分为最小周期为Tc的N个连续锯齿波;第二部分为一个持续时间为Tcheck的锯齿波;其中,Tc小于Tcheck;
将发射信号与发射信号对应的接收信号进行混频、低通滤波再进行A/D转换;
对A/D转换后的数字信号按照第一部分时间信号和第二部分时间信号分类,对第一部分时间信号进行至少两次FFT变换,求解出目标对应的距离R以及模糊速度va;对第二部分时间信号进行FFT变换求解出目距离和速度的第一耦合频率f;
将模糊速度va代入公式vr=va+N·vmax,N=0,±2,±4…计算出所有可能的实际速度vr,其中,vmax为最高检测速度;
当判定第一耦合频率f与第二耦合频率fb之间的差值小于预设门限值,将第二耦合频率fb对应的实际速度vr选定为汽车的运行速度。
进一步地,所述根据FMCW雷达芯片最高中频频率、最大检测距离、距离分辨率以及速度分辨率确定实现目标速度扩展的发射信号波形具体包括:
进一步地,所述实际速度vr可选取值为va-2vmax、va以及va+2vmax;当实际速度vr为正值;实际速度vr取值选用va-2vmax及va;当实际速度vr为正值;实际速度vr取值选用va及va+2vmax。
进一步地,所述将发射信号与发射信号对应的接收信号进行混频、低通滤波具体表示如下:
发射信号的锯齿波用公式表示如下:
若目标距离为R,速度为v,则时延τ用公式表示如下:
则发射信号对应的接收信号用公式表示如下:
将发射信号与接收信号进行混频和滤除高频分量得到中频信号用公式表示如下:
将时延表达式代入上式整理得到:
进一步地,对于所述第一部分锯齿波,中频信号可简化表示如下:
对于所述第二部分锯齿波,中频信号可简化表示如下:
基于本发明提供的上述方案,能够有效处理目标速度模糊问题,同时兼具适用面广、运算简单、更新速度和处理速度快的特点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种适用于FMCW汽车雷达目标速度扩展方法的流程示意图;
图2为本发明实施例提供的发射信号波形图;
图3为本发明实施例提供的单目标仿真测量结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明实施例提供了一种适用于FMCW汽车雷达目标速度扩展方法,包括:
S101、根据FMCW雷达芯片最高中频频率、最大检测距离、距离分辨率以及速度分辨率确定实现目标速度扩展的发射信号波形。
如图2所示,发射信号波形由两部分组成,第一部分为最小周期为Tc的N个连续锯齿波;第二部分为一个持续时间为Tcheck的锯齿波;其中,Tc小于Tcheck;优选地,Tc为的范围为10-100微秒,Tcheck的范围为1-10毫秒;第一部分锯齿波的周期较短,用以测量反射目标的位置、已模糊速度和方向;第二部分锯齿波用以校验反射目标的位置和速度,从而解速度模糊问题,实现速度扩展。
在发射信号波形的确定过程中,将距离分辨率ΔR代入公式(1)中计算出锯齿波的带宽B,
将锯齿波的带宽B、最大检测距离R以及雷达芯片最高中频频率FIF带入公式(2)计算出第一部分锯齿波的最小周期TC,
将速度分辨率Δv代入公式(3)计算出第一部分锯齿波的最小周期TC,
S102、将发射信号与发射信号对应的接收信号进行混频、低通滤波再进行A/D转换。
将发射信号与接收信后进行混频、低通滤波再进行A/D转换,将对应的数字信号送入DSP中,为进一步处理做好准备。具体如下:
图2中的锯齿波发射信号用公式(4)表示如下:
公式(4)仅仅描述了第一个锯齿波的表达式,想要表述第n个锯齿波,可以将t定义为t=nTC+ts的形式,这样第n个锯齿波发送信号可以使用公式(5)表示
若目标距离为R,速度为v,则时延τ用公式(6)表示如下:
则发射信号对应的接收信号用公式表示如下:
将发射信号与接收信号进行混频和滤除高频分量得到中频信号用公式(8)表示如下:
将时延表达式代入上式整理得到表达式(9):
中频信号xm(t)进行A/D转换,进入下一级的DSP处理。
S103、对A/D转换后的数字信号按照第一部分时间信号和第二部分时间信号分类,对第一部分时间信号进行至少两次FFT变换,求解出目标对应的距离R以及模糊速度va;对第二部分时间信号进行FFT变换求解出目距离和速度的第一耦合频率f。
对第一部分时间信号进行第一次FFT变换可以求解出距离;对第一部分时间信号进行第二次FFT变换可以求解出速度;对第一部分时间信号进行第三次FFT变换可以求解出方位信息。
按照图2所示的雷达波形,其锯齿波分为第一部分的短周期锯齿波和第二部分的长周期锯齿波,由于第一部分锯齿波周期为微秒级,而第二部分锯齿波周期为毫秒级,因此,第一部分锯齿波斜率远大于第二部分锯齿波斜率。这样在第一部分锯齿波处理中可以忽略距离和速度的耦合关系,公式(9)可以化简为公式(10):
在公式(10)中,相位差应该小于π,否则就会产生相位模糊,从而导致速度模糊,使用该方法可以求解的最高检测速度vmax由公式(13)决定,
对于所述第二部分锯齿波,对应于斜率α比较小,中频信号可简化为式(14):
对其进行FFT变换时,距离R和速度v就会产生耦合,耦合频率表示如下:
其中,α为第二部分锯齿波的斜率;fc为雷达载波频率;c为光速。
S104、将模糊速度va代入公式vr=va+N·vmax,N=0,±2,±4…计算出所有可能的实际速度vr,其中,vmax为最高检测速度。
使用模糊速度和最高检测速度即可计算出所有可能的实际速度。
在第一部分信号中当一个目标的实际速度vr大于雷达的最高检测速度vmax时,通过第二次FFT得到的速度将会产生速度模糊现象,定义这个模糊速度为va,该值可以表达式(12)得到,此时模糊速度va和实际速度vr之间满足表达式(16)
va=vr mod vmax (16)
即模糊速度va是实际速度vr对最高检测速度vmax的余数,因此我们可以认为实际速度的可能值如公式(17)所示:
vr=va+N·vmax N=0,±2,±4… (17)
实际中因为汽车速度的限制,一般N=2就已经足够,即支持va-2vmax、va以及va+2vmax扩展。当实际速度vr为正值;实际速度vr取值选用va-2vmax及va;当实际速度vr为正值;实际速度vr取值选用va及va+2vmax。
S106、当判定第一耦合频率f与第二耦合频率fb之间的差值小于预设门限值,将第二耦合频率fb对应的实际速度vr选定为汽车的运行速度。
通过上述方法有效解决速度模糊问题,实现速度扩展,达到车厂要求。
以上总体实现过程可参见如下仿真实例:
某芯片,其中频信号频率FIF应该小于5MHz,当中频信号中5MHz以上的目标反射信号会被低通滤波器抑制。目前一个典型的汽车盲点检测雷达的系统参数要求:最大检测距离为80m、距离分辨率为0.375m、最高检测速度为±30m/s、速度分辨率为0.52m/s。选用77GHz雷达波。
根据公式(1)可知扫描带宽B不能小于3*108/(2*0.375)=400MHz,根据公式(2)可知锯齿波最小周期Tc为2*80*400*109/(5*109*3*108)=42.66us;结合扫描带宽B和最小周期Tc可以得出该芯片的波形斜率不能超过9.375MHz/us,根据公式13可知此时能够实现的最高检测速度为3*108/(4*77*109*42.66*10-6)=22.86m/s,无达到±30m/s的设计指标,因此需要进行速度扩展。
仿真时目标参数设置为目标距离R=60m,目标速度为35m/s,载波中心频率为77GHz。仿真结果如图3所示。
如图3所示,左边两幅图为使用第一部分信号求解的目标距离和模糊速度,分别为60m和-9.389m/s;最右图为使用第二部分信号求解的目标距离与速度耦合频率,其值为98.5kHz。按照前面所述当探测速度为负值时,真实速度的可能值为模糊速度va或者va+2vmax,分别将(R,va)和(R,va+2vmax)带入公式(15)进行计算。此时α为第二部分信号锯齿波频率斜率,具体值为2*1011,fc为载波频率,具体值为77*109。计算出来的频率分别为74.83kHz和99.15kHz,预设门限值可以设置为5kHz。差值(99.15-98.5)kHz小于5kHz,因此应该选择(R,va+2vmax)作为最终目标距离和速度,从而完成了速度扩展。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种适用于FMCW汽车雷达目标速度扩展方法,其特征在于,包括:
根据FMCW雷达芯片最高中频频率、最大检测距离、距离分辨率以及速度分辨率确定实现目标速度扩展的发射信号波形;所述发射信号波形由两部分组成,第一部分为最小周期为Tc的N个连续锯齿波;第二部分为一个持续时间为Tcheck的锯齿波;其中,Tc小于Tcheck;
将发射信号与发射信号对应的接收信号进行混频、低通滤波再进行A/D转换;
对A/D转换后的数字信号按照第一部分时间信号和第二部分时间信号分类,对第一部分时间信号进行至少两次FFT变换,求解出目标对应的距离R以及模糊速度va;对第二部分时间信号进行FFT变换求解出目距离和速度的第一耦合频率f;
将模糊速度va代入公式vr=va+N·vmax,N=0,±2,±4…计算出所有可能的实际速度vr,其中,vmax为最高检测速度;
当判定第一耦合频率f与第二耦合频率fb之间的差值小于预设门限值,将第二耦合频率fb对应的实际速度vr选定为汽车的运行速度。
4.根据权利要求1所述的方法,其特征在于,所述实际速度vr可选取值为va-2vmax、va以及va+2vmax;当实际速度vr为正值;实际速度vr取值选用va-2vmax及va;当实际速度vr为正值;实际速度vr取值选用va及va+2vmax。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911278548.1A CN112965067B (zh) | 2019-12-13 | 2019-12-13 | 一种适用于fmcw汽车雷达目标速度扩展方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911278548.1A CN112965067B (zh) | 2019-12-13 | 2019-12-13 | 一种适用于fmcw汽车雷达目标速度扩展方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112965067A true CN112965067A (zh) | 2021-06-15 |
CN112965067B CN112965067B (zh) | 2024-09-06 |
Family
ID=76270944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911278548.1A Active CN112965067B (zh) | 2019-12-13 | 2019-12-13 | 一种适用于fmcw汽车雷达目标速度扩展方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112965067B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI802994B (zh) * | 2021-09-17 | 2023-05-21 | 為昇科科技股份有限公司 | 雷達測速系統、方法及雷達裝置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19980074938A (ko) * | 1997-03-27 | 1998-11-05 | 김영환 | 불연속 톱니파 파형을 사용한 fmcw 레이다 구동장치 및 방법 |
CN106842166A (zh) * | 2017-03-31 | 2017-06-13 | 东南大学 | 一种适用于线性调频连续波雷达系统的解速度模糊方法 |
US20180172816A1 (en) * | 2016-12-20 | 2018-06-21 | National Chung-Shan Institute Of Science & Technology | Signal processing apparatus of a continuous-wave (cw) radar sensing system |
-
2019
- 2019-12-13 CN CN201911278548.1A patent/CN112965067B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19980074938A (ko) * | 1997-03-27 | 1998-11-05 | 김영환 | 불연속 톱니파 파형을 사용한 fmcw 레이다 구동장치 및 방법 |
US20180172816A1 (en) * | 2016-12-20 | 2018-06-21 | National Chung-Shan Institute Of Science & Technology | Signal processing apparatus of a continuous-wave (cw) radar sensing system |
CN106842166A (zh) * | 2017-03-31 | 2017-06-13 | 东南大学 | 一种适用于线性调频连续波雷达系统的解速度模糊方法 |
Non-Patent Citations (1)
Title |
---|
完诚;: "一种毫米波调频连续波雷达的目标速度测量方法", 舰船电子对抗, no. 04, 25 August 2016 (2016-08-25) * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI802994B (zh) * | 2021-09-17 | 2023-05-21 | 為昇科科技股份有限公司 | 雷達測速系統、方法及雷達裝置 |
Also Published As
Publication number | Publication date |
---|---|
CN112965067B (zh) | 2024-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106371097B (zh) | 雷达系统 | |
JP4994023B2 (ja) | パルスレーダ、車載レーダおよび着陸アシストレーダ | |
JP3606257B2 (ja) | ドップラーレーダー装置 | |
JP5595496B2 (ja) | レーダ装置 | |
JP4448184B2 (ja) | 信号処理方法 | |
JP2008232832A (ja) | 干渉判定方法,fmcwレーダ | |
JPH10282227A (ja) | 周波数シフトキーイング連続波レーダに適用される距離のアンビギティーの除去方法及びその装置 | |
CN112764036A (zh) | 一种自适应多模式车载雷达系统及设计方法 | |
CN110161472B (zh) | 一种基于信号复用的宽带车载毫米波雷达解速度模糊方法 | |
CN214041733U (zh) | 一种自适应多模式车载雷达系统 | |
EP1635192B1 (en) | Radar apparatus with DC offset correction | |
CN113805166A (zh) | 一种雷达物位计的目标跟踪测距方法及系统 | |
JP2009014405A (ja) | 車載用レーダ装置 | |
CN107783132B (zh) | 自动驾驶汽车防撞毫米波雷达系统及信号处理方法 | |
CN109085568B (zh) | 一种基于二次混频的调频连续波多目标检测方法 | |
CN110531358B (zh) | 信息测量方法及信息测量装置 | |
CN112444802A (zh) | 车载雷达的目标检测方法和检测装置、车辆 | |
CN112965067B (zh) | 一种适用于fmcw汽车雷达目标速度扩展方法 | |
CN111948641B (zh) | 基于时间调制技术的调频连续波探测系统及方法 | |
CN117008130A (zh) | 一种基于毫米波雷达的车辆目标速度变化测量方法及系统 | |
TW200819775A (en) | Tracking target method of radar with frequency modulation continuous wave | |
CN105227400A (zh) | 一种针对时变ads-b报文时延的自适应估计方法 | |
WO2012028938A1 (en) | Signal processor, signal processing method and radar apparatus | |
US20240004051A1 (en) | Radar ranging method, ranging radar, and unmanned aerial vehicle thereof | |
Wang | An approach for multiple moving targets detection and velocity estimation |
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 |