数值第四章_线性方程组的迭代解法导论.pptVIP

当ω＝1时，SOR法化为 G－S迭代法 G－S法为SOR法的特例, SOR法为G-S法的加速。 例1. 用G－S法和SOR法求下列方程组的解: 要求精度1e-6 解: (1)G-S迭代法 x1 x2 x3 1 1 1 0.7500000 0.3750000 1.5000000 0.5625000 0.5312500 1.5416667 0.6510417 0.5963542 1.6145833 0.7018229 0.6582031 1.6727431 ………………………………………. 0.9999933 0.9999923 1.9999926 0.9999943 0.9999935 1.9999937 0.9999952 0.9999944 1.9999946 k = 71 x= 0.999995 0.999994 1.999995 满足精度的解 迭代次数为71次 (2)SOR迭代法 x1 x2 x3 1 1 1 0.6375000 0.0121875 1.3199063 0.2004270 0.3717572 1.3122805 0.6550335 0.5340119 1.6922848 0.7058468 0.7733401 1.7771932 ……………………………………….. 0.9999990 0.9999976 1.9999991 0.9999984 0.9999993 1.9999989 0.9999998 0.9999994 1.9999998 0.9999996 0.9999998 1.9999997 k = 24 x= 1.000000 1.000000 2.000000 满足精度的解 迭代次数为24次 选取适当的ω，SOR法的收敛速度比G-S法要快得多。 4.6 误差分析 结论：（1）的常数项b的第二个分量只有1/1000的微小变化，方程组的解变化却很大。 例 记方程组（1） 为Ax＝b，其精确解为：x1*=2，x2*=0 现考察方程组（2） 可将其表示为：A(x+?x)=b+?b 其中 ?b= (0,0.0001)T 设x为(1)的解，显然(2)的解为：x+?x= (1, 1)T 设Ax＝b的扰动方程组为(A+ ?A)(x+?x)=b+?b，其中?A叫A的扰动矩阵， ?x和?b叫x和b的扰动向量。 定义 若矩阵A或常数项b的微小变化引起方程组Ax＝b的解的巨大变化，则称此方程组为病态方程组，A为病态矩阵（相对方程组而言）；否则称方程组为良态方程组，A为良态矩阵。 研究方程组中A或b的微小误差对解的影响的分析称“扰动分析”。 (1) ?A=0，则A (x+?x)= b+?b，减去Ax=b，得 设Ax=b的扰动方程组为(A+?A) (x+?x)= b+?b，下面进行扰动分析： (2) ?b=0，则(A+?A) (x+?x)= b，同理可得 A ?x= ?b， 故 ?x= A-1 ?b， 即||?x||? ||A-1 || || ?b || ， 又由Ax=b，有||b||? ||A || || x || ， 所以 定义 设A非奇异，称数 cond(A)= ||A-1 || ||A||为矩阵A的条件数。 说明： （1）条件数小，扰动引起的解的相对误差一定小； 条件数大，扰动引起的解的误差可能很大。（条件数与所取的范数有关，最常用的是||·||∞和||·||2） （2）由于cond(A)= ||A|| ||A-1 ||≥||A A-1 ||=||I||=1， 故条件数是一放大的倍数，且总以1为下界。 （3）当A为正交矩阵（A-1=AT）时，有cond(A)2=1，所以正交矩阵的方程组是良态的。 可得x= (27.0000 -192.0000 210.0000)T x+?x= (30.0000 -210.0000 228.0000)T 例：已知Hilbert矩阵 cond(H3)=748 cond(H6)=2.9×107 cond(H7)=9.85×108 计算可得 n越大，Hn的条件数越大，故Hilbert矩阵是一个典型的病态矩阵。 如 定理 （事后误差估计）设A非奇异，x*和x分别是Ax=b的精确解和近似解，r=b-Ax为残差向