7 G* d2 \# l6 v! n4 _1 M
Radar测距测速原理介绍2 ?/ h2 v) _5 n: r
通常来说,我们可以通过确定电磁波传播的时间来估算距离,通过多普勒效应产生的频率偏移来估算相对速度。但实际使用中我们会发现,对于没有经过调制的连续波谱(Continuous Wave),我们只发出相同的正弦或余弦信号,接收到的波也是呈周期性变化的,所以我们无从知道从发出到接受到底经过了多少个波长,从而无法确定信号从发出到接受所经历的时间(time of flight),也就无法测算距离。
- S- t4 e4 I$ `8 E 于是,我们引入了调频/调幅(FMCW/AMCW),使发出信号的频率或者幅度处于不断的变化中。在汽车领域,我们一般使用调频技术。下图为脉冲FMCW和正弦FMCW时间域变化图。 8 P+ \* N. j f1 N$ Y% o
% \+ L: j/ }/ Z( l# N) |
FMCW雷达在扫频周期内发射频率变化的连续波,被物体反射后的回波与发射信号有一定的频率差,通过测量频率差可以获得目标与雷达之间的距离信息。 1 S: Y% V) f6 m. t1 f
FMCW雷达的测距/测速原理
/ o0 t, A2 _+ k+ f" H$ g7 [ 下图中,蓝色实线是发出的信号频率,红色虚线是接受到的信号频率, tct_{c} 为扫频周期的一半, fcf_{c} 为扫频带宽, τ\tau 为从信号发射到接受所经历的时间。
$ `( j, X+ U" \! c# h. B ] ) x6 p* \, g3 p" _+ x( u5 H8 n( f6 I* v
设 ,fs(t),fe(t)f_{s}(t),f_e(t)分别为发出所发出和接受信号的频率变化函数,并假设相对速度 vrv_r 为0,则在信号的上升沿有如下关系:
" D G% m3 i: N: o X fS(t)=f0+fctc⋅tf_{S}(t) = f_0+\frac{f_c}{t_c}\cdot{t} ,
7 z U. s" c+ |2 Z4 A* f fE=fS(t−τ)f_E=f_S(t-\tau) / o2 h. ~1 }5 ~% ^. X* q
且存在差频函数(beat frequency): e5 M5 z, B4 P6 m; O0 x
fb(t)=fS(t)−fE(t)f_b(t)=f_S(t)-f_E(t)
' V+ `6 a1 D! n( H# f/ f 又因为 R=τ⋅c02R=\frac{\tau\cdot{c_0}}{2} ,从上图的几何关系 fbτ=fctc\frac{f_b}{\tau}=\frac{f_c}{t_c} 中可以得出: 0 I/ K9 \9 l: `, X- b# N, k
R=c02⋅tcfc⋅fbR=\frac{c_0}{2}\cdot\frac{t_c}{f_c}\cdot{f_b} ( i$ \3 P, v$ f7 ]9 l0 m- k
从而得出:R和fbf_b 成正比 - n7 z- M& w* i6 a6 Y
设发出的信号为一余弦波,则其时间域的变化如下: ) H1 F6 }( q' r2 K5 E" D- V
; e8 W7 s' k( y, {& \ 则接收到的信号在时间域变化如下:
) N [5 F. ^: F8 g0 a+ u ) W: E% X$ _: r4 |* J
将 fEf_E 带入,得到: 7 W( ~! }& s d9 `
% @9 l# _. c7 t2 x- l* x& I 现在将 uEu_E 与 uSu_S 相乘,得到混合波 uMu_M : 3 m2 ^8 |* K' V+ ?3 ]0 e9 \
g Z& F/ v& v3 |5 Z
从上面的推导可以看出,最终得到的混合波,其频率只与 fbf_b 有关,而所要确定的相对距离R与fbf_b 成正比,所以只要得到混合波,就能确定相对距离R!
: ]9 v# N% I8 u2 y$ A 上面所讨论的是相对速度 vrv_r 为0的情况,R恒为定值。而实际情况中,我们要考虑有相对速度的情况。因为存在相对速度,距离R会随着时间进行变化,当相对速度 0">vr>0v_r>0 ,R不断增大,故 fbf_b 也不断变大;当相对速度 vr>0v_r>0 时,传递的电磁波被拉伸,频率有所减小;当 7 Q3 r% U. P% N: l% ]$ H, q
为了简化模型,且考虑到信号从发出到接收所经历的时间极短,加速度不变,汽车相对速度几乎不变,信号频率变化如下所示: ) d7 L4 h! }5 I! W
* H9 N0 u Q; U/ n# S 从上图关系可以看出,差频 fbf_b 从原来的 fb=fctc⋅2Rc0f_b=\frac{f_c}{t_c}\cdot\frac{2R}{c_0} 增加了 vrv_r 影响项(不考虑相对论的多普勒效应),转变为如下关系式: 1 @0 m; ~: |9 @% F9 L: U
, n( ]1 y! ?" c3 } 从而得到了相对速度 vrv_r 和相对距离 RR 的关系式。而为了解得 vrv_r 和 RR 的具体数值,只有上面一个式子是不够的,于是我们在信号频率的下降沿进行重复的操作,从而得到如下方程组: : {/ \+ V p" b. M3 n
3 e) B: v5 c8 `" j9 r* j" c
其交点即为所求:
: p2 f( Z# I% g7 x" O
. a1 G s# W' M' _ A. }. A 可以解得相对距离和相对速度: , A0 Q+ d" H1 w& V/ Q) z2 {4 W
& b$ z" p% J/ o1 v4 A# R8 b 如果希望测得多个物体的距离和速度,就必须对发射信号的频率进行一定的处理,如下图所示(以两个物体为例): 4 X1 w& J( _- n; J: C8 f! u
. G, w* F6 O0 b3 F 再进行上述一系列操作可以得到交点坐标图:
. {. v8 N( r/ U5 Q- F
, j. S3 U" R8 B1 L! ~6 |: C 如果按照原来的频率变换进行处理,则会得到四个交点(图中实线或者虚线各有四个交点),即得到四个距离和速度,从而无法对所探测物体进行判断。而改变第二个周期的扫频带宽后,将会得到不同扫频带宽发出信号的两个共同交点,从而确定所测两个物体与自己的相对距离和速度。实际中也可以改变扫频周期,来获得不同的扫频信号,如下图所示:
$ [, K9 u1 k0 Y4 p/ s: l . e( M7 g5 g& i8 L
实际使用中,为了避免大量的计算,提高检测速度和精度,我们又引入了快速Chirp序列方法(fast cirp sequence),下一篇对此方法进行了详细介绍,并结合了实例进行分析,放个链接,大家有兴趣可点击阅读。 # R9 W) L8 f! o: F8 @! P
$ D$ }4 m/ }7 K+ ?0 k
目前笔者正在学习自动驾驶中的信号处理相关知识,如有错误,望及时纠正,不吝赐教。 1 ]; G- k% j; z1 z$ h* v
1 u9 Z2 Q5 W6 w; F7 w
2 L4 f$ f. e% T" B$ q
# e5 g* J5 V- d" u& u# `, Z$ K5 j2 u
|