采样定理和狄拉克梳

信息处理 采样 卷积
2021-12-21 05:40:49

我正在阅读“数字信号处理的科学家和工程师指南”并试图理解下面的图 3.5,它是关于采样定理和混叠的。

我不明白频域的图片。信号可以表示为 $$f(t)=\frac{1}{2\pi}\int_{-\infty}^{\infty}F(x)e^{ixt}dx$$ 其中 $F (x)$ 是 $f$ 的傅里叶变换。这表示具有“频率”$e^{ixt}$ 的 $f$,这些可以有 $x$ 负数。那么我想如何理解图 3-5,它显示了仅大于零且实数且不复数的频域。

另外,为什么脉冲序列的频域会出现原始频谱翻转一半的情况?我目前的理解是重复来自 $$\mathcal{F^{-1}} \{f(t)\Delta(t)\} = F(x)*\Delta(x)$$ 之类的东西$\Delta$ 是狄拉克梳子,可能在某处需要某种常数,如 $2\pi$。$F(x)*\Delta(x)$ 是如何产生频域中一半时间的翻转的。我认为卷积的图像是 $\Delta$ 被翻转,并且当它向右移动 $s$ 时,$F(x)\Delta(sx)$ 的“面积”是 $\Delta(sx)$ 处的卷积值新元。这似乎多次复制原始光谱,但没有翻转一半的时间,如图 3.5 所示。

在此处输入图像描述

2个回答

令 $X(f)$ 表示 $x(t)$ 的傅里叶变换,其中 $$\begin{align*} X(f) &= \int_{-\infty}^{\infty} x(t) \ exp(-j2\pi ft) \mathrm dt\\ x(t) &= \int_{-\infty}^{\infty} X(f) \exp(+j2\pi ft) \mathrm df \end{ align*}$$ 我将通过 $x(t) \leftrightarrow X(f)$ 表示。在接下来的内容中将需要以下变换对。$$\begin{align*} \delta(t) &\leftrightarrow 1\\ \delta(t-t_0) &\leftrightarrow \exp(-j2\pi f \,t_0)\\ \sum_{n=-\ infty}^{\infty}\delta(t-nT) &\leftrightarrow \frac{1}{T}\sum_{k=-\infty}^{\infty}\delta\left(f-\frac{k }{T}\right)\\ \end{对齐*}$$

\leftrightarrow \frac{1}{T}\sum_{k=-\infty}^{\infty} X\left(f-\frac{k}{T}\right) = \hat{X}(f) \end{align*}$$ 因此,以 $T$ 秒间隔对 $x(t)$ 进行采样形成的脉冲序列的傅里叶变换是 $\hat{X}(f)$,它是通过重复 $X 获得的(f)$ 沿 $f$ 轴以 $T^{-1}$ Hz 的间隔对结果求和。此外,$\hat{X}(f)$ 是频率变量 $f$ 的周期函数,周期为 $T^{-1}$ Hz。也就是说,对于所有 $f$,$$\hat{X}\left(f + \frac{1}{T}\right) = \hat{X}(f).$$ 注意所有这些都成立不管 $X(f)$ 是什么:对于所有 $f$,$X(f)$ 都可以是非零的,并且结果仍然有效。

现在假设如果 $f < a$ 或 $f > b$,则 $X(f)$ 为零。然后,通过沿 $f$ 轴周期性重复 $X(f)$ 获得的 $\hat{X}(f)$ 仅在区间 $$\ldots, \left[a-\frac{2 }{T}, b-\frac{2}{T}\right], \left[a-\frac{1}{T}, b-\frac{1}{T}\right], [a, b] \left[a+\frac{1}{T}, b+\frac{1}{T}\right], \left[a+\frac{2}{T}, b+\frac{2}{T} \right], \ldots $$ 等如果 $$b -\frac{1}{T} < a \Rightarrow ba < \frac{1}{T},$$ 即支持$X(f)$的$ba$小于重复间隔$T^{-1}$Hz,则$X(f)$的重复不重叠。实际上,正如 OP 的数据所示,对于频谱从 $a = -f_0$ 到 $b = f_0$ 的实值信号($X(f)$ 的支持长度为 $2f_0$),采样 $x( t)$ 以 $T = (3f_0)^{-1}$ 的间隔(因此在频率轴上以 $3f_0$ Hz 的间隔重复 $X(f)$)导致在采样 $x(t )$ 以 $T = (1.5f_0)^{-1}$ 的间隔导致 $X(f)$ 的重复重叠。

最后,OP询问负频率和仅显示正频率的图片的解释。对于实值信号,$X(f)$ 具有共轭对称性,这意味着 $X(-f) = X^*(f)$,因此为 $f$ 的正值指定 $X(f)$就够了。无论如何,他正在查看的图片是 $|X(f)|$ 和 $|\hat{X}(f)|$ 的偶数函数,因此仅显示正轴可以节省空间,虽然它确实使图片看起来有点不平衡,因为只有一半的波瓣在低频下显示。对于复值信号,$X(f)$ 不具有共轭对称性并且 $|X(f)|$ 不必是 $f$ 的偶函数,因此需要显示整个轴。但是,上面的一般发展仍然适用,我们仍然需要以超过 $ba$ Hz 的速率进行采样,这有助于记住,在这种一般情况下,每个样本实际上是两个实数,而不是一个,因为我们正在对复值信号进行采样。

对于实值信号 $f(t)$,您显示的 $|\mathcal{F}\{f\}(x)|$ 的图在实轴上是对称的,因为 $\mathcal{F}\{ f\}(-x) = \overline{\mathcal{F}\{f\}(x)}$(等式右边的上划线表示复共轭)。

我认为,这也解释了您观察到的“翻转”。想象一下图中右上方的子图反映在 x 轴上;是由于有限采样而被复制的频谱

编辑:

$\begin{align*} \mathcal{F}\{f\}(-x) &= \int_{t=-\infty}^{t=\infty} f(t) e^{-2\pi它(-x)} dt \\ &= \int_{t=-\infty}^{t=\infty} f(t) \cos(2\pi t(-x)) dt - i\int_{t =-\infty}^{t=\infty} f(t) \sin(2\pi t(-x)) dt \\ &= \overline{\int_{t=-\infty}^{t=\ infty} f(t) \cos(2\pi t(-x)) dt + i\int_{t=-\infty}^{t=\infty} f(t) \sin(2\pi t(- x)) dt} \\ &= \overline{\int_{t=-\infty}^{t=\infty} f(t) e^{2\pi it(-x)} dt} \\ &= \overline{\int_{t=-\infty}^{t=\infty} f(t) e^{-2\pi i tx} dt} \\ &= \overline{\mathcal{F}\{f \}(x)} \end{对齐*}$