第5章 线性代数方程组的直接解法.pptVIP

线性方程组的直接解法 第5章 解线性方程组的直接方法 又 由分解的惟一性, 得 代入(3.6)得 定理9 (对称阵的三角分解定理) 设 为 阶对称阵，且 的所有顺序主子式均不为零，则 可唯一分解为 其中 为对角阵， 为单位下三角阵. 若 为对称正定矩阵，则在分解式 中， 的对角元素 均为正数. 事实上，由 的对称正定性，有 于是 此时 其中 为下三角矩阵.  定理10 (对称正定矩阵的三角分解或Cholesky分解) 如果 为 阶对称正定矩阵，则存在一个实的非奇异下三角阵 使 当限定 的对角元素为正时，这种 分解是唯一的. 计算 元素的递推公式. 其中 又 由矩阵乘法 解对称正定方程组 的平方根法计算公式： 对于 l. （3.7） 2. 求解 即求解两个三角形方程组： 4. 3. （3.8） 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上 解 例9 解线性方程组 于是 这个结果说明，分解过程中元素 的数量级不会增长 且对角元素 恒为正数. 当求出 的第 列元素时， 的第 行元素亦算出. 所以平方根法约需 次乘除法。 数值稳定——不选主元素的平方根法 所以 为对称阵，在计算机实现时只需存储 的下三角部分. 由 计算量——不选主元素的平方根法 存储 下三角部分共需存储 个元素，可按行主序用一维数组存放， 即 矩阵元素 在一维数组中表示为 的 元素存放在 的相应位置. 平方根法的改进——避免开方运算 即 l. （3.9） 对于 2. 引进 按行计算 元素的公式： 1. 2. 3. （3.10） 对于 Step1 Step2 Step3 Stepn 元素的储存比消去法节省一半，但比平方根法多用一个单元内存， 因此透明度很差。 5. 4. （3.11） 改进的平方根法 求解 计算公式 1. 2. 3. （3.10） 例10 用LDLT解方程组 解 例10 用LDLT解方程组 解 系数矩阵为对角占优的三对角线方程组 （3.12） 简记为 其中，当 5.3.3 追赶法——三对角线方程组 LU分解会是什么样子？ 其中 为下三角矩阵， 为单位上三角矩阵. 由系数阵 的特点，可以将 分解为两个三角阵的乘积， 即 设 （3.13） 其中 为待定系数. 由矩阵乘法 （3.14） （3.15） 即 由 得 时由对角占优的条件，(3.15)是成立的. 现设(3.15)对 成立，求证对 亦成立. 由归纳法假设 又由(3.15)及 的对角占 优条件，有 U的特点 下面介绍高斯消去法的矩阵分析. 设(2.1)的系数矩阵 的各顺序主子式均不为零. （2.1） 5.2.2 矩阵的三角分解 施行行的初等变换相当于用初等矩阵左乘 第1步消元 其中 第2步消元 一般第 步消元， 化为 , 化为 ， 相当于 其中 重复这过程，最后 （2.14） 记上三角矩阵 为 ，由(2.14)得到 单位下三角矩阵 其中 定理7（矩阵的LU分解） 设 为 阶矩阵，如果 的顺序主子式 则 可分解为一个单位下三角矩阵 和一个上三角矩阵 的乘积，且这种分解是唯一的. 证明 为非奇异矩阵 存在性已证。 下证明唯一性 设 其中 为单位下三角矩阵， 为上三角矩阵. 由于 存在，故 唯一性得证. 故 例3 对于例2，系数矩阵 由高斯消去法， 故 解： 由高斯消去法知道，在消元过程中可能出现 即使主元素 但很小时，用其作除数，会导致其他元素数量级的严重增长和舍入误差的扩散，最后也使得计算解不可靠. 这时消去法将无法进行； ? 求矩阵