glever

引数

E, A

同じ次元の正方実行列

s

文字列 (デフォルト値 's')

Bfs,Bis

多項式行列

chis

多項式

説明

一般化したLeverrierのアルゴリズムによりある行列ペンシルについて

(s*E-A)^-1

を計算します.

(s*E-A)^-1 = (Bfs/chis) - Bis.

chis = 特性多項式 (乗法定数まで).

Bfs = 分子の特性多項式行列.

Bis = 多項式行列 ( - (s*E-A)^-1 の無限大までの級数展開).

Bisの前に - 符号があることに注意してください.

注意

この関数は,Bfs,Bis および chisを 簡単化するために cleanpを使用します.

例

s=%s;F=[-1,s,0,0;0,-1,0,0;0,0,s-2,0;0,0,0,s-1];
[Bfs,Bis,chis]=glever(F)
inv(F)-((Bfs/chis) - Bis)

参照

• rowshuff — シャッフルアルゴリズム
• det — determinant of a square matrix
• invr — (有理) 行列の逆
• coffg — 多項式行列の逆
• pencan — 行列ペンシルの正準形
• penlaur — 行列ペンシルのローラン係数