一般に、すべてのクリロフ法は、行列のスペクトルで評価したときに小さい多項式を本質的に探します。特に、クリロフ法の番目の残差(初期推定がゼロ)は、次の形式で記述できます。n
rn= Pn(A )b
ここで、は次数モニック多項式です。 nPnn
場合有する、対角化可能である、我々はA = V Λ V - 1AA = VΛ V− 1
∥ Rn∥≤=∥ V∥は⋅ ∥はPn(Λ )∥ ⋅ ∥ V− 1∥ ⋅ ∥ B ∥κ (V)⋅ ∥ Pn(Λ )∥ ⋅ ∥ B ∥ 。
た場合に、(例えば、対称またはユニタリ)正常で我々は知っている CGが異なる内積を使用して多項式を構築しながら(参照、アーノルディ反復介してGMRES構築ような多項式この回答を詳細)。同様に、BiCGは非対称Lanczosプロセスを介して多項式を構築しますが、チェビシェフ反復はスペクトルに関する事前情報(通常、対称定行列の最大および最小の固有値の推定)を使用します。κ (V )= 1。Aκ (V)= 1。
クールな例(Trefethen + Bauが動機)として、スペクトルが次のようなマトリックスを考えます。
MATLABでは、次を使用してこれを構築しました。
A = rand(200,200);
[Q R] = qr(A);
A = (1/2)*Q + eye(200,200);
次数すべてのモニック多項式にわたって実際に残差を最小化する多項式を構築するGMRESを考慮すると、候補多項式を調べることで簡単に残差履歴を予測できますn
Pn( z)= (1 − z)n
私たちの場合、それは与える
|Pn( z)| = 12n
以下のためののスペクトルにおける。AzA
ランダムなRHSでGMRESを実行し、残差履歴をこの多項式と比較すると、それらは非常によく似ているはずです(ため、多項式の候補値はGMRES残差よりも小さくなります)。∥ B ∥2> 1