第四章线性方程组及求解.pptVIP

②对S-格式： 4. 充要条件.设系数矩阵A为对称矩阵，且 5.充分条件： NOTES:① ②改变方程组中方程次序不会改变方程组的解，但有时会改变格式 的收敛性。 ③SOR格式收敛的必要条件: 0ω2 四、误差估计 五、迭代法的收敛速度： 例题 解得x1=2，x2=0。 扰动b 解得x1=1，x2=1。 值得指出： ①A、b的扰动是一种客观存在； ②A、b的扰动对Ax=b的解有时会产生破坏作用。 此即为解对系数的敏感性，显然，解对A 、b的敏感性越大，对解的破坏性越大。反之亦然。 定义：对给定的线性方程组Ax=b，若A、b的微小扰动对Ax=b的解产生 很大的变化，称为病态方 程组，A为病态矩阵；反之，称Ax=b为良性方程组，A称为良态矩阵。 二、病态性分析 给定Ax=b（|A| ≠ 0），则Ax*=b，x*唯一存在。 分三种情形： 1.b扰动，A不扰动 2.A扰动，b不扰动 3.A以及b同时扰动 1.b扰动，A不扰动 b→b+δb ，则此时 x* →x*+δx A（x*+δx）=b+δb ∴ Aδx=δb ∴ δx=A-1δb ∴‖δx‖ ≤ ‖A-1‖·‖δb‖ 又‖b‖ =‖Ax*‖ ≤ ‖A‖·‖x*‖ 1/ ‖x*‖ ≤ ‖A‖/‖b‖ ∴ ≤ ‖A-1‖‖A‖‖δb‖·1/‖b‖ 易知‖A-1‖‖A‖越大,对解得变化影响越大。反之亦然。 2.A扰动，b不扰动 A→A+δA ，则此时 x*→x*+δx （A+δA）（x*+δx）=b 展开得 （A+δA）·δx=-δAx* 又 A+δA=A（E+A-1δA） 另设||A-1δA||≤||A-1||·||δA||＜1 有E+A-1δA可逆， 则‖(E+A-1δA)-1‖≤1/(1-‖A-1δA‖ ) ≤ 1/(1-‖A-1‖‖δA‖) 进一步有： δx =-(A+δA)-1δA·x*=-(E+A-1δA) -1A-1δAx* ‖δx‖ ≤ ‖(E+A-1δA)-1 ‖‖A-1 ‖‖δA ‖‖x*‖ ≤ 1/(1-‖A-1‖‖δA‖)· ‖A-1 ‖ ‖δA ‖‖x*‖ ∴ ‖δx‖/ ‖x*‖≤‖A-1‖‖δA ‖/(1-‖A-1‖‖δA‖) ≤ {‖A-1‖‖A‖ ‖δA ‖/‖A‖}/{1-‖A-1‖‖A‖‖δA‖/‖A‖)} 易知 ‖A-1‖‖A‖越大，对系数的扰动起的扩张作用越大，从而解的变化越大。反之亦然。 3.A以及b同时扰动 b→b+δb，A→A+δA ，x*→x*+δx 则（A+δA）（x*+δx）=b+δb 同2.知得： 同样可看出‖A-1‖‖A‖对解的影响一样的。 由此引出下面的条件数的概念 条件数 定义：对任意的n阶可逆矩阵A，称 Cond(A)= ‖A-1‖‖A‖ 为A的条件数。 易知，Cond(A)越大，Ax=b病态程度越大。 ① 常用的条件数 Cond(A)∞= ‖A-1‖∞ ‖A‖∞ Cond(A)2= ‖A-1‖2‖A‖2 ②条件数性质 1. Cond(A) ≥1 Cond(A)= ‖A-1‖‖ A‖ ≥ ‖A-1A‖=1 2.对常数c ≠ 0 ， Cond(cA)= Cond(A) 3.若R为正交矩阵，则Cond(R)2=1 4.对任意的n阶可逆阵A及正交阵R，有 Cond(A)2= Cond(RA)2 = Cond(AR)2 5.对任意的n阶矩阵A、B有 Cond(AB) ≤ Cond(A) · Cond(B) 三、病态性诊断及处理 已知对任意的n阶矩阵A,Cond(A)= ‖A-1‖ ‖A‖ 但实际上的计算很困难 例： 书中P170，希尔伯特矩阵Hn Cond(H3) ∞=748 诊断及处理方法 诊断 1.出现小主元 2.系数阵的行列式很小或某些行近似相关 3.系