第九章-离散傅里叶变换.ppt

WN0WN0WN0W0N-1-1-1-1P(0)P(1)Q(0)Q(1)Y(0)Y(1)Z(0)Z(1)WN0WN2WN0WN2-1-1-1-1G(0)G(1)G(2)G(3)H(0)H(1)H(2)H(3)WWWWN0N1N2N3-1-1-1-1X(0)X(1)X(2)X(3)X(4)X(5)X(6)X(7)8点DFT的FFT的运算流图如下1.蝶形运算蝶形单元蝶距:2m-1（第m级蝶形运算）五.FFT算法的特点§9.4离散傅里叶变换的性质1、线性这里，序列长度及DFT点数均为N若不等，分别为N1，N2，则需补零使两序列长度相等，均为N，且若则有限长序列的圆周移位导致频谱线性相移，而对频谱幅度无影响。(1)圆周移位2、时移特性(2)时移特性从图中两虚线之间的主值序列的移位情况可以看出：当主值序列左移m个样本时，从右边会同时移进m个样本好像是刚向左边移出的那些样本又从右边循环移了进来因此取名“循环移位”。显然，循环移位不同于线性移位3、频移特性时域序列的调制等效于频域的圆周移位eg：4、时域圆周卷积设则(1)定理圆周卷积步骤：1）翻褶2）圆周移位3）相乘4）相加(2)圆周卷积的计算Eg1.N-10nN-10n0m0m0m0m0233211N-1n最后结果：1）.线性卷积的长度为的长度为2）圆周卷积（3）线性卷积与圆周卷积的关系3)关系补‘0’加长x(n)的N点DFT是x(n)的z变换在单位圆上的N点等间隔抽样；x(n)的DTFT在区间[0,2π]上的N点等间隔抽样。§9-5一.DFT的计算工作量两者的差别仅在指数的符号和因子1/N.§9-6FFT的应用DFT与IDFT运算特点复数乘法复数加法一个X(k)NN–1N个X(k)(N点DFT)N2N(N–1)当N很大时,运算量将是惊人的,如N=1024,要完成1048576次(一百多万次)运算.这样,难以做到实时处理.二.改进的途径1.的对称性和周期性对称性:周期性:2、N点DFT2个N/2点DFT利用上述特性,可以将有些项合并,并将DFT分解为短序列,从而降低运算次数,提高运算速度.1965年,库利(cooley)和图基(Tukey)首先提出FFT算法.对于N点DFT,仅需(N/2)log2N次复数乘法运算.例如N=1024=210时，需要(1024/2)log2210=512*10=5120次。5120/1048576=4.88%,速度提高20倍三、基2FFT算法的思路把一个序列分为长度减半的偶序列和奇序列,原序列的DFT就由这两个N/2序列求得.进一步把N/2序列分解成两个N/4序列,一直分解到2点序列.Eg.N=8的FFT1,3,5,70,2,4,60,42,61,53,7x(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)x(0)x(1)x(2)x(3)x(4)x(5)x(6)x(7)第一级第二级第三级四.算法原理(基2FFT)1.先将按n的奇偶分为两组N/2点作DFT,设N=2M,不足时,可补些零。n为偶数时:n为奇数时:因此其中,2.两点结论:(1)G(k),H(k)均为N/2点的DFT。(2)X(k)=G(k)+WH(k)只能确定出X(k)的k=0,1,…(N/2-1)时的值；即前N/2点的结果。kN同理,3.X(k)后N/2个点的确定*X(k)的后N/2个点，也完全由G(k),H(k)确定*N点的DFT可由两个N/2点的DFT来计算。实现上式运算的流图称作蝶形运算k=0,1,…(N/2-1)4.蝶形运算（N/2个蝶形)11-1由X1(k)、X2(k)表示X(k)的运算是一种特殊的运算-碟形运算k=0,1,…(N/2-1)例N=8时的FFT: 1、2个N/2点DFT（1）将分成奇偶序列，并计算X(k)其中