シルベスター方程式

数学、制御理論におけるシルベスター方程式（シルベスターほうていしき、）とは、次の形式の行列方程式であるこの方程式は等価な次の形式で書かれることも多い： AX − XB = C。

ここで行列 A,B,C は与えられており、等式を満たすような行列 X が求めるべき解である。全ての行列の成分は複素数であるとする。方程式が意味を持つために、行列は適当なサイズでなければならない。例えば、全ての行列が同一サイズの正方行列である等である。より一般には、A と B がそれぞれサイズ n , m の正方行列であれば、X と C はいずれも n 行 m 列の行列でなければいけない。

シルベスター方程式が一意的な解 X を持つのは、行列 A と −B が共通する固有値を持たないとき、またそのときに限る。

より一般的には、方程式 AX + XB = C は（次元が有限とは限らない）バナッハ空間上の有界作用素間の等式であるとして考察されてきた。この場合も解についての条件はほとんど同じである：方程式が一意的な解 X を持つのは、A と −B のスペクトルが交わりを持たないとき、またそのときに限るBhatia and Rosenthal, 1997。

解の存在と一意性クロネッカー積および vec作用素の記法を用いると、シルベスター方程式は次のように書き直せる。

ここでは行列、は行列、は行列、はの単位行列である。この形に書くと、方程式は係数行列による線型方程式系と見ることができるしかし、この形に書き直した方程式は数値計算には適していない。計算量が増大する上、悪条件(ill-conditioned)となる可能性があるからである。。

命題複素数成分の行列とが与えられたとき、シルベスター方程式が任意のに対して一意的な解を持つための必要十分条件はとが共通の固有値を持たないことである。

証明で定まる線型写像を考える。

(i) とが共通の固有値を持たないとする。このときそれらの固有多項式との最大公約数は定数である。よって複素係数多項式とを、が成り立つようにとることができる（ベズーの補題）。ケイリー・ハミルトンの定理より、行列多項式としてであるので、。

をの任意の解とする。このときであり、この等式をの右辺に繰り返し適用して。よって写像の核の次元は0であり、階数・退化次数の定理よりは可逆となるから、任意のに対し一意的な解が存在する。

(ii) 逆に、が行列との共通の固有値であるとする。は転置行列の固有値でもあることに注意すると、零ベクトルでないベクトル , で , を満たすものが存在する。行列をとなるよう選ぶ。右辺はの複素共役である。

このときには解が存在しない。なぜなら、複素数体上の双線型形式（半双線型形式ではない）をと定めてを考えると、この最左辺は

となって矛盾するからである。

ロスの除去法則

複素数成分行列（サイズはそれぞれ）が与えられたとき、次の2つの n + m 次正方行列 , が互いに相似なのはどのようなときか問うことができる。この必要十分条件は AX − XB = C を満たす行列 X、言い換えるとシルベスター方程式の解が存在することである。これはロスの除去法則（Roth's removal rule）として知られている。

次のことは簡単に確認できる：もし AX − XB = C であれば、

ロスの除去法則はバナッハ空間上の無限階有界作用素へ一般化することはできないBhatia and Rosenthal, p.3。

数値計算シルベスター方程式の解を求める古典的なアルゴリズムにがある。これはとをQR法によってシューア標準形に変形し、三角行列に対する後退代入によって解を得るものである。このアルゴリズムの算術処理には（ランダウのO-記法）の計算量が必要であり、GNU OctaveにおけるLAPACKの lyap 関数でも計算量は同じであるhttps://octave.sourceforge.io/control/function/lyap.html（この言語での sylvester 関数も参照のことhttps://www.gnu.org/software/octave/doc/interpreter/Functions-of-a-Matrix.html syl 関数はGNU Octave Version 4.0以降非推奨となった。）。いくつかの特殊な画像処理への応用においては、導出されるシルベスター方程式は閉じた形で解が書ける。

シルベスター方程式

Źródło i licencja