傅里叶变换:
$$ X(j 2 \pi f) = \mathscr{F}\Big\{x(t)\Big\} \triangleq \int\limits_{-\infty}^{+\infty} x( t) \ e^{-j 2 \pi ft} \ \text{d}t $$
傅里叶逆变换:
$$ x(t) = \mathscr{F}^{-1}\Big\{X(j 2 \pi f)\Big\} = \int\limits_{-\infty}^{+ \infty} X(j 2 \pi f) \ e^{j 2 \pi ft} \ \text{d}f $$
矩形脉冲函数:
$$ \operatorname{rect}(u) \triangleq \begin{cases} 0 & \mbox{if } |u| > \frac{1}{2} \\ 1 & \mbox{if } |u| < \frac{1}{2} \\ \end{cases} $$
"Sinc" 函数 ("sinus cardinalis") :
$$ \operatorname{sinc}(v) \triangleq \begin{cases} 1 & \mbox{if } v = 0 \\ \frac{\sin(\pi v) }{\pi v} & \mbox{if } v \ne 0 \\ \end{cases} $$
定义采样频率\$ f_\text{s} \triangleq \frac{1}{T} \$作为采样周期\$T\$的倒数。
注意:
$$ \mathscr{F}\Big\{\operatorname{rect}\left( \tfrac{t}{T} \right) \Big\} = T \ \operatorname{sinc}(fT) = \ frac{1}{f_\text{s}} \ \operatorname{sinc}\left( \frac{f}{f_\text{s}} \right)$$
狄拉克梳子(又名“采样函数”又名“沙函数”):
$$ \operatorname{III}_T(t) \triangleq \sum\limits_{n=-\infty}^{+\infty} \delta(t - nT) $$
狄拉克梳是周期性的,周期为 \$T\$。傅里叶级数:
$$ \operatorname{III}_T(t) = \sum\limits_{k=-\infty}^{+\infty} \frac{1}{T} e^{j 2 \pi k f_\text{s } t} $$
采样的连续时间信号:
$$ \begin{align} x_\text{s}(t) & = x(t) \cdot \left( T \cdot \operatorname{III}_T(t) \right) \\ & = x(t) \cdot \left( T \cdot \sum\limits_{n=-\infty}^{+\infty} \delta(t - nT) \right) \\ & = T \ \sum\limits_{n=-\ infty}^{+\infty} x(t) \ \delta(t - nT) \\ & = T \ \sum\limits_{n=-\infty}^{+\infty} x(nT) \ \delta (t - nT) \\ & = T \ \sum\limits_{n=-\infty}^{+\infty} x[n] \ \delta(t - nT) \\ \end{align} $$
其中\$ x[n] \triangleq x(nT) \$。
这意味着\$x_\text{s}(t)\$仅由样本\$x[n]\$和采样周期\$T\$ 定义,并且完全丢失了\$值的任何信息x(t)\$表示采样实例之间的时间。 \$x[n]\$是一个离散的数字序列,是\$x_n\$的一种 DSP 速记符号。虽然对于\$ nT < t < (n+1)T \$确实\$x_\text{s}(t) = 0\ $ ,但对于任何\ $x[n]\$的值$n\$不是整数是未定义的。
注意:离散信号\$x[n]\$及其上的所有离散时间运算,如\$\mathcal{Z}\$ -Transform、离散时间傅里叶变换 (DTFT)、离散傅里叶变换(DFT),对于采样频率或采样周期\$T\$是“不可知的”。一旦您进入离散时间\$x[n]\$域,您就不会知道(或关心)\$T\$。只有使用Nyquist-Shannon 采样和重建定理,才能将\ $x[n]\$和\$T\$放在一起。
\$x_\text{s}(t)\$的傅里叶变换为
$$ \begin{align} X_\text{s}(j 2 \pi f) \triangleq \mathscr{F}\{ x_\text{s}(t) \} & = \mathscr{F}\left\ {x(t) \cdot \left( T \cdot \operatorname{III}_T(t) \right) \right\} \\ & = \mathscr{F}\left\{x(t) \cdot \left ( T \cdot \sum\limits_{k=-\infty}^{+\infty} \frac{1}{T} e^{j 2 \pi k f_\text{s} t} \right) \right \} \\ & = \mathscr{F}\left\{ \sum\limits_{k=-\infty}^{+\infty} x(t) \ e^{j 2 \pi k f_\text{s } t} \right\} \\ & = \sum\limits_{k=-\infty}^{+\infty} \mathscr{F}\left\{ x(t) \ e^{j 2 \pi k f_\text{s} t} \right\} \\ & = \sum\limits_{k=-\infty}^{+\infty} X\left(j 2 \pi (f - k f_\text{s })\right) \\ \end{对齐} $$
关于缩放 的重要说明:采样函数\$ T \cdot \operatorname{III}_T(t) \$和采样信号\$x_\text{s}(t)\$的因子为\$T\$你几乎不会在所有教科书中看到。这是这些教科书作者的一个教学错误,原因有多种(相关):
- 首先,省略\$T\$会改变采样信号\$x_\text{s}(t)\$的维度,而不是采样信号的维度\$x(t)\$。
- 信号链中的某处将需要该\$T\$因子。这些将其排除在采样函数之外的教科书最终将其放入采样定理的重建部分,通常作为重建滤波器的通带增益。 这在维度上令人困惑。有人可能会问:“我如何设计一个通带增益为\$T\$的砖墙 LPF ?”
- 如下所示,将\$T\$留在这里会导致零阶保持 (ZOH) 的净传递函数和净频率响应出现类似的缩放误差。 我见过的所有关于数字(和混合)控制系统的教科书都犯了这个错误,这是一个严重的教学错误。
请注意,\$x[n]\$的 DTFT 和采样信号\$x_\text{s}(t)\$的傅里叶变换在适当缩放后实际上是相同的:
DTFT:
$$ \begin{align} X_\mathsf{DTFT}(\omega) & \triangleq \mathcal{Z}\{x[n]\} \Bigg|_{z=e^{j\omega}} \\ & = X_\mathcal{Z}(e^{j\omega}) \\ & = \sum\limits_{n=-\infty}^{+\infty} x[n] \ e^{-j \omega n} \\ \end{对齐} $$
可以证明
$$ X_\mathsf{DTFT}(\omega) = X_\mathcal{Z}(e^{j\omega}) = \frac{1}{T} X_\text{s}(j 2 \pi f) \Bigg|_{f=\frac{\omega}{2 \pi T}} $$
无论\$x(t)\$是否“正确采样” ,上述数学都是正确的。如果\$x(t)\$可以从样本 \$x[n]\$和采样率或采样周期的知识中完全恢复,则\$x(t)\$是“正确采样的” 。采样定理告诉我们从\$x[n]\$和\$T\$恢复或重建\$x(t)\$需要什么。
如果\$x(t)\$的带宽限制为某个带宽限制\$B\$,这意味着
$$ X(j 2 \pi f) = 0 \quad \quad \text{for all} \quad |f| > 乙$$
考虑由原始偏移图像组成的采样信号的频谱:
$$ X_\text{s}(j 2 \pi f) = \sum\limits_{k=-\infty}^{+\infty} X\left(j 2 \pi (f - k f_\text{s })\右) $$
原始光谱\$X(j 2 \pi f)\$可以从采样光谱\$X_\text{s}(j 2 \pi f)\ $ 中恢复,如果没有移位的图像,\$X\ left(j 2 \pi (f - k f_\text{s})\right)\$,重叠它们的相邻邻居。这意味着第\$k\$图像的右边缘(即\$X\left(j 2 \pi (f - k f_\text{s})\right)\$)必须完全第 ( \$k+1\$ ) 图像的左边缘的左侧(即\$X\left(j 2 \pi (f - (k+1) f_\text{s})\right )\$ )。数学重述,
$$ k f_\text{s} + B < (k+1) f_\text{s} - B $$
这相当于
$$ f_\text{s} > 2B $$
如果我们以超过两倍带宽的采样率进行采样,则没有任何图像重叠,原始光谱\$X(j 2 \pi f)\$,也就是\$k=0\$的图像从\$X_\text{s}(j 2 \pi f)\$中提取,使用砖墙低通滤波器保持原始图像(其中\$k=0\$)未缩放并丢弃所有其他图像。这意味着它将原始图像乘以 1,并将所有其他图像乘以 0。
$$ \begin{align} X(j 2 \pi f) & = \operatorname{rect}\left( \frac{f}{f_\text{s}} \right) \cdot X_\text{s}( j 2 \pi f) \\ & = H(j 2 \pi f) \ X_\text{s}(j 2 \pi f) \\ \end{align} $$
重构滤波器是
$$ H(j 2 \pi f) = \operatorname{rect}\left( \frac{f}{f_\text{s}} \right) $$
并具有非因果 脉冲响应:
$$ h(t) = \mathscr{F}^{-1} \{H(j 2 \pi f)\} = f_\text{s} \operatorname{sinc}(f_\text{s}t) $$
这种滤波操作,表示为频域中的乘法,相当于时域中的卷积:
$$ \begin{align} x(t) & = h(t) \circledast x_\text{s}(t) \\ & = h(t) \circledast T \ \sum\limits_{n=-\infty }^{+\infty} x[n] \ \delta(t-nT) \\ & = T \ \sum\limits_{n=-\infty}^{+\infty} x[n] \ (h( t) \circledast \delta(t-nT) ) \\ & = T \ \sum\limits_{n=-\infty}^{+\infty} x[n] \ h(t-nT)) \\ & = T \ \sum\limits_{n=-\infty}^{+\infty} x[n] \ \left(f_\text{s} \operatorname{sinc}(f_\text{s}(t-nT )) \right) \\ & = \sum\limits_{n=-\infty}^{+\infty} x[n] \ \operatorname{sinc}(f_\text{s}(t-nT)) \ \ & = \sum\limits_{n=-\infty}^{+\infty} x[n] \ \operatorname{sinc}\left( \frac{t-nT}{T}\right) \\ \end {对齐} $$
这明确说明了原始的\$x(t)\$是如何从样本\$x[n]\$和采样率或采样周期的知识中重建的。
因此,实际数模转换器 (DAC)的输出既不是
$$ \sum\limits_{n=-\infty}^{+\infty} x[n] \ \operatorname{sinc}\left( \frac{t-nT}{T}\right) $$
不需要额外的治疗来恢复\$x(t)\$,也
$$ x_\text{s}(t) = \sum\limits_{n=-\infty}^{+\infty} x[n] \ T \delta(t-nT) $$
其中,使用理想的砖墙 LPF通过隔离和保留基带图像并丢弃所有其他图像来恢复\$x(t)\$ 。
如果没有对数字化信号进行处理或缩放,传统 DAC 的结果是值\$x[n]\$保持在一个恒定值,直到输出下一个样本。这会产生一个分段常数函数:
$$ x_\text{DAC}(t) = \sum\limits_{n=-\infty}^{+\infty} x[n] \ \运算符名{rect}\left(\frac{t-nT - \ frac{T}{2}}{T} \right) $$
请注意应用于\$\operatorname{rect}(\cdot)\$函数的\$\frac{1}{2}\$样本周期的延迟。这使它成为因果关系。它的意思很简单
$$ x_\text{DAC}(t) = x[n] = x(nT) \quad \quad \text{when} \quad nT \le t < (n+1)T $$
不同的说法
$$ x_\text{DAC}(t) = x[n] = x(nT) \quad \quad \text{for} \quad n = \operatorname{floor}\left( \frac{t}{T} \右)$$
其中\$\operatorname{floor}(u) = \lfloor u \rfloor\$是楼层函数,定义为不超过\$u\$的最大整数。
此 DAC 输出直接建模为接受理想采样信号\$x_\text{s}(t)\$的线性时不变系统 (LTI)或滤波器,并且对于理想采样信号中的每个脉冲,输出此脉冲回复:
$$ h_\text{ZOH}(t) = \frac{1}{T} \operatorname{rect}\left(\frac{t - \frac{T}{2}}{T} \right) $$
插入以检查此...
$$ \begin{align} x_\text{DAC}(t) & = h_\text{ZOH}(t) \circledast x_\text{s}(t) \\ & = h_\text{ZOH}(t ) \circledast T \ \sum\limits_{n=-\infty}^{+\infty} x[n] \ \delta(t-nT) \\ & = T \ \sum\limits_{n=-\infty }^{+\infty} x[n] \ (h_\text{ZOH}(t) \circledast \delta(t-nT) ) \\ & = T \ \sum\limits_{n=-\infty}^ {+\infty} x[n] \ h_\text{ZOH}(t-nT)) \\ & = T \ \sum\limits_{n=-\infty}^{+\infty} x[n] \ \frac{1}{T} \operatorname{rect}\left(\frac{t - nT - \frac{T}{2}}{T} \right) \\ & = \sum\limits_{n=- \infty}^{+\infty} x[n] \ \operatorname{rect}\left(\frac{t - nT - \frac{T}{2}}{T} \right) \\ \end{align } $$
DAC 输出\$x_\text{DAC}(t)\$作为具有脉冲响应\$h_\text{ZOH}(t)\$的 LTI 系统的输出,与上面的分段常数构造一致。这个 LTI 系统的输入是采样信号\$x_\text{s}(t)\$明智地缩放,使得 \$x_\text{s}(t)\$ 的基带图像与被采样的原始信号的频谱\$x(t)\$。那是
$$ X(j 2 \pi f) = X_\text{s}(j 2 \pi f) \quad \quad \text{for} \quad -\frac{f_\text{s}}{2} < f < +\frac{f_\text{s}}{2} $$
原始信号频谱与采样频谱相同,但丢弃所有因采样而出现的图像。
这个 LTI 系统的传递函数,我们称之为零阶保持(ZOH),是脉冲响应的拉普拉斯变换:
$$ \begin{align} H_\text{ZOH}(s) & = \mathscr{L} \{ h_\text{ZOH}(t) \} \\ & \triangleq \int\limits_{-\infty} ^{+\infty} h_\text{ZOH}(t) \ e^{-st} \ \text{d}t \\ & = \int\limits_{-\infty}^{+\infty} \frac {1}{T} \operatorname{rect}\left(\frac{t - \frac{T}{2}}{T} \right) \ e^{-st} \ \text{d}t \\ & = \int\limits_0^T \frac{1}{T} \ e^{-st} \ \text{d}t \\ & = \frac{1}{T} \quad \frac{1}{ -s}e^{-st}\Bigg|_0^T \\ & = \frac{1-e^{-sT}}{sT} \\ \end{align}$$
通过代入\$ j 2 \pi f \rightarrow s \$得到频率响应
$$ \begin{align} H_\text{ZOH}(j 2 \pi f) & = \frac{1-e^{-j2\pi fT}}{j2\pi fT} \\ & = e^{ -j\pi fT} \frac{e^{j\pi fT}-e^{-j\pi fT}}{j2\pi fT} \\ & = e^{-j\pi fT} \frac{ \sin(\pi fT)}{\pi fT} \\ & = e^{-j\pi fT} \operatorname{sinc}(fT) \\ & = e^{-j\pi fT} \operatorname{ sinc}\left(\frac{f}{f_\text{s}}\right) \\\end{align}$$
这表明线性相位滤波器具有二分之一采样周期的恒定延迟\$\frac{T}{2}\$,并且增益随着频率\$f\$的增加而减小。这是一种温和的低通滤波器效果。在 DC,\$f=0\$,增益为 0 dB,而在 Nyquist,\$f=\frac{f_\text{s}}{2}\$,增益为 -3.9224 dB。所以基带图像的一些高频分量减少了一点。
与采样信号\$x_\text{s}(t)\$一样,采样信号\$x_\text{DAC}(t)\$中存在采样频率整数倍的图像,但这些图像是幅度显着降低(与基带图像相比),因为\$|H_\text{ZOH}(j 2 \pi f)|\$当\$f = k\cdot f_\text{s}\$时通过零对于不为 0 的整数\$k\$,它位于这些图像的中间。
结论:
零阶保持 (ZOH) 是信号重建的线性时不变模型,由实际的数模转换器 (DAC) 完成,它将输出恒定保持在样本值\$x[n]\$ , 直到被下一个样本\$x[n+1]\$更新。
与常见的误解相反,ZOH与模数转换器 (ADC)之前的采样保持电路 (S/H)无关。只要 DAC 在每个采样周期内将输出保持在恒定值,无论 ADC 是否具有 S/H,ZOH 效应都会保持。如果 DAC 输出的不是上面描述为\$x_\text{DAC}(t)\$的分段常数输出(例如旨在近似狄拉克脉冲的窄脉冲序列),则不存在ZOH 效应(取而代之的是)在 ADC 之前是否有 S/H 电路。
ZOH 的净传递函数为$$ H_\text{ZOH}(s) = \frac{1-e^{-sT}}{sT} $$,ZOH 的净频率响应为$$ H_\ text{ZOH}(j 2 \pi f) = e^{-j\pi fT} \operatorname{sinc}(fT) $$ 许多教科书忽略了传递函数分母中的\$T\$因子,而这是一个错误。
ZOH显着减少了采样信号\$x_\text{s}(t)\$的图像,但并未消除它们。为了消除图像,需要像以前一样使用良好的低通滤波器。Brickwall LPF 是一种理想化。实际的 LPF 也可能会衰减高频下的基带图像(我们想要保留),并且必须将这种衰减与 ZOH 产生的衰减(小于 3.9224 dB 衰减)一样考虑在内。ZOH 还会将信号延迟半个采样周期,这可能需要考虑(连同抗成像 LPF 的延迟),尤其是在 ZOH 处于反馈环路中时。