离散信号内插和抽取Matlab,信号的采样与重构__上机实验

上机实验指导

实验一 信号的采样与重构

连续时间信号采样是获得离散时间信号的一种重要方式，但是时域上的离散化会带来信号在频域上发生相应的变化。在本实验中，我们将分别看到低通信号和带通信号在不同的采样率下得到的离散信号波形与连续信号波形在时域和频域上的对应关系。同时，离散信号的二次采样在实际的应用中可能是必须的，有时甚至是非常重要的。在实验的最后，我们也会看到离散信号的抽取和内插所带来的频谱变化。

由于matlab语言无法表达连续信号，实验中我们采用足够密的采样点来模拟连续信号(远大于奈奎斯特采样的要求)，即：

t=0:Ts:T(Ts=1/fs<

实验中，为了分析离散信号与连续信号之间的频谱关系，加深对采样定理的理解，了解模拟频谱、数字频谱、以及离散信号被加窗后各自的频谱，从而直观的理解采样频率对频谱的影响和加窗后对频谱的影响。由此可以掌握数字处理方法对模拟信号进行频谱分析的基本原则，即：如何选择合适的信号长度、采样周期以使得对模拟信号的频谱分析的误差达到分析的要求。

在该实验中，用到的Matlab函数有:

plot(x,y)，其作用是在坐标中以x为横坐标、y为纵坐标的曲线，注意x和y都是长度相同的离散向量；

xlabel(‘xxx’)，其作用是对x轴加上坐标轴说明“xxx”；

ylabel(‘yyy’)，其作用是对y轴加上坐标轴说明“yyy”；

title(‘ttt’)，其作用是对坐标系加上坐标轴说明“ttt”；

subplot(m,n,w)，其作用是当需要在同一显示面板中显示多个不同的坐标系时，m、n分别指明每行和每列的坐标系个数，w为当前显示坐标系的流水号(1到m*n之间)。

在实验中我们需要画出信号的频谱，对于连续信号频谱的逼近需要你自己编写，原理如下：

连续时间非周期信号x(t)的傅里叶变换对为：

X(j ) x(t)e j tdt 用DFT 方法对该变换逼近的方法如下：

1、将x(t)在t轴上等间隔(宽度为T)分段，每一段用一个矩形脉冲代替，脉冲的幅度为其起始点的抽样值x(t)t nT x(nT) x(n)，然后把所有矩形脉冲的面积相加。该方

法实际为平顶处理，利用采样和零阶保持器就可以完成,则有：

X(j )

n x(nT)e j nT T