小明在解方程组.doc

小明在解方程组 ?　　 ? 　　n? 　　*** 　　?xk=(bk?∑ak,j?xj)/ak,k,(k=n?1,L,2,1).?j=k+1? 　　(1)(1))(1)(1))?a11?a11a12La1(1|b1(1)?a12La1(1|b1(1)?nn?(1)(1)(1)(1)?a(1)≠0?(2)(2)(2)?11aaLa|b0aLa|b222n2?222n2?第一步消元：?21 → ? ?M?MMLM|M?MLM|M??(1)?(1)(1)(1)?(2)(2)(2)?aaLa|b0aLa|b???n1?n2nnn?n2nnn??? 　　其中 　　(1) 　　?(2)a1j(1)(1)(1)(1) 　　?aij=aij?ai1?(1)=aij?mi1?a1j,(i,j=2,3,L,n);?a11 　　?(1) 　　b?b(2)=b(1)?a(1)?1=b(1)?m?b(1),(i=2,3,L,n); 　　iii1ii11(1)?a11? 　　1)ai(1 　　和 mi1=(1) 　　a11 　　即　　A(1)x=b(1)与 A(2)x=b(2) 同解LL, 　　(1)(1))(1)?a11?a12La1(1|bn1?(2)(2)(2)?0aLa|b222n2?　　　　　　 （*） ???→??MMLM|M??(n)(n)?0Lann|bn???0? 　　(k)(k) 　　上述消去法能够顺利进行下去要求akk≠0,(k=1,2,L,n)其中akk称之为约化主元素 　　其中第k步消去公式为： 　　(k) 　　?(k+1)akj 　　(k)(k)(k)(k) 　　?aij=aij?aik?(k)=aij?mik?akj,(i,j=k+1,L,n);?akk 　　?(k) 　　?b(k+1)=b(k)?a(k)?bk=b(k)?m?b(k),(i=k+1,L,n).iiikiikk(k)?akk?(k)aik 　　其中mik=(k) 　　akk 　　(1)(1)(1) 　　定理1. 设 A=[aij]∈Rn×n,b=[b1(1),b2,L,bn]∈Rn. 若约化主元素 　　(k)akk≠0,(k=1,2,L,n),则可通过高斯（顺序）消去法将方程组Ax=b化为三角形方 　　程组（*） 　　(k) 　　注：akk≠0,(k=1,2,L,n) ? A 的所有顺序主子式不等于零 　　n3n 　　计算量（乘除法次数）：（消元和回代） +n2?=O(n3) 　　33 　　比较： 以10阶为例，运用克莱姆法则需×108次乘除法，高斯顺序消去法需 　　430次乘除法 　　可见高斯顺序消去法的计算量极少 　　(2)列主元消去法 　　(k) 　　① 当某个akk=0时，高斯消去法无法进行下去； (k)(k) ② 即使akk≠0，但|akk|很小；此时解也会失真 　　?+2x2=1 　　（用具有舍入的4位浮点数进行运算） 例 ? 　　2x+3x=212? 　　?x1= 　　； 精确解为 ? 　　x=?2?x=?x= 　　不选主元 ? ?1 ；选主元 ? ?1 　　?x2=?x2= 　　优点：提高数值计算稳定性，且|mik|≤1，增长因子为2 　　§2. （矩阵的）三角分解法 　　1．理论 　　?Ly=b 　　基本思想： Ax=b ? L?Ux=b ? ? 其中L,U皆为三角阵 　　Ux=y? 　　杜利特尔分解：A=L?U其中L为单位下三角阵，U为上三角阵； 克劳特分解： A=L?U其中L为下三角阵，U为单位上三角阵 　　(k) 　　存在性： 对方程组 Ax=b （记为A(1)x=b(1)），设akk≠0,(k=1,2,L,n)，则存在唯一的杜利特尔分解 　　引入 　　0L0??10L0 　　?0O?ML0???00L010L0? ，(k=1,2,L,n?1) 称之为初等下三角阵 　　Lk=? 　　?mk+1,k1L0?00L0?MOM?MM?? 　　? 　　?00L0?mn,k0L1??[A 　　(1) 　　|b]→[A(2)|b(2)],L, A(k+1)=LkA(k),b(k+1)=Lkb(k),(k=1,2,L,n?1)； 　　(1) 　　L1 　　?A(n)=Ln?1Ln?2LL2L1A(1)?(n)(1)?b=Ln?1Ln?2LL2L1b 　　?1?0??01 　　；而L?k=? 　　?0M??0? 　　0L00OM0L010L0mk+1,kMM0L0mn,k 　　LL0L1LO0L 　　0?0?? 　　0? 0?M??1?? 　　?1?11(n) 　　若记L=L1L2LL? ，则 A=A(1)=L?U n?1 ，U=A 　　~~~~~~ 　　?L?U=L?U?L?1?L=U?U?1,