[论文解读] Algebraic Signal Processing Theory: Cooley-Tukey Type Algorithms for Polynomial Transforms Based on Induction

Aliaksei Sandryhaila, Jelena Kovačević+1 more|arXiv (Cornell University)|Aug 17, 2010

Digital Filter Design and Implementation参考文献 22被引用 1

一句话总结

本文提出了一种新颖的代数方法，通过选取子代数并逐步归纳，对多项式代数的正则模进行分解，从而推导出多项式变换的快速 O(n log n) 算法，包括 DFT 和 DCT-4。该方法运用表示理论推广了 Cooley-Tukey FFT 框架，得到计算开销最小的新通用基底算法与结构化矩阵分解。

ABSTRACT

A <em>polynomial transform</em> is the multiplication of an input vector x∈Cⁿby a matrix P_b,_α∈C^n×n, whose (k,ℓ)th element is defined as pℓ(α_k) for polynomials pℓ(x)∈C[x] from a list b={p0(x),…,pn-1(x)} and sample points α_k∈C from a list α={α0,…,αn-1}. Such transforms find applications in the areas of signal processing, data compression, and function interpolation. An important example includes the discrete Fourier transform. In this paper we introduce a novel technique to derive fast algorithms for polynomial transforms. The technique uses the relationship between polynomial transforms and the representation theory of polynomial algebras. Specifically, we derive algorithms by decomposing the regular modules of these algebras as a stepwise induction. As an application, we derive novel O(n log n) general-radix algorithms for the discrete Fourier transform and the discrete cosine transform of type 4.<br><br>

研究动机与目标

• 开发一种通用的代数框架，用于推导标准 FFT 之外的多项式变换快速算法。
• 通过引入基于模归纳的系统化方法，扩展代数信号处理理论，实现变换因子分解。
• 推导离散傅里叶变换（DFT）和类型 4 离散余弦变换（DCT-4）的新型通用基底快速算法。
• 提供一种统一且理论基础坚实的途径，揭示快速算法的结构根源，超越临时的系数操作。

提出的方法

• 该方法利用多项式代数的表示理论，将多项式变换解释为在采样点处对基多项式求值的矩阵乘法。
• 通过选定的子代数，将多项式代数的正则模分解为一系列诱导模的链，实现逐步因子分解。
• 利用子代数的横截面对分解进行形式化，导出将变换分解为结构化、低成本矩阵的矩阵因子分解。
• 关键组件包括基变换矩阵（B）、对角缩放矩阵（D）和类似置换的矩阵（L），并为 DFT 和 DCT-4 提供了显式构造。
• 该框架通过将 Cooley-Tukey 方法嵌入代数模归纳中，实现了通用化，允许推导非 2 的幂次大小的算法。
• 通过切比雪夫多项式与三角恒等式，为 DFT 和 DCT-4 推导出具体构造，显式因子分解涉及 DCT-3、DST-3、DCT-4 和 DST-4 变换。

实验结果

研究问题

• RQ1 如何利用多项式代数的代数结构推导出多项式变换的快速算法？
• RQ2 能否通过代数归纳法将 Cooley-Tukey FFT 原理推广至非 2 的幂次大小？
• RQ3 模分解与子代数归纳在生成结构化、低复杂度变换因子分解中起什么作用？
• RQ4 如何利用该代数框架将 DFT 和 DCT-4 因子分解为 O(n log n) 次操作？
• RQ5 此类快速算法存在的最小代数条件是什么？

主要发现

• 本文推导出 DFT 和 DCT-4 的新型 O(n log n) 通用基底算法，突破了标准 2 的幂次 FFT 的限制。
• 该方法为 DFT2km 和 DCT-42km 提供了显式矩阵因子分解，涉及 DFTk、DCT-3m、DST-3m、DCT-4k 和 DST-4k 组件，以及结构化的对角矩阵与置换矩阵。
• 对于 DFT2km，因子分解包含基变换矩阵 B2km_m 和 X2km_m，对角缩放矩阵 D2km_m，并具有与 DFTk 的克罗内克积结构。
• 对于 DCT-42km，因子分解使用斜对称 DCT/DST 矩阵 X(C4)k(r) 和 X(S4)k(r)，结合对角缩放与置换矩阵，实现 O(n log n) 复杂度。
• 该框架统一并推广了先前的代数方法，表明快速算法自然源于模归纳与子代数分解。
• 所推导的算法在数学上等价于已知的快速变换，但基于表示理论从第一原理推导，揭示了其结构根源的深层洞察。