方程组的数值解法 计算方法课件及实验 培训讲解.pptVIP

计 算 方 法;§6 方程组的数值解法;§6.0 引 言;§6.0 引 言;§6.1 高斯消元法;Step k：设 ，计算因子;回代;利用高斯消元法求解方程组: ;利用;利用;§6.2 选主元素的高斯消去法;例：单精度解方程组;全主元消去法;算法：1. 消元过程，对 (1) 选主元，找 使得 (2) 若 ，则停止，推出 (3) 若 ，则换行， (4) 消元，对 有 ;回代过程： (1)若 ，则停止 (2)对 ;例：;例：用列主元法解;第二列的后两个数中选出主元 2.5;§6.3 矩阵的三角分解;显然若令;Step n ? 1:; 由上述讨论可知，高斯消去法实质上产生了一个将系数 矩阵A分解为上三角阵与下三角阵相乘的因式分解。;矩阵的三角分解;用三角分解法解方程组;由Doolittle分解得：;依次计算，原方程组可表为：; 与解f (x)=0 的不动点迭代相类似，将AX=b改写为X=BX+f 的形式，建立雅可比方法的迭代格式： 其中，B称为迭代矩阵。其计算精度可控，特别适用于求解系数为大型稀疏矩阵(sparse matrices)的方程组。;问题: (a) 如何建立迭代格式？ (b) 向量序列{ x(k) }是否收敛以及收敛条件?;2 例题分析:;2 例题分析:;2 例题分析:;迭代结果如下表：; 设方程组 AX=b , 通过分离变量的过程建立Jacobi迭代公式，即 ;? 雅可比迭代法的矩阵表示 ;例1.;依此类推,得方程组满足精度的解为x12;;… … … …;例2.;x1 =2.5000 2.0909 1.2273 d =3.4825 x2=2.9773 2.0289 1.0041 d =0.5305 x3 =3.0098 1.9968 0.9959 d =0.0465 x4 =2.9998 1.9997 1.0002 d =0.0112 x5 =2.9998 2.0001 1.0001 d =3.9735e-004 x6 =3.0000 2.0000 1.0000 d =1.9555e-004 x7 =3.0000 2.0000 1.0000 d =1.1576e-005;Seidel迭代格式为 ; （迭代法基本定理） 设有方程组 ，对于任意的初始向 量 ，迭代公式 收敛的充要条件是迭 代矩阵 的谱半径 . ; （迭代法收敛的充分条件） 设方程组 迭代法为 如果有 的某种算子范数满足 ，则 ;关于解某些特殊方程组迭代法的收敛性 定义：(对角占优阵) 设 (1) 如果 元素满足 称 为严格对角占优阵 (2) 如果 元素满足 且上式至少有一个不等式严格成立， 称 为弱对角占优阵。;例：;所以;所以Gauss-Seidel迭代法发散; 6.7 向量及矩阵范数;定义1. ;--------(1);例1 求下列向量的各种常用范数;常用的矩阵范数;例2;求矩阵A的各种常用范数;求矩阵A的各种常用范数;容易计算;谱半径定义 ;§6.8 误差分析;例: ;在解中放大了 倍。;3 系数矩阵元素a ij的误差对解的影响;称cond (A)=;例题;§6.9 总结;迭代法 是一种逐次逼近方法，注意到在使用迭代法时， X（k+1）=BX(k）+f ，其迭代矩阵B和迭代向量f在计算过程中始终不变，迭代法具有循环的计算公式、方法简单。此外，应注意收敛性与收敛速度问题.收敛性是迭代法的前提，针对不同的问题，分析并采用适当的数值算法，如Jacobi方法、Guass-Seidel方法、SOR方法等。; 复习要求