实验四FFT的计算机实现及谱分析.docxVIP

实验四FFT的计算机实现及谱分析 一、实验目的与要求 了解DFT的性质及应用。 熟悉MATLAB编程的特点。 常握有限离散傅里叶变换的计算机实现方法。 (4 )进一步加深DFT算法原理和基木性质的理解(因为FFT只是DFT的一种 快速算法，所以FPT的运算结果必然满足DFT的基本性质)o (5 )熟悉FFT算法原理和FFT子程序的应用。 (6)学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出 现的分析误差及其原因，以便在实际屮正确应用FFT。 二、实验原理与方法 1 .F F T的计算原理 设 N = 2k , 有限离散信号心…,兀心)，则它的有限离散频谱为 Xn=yxe N ^mN-\ f/l n 7 =0 令 WN =e n 则上式为 N-\ Xn=yxWfUU^mN-\ /t=0 把信号暫=牝叭,…,％J按下标偶数项和奇数项分成两部分，即令 = %, hl =兀2/+1, 它们的冇限离散频谱分别为 TOC \o 1-5 \h \z N/2-1 N/2-1 q = X gz泸=X切吒加/=() /=() N/2-1 N/2-1 H,n = E = ￡ 仏府” /=0 /=0 经分析，我们可得如下的时域分解FFT算法 Xl=Gl^Hl XN/2+l=Gl-W^Hl 若我们把X”,按偶奇序号分成两部分，则可得如下的频域分解FFT算法 7V/2-1 X2,工(￡+兀朴/2)吒2 X2, n=0 /V/2-1 E [(捡 ~ Xn+N/2 )必]叱2 ”二0 2、MATLAB提供的快速傅立叶变换函数 fft fft函数用于计算快速傅立叶变换，其语法格式为： B 二 fft(T) B - fft⑴返回信号I的fft变换矩阵，输入信号I和输出信号B大小相同。 fftshift MATLAB提供的fftshift函数用于将变换后的信号频谱中心从矩阵的原点移到矩 阵的屮心，其语法格式为： B = fftshift(1) 对于矩阵I, B二fftshift (I)将I的一、三彖限和二、四彖限进行互换。 ifft if ft函数用于计算信号的傅立叶反变换，其语法格式为： B = ifft(I) B二ifft(A)返回信号T的傅立叶反变换矩阵，输入信号T和输出信号B大小相 同。其语法格式含义与fft函数的语法格式相同，可以参考ffl函数的说明。 吋域分解FFT算法的几点说明： 1 对一般的k ,如何求信号￡ =(^0,x1,...,xv_1)的一项频谱排列 (兀0, , ? ? ?,兀N-1) ? 我们可通过二进制的逆序得到：假设J的二进制为U—2…h贝〔J X ? X. X. ? J JZ—2 …Jo JoJi…心 2、碟形运算两节点的距离： 第一级(第一列)两节点间的距离为仁 第二级(第二列)两节点间的距离为2； 第m级(第m列)两节点间的距离为2-1 ； 3、W；；的确定 因子最后一列有N/2种，顺序为为W；,吒，…,耐3 ； 因子倒数第二列有N/4种,为叫,嘴,…，昭心； 略因子倒数第三列有N/8种，为W：,叱严-4. 因子第一列有一种，为W； = 1 频域分解FFT算法的几点说明： 1＞对一般的最后得到的逆序频谱文j如何排列才能得到原信号的频谱? 我们可通过二进制的逆序得到：假设丿?的二进制为Ud??h，贝I」 Y — V — V ?一 …小一 A(册…爪］) 2、碟形运算两节点的距离： 第一级(第一列，m = 1)两节点间的距离为2k~l ； 第m级(第m列)两节点间的距离为 3、的确定 因了第一列有N /2种，顺序为为，吒,???,碎3 . 必因子第二列冇N/4种，为叫陀…，昭心； 叫因子第三列有N/8利「，为吧,呛,…,比严-4； 因子最后一列有一种，为= 1 三实验内容及步骤 1、用三种不同的DFT程序计算班力=R,(n)的傅里叶变换X”,，并比较三种程序计算 机运行时间。 用for语句的M函数文件dftl. m,用循环变量逐点计算X(k)； 编写时域分解FFT算法或频域分解FFT算法的M函数文件dft2. m； 调用FFT库函数,直接计算X(k)； 分别利用上述三种不同方式编写的DFT程序计算序列x(n)的傅立叶变换 心，并画出相应的幅频和相频特性，再比较各个程序的计算机运行时间。 (5 )复习DFT的定义、性质和用DFT作谱分析的有关内容。 (6 )复习FFT算法原理与编程思想，并对照DIT-FFT运算流图和程序框图, 读懂木实验提供的FFT子程序。 (7 )编制信号产生子程序，产生以下典型信号供谱分析用： 兀[(n) = R4 (n) n +1 0/i3 x2(n) = 8-72 4 n 7 0 -其它乃 4 -n 0 72 3 兀3(〃)= n-3 4 n 7 0 其它” x4(/t