【FMCW 02】测距

承接上篇博文 中频IF信号 ,我们已经知道得到的中频IF信号的形式为:

x I F ( t ) = A ′ ′ cos ⁡ ( 2 π K τ t + 2 π f o τ ) x_{\tiny{IF}}(t) = A^{\prime \prime} \cos(2\pi K\tau t+2\pi f_o \tau ) xIF(t)=A′′cos(2πKτt+2πfoτ)


在这里插入图片描述

其中时延 τ = 2 d c \tau = \frac{2d}{c} τ=c2d d d d 为探测目标的距离。

在实际中,我们对得到的IF信号进行ADC采样,即将模拟信号(analog signal)转换为数字信号(digital signal),然后对采样后的信号进行数字信号处理

最大探测距离

假设ADC的采样率为 f s f_s fs,采样点数为M,FFT点数为N(N≥M)。根据FFT频谱分析的知识,我们知道此时最大能分辨的频率分量为 f s 2 \frac{f_s}{2} 2fs

对采样后的IF信号作FFT,即可完成对其的频谱分析,假设我们测得在频率 f f f 处有一个峰值(peak),其就对应着在距离 d d d 处有一个目标。即:

f p e a k = K τ = K 2 d t a r g e t c f_{peak}=K\tau=K\frac{2d_{target}}{c} fpeak=Kτ=Kc2dtarget

进行距离的换算,我们就有:

d t a r g e t = c 2 K f p e a k = c ⋅ C R T 2 B W f p e a k d_{target} = \frac{c}{2K}f_{peak}=\frac{c \cdot CRT }{2BW}f_{peak} dtarget=2Kcfpeak=2BWcCRTfpeak

其中,CRT指的chirp重复的周期,即chirp repetition time,BW指的是chirp的扫频带宽,即band width。


现在我们看 f p e a k f_{peak} fpeak的取值范围。由于ADC的采样率为 f s f_s fs,FFT点数为N,即是在频率轴上( f s 2   f s 2 \frac{f_s}{2}~ \frac{f_s}{2} 2fs 2fs)分成了N个 frequency bin,目标对应的 f p e a k f_{peak} fpeak就落在这其中的某个bin中。考虑实信号有镜像频率,如果我们只取正频率部分

f p e a k = f s N k ,  0 ≤ k < N 2 f_{peak} = \frac{f_s}{N}k, \ 0≤k<\frac{N}{2} fpeak=Nfsk 0k<2N

从之前的推导很容易看出, d t a r g e t d_{target} dtarget f p e a k f_{peak} fpeak 成线性关系,即:
d m a x = c 2 K f m a x d_{max} = \frac{c}{2K} f_{max} dmax=2Kcfmax

由于在正频率部分,最大的 f p e a k f_{peak} fpeak即是 f s 2 \frac{f_s}{2} 2fs,所以
d m a x = c 4 K f s d_{max} = \frac{c}{4K}f_s dmax=4Kcfs


距离分辨率

两个可分辨的 range bin 对应着两个可分辨的 frequency bin,即
Δ d = d 1 − d 2 = c 2 K ( f 1 − f 2 ) = c 2 K Δ f = c 2 K f s N = c ⋅ C R T 2 B W f s N \Delta d = d_1 -d _2 = \frac{c}{2K}(f_1-f_2)=\frac{c}{2K}\Delta f=\frac{c}{2K}\frac{f_s}{N}=\frac{c \cdot CRT }{2BW}\frac{f_s}{N} Δd=d1d2=2Kc(f1f2)=2KcΔf=2KcNfs=2BWcCRTNfs

在连续采样的背景下(此处N表示一个CRT中的采样点数),我们有:
N f s = C R T \frac{N}{f_s} = CRT fsN=CRT
于是 Δ d \Delta d Δd进一步可以化简为:
Δ d = c 2 B W \Delta d = \frac{c}{2 BW} Δd=2BWc

一个直观的结论是,能达到的距离分辨率只与 B W BW BW,即扫频带宽有关


无法整周期采样

但在实际上,由于一些装置的设置,往往我们并不会完成整个CRT周期内的采样(此处N表示一个CRT中的采样点数),即
N f s = α C P T ≤ C R T ,其中 α ≤ 1 \frac{N}{f_s}=\alpha CPT \le CRT ,其中 \alpha \le 1 fsN=αCPTCRT,其中α1
于是在这样的情况下,实际得到的距离分辨率并没有预想的那么精细,而要略微粗些。
Δ d = c 2 B W 1 α ≥ c 2 B W \Delta d = \frac{c}{2BW} \frac{1}{\alpha} \ge \frac{c}{2BW} Δd=2BWcα12BWc

虽然这里就有聪明的小朋友说了,我可以时域补零啊,那不是看得更仔细了。但遗憾的是“时域补零相当于频域插值”,这改善的只是FFT的是视觉分辨率,而实际上的波长分辨率只由采样数决定(即采样点数M)。假如现在有两个物体落在了同一个原始采样点数换算得到的某个 range bin 中,那么这两个物体在我们最终得到的FFT频谱图中就是无法完成区分的。


我们再对这个分辨率有一个感性的认识,假设扫频带宽为3GHz,则:
Δ d = c 2 B W = 3 × 1 0 8 2 × 3 × 1 0 9 m = 5 m m \Delta d= \frac{c}{2 BW}=\frac{3 \times 10^8}{2 \times 3 \times 10^9}m=5mm Δd=2BWc=2×3×1093×108m=5mm

即一个毫米级的距离分辨率。


CRT的影响

在 TI 关于FMCW雷达的教程视频中,有提出这样一个有趣的问题,即关于 chirp 1 和 chirp 2的讨论,它们具有相同的扫频带宽,但是chirp 1的周期是chirp 2的两倍。

在这里插入图片描述
根据之前的推导,我们知道这两个chirp对应的距离分辨率应该是一致的

但假如我们现在用同一台ADC进行采样:在采样频率相同时,chirp 1能获得的采样点数更多,其对应的波长分辨率也就更高。在采样点数相同的情况下,chirp 1 会比 chirp 2需要的采样率更低,于是对ADC的容忍度提高了;但同时,chirp 1也需要花费更长的时间完成采样。那么,这里呈现的就将是一个对于ADC的 trade-off

最后,用一张图总结 range FFT的内容:

在这里插入图片描述

  • 1
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 3
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值