上式表明:X(k)可以用两个N/2点的DFT G(k)和H(K)重建出来。其中含有N次乘法 真的?那,怎么做呢? 具体的做法是:将范围
于是,关于X(k)的N个等式就可以写成下面的两组N/2个等式: 啊! 看起来好象更复杂了。不过,对上面这个计算式子,我们可以进一步简化。 我们注意到G(k)与H(k)实际上是周期的,它们的周期为N/2,因此 于是 显然,上面的这种奇偶分解过程可以不断进行下去,分解到最后,运算就成了两个原序列值进行加减运算了---因为这个时候N=2,而旋转因子的指数为零所以等于1。通过这种分解,DFT的总的运算量就减少了很多。对此,我们可以进行一下简单的分析。 后面我们将从程序设计的角度出发,来考虑如何编程实现上面的快速算法,即讨论快速傅里叶变换的算法实现流程。 |