实验三离散信号的DTFT和DFT.doc

实验三 离散信号的DTFT和DFT 一、实验目的： 加深对离散信号的DTFT和DFT的及其相互关系的理解。 DFT-Discrete Fourier Transform??离散傅立叶变换 DTFT- Discrete Time Fourier Transform??离散时间傅立叶变换 二、实验原理： 1、傅里叶级数 任何周期函数都可以用正弦函数的级数表示， 有周期函数： 那么： 根据欧拉公式： ， 可以得： 将看成，， 或者： 复指数序列定义，角频率。 是离散时间信号，得到的频谱是连续的。 3、离散傅里叶变换DFT 对第2点的频谱进行采样，频域采样，即频域相乘，等于时域卷积，DFT只是为了计算机处理方便，在频域对DTFT进行的采样并截取主值。 N点序列 的DFT定义： 在MATLAB中，对形式为的DTFT可以用函数H=freqz（num，den，w）计算；可以用函数U=fft（u，N）和u=ifft（U，N）计算N点序列的DFT正、反变换。 三、实验内容： 分别计算16点序列 的16点和32点DFT，绘出幅度谱图形，并绘出该序列的DTFT图形。 讨论DTFT和DFT之间的相互关系。说明实验产生的现象的原因。 四、实验结果： 内容提要 实验结果 （1）零极点 1程序； n=[0:15]; x=cos(5*n*pi/16); X1=fft(x,16); X2=fft(x,32); k=0:15;wk=2*k/16; subplot(3,1,1); stem(wk,abs(X1)); title(16点DFT的幅频特性); xlabel(w/pi);ylabel(幅度); k=0:31;wk=2*k/32; subplot(3,1,2); stem(wk,abs(X2)); title(32点DFT的幅频特性); xlabel(w/pi);ylabel(幅度); w=2*pi*(0:199)/200; Xw=((1-exp(j*5*pi-j*16*w))./(1-exp(j*5*pi/16-j*w)))/2+((1-exp(-j*5*pi-j*16*w))./(1-exp(-j*5*pi/16-j*w)))/2; subplot(3,1,3); plot(w/pi,abs(Xw)); title(DTFT的幅度曲线); xlabel(w/pi);ylabel(幅度); 2绘图结果； DFT和DTFT之间的相互关系：DFT是在频域上对时域DTFT进行采样的主值区间。 原因 ：DTFT是以2pi为周期的。而DFT的序列X（k）是有限长的。 参考程序：? close all;clear all;clc; k=0:15; x=cos(5*pi*k/16); k3=linspace(0,2*pi,1024); h3=freqz(x,[1],k3); 2015-2016下期数字信号处理实验指导书 3