第五讲：线性方程组的计算.pptVIP

线性方程组求解 线性方程组求解 §5.3 线性方程组的其他数值方法 1、 LU分解法 令 线性方程组求解 记为 代入方程组 则有 令 由 解得 解得 再由 例 线性方程组求解 令 得 由 得 得 由 得 由 由 得 线性方程组求解 由 得 再由 得 线性方程组求解 算法 Step1: Step2: Step3: 线性方程组求解 Step4: Step5: 线性方程组求解 2、 追赶法 · · · · · 考查 线性方程组求解 记为 令 线性方程组求解 参照 追赶过程 线性方程组求解 令 则有 线性方程组求解 再令 则有 线性方程组求解 例 令 线性方程组求解 由 得 得 由 由 得 直接得 线性方程组求解 即 由 得 线性方程组求解 由 得 线性方程组求解 线性方程组求解 §5.4 线性方程组迭代解法 1、Jacobi迭代法 考虑 令 其中 算法 则有 迭代格式 如 迭代格式 线性方程组求解 取 得 线性方程组求解 如 迭代格式 Gauss-Seidel迭代格式 线性方程组求解 x3 x2 x1 x3 x2 x1 K 1 1 1 1.0001 1.0002 0.9959 1.0041 1.2273 0 2 2 2 2.0001 1.9997 1.9968 2.0289 2.0909 0 3 3 3 2.9998 2.9998 3.0098 2.9773 2.5 0 Gauss-Seidel 1.0031 2.0026 2.9990 7 0.9998 2.0006 3.0002 8 1.0038 2.0000 2.9938 6 0.9205 1.9478 3.0241 4 1.0010 1.9840 3.0003 5 0.9716 2.0455 3.1364 3 1 2.3636 2.8750 2 3 3 2.5 1 0 0 0 0 Jacobi 线性方程组求解 x1(1)=0;x2(1)=0;x3(1)=0; for k=1:8 x1(k+1)=(20+3*x2(k)-2*x3(k))/8; x2(k+1)=(33-4*x1(k+1)+x3(k))/11; x3(k+1)=(36-6*x1(k+1)-3*x2(k+1))/12; end x1(1)=0;x2(1)=0;x3(1)=0; for k=1:8 x1(k+1)=(20+3*x2(k)-2*x3(k))/8; x2(k+1)=(33-4*x1(k)+x3(k))/11; x3(k+1)=(36-6*x1(k)-3*x2(k))/12; end Jacobi Gauss-Seidel 线性方程组求解 §5.5 迭代法的收敛性 如 定义1 记为 收敛于 则称 成立， 如果 及矩阵 设矩阵序列 线性方程组求解 迭代法的基本定理 如 860 645 6 320 285 5 20 15 2 50 45 3 140 105 4 5 5 1 0 x2 0 x1 0 k 由 得 设有方程组 对于 及 迭代格式 收敛 其中， 线性方程组求解 又如 得 -3.9815 -2.9861 6 -3.9583 -2.963 5 -3.3333 -2.5 2 -3.75 -2.7778 3 -3.8889 -2.9167 4 -2.5 -1.6667 1 0 x2 0 x1 0 k 由 线性方程组求解 其中 定理：如果方程组系数矩阵严格对角占优或不可约弱对角占优 ，则对任意初值，Jacobi迭代法，Gauss-Seidel迭代法均收敛。 1、严格对角占优 3、不可约 2、弱对角占优 如果存在置换矩阵 使得 称矩阵A可约，否则称矩阵A不可约。 线性方程组求解 如 -3.9815 -2.9861 6 -3.9583 -2.963 5 -3.3333 -2.5 2 -3.75 -2.7778 3 -3.8889 -2.9167 4 -2.5 -1.6667 1 0 x2 0 x1 0 k 860 645 6 320 285 5 20 15 2 50 45 3 140 105 4 5 5 1 0 x2 0 x1 0 k 线性方程组求解 * 第五讲 方程组求解 §5.1 Gauss消去法 线性方程组求解 记为 1、 Gauss消去过程 想法 线性方程组求解 线性方程组求解 0 . . . 0 算法 计算格式 线性方程组求解 回代 线性方程组求解 Gauss消去算法 消元 回代 线性方程组求解 §5.2 Gauss主元消去法 1、 问题 例：单精度解方程组 精确解 Gauss消去法 小主元 导致计算失败。 线性方程组求解 0 0 Gauss消去法 线性方程组求解 2、 全主元消去法 Step 1: 定义一个排序数组 Step 2: 在A(