数字信号处理 实验六 离散傅立叶变换及其快速算法.docVIP

实验六 离散傅立叶变换及其快速算法 实验目的：掌握快速傅立叶变换的应用方法；掌握离散余弦变换的应用方法；掌握Z变换的应用方法；了解Chip z变换的基本概念；掌握Hilbeit变换的初步应用；了解倒谱变换的基本概念。 实验仪器：电脑一台，MATLAB6.5或更高级版本软件一套。 实验内容： 快速傅立叶变换（FFT） 在信号处理中，DFT的计算具有举足轻重的地位，，信号的相关、滤波、谱估计等都要通过DFT来实现。然而，当很大的时候，求一个点的DFT要完成次复数乘法和次复数加法，其计算量相当大。1965年J.W.Cooley和J.W.Tukey巧妙地利用因子的周期性和对称性，构造了一个DFT快速算法，即快速傅立叶变换(FFT)。概念 通过前面的知识，已经知道有限列长为的序列的变换为 其逆变换为 由于MATLAB软件本身的特点，序列或向量元素下标从1开始记录，而不是从0开始。因此，上述两式在MATLAB中相应的表达式为 而下面所讨论使用的快速傅立叶变换并不是与不同的另外一种变换，而是为减少计算次数的一种快速有效的算法。这种快速算法，主要是利用了下面两个特性使长序列的分解为更小点数的所实现的。 利用的对称性使运算中有些项合并 利用的周期性和对称性使长序列的分解为更小点数的 快速傅立叶变换算法正是基于这一基本思想而发展起来的。快速傅立叶变换算法形式很多，但是基本上可以分为两大类，即按时间抽取（Decimation－In－Time，简称）法和按频率抽取（Decimation－In－Frequency）法。在这里，以时间抽取的算法（库利－图算法）为例，简单说明一下算法的算法原理。 为了讨论方便，设，其中为整数。如果不满足这个条件，可以认为得加上若干零点来达到。由的定义知 其中是列长为的输入序列，把它按的奇偶分成两个序列 又由于， 则 上式表明了一个N点的DFT可以被分解为两个N/2点的DFT。同时，这两个N/2点的DFT按照上式又可以合成为一个N点的DFT。 为了要用点数为N/2点的X(k)、X(k)来表达N点的X(k)值还必须要用W系数的周期性，即 这样可得 即 X（ 同理可得 X（ 另外再加上W的对称性 就可以将X(k)的表达式分为前后两个部分： 前半部分 后半部分 由以上分析可见，只要求出区间内各个整数k值所对应的X(k)、X(k)的值，即可求出区间内的全部X（k）值，这一点恰恰是FFT能大量节省计算的关键所在。 函数应用 MATLAB为计算数据的离散快速傅立叶变换，提供了一系列丰富的数学函数，主要有Fft、Ifft、Fft2 、Ifft2, Fftn、ifftn和Fftshift、Ifftshift等。当所处理的数据的长度为2的幂次时，采用基-2算法进行计算，计算速度会显著增加。所以，要尽可能使所要处理的数据长度为2的幂次或者用添零的方式来添补数据使之成为2的幂次。 Fft和Ifft函数 调用方式 Y=fft(X) 参数说明 ·如果X是向量，则采用傅立叶变换来求解X的离散傅立叶变换； ·如果X是矩阵，则计算该矩阵每一列的离散傅立叶变换； ·如果X是（N维数组，则是对第一个非单元素的维进行离散傅立叶变换； Y＝fft(X,N) 参数说明 N是进行离散傅立叶变换的X的数据长度，可以通过对X进行补零或截取来实现。 （3）Y＝fft(X,[],dim) 或Y＝fft(X,N,dim) 参数说明 ·在参数dim指定的维上进行离散傅立叶变换； ·当X为矩阵时，dim用来指定变换的实施方向：dim=1,表明变换按列进行；dim=2表明变换按行进行。 函数Ifft的参数应用与函数Fft完全相同。 应用说明 【实例6-1】fft的应用 X=[2 1 2 8]; Y=fft(X) 运行结果 Y＝ 13.0000 0+7.0000i -5.0000 0-7.0000 【实例6-2】fft(X,N,dim)的应用 A=[2 5 7 8; 1 4 0 5; 3 8 5 1; 9 1 2 7]; Z=fft(A,[],1) 【实例6-3】fft在信号分析中的应用 使用频率分析方法从受噪声污