数字信号处理实验报告6_离散傅里叶变换及快速算法.pdf

数字信号处理实验报告6_离散傅里叶变换及快速算法--第1页 实验六 离散傅立叶变换及其快速 算法 一、实验目的： 掌握快速傅立叶变换的应用方法；掌握离散余弦变换的应用 方法；掌握Z 变换的应用方法；了解Chip z 变换的基本概念； 掌握Hilbeit 变换的初步应用；了解倒谱变换的基本概念。 二、实验仪器： 电脑一台，MATLAB6.5 或更高级版本软件一套。 三、实验内容： 3.1快速傅立叶变换（FFT） DFT 是信号分析与处理中的一种重要变换。但是直接计算 DFT 的运算量与变换的长度N 的平方成正比，当N 较大时， 计算量太大。在快速傅里叶变换（简称FFT ）出现之前，直 接用DFT 算法进行谱分析和信号的实时处理是不实际的。 FFT 使得DFT 的运算效率大大提高，为数字信号处理技术应 用于各种信号的实时处理创造了条件，推动了数字信号处理 技术的发展。 数字信号处理实验报告6_离散傅里叶变换及快速算法--第1页 数字信号处理实验报告6_离散傅里叶变换及快速算法--第2页 FFT算法的MATLAB实现 MATLAB提供fft函数来计算x(n)的DFT,fft函数是有机器语 言，而不是以MATLAB指令格式写成的，因此它的执行速 度很快。 格式：y=fft(x)，计算信号x的快速傅里叶变换y。当y为矩 阵（多通道信号）时，计算x 中每一列信号的离散傅里叶变 换。当x的长度为2的幂时，有基2算法，否则采用较慢的 分裂基算法。 Y=fft(x,n)，计算 点FFT，当x的长度大于 时，截断x否 则补零。 Y=fft(x,n),计算 点FFT，当x的长度大于 时，截断x，否 则补零。IFFT可由ifft 函数来计算。 在信号处理中，DFT的计算具有举足轻重的地位，，信 号的相关、滤波、谱估计等都要通过DFT来实现。然而，当 很大的时候，求一个 点的DFT要完成N N 次复数乘法和 N N N(N 1)次复数加法，其计算量相当大。1965年J.W.Cooley和 W J.W.Tukey巧妙地利用 因子的周期性和对称性，构造了一 N 个DFT快速算法，即快速傅立叶变换(FFT)。概念 x(n) 通过前面的知识，已经知道有限列长为 的序列 的 变 N DFT 换为 N 1  nk X (k)  x(n)W n k 0,1,2 N 1 n0 其逆变换为 数字信号处理实验报告6_离散傅里叶变换及快速算法--第2页 数字信号处理实验报告6_离散傅里叶变换及快速算法--第3页 1N1 x(n)  X(k)W nk N N n0,1,N 1 k0 由于MATLAB软件本身的特点，序