毕业设计（论文）-矩阵方程AX+XB=D的极小范数最小二乘解.docVIP

题目：矩阵方程题目：矩阵方程ax+xb=dax+xb=d的极小范数最小二乘解的极小范数最小二乘解 姓 名 学 号 指导教师 专 业 学 院 河北工程大学毕业设计（论文） 1 摘摘 要要 矩阵理论既是学习经典数学的基础，又是一门最有实用价值的数学理论。它不仅 是数学的一个重要的分支，而且也已经成为现代各科技领域处理大量有限维空间形式 与数量关系的强有力的工具。特别是计算机的广泛应用，为矩阵论的应用开辟了广阔 的前景。例如，系统工程、优化方法以及稳定性理论等，都与矩阵论有着密切的联系。 当前，在矩阵理论领域，对矩阵方程解的研究一直是最热点的问题之一，矩阵方 程及其解的问题在生物学、电学、光子光谱学、振动理论、线性最优控制等很多领域 都有重要的应用。矩阵方程解的研究与探索更是没有间断过，可见该方dxbax 程的求解问题确实是一个非常重要的课题。 本文将围绕该命题展开讨论，利用矩阵的直积（kronecker 积） 、按行拉直和矩阵 的满秩分解等算法解出矩阵方程的极小范数最小二乘解，并且它可由dxbax moore-penrose 逆表出。 a 关键词：关键词：矩阵方程 极小范数最小二乘解 矩阵的直积（kroneck 积） 满秩分解算 法 moore-penrose 逆 a 河北工程大学毕业设计（论文） 2 abstract matrix theory is the foundation of learning classical mathematics, as well as one of the most meaningful mathematical theories. it is not only an important branch of mathematics, but also has become a powerful tool of handling many relationships between the finite dimension space structures and quantities in varity fields of modern technology. the extensive use of computer has brought a bright prospect to the application of matrix theory. many problem have close relation with matrix theory, such as system engineering, optimization method, stability theory, and so on. nowadays, the research on matrix equations has been turning into one of the hottest topics in matrix theory . it is widely used in different areas, for example, biology, electrics, spectroscopy, vibration theory, linear optimal control, etc. and people have researched and explored the solution of matrix equation all the time, which means the subject dxbax study is absolutely crucial. this paper discusses the available solutions of matrix equation, using kronecker product, full rank decomposition of matrix to obtain the minimal norm least squares solution of matrix equation . and it can be expressed by moore-penrose inverse.dxbax a key words: matrix equation; the minimum norm least squares solution; kronecker product; full rank factorization algorithm; moore-penrose inverse a 河北工程大学毕业设计（论文） 3 目目 录录 摘 要1 abstract.2 1 绪论4 1.1 李雅普诺夫矩阵方程dxbax应用背景及研究现状.5 1.2 本文的主要工作5 1.3 符号说明6 2 预备知识7 2.1 矩阵的范数7 2.1.1 矩阵范数的定义8 2.1.2 矩阵范数的性质8 2.2 矩阵的直积及其应用10 2.2.1 直积的概念10 2.2.2 矩阵直积的性质11 2.2.3 线性矩阵方程的可解性11 2.3 矩阵的满秩分解.12 2.3.1 矩阵满秩分解的定义12 2.3.2 矩阵满秩分解的性质12 2.4 广义逆矩阵的存在和性质.14 2.4.1 penrose 的广义逆矩阵定义.14 2.4.2 广义逆矩阵的性质14 2.4.3 moore-penrose 逆矩阵的计算.15 3 矩阵方程 dxbax 的极小范数最小二乘解16 3.1 广义逆矩阵与线性方程组的求解16 3.1.1 线性方程组极小范数最小二乘解的定义16 3.1.2 矩阵方程转化为线性代数方程组17 3.2 矩阵方程dxbax的求解.18 3.2.1 矩阵方程 dxbax 在相容条件下解的情况19 3.2.2 矩阵方程 dxbax 不相容时解的情况.21 结论22 致谢23 参考文献24 河北工程大学毕业设计（论文） 4 1 绪论绪论 矩阵理论既是学习经典数学的基础，又是一门最有实用价值的数学理论。它不仅 是数学的一个重要的分支，而且也已经成为现代各科技领域处理大量有限维空间形式 与数量关系的强有力的工具。在矩阵理论领域，对矩阵方程解的研究一直是最热点的 问题之一，而矩阵方程解的研究与探索更是没有间断过，但几乎所有的dxbax 结论要么基于矩阵为特殊矩阵（如对称矩阵等）的情形，要么解的表达形式过于ba, 繁琐。此外，还有一些文献仅仅只是讨论了解的存在性。而在力学以及控制论等诸多 领域利亚普诺夫方程都有着很重要的应用。可见该方程的求解问题确实是一个非常重 要的课题，所以，考虑给出该方程的极小范数最小二乘解的表达式也是很有必要的。 本文将围绕该命题展开讨论，并且最终给出它的重要应用矩阵方程 的最小二乘解以及极小范数最小二乘解。dxbax 1.1 李雅普诺夫矩阵方程李雅普诺夫矩阵方程应用背景及研究现状应用背景及研究现状dxbax 在科学与工程中，经常会遇到求解线性方程组的问题。矩阵是描述和求解线性方 程组最基本和最有用的数学工具。矩阵有很多基本的数学运算，如转置、内积、外积、 逆矩阵、广义逆矩阵等。 矩阵方程及其解的问题在生物学、电学、光子光谱学、振动理论、有限元、结构 设计、参数识别、自动控制理论、线性最优控制等很多领域都有重要的应用，正是这 些领域提出了许多不同类型的线性矩阵方程的模型问题刺激了理论的快速发展,使得线 性矩阵方程的求解问题成为当今计算数学领域的热门研究课题之一，经过国内外的专 家和学者的不断探索, 迄今为止，线性矩阵方程问题的研究已取得了一系列丰硕的成果。 对于矩阵方程, 1955 年, penrose 得到了它有一般解的充要条件和通解dxbax 表达式；1994 年，黄礼平博士研究了四元数体上方阵的标准形与矩阵方程 的问题； 2002 年，马飞博士研究了矩阵方程的迭代解法问dxbaxdxbax 题, 得出方程在不同精度下的解并给出迭代所需步数；2004 年, 袁永新教授研究了矩阵 方程的最优解问题。dxbax 1.2 本文的主要工作本文的主要工作 论文第一章简要的介绍了利亚普诺夫方程研究的实际背景及研究现状；第二章介 绍了已有的主要结果、范数理论、矩阵的直积（kronecker 积） 、按行拉直和矩阵的满 秩分解算法等。第三章为本文的主要工作，本文根据直积的定义及性质可以将矩阵方 程拉直为一般的线性方程组，以此为基础，求解经过按行向量拉直转化dxbax 得到的线性方程组，解出该方程的极小范数解和最小二乘解，虽然最小二乘解一般不 河北工程大学毕业设计（论文） 5 是唯一的，但是极小范数最小二乘解确实唯一的，并且它可由 moore-penrose 逆表 a 出。若矩阵方程不相容，则它的极小范数最小二乘解需满足dxbax |min |min x dxbax 求出的即为矩阵方程极小范数最小二乘解。第四章给出了结论，最后xdxbax 是本文的参考文献与致谢。 1.3 符号说明符号说明 实数域 实维向量空间n 实矩阵空间nm 秩为 的实矩阵的集合rnm 复数域 复维向量空间n 复矩阵空间nm 秩为 的复矩阵的集合rnm 矩阵的行列式a 维线性空间n 矩阵的值域，的列空间aa 矩阵的秩a 单位变换 向量与向量的内积 x y 向量生成的子空间 n xxx, 21 矩阵的转置a 矩阵的共轭转置a 矩阵的 行列处的元素aij 矩阵的任意范数a 矩阵的 frobenius 范数a 矩阵的条件数a 向量的-范数xp 矩阵的第 个特征值ai r n r nm r nm r r c n c nm c nm r c adet n v )(ar ranka e t ),(yx ),( 21n xxxl t a h a ij a)( a f a )(acond p x )(a ii 或 河北工程大学毕业设计（论文） 6 矩阵与矩阵的直积ab 矩阵按行拉直所得到的列向量a 矩阵的 moore-penrose 逆 a 矩阵的群逆 a ba )(avec a # a 河北工程大学毕业设计（论文） 7 2 预备知识预备知识 2.1 矩阵的范数矩阵的范数 在计算数学中，特别是在数值代数中，研究数值方法的收敛性、稳定性及误差分 析等问题时，范数理论都显得十分重要。本章主要讨论矩阵空间中的矩阵范数的 nm c 理论及其性质。 矩阵空间是一个维的线性空间，将矩阵 a 看做线性空间中的 nm c mnnm nm c “向量” ，可以按照定义-范数的方式定义 a 的范数。但是，矩阵之间还有乘法运算，p 它应该在定义范数时予以体现。 2.1.1 矩阵范数的定义矩阵范数的定义 定义定义 2.1 设，定义一个实值函数, 它满足以下三个条件a nm c a （1）非负性：当时，0; 当时，;0aa0a0a （2）齐次性：;caaaaa, （3）三角不等式：, ,则称为 a 的广义矩阵范数。babab nm c a 若对，及上的同类广义矩阵范数，有 nm c ln c lm c * . （4）相容性： , baabb ln c 则称为的矩阵范数。aa 2.1.2 矩阵范数的性质矩阵范数的性质 同向量的情况一样，对于矩阵序列也有极限的概念：设有一个矩阵序列，其 )(k a 中，。用记的第 行第列的元素，且都有极限，则 )(k a nm c ,.2 , 1k )(k ij a )(k aij )(k ij a ij a 称有极限，或称收敛于矩阵,记为 )(k a)( ij aa )(k aa = 或 lim x )(k aa )(k aa 不收敛的矩阵序列称为发散的。于是可以证明：的充要条件是。 )(k aa0 )( aa k 由定义 2.1 的条件（3） ，可以证明下列不等式 | |ba ba 由此可以证明矩阵范数的连续性，即由 )(k aa 可以推出 aa k )( 事实上，由上面的论述知，当时，,但是 )(k aa0 )( aa k 河北工程大学毕业设计（论文） 8 |aa k )( aa k )( 于是当时，便有。0 )( aa k aa k )( 推论推论 2.1 已知 a=（，可证明下面二函数 nn ij ca ) ， n ji ij m aa 1, 1 ij ji m ana , max 都是上的矩阵范数。 nn c 证 对于函数而言，它显然具有非负性与齐次性，现仅就三角不等式与相容 1 m a 性加以验证于下 )( 1,1, 1 ij n ji ij n ji ijij m bababa + n ji ij n ji ij ba 1,1, 1 m a 1 m b n ji njinjiji m bababaab 1, 2211 1 )( 1 1, 1njinj n ji i baba n j njj n i ini bbaa 1 1 1 1 n ij ij n ji ij ba 1,1, 1 m a 1 m b 因此，是 a 的矩阵范数。 1 m a 同理可证，也是 a 的矩阵范数。 m a 如同向量范数一样，矩阵范数也是多种多样的。但是，在数值方法中进行某种估 计时，遇到的多数情况是，矩阵范数常与向量范数混合在一起使用，而矩阵经常是作 为两个线性空间上的线性映射（变换）出现的。因此，考虑一些矩阵范数时，应该使 他能与向量范数联系起来。这可由矩阵范数与向量范数相容的概念来实现。下面引入 这一概念。 定义定义 2.2 对于上的矩阵范数和与上的同类向量范数, 如果 nm c m m c n c v , vmv xaax xca nm , n c （2.1） 则称矩阵范数与向量范数是相容的。 m v 定理定理 2.1 设 a，且与都是酉矩阵，则 nm c p nm c nm cq fff aqapa 即给 a 左乘或右乘以酉矩阵后其值不变（在时，和 q 都是正交矩阵） 。 f nm ra p 河北工程大学毕业设计（论文） 9 证 若记 a 的第 j 列为，则有), 3 , 2 , 1(njaj 2 21 2 ),( f n f aaappa 2 21 ),( f n papapa 2 2 2 1 2 2 1 f n j j n j j aapa 即。于是 ff apa f f h f hh f h f aaaqaqaq)( 推论推论 2.2 和 a 酉（或正交）相似的矩阵的 f-范数是相同的，即若,则aqqb h , 其中 q 是酉矩阵。 ff ab 2.2 矩阵的直积及其应用矩阵的直积及其应用 矩阵的直积（kronecker 积）在矩阵的理论研究和计算方法中都有十分重要的应用。 特别地，运用矩阵的直积运算，能够将线性矩阵方程转化为线性代数方程组进行讨论 或计算。 2.2.1 直积的概念直积的概念 首先从简单的例子开始。设有二元向量和三元向量,他们分别 t ),( 21 t ),( 321 经过二阶矩阵和三阶矩阵，即 ， 2221 1211 aa aa a 333231 232221 131211 bbb bbb bbb b 的变换，变成向量和，即有 t ),( 2 1 t ),( 3 2 1 ， 2 1 2 1 a 3 2 1 3 2 1 b 现在考虑以这两个向量的分量乘积为分量的六元向量 t u),( 322212312111 经过怎样的线性变换可以变成六元向量 ),( 3 2 2 2 1 2 3 1 2 1 1 1 v 由假设 2 , 1, 2211 iaa iii 3 , 2 , 1, 332211 jbbb jjjj 故有 22221212313121211111 jijijijijiji bababababa 河北工程大学毕业设计（论文） 10 + （） 3232 ji ba3 , 2 , 1, 2 , 1ji 于是所求变换的矩阵为六阶矩阵 baba baba 2221 1211 一般地，引进以下的定义。 定义定义 2.3 设则称如下的分块矩阵,)(,)( qp ij nm ij cbbcaa （2.2） 为的直积（kronecker 积） 。ba与 是一个块的分块矩阵，所以式（2.2.1）还可简写为banm = ba)(baij （2.3） 2.2.2 矩阵直积的性质矩阵直积的性质 矩阵的直积具有以下性质。 （1）)()()(kbabkabak （2）设为同阶矩阵，则 21 aa与 bababaa 2121 )( 2121 )(ababaab （3）)()(cbacba （4）设， 2211 )(,)( )2( 2 )1( 1nmijnmij aaaa 2211 )(,)( )2( 2 )1( 1qpijqpij bbbb 且，则 2121 ,pqmn 。)()()( 21212211 bbaababa （5）设都可逆，则 nnmm ba 与 111 ) baba（ （6）设都是上三角（下三角）矩阵，则也是上三角（下三角） nnmm ba 与ba 矩阵。 （7） hhh baba)( （8）设与都是正交（酉）矩阵，则也是正交（酉）矩阵。 nm a nm b ba 2.2.3 线性矩阵方程的可解性线性矩阵方程的可解性 在系统控制等工程领域，经常遇到矩阵方程 ax+xb=d （2.4） nqmp mnmm n n c bababa bababa bababa ba 21 22221 11211 河北工程大学毕业设计（论文） 11 的求解问题，其中，为以知矩阵，而,为未知矩a mm c nn cb nn cd nm cx 阵，一般地线性矩阵方程可表示为 fxba i l i i 1 （2.5） 其中为已知矩阵，而为未知矩阵。 nmnq i pm i cfcbca , qp cx 下面利用矩阵直积的性质，研究矩阵方程（2.5）的可解性问题。 设，并记的第 行为及 nq pmij cbaa ,)( qp cx i),.,2 , 1(pix t i t t p tt xxxxvec),()( 21 则有 bxaxaxaaxb t pip t i t i )( 2211 从而 )()( 2211pipii t xaxaxabaxbvec = p t mp t m t p t x x baba baba 1 1 111 =)()(xvecba t 定理定理 2.2 方程（2.5）有解的充要条件是，这里，表)()( 1 l i t ii barfvec)(ar 示矩阵a的列空间。 定理定理 2.3 设的特征值为，的特征值为，则方程 mm a m , 21 nn b n , 21 有惟一解的充要条件。fxbax), 2 , 1, 2 , 1(0njmi ji 2.3 矩阵的满秩分解矩阵的满秩分解 本节论述将非零矩阵分解为列满秩矩阵与行满秩矩阵的乘积问题。 2.3.1 矩阵满秩分解的定义矩阵满秩分解的定义 定义定义 2.4 设，如果存在矩阵和，使得)0( rca nm r rm r cf nr r cg (2.6)fga 则称式（2.6）为矩阵的满秩分解。a 当是满秩（列满秩或行满秩）矩阵时，可分解为一个因子是单位矩阵，另一aa 个因子是本身，称此满秩矩阵分解为平凡分解。a 河北工程大学毕业设计（论文） 12 2.3.2 矩阵满秩分解的性质矩阵满秩分解的性质 矩阵的满秩分解有下列性质： 定理定理 2.4 设，则有满秩分解（2.6） 。)0( rca nm r a 证 当时，根据矩阵的初等变换理论，对进行初等变换，可将化为rranka aa 阶梯形矩阵，即b nr cg g ba r , 0 行 于是存在有限个阶初等矩阵的乘积，记作,使得mp 或者 bpa bpa 1 将分块为 1 p , , 1 sfp rm r cf )(rnm rn cs 则有 1 0 gf g sfbpa 其中是列满秩矩阵，是行满秩矩阵。 证fg 毕 需要指出的是，矩阵的满秩分解（2.6）不是唯一的。这是因为若取是任意一ad 个 阶非奇异矩阵，则式（2.6）可改写为r 1 )(gfgdfda 这是的另一个满秩分解。a 在求列满秩矩阵时，需要求出矩阵及其逆矩阵，这是十分麻烦的。为了避fp 1 p 免这些运算，引入下面的定义。 定义定义 2.5 设, 且满足)0( rcb nm r （1）的前 行中每一行至少含一个非零元素，且第一个非零元素是 1，而后br 行元素均为零;rm （2）若中第 行的第一个非零元素 1 在第列（） ，则bi i jri, 2 , 1 ；r jjj 21 （3）中的列为单位矩阵的前 列，那么就称为 hermite 标准形。br jjj, 21 m irb 定义定义 2.6 以阶单位矩阵的个列向量为列构成的阶矩阵n n in n eee, 21 n ),( 21jnjj eeep （2.7） 称为置换矩阵，这里是 1,2，的一个排列。 n jjj, 21 n, 定理定理 2.5 设的 hermite 标准形为（如定义 2.5） ，那么，在的)0( rca nm r ba 河北工程大学毕业设计（论文） 13 满秩分解（2.6）中，可取为的列构成的矩阵，为的前 行构far jjj, 21 rmgbr 成的矩阵。nr 证 由知，存在阶可逆矩阵，使得ba 行 mp 或者 bpa bpa 1 根据定理，将分块为 1 p , , 1 sfp rm r cf )(rnm rn cs 可得满秩分解，其中为的前 行构成的矩阵。fga gbrnr 下面确定列满秩矩阵。参照的 hermite 标准形，作阶置换矩阵fabn ),( 11 1 nrr jjjj eeeep 划分 ， ，则有 ),( 11 1 nrr jjjj aaaaap oo bi bbbbbp r jjjj nrr 12 1 ),( 11 其中。再由可得 )( 12 rnr cb bpa 1 121 1112 () r ib appbpf sf fb oo 即为的前 列构成的矩阵，也就是的列构成的矩阵。 证f 1 aprar jjj, 21 毕 2.4 广义逆矩阵的存在和性质广义逆矩阵的存在和性质 2.4.1 penrose 的广义逆矩阵定义的广义逆矩阵定义 定义定义 2.7 设矩阵，若矩阵满足如下四个 penrose 方程 nm ca mn cx (i) aaxa (ii) xxax (iii) axax h )( (iv) xaxa h )( 则称为的 moore-penrose 逆，记为。xa a 定理定理 2.6 对任意，存在且唯一。 nm ca a 证 先证存在性 设。若，则是零矩阵，可以验证零rranka 0ranmmn 矩阵满足四个 penrose 方程。若，由定理知，可进行奇异值分解0ra ),( 21n aaaa),( 21n bbbb 河北工程大学毕业设计（论文） 14 , h udva oo o d r 1 其中是的奇异值，和分别是阶和阶酉矩阵。容易验证), 1(0ri i auvmn h mn r u oo o vx 1 1 1 满足四个 penrose 方程。可见，总是存在的。 a 唯一性 设均满足方程(i)(iv),则yx, hhhhh ayaxxaxxaxxx)()()( yyayyaxaxay hhhh )()( yyya h )( 证毕 2.4.2 广义逆矩阵的性质广义逆矩阵的性质 广义逆矩阵有以下性质 定理定理 2.7 设 ，则 nm ca pn cb c （1）;1)( )1(hh aa （2）;1 )1(h aa （3）若和非奇异，则;st1)( 1)1(1 satsat （4）;rankaranka )1( （5）和均为幂等矩阵且与同秩； )1( aaaa )1( a （6） 定理定理 2.8 给定矩阵和,则的充要条件是a1ax 2 , 1ax rankarankx 证 充分性 显然, 若，由定理 2.6 之（5）得)()(xrxarrankarankx rankarankxxarank)( 所以，即存在矩阵使得)()(xrxary xayx 从而 xxayxayxaxax)( 必要性 因为和互为1,2-逆，有定理（2.4.2）之（4）即得。 证毕ax )()( ),()(),( )1( )1()1( hh araar anaanaraar 河北工程大学毕业设计（论文） 15 定理定理 2.9 设矩阵给定，则a （1）rankaranka （2）aa ) ( （3）; tthh aaaa)()( ,)()( （4）; aaaaaaaa hhhh )()( ,)()( （5）; )()( hhhh aaaaaaa （6）)()(),()( hh ananarar 推论推论 2.2 若，则 nm ca (2.8) hh aaaa 1 )( 若，则 nm m ca (2.9) 1 )( hh aaaa 推论推论 2.3 设为维列向量，且，则n0 (2.10) hh 1 )( 而 (2.11) 1 )()()( hhh 2.4.3 moore-penrose 逆矩阵的计算逆矩阵的计算 假定矩阵的秩为 ，其中，下面介绍求 moore-penrose 逆矩阵nmar),min(nmr 的一种方法，满秩分解法。 a 前面已经介绍过矩阵的满秩分解的概念，并给出了用初等变换进行满秩 nm r ca 分解的方法。设的满秩分解为)0( rca nm r (2.12)fga 其中，那么可以按照下述定理给出的结论计算广义逆矩阵。 rm r cf nr r cg 定理定理 2.10 设的满秩分解为式（2.12） ，则)0( rca nm r （1）; 4 , 2 , 1, )1()( iiafg i （2）; 3 , 2 , 1, )()1( iiafg i （3）,3 , 2 , 1 )1( afg ;4 , 2 , 1 )1( afg （4） fgfga )4, 1()3 , 1( （5） hhhh fffgggfga 11 )()( hhhh fagfg 1 )( 证 由定理知 r iggff )1()1( 根据定义容易验证（1） ， （2）成立。 河北工程大学毕业设计（论文） 16 （3）可由（1） ， （2）直接得到。 （其中分别换成时，结论 gf 与 ）（）（ 与 4, 2, 13 , 2, 1 gf 仍成立。 ） （4）由（3）知，,又因为4 , 2 , 1 )3 , 1( afgx axffffffggax hhh )3 , 1()3 , 1()3 , 1( )()()( 所以4 , 3 , 2 , 1ax 另一式可类似的证明。 （5）由（4）知，,载根据式和可推 afg hh aaaa 1 )( 1 )( hh aaaa 得其他结论。 3 矩阵方程矩阵方程 dxbax 的极小范数最小二乘解的极小范数最小二乘解 3.1 广义逆矩阵与线性方程组的求解广义逆矩阵与线性方程组的求解 3.1.1 线性方程组极小范数最小二乘解的定义线性方程组极小范数最小二乘解的定义 考虑非齐次线性方程组 bax （3.1） 其中给定，而为待定向量。如果存在向量使方程组（3.1）成 nnm cbca , n cxx 立，则称方程组相容，否则称为不相容或矛盾方程组。 关于线性方程组的求解问题，常见的有以下几种情形。 （1）方程组（3.1）相容的条件是什么？在相容时求出其通解（若解不唯一的话） 。 （2）如果方程组（3.1）相容，其解可能有无穷多个，求出具有极小范数的解，即 x bax min （3.2） 其中是欧氏范数。可以证明，满足该条件的解是唯一的，称之为极小范数解。 （3）如果方程组（3.1）不相容，则不存在通常意义下的解。但在许多实际问题中， 需要求出极值问题 bax n cx min （3.3） 河北工程大学毕业设计（论文） 17 的解,其中是欧氏范数。称这个极值问题为求矛盾方程的最小二乘问题，相应的x 称为矛盾方程的最小二乘解。x （4）一般说来，矛盾方程的最小二乘解是不唯一的，但在最小二乘解得集合中， 具有极小范数的解 x baxmin min （3.4） 是唯一的，称之为极小范数最小二乘解。 广义逆矩阵与线性方程的求解有着极为密切的关系。利用广义逆矩阵可以给出上 述诸多问题的解，反之由线性方程的解又可以确定广义逆矩阵。 3.1.2 矩阵方程转化为线性代数方程组矩阵方程转化为线性代数方程组 对比较复杂的矩阵方程的求解，我们通常要用按行向量拉直的方法将矩阵方程转 化为线性方程组，把矩阵方程化成的形式，然后再讨论矩阵方程解的情况。bax 设，的第 行为,则可表示为 nqqppm cbcxca ,xi t i xx t p t t x x x x 2 1 t t p tt p xxx x x x xvec 21 2 1 )( 由此可得 b x x x aaa aaa aaa axb t p t t mpmm p p 2 1 21 22221 11211 b xaxaxa xaxaxa t pmp t m t m t pp tt 2211 1212111 bxaxaxa bxaxaxa t pmp t m t m t pp tt )( )( 2211 1212111 所以 )( )( )( 2211 1212111 pmpmm t pp t xaxaxab xaxaxab axbvec pmp t m t m t pp ttt xabxabxab xabxabxab 2211 1212111 河北工程大学毕业设计（论文） 18 p t mp t m t m p t p tt xbaxbaxba xbaxbaxba 2211 1212111 )()( 1 1 111 xvecba x x baba baba t o t mp t m t p t 由上面推理可得 )()()(xvecbiiaxbaxvec t 即矩阵方程按行向量拉直后为dxbax )()()(dvecxvecbiia t 令为, ，则矩阵方程可写成)( tt biiagxxvec为)(cdvec为)(dxbax 等价形式 cgx 3.2 矩阵方程矩阵方程的求解的求解dxbax 考虑将矩阵方程（其中） ，按行向量拉直后有dxbax nn cdxba , ，即)()()(dvecxvecbiia t nn t t t nnn n d d d x x x b b b iaia iaia 2 1 2 1 1 111 00 0 0 00 这里表示的第 行，表示的第 行。 t i xxi t i ddi 经整理得 nn t nnnnn nn t n t d d d x x x biaiaia ia biaia iaiabia 2 1 2 1 1,1 , 1 2221 11211 3.2.1 矩阵方程矩阵方程 dxbax 在相容条件下解的情况在相容条件下解的情况 在本节主要考虑矩阵方程在相容条件下解的情况，并且其中dxbax 与。由前面分析知道，这时矩阵方程成为 mm ca nn cb )()()(dvecxvecbiia t mn 可知该矩阵相容的条件是,它有唯一解)(,()(dvecbiiarbiiar mn t mn 当且仅当矩阵为非奇异矩阵。)( t mn biia 河北工程大学毕业设计（论文） 19 定理定理 3.1 设与，则当且仅当,其中分别是 mm ca nn cb 0 ji ji 与 的特征值，矩阵方程有唯一解。ba和dxbax 证明 矩阵方程有唯一解dxbax 有唯一解)()(dvecxbaxvec 有唯一解)()()(dvecxvecbiia t mn )(xvec 0)det( t mn biia 令，并注意与的特征值完全相同，由 1001 ),(yxyxyxf t bb n t iabaf),( 的特征值为 ，故 t m bi jiji f),( 0)det( t mn biia0 ji 其中。 证njmi, 2 , 1;, 2 , 1 毕 推论推论 3.1 设的特征值为，的特征值为，则齐次方 mm a m , 21 nn b n , 21 程有非零解的充要条件是存在使。0 xbax 00 ji 与0 00 ji 虽然上述定理给出了矩阵方程有唯一解的充分与必要条件，但要写dxbax 出它们解矩阵的明显表达式一般是不容易的。若与为稳定矩阵，x mm ca nn cb 即所有特征值的实部小于零的矩阵，则的解有明显的表达式。dxbaxx 定理定理 3.2 设与为稳定矩阵，且给定。则线性矩阵方程 mm ca nn cb nm cd 有唯一解，并且，dxbaxx dtdeex btat 0 （3.5） 证 本定理的假定保证了与-没有相同的特征值，因而按定理 3.1，ab 有唯一解。现考虑初值问题dxbaxx ,)()()(btztaztzdz)0( （3.6） 式中，为定义在上的矩阵值函数。直接验证可得，为 nm ctz )(), 0 btatce etz)( 问题（3.6）的解。自对方程（3.6）的两端求积分便有tt到0 , 00 )()()0()(bdttzdttzazz 即有 .)()()( 00 zdbdtceedtceea btatbtat 因此，为了证明由（3.5）式确定的为的解，只需验证上式中xdxbax 。为此只要证明：对任意稳定矩阵成立。但按普分0lim)( btat t ceez0lim at t ea 确定理，我们有 ， ， kj t m j j s k at eete k k 1 01 0t 河北工程大学毕业设计（论文） 20 式中，为不同的特征值， 分别为它们的指标，为的分 s , 21 a s mmm, 21 kj ea 量。现令，分别为的实部和虚部。按为稳定矩阵的假设， kkk ia kk a与 k a ，因而skak, 2 , 1, 0 .tetet tajtj kk , 0 于是，。所以结论成立。0lim at t e 前面我们讨论了有解与有唯一解的条件，现在讨论当此方程有解时，dxbax 解的结构形式。具体步骤如下。首先取,我们得到如下特殊的矩阵方程：ocab, , xaax nn ca （3.7） 求解方程（3.7）等价于寻求与可交换的所有矩阵。显然，方程（3.7）有无限 nn ca 多个解。事实上，为复系数多项式）总是它的解。在下面定理 3.3 中将)()(apapx 给出方程（3.7）的解的表达式，接着，应用这个结果讨论解的表达式与oxbax 它的解空间的构造。由于的解可以表示为齐次方程的同解与dxbaxoxbax 非齐次方程某个特解之和，故最后我们只要讨论有解的条dxbaxdxbax 件即可。 定理定理 3.3 设有形式,这里为的 jordan nn ca 1 pjpa, 21p jjjdiagja 正规形式，分别是对应于的特征值的 jordan 块。则 p jjj, 21 a p , 21 为方程（3.7）解的充分与必要条件为,其中，是与有相 nn cx 1 pypx st yy j 同分块形式的分块矩阵。 接着考虑齐次线性矩阵方程 oxbax （3.8） 的解的形式。主要的技巧是将它归结为前面讨论过的问题。显然，方程（3.8）有解。 假定为它的解，则可以验证x bo oa io xi n m 与 （3.9） 可交换，反之，若（3.9）式中两个阶矩阵可交换，则必定为方程（3.8）的)nm（x 解。因此，对（3.9）式中这两个矩阵应用定理 3.3 的结果，我们有 ，, 11 1 1 p pk kpkka ninidiagappj ，, 11 1 1 q qr rqrrb ninidiagaqqj 则矩阵方程（3.8）的任意解有形式,其中（）是x 1 pyqx qp tsst yy , 1, )( , cmy tr k 矩阵方程的一般解。oyjyj ba 类似地，将推论应用到矩阵方程 河北工程大学毕业设计（论文） 21 n m io xi bo oa bo oa io xi n m （3.10） 我们有如下结果 推论推论 3.2 矩阵方程（3.8）的一般解呈如下形式： ， i a i ix vx 1 （3.11） 式中，为方程（3.8）的线性无关解，而 a xx, 1 ， q t st p s aa 11 （3.12） 这里，为的初等因子的初等因子的最)1 ,1 (qtpsastab s k s 与)( tr t) ( 大公因式的次数。 我们顺便指出， （3.12）式中恰为