测试信号分析与处理第4章 离散傅里叶变换及其快速算法.pptVIP

第四节 离散傅里叶变换的性质 帕斯瓦尔定理 若 ，则 式左端代表离散信号在时域中的能量，右端代表 在频域中的能量，该式表明变换过程中能量是守恒 的。 第五节 快速傅里叶变换 一、DFT运算的特点 1.DFT直接运算的工作量 N点序列x(n)的DFT为 按定义计算DFT时，需作 次复数乘和 次复数加。 由于直接计算量非常大，不可能实现信号的实时处理，因此必须改进DFT的算法。 第五节 快速傅里叶变换 2.DFT运算特点 1) 的周期性 2) 的对称性 因为 故 第五节 快速傅里叶变换 3) 的可约性 第五节 快速傅里叶变换 利用上述结果，如对应于N=4的矩阵W可以简化为 上式右端的矩阵中的元素有许多是相等的，计算量明显减少。由原来的16个元素变成只有两个独立元素需要计算。 第五节 快速傅里叶变换 二、基2时析型FFT算法（时间抽取法） 1.算法原理 对长度为 （L为正整数，若原序列的长度不满足此条件，则可用零补足）的序列x(n),按序列各项序号的奇偶将序列分成两个子序列（大点数化为小点数），有 偶序号序列 奇序号序列 其中 ，两序列的长度均为N/2。且其DFT 分别为 第五节 快速傅里叶变换 上式中，右边即为x(n)的前一半DFT值 第五节 快速傅里叶变换 对于另外N/2个点的DFT,即 的 ，可表示为 根据周期性，有 由对称性，有 第五节 快速傅里叶变换 x(n)的DFT最后结果 x(n)的DFT可由奇偶对分序列的DFT合成而获得。该计算关系可用右图来表示。由于图形酷似蝴蝶，所以称之为蝶形图（或蝶形单元）。上式也因而称为蝶形算法。 第五节 快速傅里叶变换 2.算法的具体实现 对于长度为 的序列逐次奇偶对分，则最后一定能得到N个单项序列（序列的长度为1），而单项序列的DFT就是其本身。因此根据蝶形图，计算N项序列的DFT，只需要按照蝶形算法逐次合成，即由N个1点长序列的DFT合成N/2个2点长序列的DFT，再由N/2个2点长序列的DFT合成N/4个4点长序列的DFT，如此继续下去，最后由两个N/2点长序列的DFT合成N点长序列的DFT。 第五节 快速傅里叶变换 第五节 快速傅里叶变换 当序列的长度为 时，根据基2时析型FFT的 算法，可画出如下蝶形图，求出序列的DFT值。 第五节 快速傅里叶变换 3.流程图规律 1)蝶形图参数 蝶群序号 蝶距（序号差） 蝶群宽（点数） 蝶群数 第一级（2点DFT） 第i级（ 点DFT) 第L级（ 点DFT) 第五节 快速傅里叶变换 2）每个蝶形单元的运算，都包括乘 ，并与相应的DFT结果加减各一次 3）同一级中， 的分布规律相同 第一级（2点DFT）： 第二级（4点DFT）： ； ； 第 i级( 点DFT): ; ;...; 第五节 快速傅里叶变换 序列输入的 自然顺序 十进制 二进制码 码位倒置结果（二进制码） 乱序十进制 序列乱序的输入顺序 4）输入重排 第五节 快速傅里叶变换 4.运算量比较 N（ ）点的FFT总运算量为 复数乘 复数加 利用基2时析型FFT求序列的DFT同直接计算序列的DFT的复数乘运算次数之比为 第六节 IDFT的快速算法（IFFT） 利用FFT的程序求IFFT的方法 先对DFT和IDFT两者的定义式进行比较，二者没有本质上的区别，可以利用FFT流图（蝶形图）来计算IFFT。 第六节 IDFT的快速算法（IFFT） 对DFT的反变换取共轭，有 与DFT的正变换式比较可知，完全可以利用FFT的计算 程序，只需将 作为输入序列，并将最后结果取共 轭，再除以N就得到了 第七节 实序列的FFT高效算法 一、同时计算两组实序列的DFT 设有两组长度均为N的实序列 和 ，构成新序列 ，根据DFT定义 因为 ， ，有 因此 第七节 实序列的FFT高效算法 解得 由此说明用FFT子程序求得 的DFT 后， 和 的DFT（ 、 ）可由 的DFT( ) 分离出来。 第七节 实序列的FFT高效算法 二、用N点序列的DFT结果获得2N点长实序列的DFT结果 将一个2