连续信号的采样与恢复实验报告.pdfVIP

实验六、连续信号的采样与恢复

一、实验目的

1.加深理解采样对信号的时域和频域特性的影响；

2.加深对采样定理的理解和掌握，以及对信号恢复的必要性；

3.掌握对连续信号在时域的采样与重构的方法。

二、实验原理

(1)信号的采样

信号的采样原理图如下图所示，其数学模型表示为：

=

其中的f(t)为原始信号，为理想的开关信号（冲激采样信号）δTs(t)=，fs(t)为采样

后得到的信号称为采样信号。由此可见，采样信号在时域的表示为无穷多冲激函数的线性组合，其权值为

原始信号在对应采样时刻的定义值。

令原始信号f(t)的傅立叶变换为F(jw)=FT(f(t))，则采样信号fs(t)的傅立叶变换

Fs(jw)=FT(fs(t))=。由此可见，采样信号fs(t)的频谱就是将原始信号f(t)的频谱在频

率轴上以采样角频率ws为周期进行周期延拓后的结果（幅度为原频谱的1/Ts）。如果原始信号为有限带宽

的信号，即当|w||wm|时，有F(jw)=0，则有：如果取样频率ws≥2wm时，频谱不发生混叠；否则会出现频

谱混叠。

(2)信号的重构

设信号f(t)被采样后形成的采样信号为fs(t)，信号的重构是指由fs(t)经过内插处理后，恢复出原来的信号f(t)

的过程。因此又称为信号恢复。

由前面的介绍可知，在采样频率ws≥2wm的条件下，采样信号的频谱Fs(jw)是以ws为周期的谱线。选择一

个理想低通滤波器，使其频率特性H(jw)满足：

H(jw)=

式中的wc称为滤波器的截止频率，满足wm≤wc≤ws/2。将采样信号通过该理想低通滤波器，输出信号的频

谱将与原信号的频谱相同。因此，经过理想滤波器还原得到的信号即为原信号本身。信号重构的原理图见

下图。

通过以上分析，得到如下的时域采样定理：一个带宽为wm的带限信号f(t)，可唯一地由它的均匀取样信号

fs(nTs)确定，其中，取样间隔Tsπ/wm,该取样间隔又称为奈奎斯特（Nyquist）间隔。

根据时域卷积定理，求出信号重构的数学表达式为：

式中的抽样函数Sa(wct)起着内插函数的作用，信号的恢复可以视为将抽样函数进行不同时刻移位后加权求

和的结果，其加权的权值为采样信号在相应时刻的定义值。利用MATLAB中的抽样函数

来表示Sa(t)，有，，于是，信

号重构的内插公式也可表示为：

(3)模拟低通滤波器的设计

在任何滤波器的设计中，第一步是确定滤波器阶数N及适当的截止频率Ωc。对于巴特沃斯滤波器，可使

用MATLAB命令buttord来确定这些参数，设计滤波器的函数为butter。其调用形式为

[N,wn]=buttord(wp,ws,rp,rs,s)

[b,a]=butter[N,wn,s]

其中，wp、ws、rp、rs为待设计滤波器的技术指标，分别代表通带截止频率、阻带截止频率、通带最大衰

减和阻带最小衰减；s表示设计滤波器的类型为模拟滤波器；N、wn为设计得到的滤波器的阶数和3dB截

止频率，b、a为滤波器系统函数的分子和分母多项式的系数矢量，假定系统函数的有理分式表示为：

三、程序示例

示例1：选取门信号f(t)=g2(t)为被采样信号。利用MATLAB实现对信号f(t)的采样，显示原信号与采样信

号的时域和频域波形。

因为门信号并非严格意义上的有限带宽信号，但是，由于其频率f1/τ的分量所具有的能量占有很少的比重，

所以一般定义fm=1/τ为门信号的截止频率。其中的τ为门信号在时域的宽度。在本例中选取fm=0.5，临界

采样频率为fs=1，过采样频率为fs1（为了保证精度，可以将其值提高到该值的50倍），欠采样频率为

fs1。

%显示原信号及其Fourier变换示例

R=0.01;%采样周期

t=-4:R:4;

f=rectpuls(t,2);

w1=2*pi*10;%显示从-20*pi到20*pi频率范围内的频谱

N=1000;%计算出2*1000+1个频率点的值

k=0:N;

wk=k*w1/N;

F