回答:
繰り返し行う場合、とからとを計算することにより、浮動ポイントエラーは蓄積によって爆破されません。これは、遷移行列が正規直交であるためです。単純な回転行列です。のcos (nはθ )の罪((N - 1 )θ )のcos ((N - 1 )θ )
一部のFFT実装でも使用されていることがわかっている、私が好む手順も3項の反復を含みますが、以前のおよびの計算よりも少し安定しています。(Bulirsch / Stoerでのこれらの再発の優れた分析があります。例2と3を参照してください。不安定性を幾何学的に見る方法は、最初はおおよその円をトレースしますが、長い目で見れば、文字通り制御不能になっているということです。)
再発(事実上、Bulirsch / Stoerの例4と同じです。詳細な分析については、次を参照してください)は次のように進行します。
事前に計算された定数があります
前述のように、特にが小さい場合は、非常に優れたパフォーマンスを発揮します。