数值代数精品课件-第3章-线性代数方程组的解法.pptVIP

第三章 线性代数方程组的解法 第一节 引言 举例（一） 一、Gauss 消去法计算过程 讨论第K次消元，得到消元公式 第K次消元目的：aik(k)=0 (i=k+1….n), 设akk(k) ≠ 0 ,为使aik(k)=0 (i=k+1….n）选取适当因子M使 aik(k)- M akk(k) =0 可求出 M = aik(k)/ akk(k) 第i个方程其他系数的变化为 （ Ri - M.Rk ) aij(k+1) = aij(k) – M. akj(k) ( i= k+1….n , j=k+1…..n+1) 所以，第K次消元时（后）（K=1….n-1) 消元因子： M = aik(k)/ akk(k) ( i= k+1,n ） 系数变化： ① aij(k+1) = aij(k) （i≤ k) ② aij(k+1) = aij(k) - M. akj(k) (ik , j=k+1…n+1) 最后得到： [A(n) | b(n)] (n-1 次消元后）其增广矩阵为： 回代： Xn= ann+1(n)/ann(n) 编程时， 为节省内存将Xn放在a nn+1(n) 同理X i 放在 a in+1(i) 第i次回代公式（i=n-1…..1) Xi(即a in+1(i))=(a in+1(i) - ）/aii(i) 二、算法描述 1、消元 对k=1…..n-1 消元因子： C= aik(k)/ akk(k) ( i= k+1….n ） 系数变化： ① aij(k+1) = aij(k) （i≤ k) ② aij(k+1) = aij(k) - C. akj(k) ( i k , j=k+1,…,n+1 ) 2、回代 第 i次回代公式 （ i=n,n-1….1) Xi(即a in+1(i))=(a in+1(i) - ）/aii(i) 三、特点 1、优点： 公式简明，容易程序化 2、缺点： 第k次消元时， 必须 a kk ≠ 0 , 且当 a kk ?0 时 ，误差很大， 数值不稳定 (p35 N-S图） 第三节 GAUSS列主元消去法 一、高斯消去法的缺点 在简单的高斯消去法中，如果遇到 a（k）kk =0,则消去过程就会中断，如果a（k）kk ≠ 0 ，但其绝对值很小时，在高斯消去法中，因为M = aik(k)/ akk(k) ，所以不是因为M数值过大，就是舍入误差过大，与实际的解相差很远。 例如：求解下列方程组 此方程组的准确解为：x1=10.00；x2=1.00 下面采用高斯消去法解方程组（取四位有效数字）， 消元后得-1043x2=-1044，则x2=1.001. 将x2代入第一个方程中，得 x1=-9.713 显而易见，利用高斯消去法得到的结果与精确解相差太悬殊了 从求解过程可以看出，在第一种解法中，乘以的系数为m=5.291/0.0030 如果将两个方程互换，再采用高斯消去法解方程组 消元后得59.14x2=59.14，则x2=1.000 将x2代入第一个方程中，得 x1=10.000 此时，利用高斯消去法得到的结果与精确解一致。 从求解过程可以看出，在第一种解法中，乘以的系数为m=5.291/0.0030，第二种解法中，乘以的系数为m=0.0030/5.291。我们希望m尽量小，希望主元尽可能大，就有了列主元消去法（按列）。 二、列主元消去法 1、选主元再消元 在每一步消元之前，即第k次消元时 ，在第k列上选取绝对值最大的元素为主元素 ai0k(最大值在i0行上） 若 | a i0k | εp 中断 否则 若i0 ≠ k 将i0 和 k 行互换 再按高斯消去法消元 2、回代 不改变方程组的解，同时又有效地克服了Gauss消元的缺陷 三、算法描述 1、选主元 （ 对 k=1…..n-1 ) (1) p(主元素） ?0 (2) 对 i=k……n 做 若 | aik | p 则 p ? | aik | i0?i 2、换行 a kj -- ai0j ( j= k……n) 3、消元： 同Gausss 4、回