基本介紹
- 中文名:西爾維斯特方程
- 外文名:Sylvester equation
- 性質:矩陣方程
- 學科:控制理論
定義,解的存在及唯一,Roth消去法則,數值解,
定義
![](/img/1/1ed/e6774de86555a284655d8f9567d9.jpg)
其中A、B及C是已知的矩陣,問題是要找出符合條件的X。其中所有矩陣的係數都是複數。為了要使方程成立,矩陣的行和列需要滿足一定條件,A和B都要是方陣,大小分別是n和m,而X和C要是n行m列的矩陣,n和m也可以相等,四個矩陣都是大小相同的方陣。
西爾維斯特方程有唯一解X的充份必要條件是A和-B沒有共同的特徵值。
解的存在及唯一
![](/img/2/442/25e37fa660c093b3d3f86991c33f.jpg)
![](/img/5/a33/aba9a26c39b71368692b22a0a4fc.jpg)
- 命題
假定複數的
矩陣
和
,西爾維斯特方程針對任意
有唯一解,若且唯 若
和
沒有共同的特徵值。
![](/img/1/18a/6387f27ec8e0736d2b2ecca25d2a.jpg)
![](/img/2/c55/52f89362b74bbd06d301d59bf879.jpg)
![](/img/9/e85/38dcd23d2260794419bbf0f8bc00.jpg)
![](/img/4/149/14e7afa211bd7d74856fa0a2acf2.jpg)
![](/img/b/9fb/75b7b170a44c62b1ba339d984dde.jpg)
![](/img/9/663/0444d2124a876f7f0c5e10dc64aa.jpg)
- 證明
考慮線性轉換
,![](/img/4/ad0/8dedcc0d07d167407c3aec3625aa.jpg)
![](/img/d/775/26cf3e2e1ebe04c80bbab4c55219.jpg)
![](/img/4/ad0/8dedcc0d07d167407c3aec3625aa.jpg)
(i)假設
和
沒有共同的特徵值,則其特徵方程式
和
的最大公因式為
,因此存在複數多項式
和
,使得
。依照Cayley–Hamilton定理,
;因此
。令
為
的解,則
,重複上述作法,可得
。因此依照秩-零化度定理,
是可逆的,因此針對所有的
都存在唯一的解
。
![](/img/1/483/3202055a658c9ea7b7845232b848.jpg)
![](/img/c/3bb/c670412c6d858ef168f5d82ec8cc.jpg)
![](/img/f/fca/5f502e5a5687a9b39550516898f8.jpg)
![](/img/a/d69/f0da48f2d5c808a942c703cbda4b.jpg)
![](/img/1/a2c/ce5e94019e0ed3d01313f625092b.jpg)
![](/img/8/79f/ac6825c20e3f4491c76667a49477.jpg)
![](/img/b/fea/02b789eb3638882e6951627c239d.jpg)
![](/img/6/987/a8e10b175cc754b722f2fced594b.jpg)
![](/img/4/5e5/6512bccf410f30c6c58de1dbe30d.jpg)
![](/img/0/849/eb0ad45335eb88e69e7d39ca90a8.jpg)
![](/img/2/5a1/08d868f9f082297eec9fcc86ff55.jpg)
![](/img/1/629/44851d50c92211be57f0bdebd9d2.jpg)
![](/img/0/e48/45572dc1bf8c1f2159bbd0fc4b03.jpg)
![](/img/e/093/1c6e692ff046d78d92eca442951a.jpg)
![](/img/d/d9c/1fdf7deef0c2a0ff7515768ff77c.jpg)
![](/img/0/4d7/b5f1167a69ac1965c3103b4fcd3e.jpg)
![](/img/c/059/6d2b335a874040b3b5675113d75a.jpg)
(ii) 相對的,若假設
是
和
的共同特徵值,則
也是
的特徵值。存在非零向量
和
使得
以及
。選擇
使得
,向量的元素是
的共軛複數,則
沒有解
,因為複數的雙線性pairing
,等號的右邊為正值,而左側為零。
![](/img/3/c42/5dffd74d12779d2fb1ec4cebe83b.jpg)
![](/img/d/38d/0baa75fc64640f0f8b6d0765e24d.jpg)
![](/img/5/20b/57270d9e22edf10ee551f51fe971.jpg)
![](/img/2/883/7ff70e030736d1f2d2a5282ab371.jpg)
![](/img/b/9ca/105b7ca48976fb48d07511db0739.jpg)
![](/img/7/265/37381eafb7c2fd7c015d0f0e2d63.jpg)
![](/img/f/bda/146e63035d972437e1f696da8447.jpg)
![](/img/4/a70/97f490f5f202cca6a8d41a5564e3.jpg)
![](/img/b/8d2/586157315b63e8c5ae3a2b52cf7c.jpg)
![](/img/6/fbe/b483bd23d409f4eefe8f56247776.jpg)
![](/img/6/105/c1a5d4a38bf28493abb4d94d95ac.jpg)
![](/img/2/cce/4003d53a4d375816dc23aab38839.jpg)
![](/img/0/ced/9b9d01bea25e0dfdefee8d29df0f.jpg)
![](/img/c/738/72964d29dc6af9f213af464da504.jpg)
![](/img/e/9aa/eba4dc8c74efe0eadedccd1ef2aa.jpg)
Roth消去法則
假設二個大小分別為n和m的方陣A和B,以及大小為n乘m的矩陣C,則可以確認以下二個大小為n+m的方陣
和![](/img/a/499/87e6c33a172a522d7337dcb136ca.jpg)
![](/img/a/e91/d1426f954427ef169d90c56b3b73.jpg)
![](/img/a/499/87e6c33a172a522d7337dcb136ca.jpg)
是否彼此相似。這二個矩陣相似的條件是存在一矩陣X使得AX-XB=C,換句話說,X為西爾維斯特方程的解,這稱為Roth消去法則(Roth's removal rule)。
可以用以下方式檢查,若AX-XB=C,則
![](/img/4/544/e5c6fbbb607674b7504b60d3de5d.jpg)
Roth消去法則無法延伸到巴拿赫空間中的無窮維有界運算元中。
數值解
西爾維斯特方程數值解的經典算法是Bartels–Stewart算法,利用QR算法將矩陣
和矩陣
轉換為舒爾形式,再用逆向取代法求解三角矩陣。此算法若用LAPACK計算,或是GNU Octave的lyap函式計算,計算複雜度是
個數學運算。也可以參考其中的sylvester函式。在一些特定的影像處理套用中,西爾維斯特方程會有解析解。
![](/img/1/d5d/f515d3fba5021b817281953edb78.jpg)
![](/img/a/00c/09789dc921b4b9a457a32a1bc267.jpg)
![](/img/2/cb0/9c082d811d8678fab7e01a061c29.jpg)