一个复数傅立叶级数之机械计算算法.docVIP

一个复数傅立叶级数之机械计算算法 Yates介绍了一个2m的析因实验的相互作用的有效的计算方法，让他的名字家喻户晓。关于3m的一般化是由Box [1]和Good [2]等人提出的。他们总结这些方法，并给一类关于傅立叶级数计算的应用的快速算法。在他们完整的概括中，Good的方法被应用到某些问题，在其中一个N维向量必须乘以一个可以分解成m个稀疏矩阵的N*N矩阵，其中m是成正比logN。这样导致一个需要一些正比与N*logN而非N运算的过程。这些方法正应用于计算复数傅立叶级数。当数据点的数量是一个高度综合的数字，这些情况下这种方法是非常有用的。该算法是在一个相当不同的的形式提出的，同时需注意N的选择。它也显示N = 2m时使用这种方法的特殊优势以及如何傅立叶系数在N个数据存储地点的数组的使用的整个计算。 考虑当计算复数傅立叶级数时，给定的傅立叶系数A（K）是复杂的和W是N次根数的问题，一个使用（1）的简单计算将需要N2次操作，“操作”是指一个复杂的乘法和一个复杂的加法将贯穿始终。 这个算法描述了对于给定的复傅里叶幅度的数组的迭代，只需要小于2N?log2?N次操作，而且比数组A需要更少的存储空间。为了推导这个算法，假设N是一个合数，即N=r1×r2。接着将这个数代入式（1）， 这样式子（1）就可以写为 由于 将用和表示，则里面的可表示成一个新的数组 这样式（1）就可