用dft(fft)对时域离散信号进行频谱分析报告.pdfVIP

电子科技大学中山学院电子工程系 学生实验报告 课程名称 数字信号处理实验 实验名称 用DFT(FFT)对时域离散信号进行频谱分析 班级，分组 实验时间 2011 年 11 月29 日 姓名，学号 指导教师 报 告 内 容 一、实验目的和任务 1.进一步加深DFT 算法原理和基本性质的理解（因为FFT 只是DFT 的一种快速算法，所以FFT 的运算结果必然满足DFT 的基本性质）。 2.掌握DFT （FFT ）对时域离散信号进行频谱分析的方法。 二、实验原理简介 1.DFT 和FFT 原理： 长度为N 的序列x(n) 的离散傅里叶变换为X （k ）： N1  nk X(k)  x(n)W ,k0,1,,N1 N n0 首先按 的奇偶把时间序列x(n)分解为两个长为N/2 点的序列 x (n)=x(2r) r=0.1, ….,N/2-1 1 x (n)=x(2r+1) r=0,1, …..,N/2-1 2 则x(n) 的DFT 为X(k) N/ 21 N/ 21  kr k  kr X(k)  x (r)W W x (r)W 1 N/ 2 N 2 N/ 2 r 0 r0 k X (k）W X (K), k0,1,, N/ 2 1 1 N 2 通过上面的推导可以看出，N 点的DFT 可以分解为两个N/2 点的DFT，每个N/2 点的DFT 可以分解为两个点的DTF 。依次类推，当N 为2 的整数次幂时（N=2M ）, 由于每分解一次降低一 阶幂次，所以通过M 次的分解，最后全都成为一系列2 点DFT 运算。以上就是按时间抽取的快 速傅里叶变换(FFT)算法。 2，MATLAB 中计算DFT(FFT)的函数 函数fft 用来求序列的DFT，调用格式为：[X ]=fft(x,N) 其中，x 为有限长序列，N 为序列x K 的长度，X 为序列x 的DFT 。 k 函数ifft 用来求IDFT，调用格式为：[x]=ifft(X,N) 其中，X 为有限长序列，N 为序列X 的 k k 长度，x 为序列X 的IDFT 。 k 三、实验内容和数据记录 (1)复习DFT 的定义，性质和用DFT作频谱分析的有关内容。 1 / 7 （2 ）用MATLAB编制程序产生以下典型信号供谱分析用： x n)=R (n) 1( 4 n1,0 n3  x (n)  8 n,