（规格为A4纸或A3纸折叠）
一、实验目的 1、了解信号的抽样方法与过程以及信号恢复的方法。 2、验证抽样定理。 二、实验原理 1.连续时间信号的离散化与抽样 从连续时间信号转化为离散时间信号，需要通过抽样来完成，以便计算机和其他系统进行处理。离散时间信号可以从离散信号源获得，也可以从连续时间信号抽样而得。抽样信号fs(t)可以看成连续信号f(t)和一组开关函数s(t)的乘积fs(t)?f(t)?s(t)。s(t)是一组周性窄脉冲，见图4-1，Ts称为抽样周期，其倒数fs?1称抽样频率。 TsF(t)E0?1/2(a)S(t)?13?1/2tAI-4TS-TS0TS4TS(b)8TS12TStFS(t)??1/20?1/2(c)?13?1/2t 图 4-1矩形抽样脉冲 对抽样信号进行傅里叶分析可知，抽样信号的频率包括了原连续信号以及无限个经过平移的原信号频率。平移的频率等于抽样频率fs及其谐波频率2fs、3fs……。当抽样信号是周性窄脉冲时，平移后的频率幅度按?sinx?规律衰减。抽样信号的频谱是原信号频谱周期的延拓，它占有的频带要比原信号频谱x宽得多。 2.信号恢复 正如测得了足够的实验数据以后，我们可以在坐标纸上把一系列数据点连起来，得到一条光滑的曲线一样，抽样信号在一定条件下也可以恢复到原信号。只要用一截止频率等于原信号频谱中最高频率fn的低通滤波器，滤除高频分量，经滤波后得到的信号包含了原信号频谱的全部内容，故在低通滤波器输出可以得到恢复后的原信号。 3.抽样定理 抽样定理指出，一个有限频宽的连续时间信号f(t)，其最高频率为?m，经过等间隔抽后只要抽样频率?s不小于信号最高频率的两倍，即满足?s?2?m，就能从抽样信号fs(t)中恢复出原信号，而fmin?2fm称为最低抽样频率又称“奈奎斯特抽样率”。抽样定理图形解释如下。 当fs?2fm时，抽样信号的频谱会发生混叠，从发生混叠后的频谱中我们无法用低通滤波器获得原信号频谱的全部内容。在实际使用中，仅包含有限频率的信号是极少的。因此即使fmin?2fm，恢复后的信号失真还是难免的。图6-2画出了当抽样频率fs?2fm（不混叠时）当抽样频率fs?2fm（混叠时）两种情况下冲激抽样信号的频谱。 三、实验内容及步骤 1.信号的抽样 1.1 正弦信号的采样 （1） 参考下面程序，得到50Hz正弦信号在采样时间间隔分别为0.01s、0.002s和0.001时的采样信号。 （2）在(1)基础上恢复正弦信号，比较那个采样间隔能较好的恢复原正弦信号。改变几个不同的采样间隔，比较恢复信号。 2.抽样信号的恢复 设信号f(t)?Sa(t)?sint，在抽样间隔分别为 t(1) Ts?0.7?(令?m?1，?c?1.1?m) (2) Ts?1.5?(令?m?1，?c?1.1?m) 的两种情况下， 对信号f(t)进行采样， 试编写MATLAB 程序代码， 并绘制出抽样信号波形、由抽样信号得到的恢复信号波形。 选取信号f(t)?Sa(t)?sint作为被抽样的信号，显然，信号的带宽?m?1。当采样频率?s?2?mt时，被称为临界抽样(取?c??m)。在临界抽样状态下实现对信号Sa(t)的抽样及由该抽样信号恢复Sa(t)的参考程序如下： 图 4-3 Sa(t)信号抽样 四、实验报告要求 1、简述实验目的和实验原理。 2、编程实现实验内容要求附上源代码及运行结果图示。 3、总结实验中的主要结论、收获和体会。 五、实验结果： 1.信号的抽样 1.1 正弦信号的采样 (1) 参考下面程序，得到50Hz正弦信号在采样时间间隔分别为0.01s、0.002s和0.001时的采样信号。 代码： fs=1000; t=0:1/fs:0.2; f0=50; x=sin(2*pi*f0*t); subplot(2,2,1); n1=0:0.01:0.2; plot(t,x); x1=sin(2*pi*f0*n1); subplot(2,2,2); stem(n1,x1); n2=0:0.005:0.2; x2=sin(2*pi*f0*n2); subplot(2,2,3); stem(n2,x2); n3=0:0.001:0.2; x3=sin(2*pi*f0*n3); subplot(2,2,4); stem(n3,x3,'.'); (2) 在(1)基础上恢复正弦信号，比较那个采样间隔能较好的恢复原正弦信号。改变几个不同的采样间隔，比较恢复信号。 代码： f0=50; n1=0:0.01:0.2; x1=sin(2*pi*f0*n1); n2=0:0.002:0.2; x2=sin(2*pi*f0*n2); n3=0:0.001:0.2; x3=sin(2*pi*f0*n3); subplot(3,3,1); stem(n1,x1); subplot(3,3,4); plot(n1,x1); subplot(3,3,2); stem(n2,x2); subplot(3,3,5); plot(n2,x2); subplot(3,3,3); stem(n3,x3); subplot(3,3,6); plot(n3,x3); 2.抽样信号的恢复 设信号f(t)?Sa(t)?sint，在抽样间隔分别为 t(1) Ts?0.7?(令?m?1，?c?1.1?m) (2) Ts?1.5?(令?m?1，?c?1.1?m) 的两种情况下， 对信号f(t)进行采样， 试编写MATLAB 程序代码， 并绘制出抽样信号波形、由抽样信号得到的恢复信号波形。 代码： （1） wm=1;%信号带宽 wc=1.1*wm;%滤波器截止频率 Ts=0.7*pi;%抽样间隔 ws=2*pi/Ts;%抽样角频率 n=-100:100;%时域抽样点数 nTs=n*Ts;%时域抽样点 f=sinc(nTs/pi); Dt=0.005;
信号与系统实验4：抽样定理与信号恢复
