QZ変換により列階段型ペンシルを計算する
[AE,EE,QE,ZE,blcks,muk,nuk,muk0,nuk0,mnei]=fstair(A,E,Q,Z,stair,rk,tol)
実数エントリを有するm x n行列.
実数の正のスカラー.
列階段型行列
m x m ユニタリ行列
n x n ユニタリ行列
添字ベクトル (ereduc参照)
整数, 行列ランクの推定値
実数エントリを有するm x n行列.
列階段型行列
m x m ユニタリ行列
n x n ユニタリ行列
行列A
で検出された
フル行ランクを有するサブ行列の数(>= 0).
次元 (n) の整数配列. ペンシルsE(eps)-A(eps)において 列フルランクを有するサブ行列の列次元 mu(k) (k=1,...,nblcks) を含みます.
次元 (m+1) の整数配列. ペンシルsE(eps)-A(eps)において 行フルランクを有するサブ行列の行次元 nu(k) (k=1,...,nblcks) を含みます.
次元 (n) の整数配列. ペンシルsE(eps,inf)-A(eps,inf)において 列フルランクを有するサブ行列の列次元 mu(k) (k=1,...,nblcks) を含みます.
次元 (m+1) の整数配列. ペンシルsE(eps,inf)-A(eps,inf)において 行フルランクを有するサブ行列の行次元 nu(k) (k=1,...,nblcks) を含みます.
次元 (4) の整数配列. mnei(1) = sE(eps)-A(eps)の行の次元
行列 E
を列階段形として,
ペンシル sE-A
を指定すると,
関数fstair
は
ユーザの指定に基づき
ユニタリ変換されたペンシルQE(sEE-AE)ZE
を
計算します.
このペンシルは, ほぼペンシルsE-A
の一般化Schur形式です.
この関数は,指定したペンシルの
クロネッカー構造の部分も出力します.
Q,Z
はユニタリ行列で,
ペンシルを計算する際に使用されます.
ただし, E は列階段形です (ereduc参照)