[論文レビュー] A Factorization Scheme for Some Discrete Hartley Transform Matrices

H. M. de Oliveira, Renato J. Cintra+1 more|arXiv (Cornell University)|Feb 1, 2015

Digital Filter Design and Implementation参考文献 9被引用数 7

ひとこと要約

本稿では、ブロック長N ∈ {3, 5, 6, 12, 24}において理論的最小乗算複雑度を達成する離散ハートリー変換（DHT）行列の因子分解方式を提示する。DHTカーネルの対称性を活用し、前加算、スケーリング、後加算に分離された階層的演算に変換することで、実数乗算の下界に一致する高速アルゴリズムを導出する。これにより、低消費電力のデジタル信号プロセッサへの効率的な実装が可能となる。

ABSTRACT

Discrete transforms such as the discrete Fourier transform (DFT) and the discrete Hartley transform (DHT) are important tools in numerical analysis. The successful application of transform techniques relies on the existence of efficient fast transforms. In this paper some fast algorithms are derived. The theoretical lower bound on the multiplicative complexity for the DFT/DHT are achieved. The approach is based on the factorization of DHT matrices. Algorithms for short blocklengths such as $N \in \{3, 5, 6, 12, 24 \}$ are presented.

研究の動機と目的

• 離散ハートリー変換（DHT）の高速アルゴリズムを導出し、乗算複雑度の理論的下界を達成すること。
• 従来の手法が最適でない可能性がある短いブロック長において、効率的なDHTアルゴリズムを設計する課題に対処すること。
• DHTカーネルの自己逆性および対称性を活用して、行列因子分解を簡素化し、計算コストを低減すること。
• 小規模なNに適用可能な低複雑度DHTアルゴリズムを構築する一般化可能なフレームワークを提供すること。
• 実数乗算と加算を最小限に抑えることで、低消費電力のデジタル信号プロセッサへの実装を可能にすること。

提案手法

• DHT行列を、前加算（pre-additions）、無理数定数（例：(√3−1)/2）によるスケーリング、後加算（post-additions）を表すスパース行列の積に因子分解する。
• DHTカーネルの対称性を活用：絶対値が等しい要素同士をハダマードに類似した変換により列ペアリング可能。
• 行列を整数部（行列Lを介して）と残りの無理数成分（行列H′′′に分離）に分解し、後者をさらに因子分解する。
• アルゴリズムは階層的構造を採用：前加算行列（A₁, A₂, ...）、スケーリング行列（B₁, B₂, ...）、後加算行列（C₁, C₂, ...）を組み合わせ、再帰的分解を実現。
• 最終的な変換は、これらの行列の合成に加えて補正行列Lを加算し、V = (C B A + L) A₁ v の形で表現される。
• 12および24点への一般化は、3-, 6-, 12点DHTアルゴリズムを再帰的に埋め込むことで実現される。

実験結果

リサーチクエスチョン

• RQ1 短いブロック長において、乗算複雑度の理論的下限に達する高速DHTアルゴリズムを構築可能か？
• RQ2 DHTカーネルが持つ固有の対称性および自己逆性を活用して、算術演算を最小限に抑える方法は何か？
• RQ3 どのような因子分解戦略により、DHT行列を低複雑度の構成要素（加算、乗算、スパース演算）に分解できるか？
• RQ4 複数のブロック長に適用可能な最小複雑度DHTアルゴリズムを構築する一般化可能なフレームワークを確立可能か？
• RQ5 N = 12およびN = 24における最小乗算DHTアルゴリズムの加算複雑度は何か？

主な発見

• 提案された3点DHTアルゴリズムは、わずか1つの実数乗算と10回の加算で実現され、乗算複雑度の理論的下限に一致する。
• 6点DHTアルゴリズムは2つの乗算と20回の加算を要し、N=6における最小乗算複雑度を達成する。
• 12点DHTアルゴリズムは4つの乗算と52回の加算を実現し、μ(N)の理論的下限を満たす。
• 24点DHTアルゴリズムは12つの乗算と138回の加算を達成し、最小乗算複雑度が達成されていることを確認した。
• 因子分解フレームワークは複数のブロック長に一般化可能であり、小規模なDHT変換の再帰的埋め込みを可能にする。
• 最小限の乗算回数と構造化された計算により、低消費電力のデジタル信号処理に適したアルゴリズムである。