数字信号处理(第四版)高西全第4章.pptVIP

　　设x(n)为N点实序列，取x(n)的偶数点和奇数点分别作为新构造序列y(n)的实部和虚部，即   对y(n)进行N/2点FFT，输出Y(k)，则 　　根据DIT-FFT的思想及式(4.2.7)和(4.2.8)，可得到 X(k)的前 个值：  (4.3.5) 式中， 。由于x(n)为实序列， 因此X(k)具有共轭对称性，X(k)的另外N/2点的值为 计算 点FFT的复乘次数为 　　　，计算式(4.3.5)的复乘次数为　　，所以用这种算法, 计算X(k)所需复数乘法次数为　　　　　　　　　　　　。相对一般的N点FFT算法，上述算法的运算效率为 　　　　　　　　　　　　　　，当N=2M=210时，η=20/11，运算速度提高近1倍。 4.4 其他快速算法简介 　　快速傅里叶变换算法是信号处理领域重要的研究课题。由于教材篇幅和教学大纲所限，本章仅介绍算法最简单、编程最容易的基2FFT算法原理及其编程思想，使读者建立快速傅里叶变换的基本概念，了解研究FFT算法的主要途径和编程思路。其他高效快速算法请读者参考文献［1］、［3］、［12］。例如，分裂基FFT算法、离散哈特莱变(换DHT)、基4FFT、基8FFT、基rFFT、混合基FFT，以及进一步减少运算量的途径等内容，对研究新的快速算法都是很有用的。本节简要介绍其他几种快速算法的运算量及其主要特点，以便读者选择快速算法时参考。 　　从理论上讲，不同基数的FFT算法的运算效率不同，实际中最常用的是基2FFT、基4 FFT、分裂基FFT和DHT ［1］。为此，下面简要介绍后三种FFT算法的特点和运算效率，以扩展读者的视野。其具体算法请参考文献［1］、［12］。 　　在基rFFT算法中，基4FFT算法运算效率与基8FFT很接近，但基4FFT算法实现程序简单，且判断开销少。可以证明，当FFT的基大于4时，不会明显降低计算量。基4FFT要求N=4M，M为自然数。其复数乘法次数为［12］ （4.4.1） CM(基4) 其中未计入乘以±j和1的计算。比较基2FFT的复数乘法次数 ，基4FFT的复数乘法次数减少25%。 　　1984年，法国的杜梅尔（P.Dohamel）和霍尔曼（H. Hollmann）将基2分解和基4分解糅合，提出了分裂基FFT算法，其复数乘法次数接近FFT理论最小值，但其运算流图却与基2FFT很相似，编程简单，运算程序也很短，是一种很实用的高效算法。分裂基FFT算法复数乘法次数为［2］  CM(分裂基)= 　　　　　　 （4.4.2） 只考虑（4.4.2）式的第一项，分裂基FFT算法的复数乘法次数就比基2FFT减少33%，比基4FFT减少11%。应当说明，在比较时，未考虑（4.4.2）式后2项减少的运算量，所以分裂基FFT算法的效率更高。 　　由以上比较可见，分裂基FFT算法的效率最高，所以得到广泛应用。但是，对实序列x(n)，上述各种FFT算法仍将其看成虚部为零的复序列存储和计算。而一次复数乘法需要四次实数乘法和二次实数加法。所以，必然浪费存储资源和增加多余的运算量。我们知道，实序列的N点DFT具有共轭对称性，　即  所以，只要计算出X（k）的前面N/2个值，则其后面的N/2个值可以由对称性求得。因此，FFT算法得到的N个 X（k）值有一半是多余的。由以上分析可见，对实序列一定存在更高效的快速算法。离散哈特莱变换（DHT）就是针对实序列的一种高效变换算法，相对一般的FFT算法，DHT的快速算法FHT可以减少近一半的计算量［1］。N点基2时域抽取快速DHT（基2DIT-FHT）算法的实数乘法次数为［1］  （4.4.3） N点基2DIT-FHT算法的实数加法次数为［1］  （4.4.4）  由式（4.4.3）可见，基2DIT-FHT算法的实数乘法次数约为基2DIT-FFT算法的一半。与前面三种FFT算法比较