数字信号处理4：采样定理

采样定理又称奈奎斯特定理，它规定了连续信号抽样成为离散信号时的采样率与信号最大频率之间的关系，那就是：

对于连续信号$x(t)$进行抽样时，抽样信号的最小频率$p(t)$的频率要大于$x(t)$的最大频率的2倍，采样得出的信号$x[n]$才能还原出原始信号$x(t)$。

这一点其实容易想明白，比如下面三个信号：

如果在T的整数倍处对三个信号采样，你会发现采样出的信号是完全一致的，这就导致你无法还原原始信号。因此，我们要明确采样频率和原信号的中包含的最高频率分量的关系。

通常来说，我们对于一个连续信号采样时，是截取$x(t)$在某一时刻的值，因此不妨将整个采样过程看做$x(t)$与$p(t)$的乘积，其中$p(t)$是梳妆函数，它是一系列冲激信号的延时叠加：
$$p(t)=\sum _{n=-\infty }^{+\infty }\delta (t-nT)$$

它的时域图像如下：

对它做傅里叶变换得到频域图像如下：

对应的表达式为

$$P(jw)=\sum _{k=-\infty }^{+\infty }\delta (w-k{w}_s)$$

而对于$x(t)$，我们需要假设它是一个带限信号，意味着它在频率上是有界的（否则就会发生频谱混叠）。它的频域图像如下：

以上就是我们对输入信号$x(t)$和采样信号$p(t)$的直观认识，接下来推导采样定理。

1. 卷积定理

之前我们有说过，时域卷积等于频域相乘，频域相乘等于时域卷积。这个结论对于我们推导采样定理以及建立对傅里叶变换后的频域的认识非常重要。因此先推导这个结论。

假设有信号
$$y(t)=x(t)\ast h(t)={\int }_{-\infty }^{+\infty }x(\tau )h(t-\tau )d\tau$$

直接对$y(t)$做傅里叶变换：
$$Y(jw)={\int }_{-\infty }^{+\infty }[{\int }_{-\infty }^{+\infty }x(\tau )h(t-\tau )d\tau ]e^{-jwt}dt$$
交换积分次序，且由于$x(\tau )$与$t$无关

$$={\int }_{-\infty }^{+\infty }x(\tau )[{\int }_{-\infty }^{+\infty }h(t-\tau )e^{-jwt}dt]d\tau$$

根据傅里叶变换的时移性质：
若$x(t)$的傅里叶变换为$X(jw)$，则$x(t-\tau )$的傅里叶变换为：
$${\int }_{-\infty }^{+\infty }x(t-\tau )e^{-jwt}dt={\int }_{-\infty }^{+\infty }x(t-\tau )e^{-jw(t-\tau )}{e}^{-jw\tau }dt=e^{-jw\tau }X(jw)$$

因此，上式中方括号内的部分就可以用时移定理处理。

$$={\int }_{-\infty }^{+\infty }x(\tau )[e^{-jw\tau }H(jw)]d\tau$$

且$H(jw)$与$\tau$无关，可移出积分外。

$$=H(jw){\int }_{-\infty }^{+\infty }x(\tau )e^{-jw\tau }d\tau =H(jw)X(jw)$$

故而
$$y(t)=x(t)\ast h(t)$$
的傅里叶变换为
$$Y(w)=X(jw)H(jw)$$

对于频域卷积的情况，由于情况是对称的，所以不再单独推导。

2. 采样定理推导

之前的结论是，对于采样过程，可以写成
$$x_p(t)=x(t)p(t)$$

那么，对于$x_p(t)$的傅里叶变换$X_p(jw)$，就是频域卷积$X(jw)\ast P(jw)$。

之前我们已经给出了$P(jw)$的图像。它在频域上也是冲激信号的延时叠加。

而对于一个信号，与冲激信号进行卷积，就相当于将该信号进行移位：
$$X(jw)\ast \delta (w-k{w}_s)=X(j(w-k{w}_s))$$
故而，$X_p(jw)$的图像就是$X(jw)$的一系列位移叠加。

图中，$w_s$是采样信号$p(t)$的频率，$w_s=2\pi /T$。$w_M$为原始信号中包含的最大频率分量。显然，从图中可以看到，仅当$w_s>2w_M$时，频谱才不会发生混叠，才能唯一还原成原信号$x(t)$。

3. 傅里叶变换的频域标度

通常在做傅里叶变换时，我们取N个点的时域信号，然后做DFT之后变成频域表示，同样有N个点。那么这N个点代表什么呢？

这里要回想一下之前的离散傅里叶变换

$$X[k]=\sum _{n=<N>}x[n]e^{-jk\frac{2\pi }{N}n}$$

显然$k$是以$N$为周期的。并且，$|X[-k]|=|X[-k]|$，这和上图中的c保持一致。

当$k=N$时，$k\frac{2\pi }{N}=2\pi$，从上面的图不难发现，这里的$2\pi$对应的就是采样率$w_s$，如果换算成Hz为单位的频率，它其实就代表的是我们通常在音频文件中所说的采样率，比如44100Hz。而$k=N/2$时，就代表采样率的一半。而由于频域图像的对称性，所以通常来说，真正表示频域信息的区域仅在$[0,\pi ]$内，这也是为什么通常我们从matlab中看到的频率分析定义域都是$[0,\pi ]$。

以通常的音频文件采样率举例，最常见的采样率就是44.1kHz，那么这个采样率能表示的信号的最高频率就是22.05kHz。而人耳能感受到的声音频率范围是20Hz~20kHz，因此，这个采样率可以说是专为人耳特定的。而由于现在存储媒介容量大幅增长，网络速度也越来越快，更高的采样率如48kHz、96kHz也逐渐开始流行。