实验14 快速傅里叶变换(FFT) (完美格式版,本人自己完成,所有语句正确,不排除极个别错误,特别适用于山大,勿用冰点等工具下载,否则下载之后的w ord 格式会让很多部分格式错误,谢谢) XXXX 学号姓名处XXXX 一、实验目的 1、加深对双线性变换法设计 IIR 数字滤波器基本方法的了解。 2、掌握用双线性变换法设计数字低通、高通、带通、带阻滤波器的方法。 3、了解 MATLAB 有关双线性变换法的子函数。 二、实验内容 1、双线性变换法的基本知识 2、用双线性变换法设计 IIR 数字低通滤波器 3、用双线性变换法设计 IIR 数字高通滤波器 4、用双线性变换法设计 IIR 数字带通滤波器 三、实验环境 MATLAB7.0 四、实验原理 1、实验涉及的MATLAB子函数 (1)fft 功能:一维快速傅里叶变换(FFT)。 调用格式: )( xffty ;利用FFT算法计算矢量 x的离散傅里叶变换,当 x为矩阵时,y为矩阵 x每一列的FFT。当 x的长度为 2的幂次方时,则fft函数采用基 2的FFT算法,否则采用稍慢的混合基算法。 ),(nxffty ;采用n点FFT。当 x的长度小于n时,fft函数在 x的尾部补零,以构成n点数据;当 x的长度大于n时,fft函数会截断序列 x。当 x为矩阵时,fft函数按类似的方式处理列长度。 (2)ifft 功能:一维快速傅里叶逆变换(IFFT)。 调用格式: )( xiffty ;用于计算矢量 x的IFFT。当 x为矩阵时,计算所得的y为矩阵 x中每一列的IFFT。 ),(nxiffty ;采用n点IFFT。当 length(x)n时,将 x截断,使 length(x)=n。 (3)fftshift 功能:对 fft的输出进行重新排列,将零频分量移到频谱的中心。 调用格式: )( xfftshifty ;对 fft的输出进行重新排列,将零频分量移到频谱的中心。当 x为向量时,)( xfftshift直接将 x中的左右两半交换而产生 y。 当x为矩阵时,)( xfftshift同时将x的左右、上下进行交换而产生 y。 2、用 MATLAB提供的子函数进行快速傅里叶变换 从理论学习可知,DFT 是唯一在时域和频域均为离散序列的变换方法,它适用于有限长序列。尽管这种变换方法是可以用于数值计算的,但如果只是简单的按照定义进行数据处理,当序列长度很大时,则将占用很大的内存空间,运算时间将很长。 快速傅里叶变换是用于 DFT 运算的高效运算方法的统称,FFT 只是其中的一种。FFT主要有时域抽取算法和频域抽取算法,基本思想是将一个长度为N 的序列...