. [9 E2 o5 I; ]0 t6 f4 Z1 f3 n Radar测距测速原理介绍2 t) n& b& J v4 F. L1 P
通常来说,我们可以通过确定电磁波传播的时间来估算距离,通过多普勒效应产生的频率偏移来估算相对速度。但实际使用中我们会发现,对于没有经过调制的连续波谱(Continuous Wave),我们只发出相同的正弦或余弦信号,接收到的波也是呈周期性变化的,所以我们无从知道从发出到接受到底经过了多少个波长,从而无法确定信号从发出到接受所经历的时间(time of flight),也就无法测算距离。 9 H( Z/ Q0 `4 S- g2 g2 R' g8 d
于是,我们引入了调频/调幅(FMCW/AMCW),使发出信号的频率或者幅度处于不断的变化中。在汽车领域,我们一般使用调频技术。下图为脉冲FMCW和正弦FMCW时间域变化图。
" Z$ g4 i; o; U ( C1 M4 f+ j3 P/ U
FMCW雷达在扫频周期内发射频率变化的连续波,被物体反射后的回波与发射信号有一定的频率差,通过测量频率差可以获得目标与雷达之间的距离信息。 " A4 W2 _) m! M/ @9 f8 H: r ]6 {
FMCW雷达的测距/测速原理
2 X& M6 V. _( a# l 下图中,蓝色实线是发出的信号频率,红色虚线是接受到的信号频率, tct_{c} 为扫频周期的一半, fcf_{c} 为扫频带宽, τ\tau 为从信号发射到接受所经历的时间。 + L1 o& v: T3 b3 W
2 W# B% h M9 u* B6 C# S; R/ a C
设 ,fs(t),fe(t)f_{s}(t),f_e(t)分别为发出所发出和接受信号的频率变化函数,并假设相对速度 vrv_r 为0,则在信号的上升沿有如下关系:
1 v1 `0 _! ]2 n+ S/ ]8 ~, o fS(t)=f0+fctc⋅tf_{S}(t) = f_0+\frac{f_c}{t_c}\cdot{t} , * ?3 C( a D2 }
fE=fS(t−τ)f_E=f_S(t-\tau) ) a% H% b/ W5 x* ~4 ^* M* f
且存在差频函数(beat frequency): + C% k) z0 h( y: d
fb(t)=fS(t)−fE(t)f_b(t)=f_S(t)-f_E(t)
- u/ U ]$ R: j( F, S [" e6 R% Q 又因为 R=τ⋅c02R=\frac{\tau\cdot{c_0}}{2} ,从上图的几何关系 fbτ=fctc\frac{f_b}{\tau}=\frac{f_c}{t_c} 中可以得出:
8 Q6 n" }+ j( v' ]& j: G3 d R=c02⋅tcfc⋅fbR=\frac{c_0}{2}\cdot\frac{t_c}{f_c}\cdot{f_b}
% P+ j3 I; j7 q, c/ }* o- E 从而得出:R和fbf_b 成正比 5 T- Q9 i: }4 f- W* C+ w
设发出的信号为一余弦波,则其时间域的变化如下:
0 U+ H, F) ]0 ?& | 8 R( g& Z" j" s- \) C1 l. W
则接收到的信号在时间域变化如下: * }; T) g$ X9 P. U# U- r
+ D/ n/ Z8 E! N v
将 fEf_E 带入,得到: " {" }4 u3 m# O' x$ e" o; H, X
4 D8 ]9 c# o: j. } 现在将 uEu_E 与 uSu_S 相乘,得到混合波 uMu_M :
1 k* C8 t9 Y$ e5 n1 ~& O
* h2 F: Z) C1 J t6 C- T8 q+ S+ G 从上面的推导可以看出,最终得到的混合波,其频率只与 fbf_b 有关,而所要确定的相对距离R与fbf_b 成正比,所以只要得到混合波,就能确定相对距离R!
1 @$ ~4 L2 k0 V6 e' v; O7 X5 W! J 上面所讨论的是相对速度 vrv_r 为0的情况,R恒为定值。而实际情况中,我们要考虑有相对速度的情况。因为存在相对速度,距离R会随着时间进行变化,当相对速度 0">vr>0v_r>0 ,R不断增大,故 fbf_b 也不断变大;当相对速度 vr>0v_r>0 时,传递的电磁波被拉伸,频率有所减小;当
, W9 p9 V* O1 N T: r+ } 为了简化模型,且考虑到信号从发出到接收所经历的时间极短,加速度不变,汽车相对速度几乎不变,信号频率变化如下所示: 4 P' ^& N0 {* Y. J$ \# s* t4 \
( J# m5 N% l4 f. J' J
从上图关系可以看出,差频 fbf_b 从原来的 fb=fctc⋅2Rc0f_b=\frac{f_c}{t_c}\cdot\frac{2R}{c_0} 增加了 vrv_r 影响项(不考虑相对论的多普勒效应),转变为如下关系式: ( c$ p F: d( t7 i
/ U8 Z+ l0 f) B* w/ {' q# G
从而得到了相对速度 vrv_r 和相对距离 RR 的关系式。而为了解得 vrv_r 和 RR 的具体数值,只有上面一个式子是不够的,于是我们在信号频率的下降沿进行重复的操作,从而得到如下方程组: c' i. p3 ]$ r" ^% h0 _
2 j: b5 o# X8 S 其交点即为所求: 0 t! [! _! i. W* f2 D5 c2 }
7 ?/ m& i" B1 L
可以解得相对距离和相对速度: ! i1 M# O) Z# J( c( u" d8 D
+ O4 |: I3 k, _9 o0 K" v' D 如果希望测得多个物体的距离和速度,就必须对发射信号的频率进行一定的处理,如下图所示(以两个物体为例):
1 T: M- Z6 h7 o4 z; t9 U7 T. t( m # I5 \3 N8 h5 k: a2 \! W
再进行上述一系列操作可以得到交点坐标图:
" [. t* H k- L+ ?% I( q7 u ( o, N; ?$ n7 H. C0 G5 ^
如果按照原来的频率变换进行处理,则会得到四个交点(图中实线或者虚线各有四个交点),即得到四个距离和速度,从而无法对所探测物体进行判断。而改变第二个周期的扫频带宽后,将会得到不同扫频带宽发出信号的两个共同交点,从而确定所测两个物体与自己的相对距离和速度。实际中也可以改变扫频周期,来获得不同的扫频信号,如下图所示:
O* o* M3 v7 O( d
! \$ F* \1 {) D! _) v6 T3 Z9 w 实际使用中,为了避免大量的计算,提高检测速度和精度,我们又引入了快速Chirp序列方法(fast cirp sequence),下一篇对此方法进行了详细介绍,并结合了实例进行分析,放个链接,大家有兴趣可点击阅读。
% t+ j0 M' J: {$ p/ A
" t1 s) a! X' A+ }" F' s; j1 M 目前笔者正在学习自动驾驶中的信号处理相关知识,如有错误,望及时纠正,不吝赐教。 9 C/ Z- k6 z* ^9 g- r2 w/ `
" C/ B0 X2 e; u6 X) b5 r2 U: Y/ v' a
0 \! ]& e1 a# Z" J" b
1 O' ], r y) d5 Y, W9 C V
/ N! q" H0 q% W t% s0 h |