🚀 go-pugleaf

テトリスのソースからデジタル信号処理まで変な分野しかしらない。
ようめい＠明治大学です。

> for( j=1 ; j<n ; j++ )
> {
> for(sigma=i=0 ; i<n ; i++) sigma += x[i] * cosTable[(( 2 * i + 1 ) *
>j)%nn] ;
>
> y[j] = nrm * sigma;
> }


式からすると離散コサイン変換（ＤＣＴ）ですよね？
これは離散フーリエ変換（ＤＣＴ）と同じで高速フーリエ変換（ＦＦＴ）を使えば高
速化できたと思います。確かＤＣＴの厚い本によく載っています。

それと、cosTable[]を求めるときも全てのcosを計算するのではなく、z変換を用いて
差分式に変換するのが良いでしょう。最初だけcos()を求めてあとは差分を計算すれ
ば高速でcosTableを求められます。（差分の場合、誤差がでますけれど。）
この辺も、僕が全く理解不能な「デジタル信号処理うんぬん」とか「Ｚ変換うんぬ
ん」とか書いてある本によく載っています。でもこれらの本は理論ばかり書いてあっ
て実際のことはあまり書いてないのでよくわからないです．