+ J7 X8 _& N! v6 w, o8 c
Radar测距测速原理介绍
& G; d! V& A' q# r# J 通常来说,我们可以通过确定电磁波传播的时间来估算距离,通过多普勒效应产生的频率偏移来估算相对速度。但实际使用中我们会发现,对于没有经过调制的连续波谱(Continuous Wave),我们只发出相同的正弦或余弦信号,接收到的波也是呈周期性变化的,所以我们无从知道从发出到接受到底经过了多少个波长,从而无法确定信号从发出到接受所经历的时间(time of flight),也就无法测算距离。 8 H, o7 y, ?. q! d& K, A4 l7 k4 S
于是,我们引入了调频/调幅(FMCW/AMCW),使发出信号的频率或者幅度处于不断的变化中。在汽车领域,我们一般使用调频技术。下图为脉冲FMCW和正弦FMCW时间域变化图。
6 I! m! n0 K6 i3 _6 I( s# @ 7 p5 Y5 r8 F/ t. H# T2 ~4 V& v
FMCW雷达在扫频周期内发射频率变化的连续波,被物体反射后的回波与发射信号有一定的频率差,通过测量频率差可以获得目标与雷达之间的距离信息。 : I. W+ ?- N7 X5 t$ `8 n
FMCW雷达的测距/测速原理" s. Z- M2 v; j( I( O1 a, f
下图中,蓝色实线是发出的信号频率,红色虚线是接受到的信号频率, tct_{c} 为扫频周期的一半, fcf_{c} 为扫频带宽, τ\tau 为从信号发射到接受所经历的时间。
" m# L4 F+ m- B, U( u( T& k2 A ( V5 W3 G8 t& S7 ?1 W$ D9 I! v/ y
设 ,fs(t),fe(t)f_{s}(t),f_e(t)分别为发出所发出和接受信号的频率变化函数,并假设相对速度 vrv_r 为0,则在信号的上升沿有如下关系: 8 M% m( O8 ~* E( w
fS(t)=f0+fctc⋅tf_{S}(t) = f_0+\frac{f_c}{t_c}\cdot{t} , * I9 I* \2 e& V" ? n
fE=fS(t−τ)f_E=f_S(t-\tau)
& M) [3 x) B7 ^9 ^ 且存在差频函数(beat frequency):
7 ]) ?$ {" y' _. U5 g3 ] fb(t)=fS(t)−fE(t)f_b(t)=f_S(t)-f_E(t) + W4 p( o* [2 d( ]: I. o
又因为 R=τ⋅c02R=\frac{\tau\cdot{c_0}}{2} ,从上图的几何关系 fbτ=fctc\frac{f_b}{\tau}=\frac{f_c}{t_c} 中可以得出:
+ _$ _ j+ Q# l% u. r2 r' d R=c02⋅tcfc⋅fbR=\frac{c_0}{2}\cdot\frac{t_c}{f_c}\cdot{f_b} & p: O N1 X4 v v) p( s
从而得出:R和fbf_b 成正比
8 ]$ C a! q3 o# p1 W( V/ I) n, E 设发出的信号为一余弦波,则其时间域的变化如下: - ]9 q3 z2 X) B) e
. P1 [# z+ }1 ~2 j3 q/ C2 k1 A 则接收到的信号在时间域变化如下:
5 Z4 ?& A0 b1 d K# ]& K & ^6 b8 S) b! ?7 e* ^& ^4 T
将 fEf_E 带入,得到:
$ j, Q! {+ `5 M9 s9 g9 B) X `
3 ?) e3 [/ Q9 ?% l" ]0 p4 c) K( t/ Z 现在将 uEu_E 与 uSu_S 相乘,得到混合波 uMu_M : 4 L3 r. B; i$ {$ w3 @7 f
6 l, J: O4 s! H+ e 从上面的推导可以看出,最终得到的混合波,其频率只与 fbf_b 有关,而所要确定的相对距离R与fbf_b 成正比,所以只要得到混合波,就能确定相对距离R! , x& \6 M2 E3 P, {; [1 u
上面所讨论的是相对速度 vrv_r 为0的情况,R恒为定值。而实际情况中,我们要考虑有相对速度的情况。因为存在相对速度,距离R会随着时间进行变化,当相对速度 0">vr>0v_r>0 ,R不断增大,故 fbf_b 也不断变大;当相对速度 vr>0v_r>0 时,传递的电磁波被拉伸,频率有所减小;当 : ?: G2 H1 V3 R6 B
为了简化模型,且考虑到信号从发出到接收所经历的时间极短,加速度不变,汽车相对速度几乎不变,信号频率变化如下所示:
% a/ u$ M2 O; H }, \ + g+ P% T. p/ x' K {$ g& C
从上图关系可以看出,差频 fbf_b 从原来的 fb=fctc⋅2Rc0f_b=\frac{f_c}{t_c}\cdot\frac{2R}{c_0} 增加了 vrv_r 影响项(不考虑相对论的多普勒效应),转变为如下关系式: . Q5 k8 p# a$ D! T- Y
$ A9 E4 q: \; ~! @: I2 I. c 从而得到了相对速度 vrv_r 和相对距离 RR 的关系式。而为了解得 vrv_r 和 RR 的具体数值,只有上面一个式子是不够的,于是我们在信号频率的下降沿进行重复的操作,从而得到如下方程组: 0 p. i$ Y& N8 _2 X. q
4 w2 Y) E3 C- M# q7 x 其交点即为所求:
6 _7 N8 g: W# a' T5 ?) M) S 6 X v( B( Y. w3 V, O; Y
可以解得相对距离和相对速度:
# P$ |& V0 x o: h
: _ A- H0 a2 ~8 J S( C7 y+ `& v 如果希望测得多个物体的距离和速度,就必须对发射信号的频率进行一定的处理,如下图所示(以两个物体为例): $ S" S& j! S- n! O6 w: {
8 L: [% E' ?# t
再进行上述一系列操作可以得到交点坐标图:
- i }' b* w' D& ] + z0 u5 H6 j9 d2 _0 h# L. o! y% j9 Q
如果按照原来的频率变换进行处理,则会得到四个交点(图中实线或者虚线各有四个交点),即得到四个距离和速度,从而无法对所探测物体进行判断。而改变第二个周期的扫频带宽后,将会得到不同扫频带宽发出信号的两个共同交点,从而确定所测两个物体与自己的相对距离和速度。实际中也可以改变扫频周期,来获得不同的扫频信号,如下图所示:
- L/ v+ d4 F0 K
+ S6 @7 q$ t! c" f% e 实际使用中,为了避免大量的计算,提高检测速度和精度,我们又引入了快速Chirp序列方法(fast cirp sequence),下一篇对此方法进行了详细介绍,并结合了实例进行分析,放个链接,大家有兴趣可点击阅读。
" ^0 f& c0 [: t: @; A! d
" R# s2 f3 c7 y: ?* d+ I 目前笔者正在学习自动驾驶中的信号处理相关知识,如有错误,望及时纠正,不吝赐教。
$ i& W# q" m+ ]1 X
$ q e3 b7 x2 t* y8 K2 J3 m6 }: @, s8 |. f" O
% L. R+ o# [7 ]1 g1 G+ u
0 \5 w0 z$ h3 Z |