利用（3.1）得到课程.ppt

5.3.2 高斯-若当消去法 5.4 矩阵三角分解法 5.4.2 平方根法 5.4.3 追赶法 定理11 如果 为 阶对称正定矩阵，则存在一个实的非奇异下三角阵 当限定 的对角元素为正时，这种 可以用直接分解方法来确定计算 元素的递推公式. (对称正定矩阵的三角分解或Cholesky分解) 使 分解是唯一的. 因为 其中 于是得到解对称正定方程组 的平方根法计算公式： 对于 l. （4.7） 由矩阵乘法及 按等式两边对应元素相等，得 2. 求解 即求解两个三角形方程组： 4. 3. （4.8） 由计算公式1知 所以 于是 这个结果说明，分解过程中元素 的数量级不会增长 且对角元素 恒为正数. 当求出 的第 列元素时， 的第 行元素亦算出. 所以平方根法约需 次乘除法，大约为一般直接分解法计算量的一半. 于是不选主元素的平方根法是一个数值稳定的方法. 由于 为对称阵，因此在计算机实现时只需存储 的 下三角部分. 下三角部分共需存储 个元素，可按行主序用一维数组存放， 即 矩阵元素 在一维数组中表示为 的 元素存放在 的相应位置. 用平方根法解对称正定方程组时，需要用到开方运算. 为了避免开方，用定理10的分解式 即 由矩阵乘法， 并注意 得 于是得到计算 的元素及 的对角元素公式： l. （4.9） 对于 2. 为避免重复计算，引进 由(4.9)得到按行计算 元素的公式： 对于 1. 2. 3. 计算出 的第 行元素 后， 存放在 的第 行相应位置， 然后再计算 的第 行元素， 存放在 的第 行. 的对角元素存放在 的相应位置. 对称正定矩阵 按  分解和按 分解计算量差 不多，但 分解不需要开方计算. 例如 5. 4. （4.11） 计算公式(4.10)，(4.11)称为改进的平方根法. 求解 计算公式 实际问题中，通常要求解系数矩阵为对角占优的三对角线方程组 （4.12） 简记为 其中，当 利用(3.1)得到 则有 若记 考虑 时的 . （3.2） 为单位下三角阵，其元素的绝对值不超过1. 记 由(3.2)得到 其中 为排列矩阵， 为单位下三角阵， 为上三角阵. 其中 这说明对(2.1)应用列主元素消去法相当于对 先 进行一系列行交换后对 再应用高斯消去法. 定理8 如果 为非奇异矩阵， 则存在排列矩阵 使 其中 为单位下三角阵， 为上三角阵. 编程时， 元素存放在数组 的下三角部分， 元素 存放在 上三角部分，由记录主行的整型数组 可知 的情况. 而在实际计算中只能在计算过程中做行的交换. (列主元素的三角分解定理) 高斯消去法中，若同时消去对角线下方和上方的元素， 设用高斯-若当消去法已完成 步， 化为等价 方程组 ，其中 这种方法称为高斯-若当(Gauss-Jordan)消去法. 在第 步计算时 对上述矩阵第 行的上、下都进行消元. 1. 按列选主元素，即确定 使 2. 换行(当 时)交换 第 行与第 行元素. 3. 计算乘数 ( 可保存在存放 的单元中). 5. 计算主行 上述过程结束后有 4. 消元计算 用高斯-若当方法将 约化为单位矩阵，计算解就在 常数项位置得到，用不着回代求解，计算量大约需要 次乘除法，比高斯消去法大，但用高斯-若当方法求矩阵 的逆矩阵还是比较合适的. 定理9 设 为非奇异矩阵， 方程组 的增广矩阵为 . 如果对 应 用高斯-若当方法化为