xie数字信号实验二重点.docVIP

电子科技大学中山学院电子信息学院 学生实验报告 课程名称 数字信号处理实验 实验名称 用DFT（FFT）对时域离散信号进行频谱分析 班级，分组 通信工程 实验时间 2016年　10 月17　 日 姓名，学号 指导教师 报 告 内 容 一、实验目的和任务 1． 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 2．掌握DFT（FFT）对时域离散信号进行频谱分析的方法。 二、实验原理简介 1、DFT和FFT原理： 长度为N的序列x(n)的离散傅立叶变换为X(k)： 首先按n的奇偶把时间序列x(n)分解为两个长为N/2点的序列 r=0,1,...,N/2-1 r=0,1,...,N/2-1 则x(n)的DFT为X(k) 由于，故有 其中X1(k) 和X2(k)分别为x1(n) 和x2(n)的N/2点DFT。因为X1(k) 和X2(k)均是以N/2为周期的，且。因此可将N点DFT X(k)分解为下面的形式 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)算法。 序列X(k)的离散傅立叶反变换为： 离散傅立叶反变换与正变换的区别在于WN变为WN-1，并多了一个1/N的运算。因为WN和WN-1对于推导按时间抽取的快速傅立叶变换算法并无实质性区别，因此可将FFT和快速傅立叶反变换（IFFT）算法合并在同一个程序中。 2、MATLAB中计算DFT(FFT)的函数 函数fft用来求序列的DFT,调用格式为：[Xk]=fft(x ,,N) 其中，x 为有限长序列，N为序列x的长度，Xk为序列xn的DFT. 函数ifft用来求IDFT,调用格式为：[x]=ifft(Xk,N) 其中，Xk为有限长序列，N为序列Xk的长度，x为序列Xk的IDFT 三、实验内容和数据记录 (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; （2）分别以变换区间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; （3）分别以变换区间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