Leverrierのアルゴリズム
[num,den]=nlev(A,z [,rmax])
実数正方行列
文字列
オプションのパラメータ (bdiag
参照)
[num,den]=nlev(A,z [,rmax])
は,
(z*eye()-A)^(-1)
を計算します.
計算の際には,Aのブロック対角化の後に Leverrierアルゴリズムを各ブロックに適用します.
このアルゴリズムは通常の Leverrierアルゴリズム より優れていますが,まだ完全ではありません!