实验六信号抽样与重建.docVIP

实验六 信号抽样与重建 1 实验目的 (1) 掌握信号的抽样及抽样定理。 (2) 掌握利用MATLAB分析抽样信号的频谱。 (3) 掌握和理解信号抽样以及信号重建的原理。 (4) 理解频率混叠的概念。 2实验原理及方法 2.1信号的抽样及抽样定理 抽样就是从连续时间信号中抽取一系列的信号样本，从而得到一个离散时间序列，这个离散序列经量化后，就成为所谓的数字信号。今天很多信号在传输与处理时，都是采用数字系统进行的，但是数字系统只能处理数字信号，不能直接处理连续时间信号或模拟信号。为了能够处理模拟信号，必须先将模拟信号进行抽样，使之成为数字信号，然后才能进行传输与处理。所以，抽样是将连续时间信号转换成离散时间信号必要过程。模拟信号经抽样、量化、传输和处理之后，其结果仍然是一个数字信号，为了恢复原始连续时间信号，还需要将数字信号经过所谓的重建和平滑滤波。 图6-1给出了信号理想抽样的原理图。 图6-1 (a) 抽样原理图 (b)带限信号的频谱 上图中，假设连续时间信号x(t)是一个带限信号，其频率范围为，抽样脉冲为理想单位冲激串，其数学表达式为： 6-1 由图可见，模拟信号x(t)经抽样后，得到已抽样信号xs(t)： 6-2 将p(t)的数学表达式代入上式得到： 6-3 显然，已抽样信号xs(t) 也是一个冲激串，只是这个冲激串的冲激强度被x(nTs) 加权了。 从频域上来看，p(t) 的频谱也是冲激序列，且为： 6-4 根据傅里叶变换的频域卷积定理，时域两个信号相乘，对应的积的傅里叶变换等于这两个信号的傅里叶变换之间的卷积。所以，已抽样信号xs(t)的傅里叶变换为： 6-5 表达式6-5告诉我们，如果信号x(t)的傅里叶变换为X(j()，则已抽样信号xs(t) 的傅里叶变换Xs(j()等于无穷多个加权的移位的X(j()之和，或者说，已抽样信号的频谱等于原连续时间信号的频谱以抽样频率(s为周期进行复制的结果。如图6-2所示： 图6-2 信号抽样及其频谱图 由图6-2可见，如果抽样频率不小于信号带宽的2倍时，xs(t) 的频谱中，X(j()的各个复制品之间没有混叠，因此，可以用一个理想低通滤波器来恢复原始信号。由抽样信号恢复原来的原始信号的过程称为信号的重建。反之，如果抽样频率小于信号带宽的2倍时，xs(t) 的频谱中，X(j()的各个复制品之间的距离(也就是(s)太近，所以必将造成频谱之间的混叠，在这种情况下，是无论如何也无法恢复出原来的连续时间信号的。 由此，我们得出下面的结论：当抽样频率 (s2(M 时，将原连续时间信号x(t)抽样而得到的离散时间序列xs(t)可以唯一地代表原连续时间信号，或者说，原连续时间信号x(t)可以完全由xs(t)唯一地恢复。 以上讨论的是理想抽样的情形，由于理想冲激串是无法实现的，因此这种理想抽样过程，只能用来在理论上进行抽样过程的分析。在实际抽样中，抽样往往是用一个A/D转换器实现的。一片A/D转换芯片包含有抽样保持电路和量化器。模拟信号经过A/D转换器后，A/D转换器的输出信号就是一个真正意义上的离散时间信号，而不再是冲激串了。 A/D转换器的示意图如图6-3所示。 图6-3 A/D转换器示意图 上述的实际抽样过程，很容易用简单的数学公式来描述。设连续时间信号用x(t)表示，抽样周期为Ts，抽样频率为(s，则已抽样信号的数学表达式为 6-6 在MATLAB中，对信号抽样的仿真，实际上就是完成式6-6的计算。下面给出一个例题和相应的范例程序，来实现信号抽样的仿真运算。 例6-1：设连续时间信号为一个正弦信号 x(t) = cos(0.5πt)，抽样周期为Ts = 1/4秒，编程序绘制信号x(t)和已抽样信号x[k]的波形图。 % Program6_1 t = 0:0.01:10;Ts = 1/4; k = 0:Ts:10; x = cos(0.5*pi*t); xk = cos(0.5*pi*k); subplot(211),plot(t,x),box off title(原信号x(t)),xlabel(t) subplot(212),stem(k,xk),box off hold on,plot(k,xk,o) title(抽样信号x[k] )),xlabel(k) 执行该程序后，得到的波形图如图6-4所示。 图6-4 连续信号及抽样后的离散序列 在这个范例程序中，先将连续时间t进行离散化，使之成为以Ts = 1/4秒的离散时间k，然后将k代入到信号