矩阵的分解.doc

9. 矩阵的分解矩阵分解是将一个矩阵分解为比较简单的或具有某种特性的若干矩阵的和或乘积，这是矩阵理论及其应用中常见的方法。由于矩阵的这些特殊的分解形式，一方面反映了原矩阵的某些数值特性，如矩阵的秩、特征值、奇异值等；另一方面矩阵分解方法与过程往往为某些有效的数值计算方法和理论分析提供了重要的依据，因而使其对分解矩阵的讨论和计算带来极大的方便，这在矩阵理论研究及其应用中都有非常重要的理论意义和应用价值。这里我们主要研究矩阵的三角分解、谱分解、奇异值分解、满秩分解及特殊矩阵的分解等。一、矩阵的三角分解是矩阵的一种有效而应用广泛的分解法。将一个矩阵分解为酉矩阵（或正交矩阵）与一个三角矩阵的乘积或者三角矩阵与三角矩阵的乘积，这对讨论矩阵的特征、性质与应用必将带来极大的方便。首先我们从满秩方阵的三角分解入手，进而讨论任意矩阵的三角分解。定义1 如果均为正实数，则上三角矩阵称为正线上三角复（实）矩阵，特别当时，称为单位上三角复（实）矩阵。定义2如果均为正实数，则下三角矩阵称为正线下三角复（实）矩阵，特别当时，称为单位下三角复（实）矩阵。定理1设则可唯一地分解为其中是酉矩阵，是正线上三角复矩阵；或者可唯一地分解为其中是酉矩阵，是正线下三角复矩阵。推论1设则可唯一地分解为其中是正交矩阵，是正线上三角实矩阵；或者可唯一地分解为其中是正交矩阵，是正线下三角实矩阵。推论2 设是实对称正交矩阵，则存在唯一的正线上三角实矩阵，使得推论3设是正定Hermite矩阵，则存在唯一的正线上三角复矩阵，使得定理2设用表示下三角复矩阵，表示单位下三角复矩阵，表示上三角复矩阵，表示单位上三角复矩阵，表示对角矩阵，则下列命题等价：（1）的各阶顺序主子式；（2）可唯一地分解为，并且的主对角线上元素不为零；（3）可唯一地分解为，并且的主对角线上元素不为零；（4）可唯一地分解为，并且的主对角线上元素不为零；说明：若是阶满秩实方阵，则对于实矩阵、，定理2 仍然成立。例1.设，求的三角分解。解. 由所以阶方阵的三角分解对求解非其次线性方程组非常方便。比如，设方程组，系数矩阵有三角分解式，则有，于是令，有先求第一个方程组中的未知向量，然后将代入第二个方程组再求解。由于它们都是以三角矩阵为系数矩阵的方程组，所以很容易求出方程组的解，并且易于利用计算机求解。例2 用三角分解求解方程组：解：系数矩阵可以分解为代入上面的新方程组中的第一式可得：，再将此结果代入新方程组中的第二式可得：，此即所求方程的解。二、任意矩阵的三角分解前面讨论的矩阵分解仅仅是对阶方阵的三角分解，而且所分解的矩阵是可逆矩阵，下面我们将以上的矩阵分解进行推广，即讨论任意矩阵的三角分解。定义3 设是阶复（实）矩阵，如果，则称是行满秩矩阵，记为；如果，则称是列满秩矩阵，记为。定理3 当是行满秩或列满秩复矩阵时，有（1）若，则存在阶正线下三角复矩阵和阶酉矩阵，使得（2）若，则存在阶酉矩阵和阶正线上三角复矩阵，使得注：该定理表明了行（列）满秩矩阵能分解为一个酉矩阵与一个长（高）三角矩阵的乘积。下面我们进一步给出行（列）满秩矩阵能分解为一个正线三角矩阵与一个长（高）酉矩阵的乘积。记表示以个两两正交的单位向量为行组成的矩阵的集合，表示个两两正交的单位向量为列组成的矩阵的集合。定理4 （1）若，则可唯一地分解为其中是阶正线下三角矩阵，。（2）若，则可唯一地分解为其中，是阶正线上三角矩阵。说明：当是行满秩或列满秩实矩阵时，亦有类似于定理3和定理4的结论。当既不是行满秩矩阵，也不是列满秩矩阵时，则有定理5 设，则存在酉矩阵、及阶正线下三角矩阵，使得推论4设，则存在酉矩阵、及阶正线上三角矩阵，使得三、矩阵的谱分解在线性代数中，已经讨论了一个方阵的特征值和特征向量的问题，已经发现特征值有着非常重要的作用。由于相似矩阵有相同的特征值，因而人们总希望在相似矩阵中找到结构最简单的矩阵，这就是对角矩阵或Jordan标准形矩阵。下面我们将根据矩阵的特征值，进一步寻求利用简单矩阵来表示已知的矩阵，即讨论矩阵的谱分解。 1.单纯矩阵的谱分解定义1若矩阵的每个特征值的代数重复度与几何重复度相等，则称矩阵为单纯矩阵。注意到“属于每个特征值的线性无关的特征向量合起来也是线性无关的”这一事实，则可知如下定理是成立的。定理1 是单纯矩阵与对角矩阵相似。下面给出单纯矩阵的谱分解定理。定理2 设是单纯矩阵，则可分解为一系列幂等矩阵的加权和，即， （3.1）其中是的特征值。定理2中的分解式称为的谱分解，分解式中的矩阵具有如下的性质：（1）幂等性：；（2）分离性：；（3）可加性：。由这些性质容易得出：。当是的多项式或是的解析函数时，容易得到： （3.2）称上式为矩阵函数的谱分解。例3 求的谱分解。解：由（3.2）式若设由HamiltonCaylay定理可知：， （3.3）则有：由此可知，对任意的，都是矩阵的线性组合。同时由（3.3）式，当时，可知可逆，且的逆矩阵为由（3.2）式容易求得的谱分解为。 （3.4）把一个单纯矩阵分解为一系列幂等矩阵的加权和，无论从代数上，还是从几何上进行研究，都有它的方便之处。特别对于（3.2）和（3.4）的分解，在自动控制中有许多应用。更一般地，单纯矩阵的谱分解定理为：定理3设，它有个相异特征值，则是单纯矩阵存在个矩阵满足：（1）（2）；（3）。其中。值得注意的是定理中的条件（1）中的矩阵是幂等矩阵，故定理中存在个矩阵，又可看作是存在个投影算子（简称为谱算子）。例4 求矩阵的谱分解。解：首先求得的特征值为，则 ，所以：因此， 2.正规矩阵及其分解引理1 设是正规矩阵，与酉相似，则也是正规矩阵。引理2设，则存在酉矩阵，使得其中是一个上三角矩阵且主对角线上的元素为的特征值。引理3 设是三角矩阵，则是正规矩阵是对角矩阵。定理4 阶复矩阵是正规矩阵与对角矩阵酉相似，即存在阶酉矩阵，使得：其中是的个特征值。 下面给出正规矩阵的谱分解定理。定理5设，它有个相异特征值，则是正规矩阵存在个矩阵使其满足：（1）；（2）；（3）；（4）。其中。 3.与Jordan标