xie数字信号实验二.doc

电子科技大学中山学院电子信息学院 学生实验报告 课程名称 数字信号处理实验 实验名称 用DFT（FFT）对时域离散信号进行频谱分析2016年　10 月17　 日 姓名，学号 指导教师 报 告 内 容 一、实验目的和任务 1． 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 2．掌握DFT（FFT）对时域离散信号进行频谱分析的方法。1、DFT和FFT原理： 长度为N的序列x(n)的离散傅立叶变换为： 首先按n的奇偶把时间序列x(n)分解为两个长为N/2点的序列 r=0,1,...,N/2-1 r=0,1,...,N/2-1 则x(n)的DFT为 由于，故有 其中分别为的N/2点DFT。因为均是以N/2为周期的，且。因此可将N点DFT分解为下面的形式 k=0,1,...,N/2-1 k=0,1,...,N/2-1 通过上面的推导可以看出，N点的DFT可以分解为两个N/2点的DFT，每个N/2点的DFT又可以分解为两个N/4点的DFT。依此类推，当N为2的整数次幂时()，由于每分解一次降低一阶幂次，所以通过M次的分解，最后全部成为一系列2点DFT运算。以上就是按时间抽取的快速傅立叶变换(FFT)算法。 序列的离散傅立叶反变换为： 离散傅立叶反变换与正变换的区别在于变为，并多了一个的运算。因为和对于推导按时间抽取的快速傅立叶变换算法并无实质性区别，因此可将FFT和快速傅立叶反变换（IFFT）算法合并在同一个程序中。 2、MATLAB中计算DFT(FFT)的函数 (1) 用MATLAB编制程序产生以下典型信号供谱分析用： n1=0:3; n2=4:7; n=0:7; m=0:8; x1n=ones(1,4); x2n=[n1+1,8-n2]; x3n=[4-n1,n2-3]; x4n=[cos(pi*n/4)]; x5m=10*0.8.^m; subplot(2,3,1);stem(n1,x1n,.); xlabel(n); ylabel(x1n); grid on; subplot(2,3,2);stem(n,x2n,*); xlabel(n); ylabel(x2n); grid on; subplot(2,3,3);stem(n,x3n,.); xlabel(n); ylabel(x3n); grid on; subplot(2,3,4);stem(n,x4n,*); xlabel(n); ylabel(x4n); grid on; subplot(2,3,5);stem(m,x5m,.); xlabel(m); ylabel(x5m); grid on; （）分别以变换区间N＝8，16，32对进行DFT(FFT)，画出相应的幅频特性曲线； x1n=ones(1,4); X1k8=fft(x1n,8); X1k16=fft(x1n,16); X1k32=fft(x1n,32); subplot(2,2,1);stem(abs(X1k8),.); xlabel(\omega/\pi); ylabel(X1k8); grid on; subplot(2,2,2);stem(abs(X1k16),*); xlabel(\omega/\pi); ylabel(X1k16); grid on; subplot(2,2,3);stem(abs(X1k32),.); xlabel(\omega/\pi); ylabel(X1k32); grid on; （）分别以变换区间N＝8，16对分别进行DFT(FFT)，画出相应的幅频特性曲线；n1=0:3; n2=4:7; x2n=[n1+1,8-n2]; x3n=[4-n1,n2-3]; X2k8=[fft(x2n,8)]; X2k16=[fft(x2n,16)]; X3k8=[fft(x3n,8)]; X3k16=[fft(x3n,16)]; subplot(2,2,1);stem(abs(X2k8),*); xlabel(\omega/\pi); ylabel(X1k8); grid on; subplot(2,2,2);stem(abs(X2k16),.); xlabel(\omega/\pi); ylabel(X1k8); grid on; subplot(2,2,3);stem(abs(X3k8),*); xlabel(\omega/\pi); ylabel(X1k8); grid on; subplot(