数字信号处理 --- 信号的采样和奇妙的混叠（Aliasing） 壹

信号的采样与混叠(时域)

     在一定条件下，一个连续的时间信号完全可以用该信号在等时间间隔上的样本来表示，并且可以用这些样本的值把该信号完全恢复出来。这一非常重要的发现被称之为采样定理。

     采样定理在连续时间信号（一个是连续的波形）和离散时间信号（一个是数组）之间架起了一座桥梁！！！（就好像微积分基本定理一样，把斜率和面积这两个毫不相干的概念紧密的联系在了一起。）

一维信号的采样：

     通过一个固定周期的冲击串来实现信号的采样（即，连续信号的离散化）。如下图，一个本来是连续的一维信号被分割成了很多小块，这非常像黎曼积分（Riemann integral），定积分作黎曼和的极限！

上图中原来是连续的数据通过等间隔（采样周期为T）的采样，被分割成了等间隔，具有不同幅度的冲击串，见下图。

附：下图为黎曼和的极限，图中矩形的宽度就好比是采样间隔T！

二维信号的采样：

     二维信号的采样常被用于图像的数字化，相对于一维采样，他只是增加了采样的维度。不仅在X方向采样还在Y方向采样，采样后的结果是一个二维矩阵，矩阵中的每一个点代表了该位置的响应灰度。   

      例如，在出版业里，一幅图是由很多很小的采样点构成的，如果这些采样点的间距很小，小到你无法用肉眼分辨的话。（或者说，即便报纸上的采样点间距不是很小，但是读者眼睛距离报纸不是很近的话。）那么，报纸上的照片在我们看来就是一幅连续的图片。但是，如果你拿着放大镜或者眼睛靠着很近观察的话，你就会发现报纸上的一幅图像其实是由一个个等间隔的离散化的点拼出来的。

        下面我们用Matlab软件来仿真信号的采样，以及如果采样频率（采样周期的倒数）不能满足一定要求是所带来的混叠。

Tips : 本文的后面会有一个录像说明，现在的MATLAB仿真和录像说明有一个很大的不同就是。这里的仿真是保持原始信号的频率Fo不变，不断地改变采样频率Fs.

---

第一种情况 ------ 采样频率高于原始信号频率最高频的两倍（Fs > 2*Fo）

        下图所绘的是一个频率为60Hz的连续的余弦波，图中的红点为采样点。当采样冲击串的频率（2000Hz）远高于原始信号时（原始信号的33.3倍），采样点在曲线上的分布非常密集，重建后的效果也很好。如果采样频率（400Hz）只是略高于原始信号频率

（原始信号的6.6倍）,则原曲线上的采样点明显不足，但不用担心，即便是这种看起来不太充足的采样还是可以重建原始的连续信号的。

（点击图片放大）

Matlab代码： 

%% signal sampling and aliasing demo for CSDN
% Created: Late spring, 2018. (2018/05/03)
% Author: Z.Zhu, zzz0029@tigermail.auburn.edu
% Copy Rights Reserved.
% May not be copied, scanned, or duplicated, in whole or in part.
% Reference of algorithm by: Steve Eddins's Aliasing and