nlev

引数

A

実数正方行列

z

文字列

rmax

オプションのパラメータ (bdiag参照)

説明

[num,den]=nlev(A,z [,rmax]) は, (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))

参照

• coff — レゾルベント (余因子法)
• coffg — 多項式行列の逆
• glever — 行列ペンシルの逆
• ss2tf — 状態空間から伝達関数への変換