matlab在DSP中的应用（四）---时域抽样与信号的重建

一、实验目的

(1)掌握用MATLAB语言进行离散时间傅里叶变换和逆变换的方法。

(2)了解用MATLAB语言进行时域抽样与信号重建的方法。

(3)进一步加深对时域信号抽样与恢复的基本原理的理解。

(4)观察信号抽样与恢复的图形，掌握采样频率的确定方法和内插公式的编程方法。

二、实验原理

1.DTFT

离散时间傅里叶变换(DTFT)是指信号在时域上为离散的，而在频域上则是连续的。

如果离散时间非周期信号为x(n)，则它的离散傅里叶变换对(DTFT)表示为

反变换为：

其中，信号的频谱为：，其中称为序列的幅度谱，称为序列的相位谱。

从离散时间傅里叶变换的定义可以看出，信号在时域上是离散的、非周期的，而在频域上则是连续的、周期性的。

实例1：求x(n)＝［0，1，2，3，4，5，6，7］，0≤n≤7的DTFT，将(－2* pi，2* pi)区间分成500份。要求：
(1)画出原信号。
(2)画出由离散时间傅里叶变换求得的幅度谱和相位谱图形。

代码：

xn=0:7;
n=0:length(xn)-1;
w=linspace(-2*pi,2*pi,500);%将[-2*pi,2*pi]频率区间分割为500份
X=xn*exp(-1j*n'*w);%离散时间傅里叶变换，这里的1j也为i或j

subplot(3,1,1),stem(n,xn);%显示原序列
ylabel('x(n)')

subplot(3,1,2),plot(w,abs(X));%显示序列幅度谱
axis([-2*pi 2*pi 1.1*min(abs(X)) 1.1*max(abs(X))])
ylabel('幅度谱')

subplot(3,1,3),plot(w,angle(X));%显示序列相位谱
axis([-2*pi 2*pi 1.1*min(angle(X)) 1.1*max(angle(X))])
ylabel('相位谱')

输出：

在这个例子中，最重要的就是离散时间傅里叶变换公式的理解，结合本题，可将公式写成如下矩阵形式：

那么，自然地，就可得到DTFT的matlab代码：X=xn*exp(-1j*n'*w);

注：代码中之所以将j写成1j，是因为matlab自动提示这么干，原理是可提到程序运行速度和鲁棒性。

2.信号采样

离散时间信号大多由连续时间信号(模拟信号)抽样获得。在模拟信号进行数字化处理的过程中，主要经过A/D转换、数字信号处理、D/A转换和低通滤波等过程，如下图所示。其中，A/D转换器的作用是将模拟信号进行抽样、量化、编码，变成数字信号。经过处理后的数字信号则由D/A转换器重新恢复成模拟信号。

如果A/D转换电路输出的信号频谱已经发生了混叠现象，则信号再经过后面的数字信号处理电路和D/A转换电路就没有实际使用的意义了。因此，信号进行A/D转换时，采样频率的确定是非常重要的。

表示了一个连续时间信号xa(t)、对应的抽样后获得的信号 以及对应的频谱。在信号进行处理的过程
中，要使有限带宽信号xa(t)被抽样后能够不失真地还原出原模拟信号，抽样信号p(t)的周期Ts及抽样频率Fs的取值必须符合奈奎斯特(Nyquist)定理。假定xa(t)的最高频率为fm，应有Fs≥2fm，即图中的。

因此，若Fs>=2fm时，只要将抽样信号通过一个低通滤波器，就可以不失真地还原出原信号。
（注：原信号频谱在f=fm这个点的值不为0时，要求：Fs>2fm，不能等于）

2.1对连续信号进行采样

在实际使用中，绝大多数信号都不是严格意义上的带限信号。为了研究问题的方便，我们选择两个正弦频率叠加的信号作为研究对象。

实例1：已知一个连续时间信号，f0＝1 Hz，取最高有限带宽频率fm＝5f0。分别显示原连续时间信号波形和Fs>2fm、Fs＝2fm、Fs<2fm三种情况下抽样信号的波形。
（为方便，可分别取Fs＝fm、Fs＝2fm和Fs＝3fm来研究问题。）

代码：

f0=1;%信号频率
fm=5*f0;%信号在频谱的最高频率
T0=1/f0;
t=-2:0.1:2;
f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t);

subplot(4,1,1),plot(t,f);
ylabel('f(t)')

for i=1:3
    fs=i*fm;%采样频率
    Ts=1/fs;%采样周期，采样步长
    n=-2:Ts:2;%采样
    f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n);
    subplot(4,1,i+1),stem(n,f);
    label=['fs=',num2str(i),'*fm'];