只要满足了采样定理,就不会有任何信息丢失,剩下的问题仅仅是如何重构而已,而重构过程跟采样过程一样同样不允许出现混叠,从频域上看,就是一个矩形低通滤器必须准准地截止在【采样频率】/2上(少了可能会产生信息损失,超过了可能会产生重构混叠),将这个频域的矩形返回到时域看,就是个SINC函数。因此重构函数的选择并非跟对未知趋势规律的实验数据做Curve Fitting一样可以比较随意,只要光滑地通过实验点就行了,而是必须基于SINC函数,这一唯一理论正确的重构函数,SINC函数的现实问题仅仅是它是无限长的而无法实现,不过好在其在无限远的地方衰减到零,因此只要有一定长度,截断造成的误差就会很小。实际上重构滤波函数的选择都是基于或逼近理想SINC函数的,一般也就是对理想SINC函数截短加窗,当然模拟滤波器更受限于其因果性局限,其他一些偏离SINC函数的滤波器都有重构原波形以外的其他目的。
为了避免将SINC插值理解为正弦波插值,下面给出一个非简单正弦波的例子。信号为20kHz内带限信号,包括20kHz, 19kHz, 7kHz和1.1kHz分量。
(1)下图的采样率为44.1kHz,采样频率与信号最高频率比为44.1/20=2.2倍,直线(插值)连接采样点,显示出来的时域波形不令人满意,如果把这些线性插值点也包括进来做FFT,将有很多带限范围外的高频。
(2)下图的采样率为44.1kHz,采样频率与信号最高频率比为44.1/20=2.2倍,SINC(插值)连接采样点,显示出来的时域波形令人满意。
(3)下图的采样率为200kHz,采样频率与信号最高频率比为200/20=10倍,直线(插值)连接采样点,显示出来的时域波形基本令人满意,但仍然可见一些突兀转折。
比较(2)和(3)的波形可见,尽管(3)的采样率是(2)的大约5倍,但由于(3)采用了Lousy的线性插值,其波形仍然不如(2)的采用SINC插值出来的波形更接近真实情况。
|