西尔维斯特方程(Sylvester equation)是控制理论中的矩阵方程。
定义西尔维斯特方程(Sylvester equation)是控制理论中的矩阵方程,形式如下:
其中A、B及C是已知的矩阵,问题是要找出符合条件的X。其中所有矩阵的系数都是复数。为了要使方程成立,矩阵的行和列需要满足一定条件,A和B都要是方阵,大小分别是n和m,而X和C要是n行m列的矩阵,n和m也可以相等,四个矩阵都是大小相同的方阵。
西尔维斯特方程有唯一解X的充份必要条件是A和-B没有共同的特征值。
AX+XB=C也可以视为是(可能无穷维中)巴拿赫空间中有界算子的方程。此情形下,唯一解X的充份必要条件几乎相同:唯一解X的充份必要条件是A和-B的谱互为不交集。1
解的存在及唯一 利用克罗内克积以及向量化量子,可以改写西尔维斯特方程为
其中 为 单位矩阵。在此形式下,可以将问题改为 维的线性系统。
命题
假定复数的 矩阵 和 ,西尔维斯特方程针对任意 有唯一解,若且唯 若 和 没有共同的特征值。
证明
考虑线性转换 ,
(i)假设 和 没有共同的特征值,则其特征方程式 和 的最大公因式为 ,因此存在复数多项式 和 ,使得 。依照Cayley–Hamilton定理, ;因此 。令 为 的解,则 ,重复上述作法,可得 。因此依照秩-零化度定理, 是可逆的,因此针对所有的 都存在唯一的解 。
(ii) 相对的,若假设 是 和 的共同特征值,则 也是 的特征值。存在非零向量 和 使得 以及 。选择 使得 ,向量的元素是 的共轭复数,则 没有解 ,因为复数的双线性pairing ,等号的右边为正值,而左侧为零。2
Roth消去法则假设二个大小分别为n和m的方阵A和B,以及大小为n乘m的矩阵C,则可以确认以下二个大小为n+m的方阵 和
是否彼此相似。这二个矩阵相似的条件是存在一矩阵X使得AX-XB=C,换句话说,X为西尔维斯特方程的解,这称为Roth消去法则(Roth's removal rule)。
可以用以下方式检查,若AX-XB=C,则
Roth消去法则无法延伸到巴拿赫空间中的无穷维有界算子中。1
数值解西尔维斯特方程数值解的经典算法是Bartels–Stewart算法,利用QR算法将矩阵 和矩阵 转换为舒尔形式,再用逆向取代法求解三角矩阵。此算法若用LAPACK计算,或是GNU Octave的lyap函数计算,计算复杂度是 个数学运算。也可以参考其中的sylvester函数。在一些特定的影像处理应用中,西尔维斯特方程会有解析解。2
本词条内容贡献者为:
李斌 - 副教授 - 西南大学