承接上篇博文 中频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=2BWc⋅CRTfpeak
其中,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, 0≤k<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=d1−d2=2Kc(f1−f2)=2KcΔf=2KcNfs=2BWc⋅CRTNfs
在连续采样的背景下,我们有:
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
f
s
=
α
C
P
T
≤
C
R
T
,其中
α
≤
1
\frac{N}{f_s}=\alpha CPT \le CRT ,其中 \alpha \le 1
fsN=αCPT≤CRT,其中α≤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α1≥2BWc
虽然这里就有聪明的小朋友说了,我可以时域补零啊,那不是看得更仔细了。但遗憾的是“时域补零相当于频域插值”,这改善的只是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的内容: