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