第5篇 解线性方程组的直接方法.ppt

§1 引言与预备知识 10) 初等置换阵 由单位矩阵I交换第i行与第j行（或交换第i列与第j列），得到的矩阵记为Iij,且 IijA=A（为交换A第i行与第j行得到的矩阵）； AIij=B（为交换A第i列与第j列得到的矩阵）。 11) 置换阵 由初等置换阵的乘积得到的矩阵. §2 高斯消去法 解 第1步.将方程（2.2）乘上-2加到方程（2.4）式中的未知数X1，得到 －4X2－X3=－11. (2.5) 第2步. 将方程（2.3）加到方程（2.5）上去，消去（2.5）中的未知数X2。得到与原方程组等价的三角形线性方程组 由此看出，用消去法解线性方程组的基本思想是用逐次消去 未知数的方法把原线性方程组AX=b化为与其等价的三角形线性 方程组，从而求解三角形线性方程组可用回代的方法求解。 下面我们讨论求解一般线性方程组的高斯消去法。 将方程组AX=b 记为A（1）X=b（1） .其中 A（1） =（aij (1)）=(aij), b(1)=b. §3 高斯主元素消去法 重复上述过程，设已完成第k-1步的选主元素，交换两行及消元计算， 约化为 (5)det←akk*det 3.如果ann=0，则计算停止（det（A）=0） 4.回代求解 (1)bn ←bn/ann (2)对于i=n-1,…,2,1 §4 矩阵的三角分解法 克劳特分解方法 设A为n×n阶非奇异矩阵,且各阶主子矩阵为非奇异,则矩 阵A的克劳特(Crout)分解为 A=LU 其中 这样,L、U中的元素都已求出。计算L的各列与U的各行的次序如图所示 。 对方程组Ax=b的系数矩阵A作出LU分解后,方程组便化为 LUx=b 则求解上列方程组就化为依次解方程组 Ly=b Ux=y 由于L为下三角矩阵,U为单位上三角矩阵,故上述方程组的求解极为方便。他们的计算公式分别为 用克劳特分解求解线性方程组Ax=b的计算过程为: ① LU分解过程:对于k=1,2,…,n依次计算 例4 用克劳特分解方法求解下列方程组 利用矩阵乘法可得到 这样原方程组就化为依次求下列两个三角形方程组 §5 向量和矩阵的范数 §6 误差分析 证明 用反证法。若det（I―B）=0,则（I―B）X=0有非零解， 即存在X0≠0使BX0=X0, 故‖B‖≧1，与假设矛盾，又由 (I－B)（I － B）-1 =I，有 （I － B）-1=I+B （I － B）-1 ， 从而 ‖（I ±B）-1 ‖≦‖I‖+‖B‖‖ （I ±B）-1 ‖ ‖（I －B）-1 ‖ ≦ 一、矩阵的条件数 考虑线性方程组 AX=b 系数矩阵A和右端b的小扰动所产生的相对误差. 例8 方程组 准确解为 常数项微小变化后 准确解 定义7 如果矩阵A或常数项b的微小变化,引起线性方程组AX=b的解的巨大变化,则称此方程组为病态方程组矩阵A称为病态矩阵,否则称方程组为良态方程组,矩阵A为良态矩阵. 条件数刻画了线性方程组AX=b的解对数据误差的灵敏程度，它只与此方程组的系数有关，反映了方程组固有的本性。故可用条件数来描述方程组的性态. 例9 求Hilbert矩阵H3的条件数. 注：一般判断矩阵是否病态，并不计算A?1，而由经验得出。 ? 行列式很大或很小（如某些行、列近似相关）； ? 元素间相差大数量级，且无规则； ? 主元消去过程中出现小主元； ? 特征值相差大数量级。 如何发现判断矩阵是病态的? 如何解决和处理? 预处理方法，即将AX=b转化为等价的方程组 例10 设 则 化为 则 二、高斯—若当消去法 算法(高斯—若当消元法). 例4 采用高斯—若当消去法求矩阵 的逆A-1 . 设有线性方程组：AX=b 一、直接三角分解法 1、不选主元三角分解算法 当A非奇异时，由不需选主元的顺序高斯消去法知 就有,不选主元的三角分解算法： 于是，可以通过求解两个三角形方程组 得到原方程