CN109946659B - 一种车载毫米波雷达线性调频连续波运动频率扩展校正方法 - Google Patents
一种车载毫米波雷达线性调频连续波运动频率扩展校正方法 Download PDFInfo
- Publication number
- CN109946659B CN109946659B CN201910170757.8A CN201910170757A CN109946659B CN 109946659 B CN109946659 B CN 109946659B CN 201910170757 A CN201910170757 A CN 201910170757A CN 109946659 B CN109946659 B CN 109946659B
- Authority
- CN
- China
- Prior art keywords
- time dimension
- slow time
- fft operation
- frequency
- slow
- 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
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012937 correction Methods 0.000 title claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 40
- 230000008707 rearrangement Effects 0.000 claims abstract description 12
- 230000003595 spectral effect Effects 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种车载毫米波雷达线性调频连续波频率扩展校正方法,该方法步骤为:首先对每根接收天线的线性调频连续波中频信号,按慢时间维补零并进行FFT运算;其次,根据系统的速度分辨率计算新的慢时间维频谱单元对应的速度刻度;接着按照每个新的慢时间维频谱单元和快时间维时域单元,根据对应的速度刻度,对每根接收天线中频信号慢时间维FFT结果进行匹配重排并进行运动相位补偿;最后对匹配重排及相位补偿后的慢时间维FFT运算结果按照快时间维进行FFT运算,得到无运动频率扩展影响的线性调频连续波二维FFT运算结果。该方法能够解决已有技术方案中车载毫米波雷达线性调频连续波的运动频率扩展问题,明显地提升系统性能,其复杂度较低、实现简单。
Description
技术领域
本发明属于雷达信号处理和汽车电子技术领域,尤其涉及一种适用于车载毫米波雷达线性调频连续波的运动频率扩展校正方法。
背景技术
车载毫米波的研究始于20世纪70年代,早期受到元器件的限制发展较慢,随着近年来集成电路的进步、人工智能以及无人驾驶研究的深入,包括电路实现、市场分析以及信号处理在内的车载毫米波雷达研究变得活跃起来。车载毫米波雷达技术的发展与演变,自上世纪70年代最先路测的35GHz频段,到近来得到广泛商业应用的24/77GHz频段。
毫米波雷达探测范围适中,距离分辨率及距离精度较高,并且对目标的速度信息比较敏感。在常见的雾、雨、雪等天气状况下,毫米波雷达均能正常工作。即使在比较恶劣的环境下,毫米波雷达也能维持一定的性能水平。同时,毫米波电路由于其在硅片中的高度集成且价格低廉,完美适应无人驾驶场景下对传感器小体积、低功率以及低成本的限制和需求。基于以上原因,在无人驾驶汽车及高级驾驶辅助系统中,毫米波雷达相比于激光雷达、超声波雷达以及摄像机等传感器有着明显的优势和不可替代的作用。
车载毫米波雷达波形包括连续波、脉冲连续波、调频连续波以及步进连续波等,不同雷达波形因其距离/速度分辨率不同而具有各自的适用场景。其中,调频连续波的频率随着时间线性增长,信号产生和处理简单,通过发送和接收信号的频率差来检测并区分目标,能同时获得较高的距离/速度分辨率,并且距离分辨率与带宽成反比。在调频连续波中,锯齿波调制方式通过2维FFT处理来解目标距离参数及速度参数的频率耦合,其检测精度高、不涉及多目标的匹配问题,能获得更低的虚警率和漏警率,因此被广泛研究和使用。
FMCW一方面通过检测目标中频信号频率计算其对应的距离参数,另一方面针对多个扫频的相位变化来检测其多普勒频率。但由于FMCM为连续波,即在整个检测时间内持续发射信号,因此目标与雷达的径向距离由于相对运动的存在不能再被认为是固定的。具体而言,目标与雷达的距离变化不仅与慢时间维有关,而且还与快时间维即扫频时间有关。因此,FMCW中频信号中会随着目标与雷达的相对运动被引入一个运动频率偏移,该频率偏移会造成目标频率的扩展继而会对系统的整体性能造成影响。
目前针对FMCW车载毫米波雷达信号处理算法的研究中,多是忽略这一运动频率偏移带来的目标频率扩展现象。当系统扫频带宽较小、目标与雷达的相对速度较低时,该运动频率偏移较小,因此这种忽略处理对系统性能的影响较小。但是在系统扫频带宽较大,目标与雷达相对速度较高的情况下,这种忽略处理会严重降低系统性能。
发明内容
发明目的:针对目前车载毫米波雷达线性调频连续波的运动频率扩展问题,提出一种实现简单且能够灵活控制复杂度与性能折中的运动频率扩展校正方法。
技术方案:为实现本发明的目的,本发明所采用的技术方案是:一种车载毫米波雷达线性调频连续波运动频率扩展校正方法,所述的方法步骤包括:
(1)对每根接收天线的中频信号按照慢时间维补零并进行FFT运算;
(2)按照雷达系统参数计算系统速度分辨率,并由此计算得到速度刻度向量,作为新的慢时间维坐标刻度对应的速度值;
(3)对每根接收天线中频信号的慢时间维FFT运算结果,按照每个快时间维时域单元和每个新的慢时间维频谱单元,根据对应的速度刻度向量进行频谱单元位置重排;
(4)对每根接收天线中频信号进行频谱单元位置重排后的慢时间维FFT运算结果,按照每个快时间维时域单元和新的慢时间维频谱单元,计算其对应的运动相位补偿因子,并根据其进行运动相位补偿;
(5)对每根接收天线中频信号进行频谱单元位置重排和相位补偿后的慢时间维FFT运算结果,按照快时间维进行FFT运算,得到无运动频率扩展的车载毫米波雷达线性调频连续波中频信号二维FFT运算结果。
进一步的,所述步骤(1)中,对每根接收天线的中频信号按照慢时间维补零并进行FFT运算,得到其计算结果为Ns=M+Nc,其中Nq为系统单个扫频率的采样点数即快时间维采样点数,Ns为系统的慢时间维FFT运算点数,M为系统的扫频重复个数即慢时间维采样点数,Nc为系统慢时间维FFT运算时的补零点数。
进一步的,所述步骤(2)中的系统速度分辨率为:
式中,c为电磁波在自由空间的传播速度,T为系统扫频重复时间,f0为载波中心频率,由式(1)与系统慢时间维采样点数M计算新的慢时间维频谱单元位置m对应的速度值为:
V[m]=(M-m+1)·vres,m=1,2,...M。 (2)
进一步的,所述步骤(3)中,对每根接收天线中频信号慢时间维FFT运算结果,按照每个新的慢时间维频谱单元和快时间维时域单元根据其对应的速度值进行频谱单元位置重排,具体方法为:
式中,μ为扫频斜率,μ=B/T,B为扫频带宽;
将搜索得到的原慢时间维频谱单元位置z处的频谱信息赋值给新的慢时间维位置m处的频谱单元内:
式中,m=1,2,...M,对重排后的慢时间维FFT运算结果按照每个新的慢时间维坐标位置m和快时间维时域单元位置根据上式中计算的运动相位补偿因子,进行运动相位补偿,得到补偿后的慢时间维FFT运算结果为:
进一步的,所述步骤(5)中,对每根接收天线中频信号频谱单元位置重排和相位补偿后的慢时间维FFT运算结果按照快时间维进行FFT运算,得到无运动频率扩展的车载毫米波雷达线性调频连续波中频信号二维FFT运算结果其中
进一步的,所述的慢时间维FFT运算补零点数Nc可以灵活设置,通过控制Nc的数量以实现方法复杂度与系统性能之间的折中。其中,Nc的值越大,方法的复杂度越高,系统性能也越高。与之相反的,Nc的值越小,方法的复杂度越低,系统性能也越低。
有益效果:与现有技术相比,本发明的技术方案具有以下有益技术效果:
本发明方法能够有效地解决已有技术方案中车载毫米波雷达线性调频连续波的运动频率扩展问题,明显地提升系统性能,同时其复杂度较低、实现简单且能够通过控制慢时间维FFT运算的补零点数灵活地控制复杂度与性能之间的折中。
附图说明
图1为本发明方法流程图;
图2为本发明具体实施方式中不同方案的系统距离参数估计性能比较;
图3为本发明具体实施方式中不同方案的系统距离参数估计性能比较;
图4为本发明具体实施方式中不同方案的系统距离参数估计性能比较。
具体实施方式
下面结合具体实施案列,进一步阐明本发明,应理解这些实施案例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本具体实施方式中,车载毫米波雷达的实际系统参数如表1所示。
表1实际系统参数
本具体实施方式公开了一种适用于短距离车载雷达的测角补偿方法,包括以下步骤:
步骤2:由表1中的实际系统参数和式(1)计算系统分辨率为v'res=0.5945m/s,并由式(2)计算系统刻度向量为V'M=[0.5945 1.1890 1.7835 … 75.5013 76.0958];
步骤3:对每根接收天线中频信号慢时间维FFT运算结果Y',按照每个慢时间维频谱单元位置ns和快时间维时域单元位置由式(3)计算其对应的速度参数为按并根据其对应的速度参数搜索其与速度刻度向量中最接近的坐标位置,并作为新的慢时间维频谱单元位置,得到重排后的慢时间维FFT运算结果为其中
步骤4:对每根接收天线中频信号频谱单元位置重排后的慢时间维FFT运算结果,按照每个速度刻度单元m和快时间维时域单元由式(5)计算其对应的运动相位补偿因子为根据上式中计算的运动相位补偿因子及式(6),进行运动相位补偿,得到补偿后的慢时间维FFT运算结果为其中
图2、图3和图4分别为不同方案下的系统距离、速度和角度参数估计的性能比较仿真图。从图中可以看出,基于本发明所提的车载毫米波雷达线性调频连续波运动频率校正法方案的系统距离、速度和角度参数估计,在整个信噪比范围内均要优于已有方案。如图2和图3所示,尤其是在系统距离和速度参数估计方面,相比于已有方案,基于本发明所提方案的系统距离和速度参数估计平均误差降低了90%以上。
Claims (3)
1.一种车载毫米波雷达线性调频连续波频率扩展校正方法,其特征在于,所述的方法步骤包括:
(1)对每根接收天线的中频信号按照慢时间维补零并进行FFT运算;
(2)按照雷达系统参数计算系统速度分辨率,并由此计算得到速度刻度向量,作为新的慢时间维坐标刻度对应的速度值;
(3)对每根接收天线中频信号的慢时间维FFT运算结果,按照每个快时间维时域单元和每个新的慢时间维频谱单元,根据对应的速度刻度向量进行频谱单元位置重排;
(4)对每根接收天线中频信号进行频谱单元位置重排后的慢时间维FFT运算结果,按照每个快时间维时域单元和新的慢时间维频谱单元,计算其对应的运动相位补偿因子,并根据其进行运动相位补偿;
(5)对每根接收天线中频信号进行频谱单元位置重排和相位补偿后的慢时间维FFT运算结果,按照快时间维进行FFT运算,得到无运动频率扩展的车载毫米波雷达线性调频连续波中频信号二维FFT运算结果;
所述步骤(1)中,对每根接收天线的中频信号按照慢时间维补零并进行FFT运算,得到其计算结果为Ns=M+Nc,其中,Nq为系统单个扫频率的采样点数,即快时间维采样点数,Ns为系统的慢时间维FFT运算点数,M为系统的扫频重复个数,即慢时间维采样点数,Nc为系统慢时间维FFT运算时的补零点数;
所述步骤(2)中的系统速度分辨率为:
式中,c为电磁波在自由空间的传播速度,T为系统扫频重复时间,f0为载波中心频率,由式(1)与系统慢时间维采样点数M计算新的慢时间维频谱单元位置m对应的速度值为:
V[m]=(M-m+1)·vres,m=1,2,...M (2)
所述步骤(3)中,对每根接收天线中频信号慢时间维FFT运算结果,按照每个新的慢时间维频谱单元和快时间维时域单元根据其对应的速度值进行频谱单元位置重排,具体方法为:
式中,μ为扫频斜率,μ=B/T,B为扫频带宽;
将搜索得到的原慢时间维频谱单元位置z处的频谱信息赋值给新的慢时间维位置m处的频谱单元内:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910170757.8A CN109946659B (zh) | 2019-03-07 | 2019-03-07 | 一种车载毫米波雷达线性调频连续波运动频率扩展校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910170757.8A CN109946659B (zh) | 2019-03-07 | 2019-03-07 | 一种车载毫米波雷达线性调频连续波运动频率扩展校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109946659A CN109946659A (zh) | 2019-06-28 |
CN109946659B true CN109946659B (zh) | 2022-11-04 |
Family
ID=67009180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910170757.8A Active CN109946659B (zh) | 2019-03-07 | 2019-03-07 | 一种车载毫米波雷达线性调频连续波运动频率扩展校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109946659B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111537966B (zh) * | 2020-04-28 | 2022-06-10 | 东南大学 | 一种适用于毫米波车载雷达领域的阵列天线误差校正方法 |
CN113267778B (zh) * | 2021-06-18 | 2024-08-30 | 珠海上富电技股份有限公司 | 一种车载毫米波雷达出厂检测装置及方法 |
US20230194705A1 (en) * | 2021-12-20 | 2023-06-22 | Aptiv Technologies Limited | Movement Compensation for Multi-Dimensional Radar Signal Processing |
CN114814815B (zh) * | 2022-03-16 | 2024-08-16 | 华域汽车系统股份有限公司 | 一种基于插值和单点dft滤波的解决信号处理跨越损失的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106814359B (zh) * | 2016-12-09 | 2019-06-04 | 北京无线电计量测试研究所 | 一种机场跑道异物检测雷达信号预处理方法 |
CN106646447B (zh) * | 2017-01-18 | 2019-03-26 | 武汉雷博合创电子科技有限公司 | 基于线性调频连续波的雷达目标长时间积累检测方法 |
CN107688178A (zh) * | 2017-08-25 | 2018-02-13 | 上海通趣科技有限公司 | 一种基于77GHz毫米波雷达的锯齿波测距测速方法 |
CN109061589B (zh) * | 2018-07-06 | 2022-08-26 | 西安电子科技大学 | 随机跳频雷达的目标运动参数估计方法 |
-
2019
- 2019-03-07 CN CN201910170757.8A patent/CN109946659B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109946659A (zh) | 2019-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109946659B (zh) | 一种车载毫米波雷达线性调频连续波运动频率扩展校正方法 | |
CN100365429C (zh) | 一种合成孔径雷达运动目标成像方法 | |
Zhou et al. | Approach for single channel SAR ground moving target imaging and motion parameter estimation | |
US9746554B2 (en) | Radar imaging system and related techniques | |
CN112098990B (zh) | 车载高分辨毫米波雷达对于中高速车辆的检测与跟踪方法 | |
US11092686B2 (en) | Method, apparatus and device for doppler compensation in a time switched MIMO radar system | |
CN108490443B (zh) | 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法 | |
CN102707269B (zh) | 一种机载雷达距离走动校正方法 | |
CN103454632A (zh) | 一站固定式调频连续波双基地sar成像方法 | |
CN111505618A (zh) | 一种适用于车载毫米波雷达领域的基于频率估计的去耦合校正测距方法 | |
CN111239721A (zh) | 车载mimo雷达求熵解速度模糊的方法 | |
CN116075745A (zh) | 低信号处理负荷情况下具有高距离分辨率的雷达调制方法 | |
CN115421134A (zh) | 一种雷达的速度解模糊的方法、装置及毫米波雷达 | |
CN114200411A (zh) | 一种mimo雷达多目标测速扩展方法 | |
CN108983193A (zh) | 快速非搜索的地面运动目标参数估计方法 | |
CN112098965A (zh) | 一种基于毫米波雷达的载体姿态角解算方法 | |
US12000950B2 (en) | Estimation of transverse velocities or cartesian velocities of point targets with a radar sensor | |
CN117008130A (zh) | 一种基于毫米波雷达的车辆目标速度变化测量方法及系统 | |
TW200819775A (en) | Tracking target method of radar with frequency modulation continuous wave | |
WO2023137682A1 (zh) | 一种数据处理的方法、装置和交通工具 | |
CN114047501B (zh) | 一种基于毫米波雷达的室内定位系统 | |
CN110988860B (zh) | 基于三角波调制的调频连续波sar运动补偿方法 | |
CN114994676A (zh) | 一种一站固定式双站低频超宽带sar运动目标成像方法 | |
CN112965067B (zh) | 一种适用于fmcw汽车雷达目标速度扩展方法 | |
US11493596B2 (en) | Estimation of cartesian velocities of extended radar objects using a radar sensor |
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 |