矩阵的分解将矩阵表示成特定类型矩阵的乘积这种表示称为.doc

第章 将矩阵表示成特定类型矩阵的乘积，这种表示称为矩阵的分解，矩阵的分解理论与方法是矩阵分析理论的重要内容，也在解线性方程组，求特征值，求广义逆矩阵的实际计算中有着重要的应用价值 4.1 矩阵的分解 定义 设，若可以表示成一个下三角矩阵与一个上三角矩阵的乘积 则称其为矩阵的分解（三角分解）。 定理 设，如果的顺序主子式 则存在惟一的主对角线上元素全为的下三角阵与惟一的上三角矩阵，使得。 证 我们对矩阵的阶数使用数学归纳法。明显地，当时，就是惟一的分解式。 现在假定对阶矩阵，定理的结论成立。我们证明对阶矩阵结论成立。 对进行分块 其中。由于阶矩阵的阶顺序主子式就是的阶顺序主子式，故它们都不为零。从而由归纳假设，有惟一的分解 其中的主对角线上的元都是。由于 所以，及是n-1阶可逆矩阵. 我们先假设已有，其中 ， 是待定向量.作乘积 则必须满足 ， 注意到及都是可逆的，则由上式可以唯一确定 ， 这就证明了的三角分解的存在性与唯一性. 矩阵的三角分解在求解线性方程组时将十分简便.如对线性方程组，设是其分解.我们先求解方程组.由于是下三角矩阵，则解向量可以通过依次求出其分量而求出，再求解方程组.解向量可以通过该方程组依次求出分量而快速得出.于是由两个方程组，的求解而给出的解. 当n阶矩阵满足定理4.1的条件时，我们可以用初等变换的方法求出和. 因为当时，由于可逆，故必存在可逆阵组 即.也就是说，我们可以先对施行行的初等变换得出上三角阵，而矩阵可以通过对单位阵进行相同的行初等变换得出.即有 于是，为保持为下三角（从而也是下三角阵），在进行行初等变换时，不能进行行的对换，上行的倍数应加到下行的对应元. 例1 试将矩阵 进行分解. 解 从而 ，， 因为 ，所以 . 定理4.2 在定理4.1的条件下，存在唯一的主对角元都是1的下三角阵L与上三角阵U以及唯一的对角阵D，使得. 证 我们首先证明这种形式的分解的存在性.由于 其中L是主对角元为1的下三角阵.作分块矩阵的乘法： 则.于是 即 时，有所以 L,就是主对角元全为1的下、上三角矩阵，D是对角阵。 假定，其中、是主对角元为1的下三角阵，，是主对角元为1的上三角阵，则 （4-1） 是A的两个LU分解。根据A的LU分解的惟一性知，。但由式（4-1）又得 （4-2） 式（4-2）左边的矩阵是下三角矩阵，且主对角线上的元素是的主对角线上的对应元，而右边的矩阵是上三角矩阵，其主对角线上的元素是的主对角线上的对应元。从而是对角阵。所以，。 令 ， ，，则由 ，得 = = 由于已有，故，均不为零，所以，有 即 。所以，。 推论4.1 假定矩阵A满足定理4.1的条件且 ，则存在主对角元为1的下三角阵L与实对角阵D,使得 称为A的Cholesky分解。 证 由定理4. 1 ,存在主对角元为1的下、上三角阵L，U及对角阵D，使 并且这样的表示是唯一的，根据 及，分别是主对角元为1的下、上三角阵，由分解的惟一性，可知 所以，D是实对角阵且 § 4.2 矩阵的QR分解 在第2章 §2.2定理2.2中，我们利用Gram-Schmidt方法将一组线性无关的向量 变为与之等价的一组单位正交向量.这一过程将给我们提供另一种矩阵的分解形式. 定理4.3 设,则必存在非奇异上三角矩阵R及矩阵Q，，使得 证 将A写成分块矩阵 则因,故线性无关，由第2章§2.2推论2.1知 令 则因Q的列向量组是单位正交组，故，又因，,故R的主对角线上的元素都不为0,从而是非奇异的上三角阵。 特别地，当时，是n阶酉矩阵。 我们将等式称为矩阵的分解。 设矩阵 利用Gram-Schmidt 方法对进行分解 解 令 即 所以 用Gram-Schmidt 方法对矩阵进行分解时，所论矩阵必须为列满秩矩阵。不是列满秩的矩阵只要是方阵，也可以做类似分解。这就是所谓的Householder 方法。 定义 4.2 设是单位向量，令 称为Householder 矩阵或Householder变换。 定理 4.4 设为阶Householder 矩阵，则 （1） （2） （3）以1为重特征值，为单特征值； （4）是阶Householder 矩阵，其中。 证 容易由定义证明（1）和（2），下证（3）（4）. 当n=1时，结论成立。假定。若为的一个特征值，为对应于的特征向