% S) c" i& e* @5 y Radar测距测速原理介绍
6 s: m. `- p* v( ^0 s5 x 通常来说,我们可以通过确定电磁波传播的时间来估算距离,通过多普勒效应产生的频率偏移来估算相对速度。但实际使用中我们会发现,对于没有经过调制的连续波谱(Continuous Wave),我们只发出相同的正弦或余弦信号,接收到的波也是呈周期性变化的,所以我们无从知道从发出到接受到底经过了多少个波长,从而无法确定信号从发出到接受所经历的时间(time of flight),也就无法测算距离。 1 n+ d8 f6 V; p0 {
于是,我们引入了调频/调幅(FMCW/AMCW),使发出信号的频率或者幅度处于不断的变化中。在汽车领域,我们一般使用调频技术。下图为脉冲FMCW和正弦FMCW时间域变化图。 ( O- i- s" a* R; |0 t' [
, I7 \" [: P7 ?6 }0 {0 ?( a FMCW雷达在扫频周期内发射频率变化的连续波,被物体反射后的回波与发射信号有一定的频率差,通过测量频率差可以获得目标与雷达之间的距离信息。 6 q1 r% @& w2 @, ]
FMCW雷达的测距/测速原理
- r9 u8 h1 S4 B 下图中,蓝色实线是发出的信号频率,红色虚线是接受到的信号频率, tct_{c} 为扫频周期的一半, fcf_{c} 为扫频带宽, τ\tau 为从信号发射到接受所经历的时间。
% R0 O. m G$ G* k5 D% z 6 I" D9 I& |3 x( R# K
设 ,fs(t),fe(t)f_{s}(t),f_e(t)分别为发出所发出和接受信号的频率变化函数,并假设相对速度 vrv_r 为0,则在信号的上升沿有如下关系: ; v [% H& d0 e1 ~
fS(t)=f0+fctc⋅tf_{S}(t) = f_0+\frac{f_c}{t_c}\cdot{t} ,
- @9 P* ?% p8 y" |( h9 Z fE=fS(t−τ)f_E=f_S(t-\tau) , _0 i% \1 d/ U: q0 t. L5 u H
且存在差频函数(beat frequency): : p9 d a- e! E3 t6 ~
fb(t)=fS(t)−fE(t)f_b(t)=f_S(t)-f_E(t)
6 ~1 G; _9 e L+ u' K! h 又因为 R=τ⋅c02R=\frac{\tau\cdot{c_0}}{2} ,从上图的几何关系 fbτ=fctc\frac{f_b}{\tau}=\frac{f_c}{t_c} 中可以得出:
4 P+ x: [+ e4 p4 O R=c02⋅tcfc⋅fbR=\frac{c_0}{2}\cdot\frac{t_c}{f_c}\cdot{f_b}
) Z+ Z. M9 q4 k* l6 X 从而得出:R和fbf_b 成正比
; N" g U7 K- U% D 设发出的信号为一余弦波,则其时间域的变化如下: C O$ c) n) n6 y3 T
; Z$ _' C& ~3 H9 n1 K/ f) M8 E 则接收到的信号在时间域变化如下:
" u( c0 ]1 M7 O( q . k u& f# i9 l5 O" t6 j1 Z
将 fEf_E 带入,得到: 9 g4 W5 i0 q& R) \! S
% Z& E) a2 ^% s7 t* m1 f0 m1 j 现在将 uEu_E 与 uSu_S 相乘,得到混合波 uMu_M : * S! N: K+ {. L. t- `
) q) Y/ I) N: m8 M7 Z 从上面的推导可以看出,最终得到的混合波,其频率只与 fbf_b 有关,而所要确定的相对距离R与fbf_b 成正比,所以只要得到混合波,就能确定相对距离R!
8 ~6 U: e R9 v: X 上面所讨论的是相对速度 vrv_r 为0的情况,R恒为定值。而实际情况中,我们要考虑有相对速度的情况。因为存在相对速度,距离R会随着时间进行变化,当相对速度 0">vr>0v_r>0 ,R不断增大,故 fbf_b 也不断变大;当相对速度 vr>0v_r>0 时,传递的电磁波被拉伸,频率有所减小;当 4 x- D P! m* x' t4 @- A
为了简化模型,且考虑到信号从发出到接收所经历的时间极短,加速度不变,汽车相对速度几乎不变,信号频率变化如下所示: ( t0 Y; J7 P1 N4 k) v8 C& D/ Y
% Q# A k S& `7 \ 从上图关系可以看出,差频 fbf_b 从原来的 fb=fctc⋅2Rc0f_b=\frac{f_c}{t_c}\cdot\frac{2R}{c_0} 增加了 vrv_r 影响项(不考虑相对论的多普勒效应),转变为如下关系式: , @0 E* g- ]" q' M5 B
+ F! \9 {, C# m2 S/ h9 _
从而得到了相对速度 vrv_r 和相对距离 RR 的关系式。而为了解得 vrv_r 和 RR 的具体数值,只有上面一个式子是不够的,于是我们在信号频率的下降沿进行重复的操作,从而得到如下方程组:
! f$ j' N; f! Y) Z% `
% {9 N& S# }0 D& X g0 g 其交点即为所求: 8 K$ {% J( @; X. ~# {" @
. ~' f: P9 I3 @7 N# u! | 可以解得相对距离和相对速度: : s% R5 U7 j6 `4 }) J7 F( }; V
# W. k) n& |. X' l- [$ b
如果希望测得多个物体的距离和速度,就必须对发射信号的频率进行一定的处理,如下图所示(以两个物体为例):
; K& @1 ~6 Q/ C5 ^ * J: F" S( H$ h
再进行上述一系列操作可以得到交点坐标图: 5 P0 L2 I# r- D: w2 q
; t* G |+ ^, a
如果按照原来的频率变换进行处理,则会得到四个交点(图中实线或者虚线各有四个交点),即得到四个距离和速度,从而无法对所探测物体进行判断。而改变第二个周期的扫频带宽后,将会得到不同扫频带宽发出信号的两个共同交点,从而确定所测两个物体与自己的相对距离和速度。实际中也可以改变扫频周期,来获得不同的扫频信号,如下图所示:
7 H2 N- O# w2 \7 M- g5 r
1 l8 I- ]# S0 [5 r 实际使用中,为了避免大量的计算,提高检测速度和精度,我们又引入了快速Chirp序列方法(fast cirp sequence),下一篇对此方法进行了详细介绍,并结合了实例进行分析,放个链接,大家有兴趣可点击阅读。
9 `! B: Z6 o2 r/ F9 t k6 c
: M, _% k- w* T( E6 E 目前笔者正在学习自动驾驶中的信号处理相关知识,如有错误,望及时纠正,不吝赐教。 $ _. A' k/ M6 w2 |1 M( ?$ @4 \: G
. C6 J* G3 L, l
7 x7 A, T8 p( _% S
1 _- ]6 c0 q1 S, K2 d% Y. R
$ y) e' p% ?9 [- ?+ h9 @1 @2 Z |