一文看懂!FMCW毫米波雷达测距原理

FMCW体制雷达

FMCW(Frequency Modulated Continuous Wave,调频连续波) Radar,顾名思义,这种雷达发射的波形为调频连续波。

信号模型

假设发射信号 s 1 ( t ) \rm{s}_1(t) s1(t)为:
s 1 ( t ) = exp ⁡ { j 2 π ( f 0 t + 0.5 K t 2 ) } , 0 ≤ t ≤ T \rm{s}_1(t)=\exp \left\{ {j2\pi ({f_0}t + 0.5K{t^2})} \right\},0 \le t \le T s1(t)=exp{j2π(f0t+0.5Kt2)},0tT
其中, T \rm{T} T是一个chirp的持续时间, K \rm{K} K是调频斜率。在一个周期T内,此信号的频率随着时间变化, f = k t f = kt f=kt
f 0 {{f_0}} f0表示载波,这个线性调频信号类似于一些调制,例如幅度调制角度调制等。
此信号也可以表示为:
s 1 ( t ) = exp ⁡ ( j 2 π f 0 t ) × exp ⁡ ( j 2 π × 0.5 K t 2 ) {{\rm{s}}_1}(t) = \exp (j2\pi {f_0}t) \times \exp (j2\pi \times 0.5K{t^2}) s1(t)=exp(j2πf0t)×exp(j2π×0.5Kt2)
exp ⁡ ( j 2 π f 0 t ) \exp (j2\pi {f_0}t) exp(j2πf0t)就是频率为 f 0 {f_0} f0的信号,成为载波。那么第二项就是调制信号(为方便理解,这里也可称之为“被调制信号”,需要注意的是所有教材中均称之为“调制信号”)。

反射信号

发射信号经过目标反射后,得到回波 s 2 ( t ) = p s 1 ( t − τ ) {{\rm{s}}_2}(t) = p{{\rm{s}}_1}(t - \tau ) s2(t)=ps1(tτ),这里假定 p = 1 p = 1 p=1

所以,回波信号 s 2 ( t ) {{\rm{s}}_2}(t) s2(t)可以表示为:
s 2 ( t ) = p s 1 ( t − τ ) = exp ⁡ { j 2 π ( f 0 ( t − τ ) + 0.5 K ( t − τ ) 2 ) } {{\rm{s}}_2}(t) = p{{\rm{s}}_1}(t - \tau ) = \exp \left\{ {j2\pi ({f_0}(t - \tau ) + 0.5K{{(t - \tau )}^2})} \right\} s2(t)=ps1(tτ)=exp{j2π(f0(tτ)+0.5K(tτ)2)}
对于FMCW体制的雷达,一般将回波信号与参考信号进行混频处理,以降低器件的成本。

混频信号

那么混频信号 s 3 ( t ) {{\rm{s}}_3}(t) s3(t)可以表示为:
s 3 ( t ) = s 1 ( t ) × s 2 ∗ ( t ) = exp ⁡ { j 2 π ( f 0 τ + K τ t − 0.5 K τ 2 ) } {{\rm{s}}_3}(t) = {{\rm{s}}_1}(t) \times {\rm{s}}_2^*(t) = \exp \{ j2\pi ({f_0}\tau + K\tau t - 0.5K{\tau ^2})\} s3(t)=s1(t)×s2(t)=exp{j2π(f0τ+Kτt0.5Kτ2)}
上式第一项是载频相位延迟,第二项 K τ K\tau Kτ保存着距离信息( τ = 2 R c \tau = \frac{{2R}}{c} τ=c2R, R R R为雷达与目标之间的距离, c c c为光速),而第三项也就是RVP(residual video phase)项可以忽略掉。

提取距离信息

保留有距离信息的频率 K τ K\tau Kτ(K的单位为Hz/s, τ \tau τ的单位为s)称之为中频或者拍频,用符号 f I F {f_{IF}} fIF表示。
对混频信号 s 3 ( t ) {{\rm{s}}_3}(t) s3(t)进行傅里叶变换即可得到距离信息。
S 3 ( f ) = F F T ( s 3 ( t ) ) = e j 2 π f 0 τ × sin ⁡ c ( π T ( f − 2 K R c ) )   {{\rm{S}}_3}{\rm{(f) = FFT(}}{{\rm{s}}_3}(t){\rm{) = }}{{\rm{e}}^{j2\pi {f_0}\tau }} \times \sin c(\pi T(f - \frac{{2KR}}{c}))\ S3(f)=FFT(s3(t))=ej2πf0τ×sinc(πT(fc2KR)) 
对上述函数进行傅里叶变换后是一个辛克函数,显然,在没有噪声的理想情况下,最终会在频率 f = K τ f=K\tau f=Kτ处得到最大值,也就是距离对应的中频频率。
f I F = K τ = K 2 R c {f_{IF}} = K\tau = K\frac{{2R}}{c} fIF=Kτ=Kc2R

仿真结果

假设有个目标距离雷达距离为10米,并进行仿真结果如下:
在这里插入图片描述

Matlab代码

作者上传资源 FMCW毫米波雷达测距仿真代码(matlab)(免费)

  • 5
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

月茗爱喝茶

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值