第五课线性方程组的直接解法.ppt

 在自然科学和工程技术中,很多问题归结为解线性方程组.有的问题的数学模型中虽不直接表现为含线性方程组,但它的数值解法中将问题“离散化”或“线性化”为线性方程组.因此线性方程组的求解是数值分析课程中最基本的内容之一. 2换行,交换增广阵第k行和第ik行 (j=k,k+1,?,n), 结束 * 3计算乘数 mik= - aik／akk(i=1,2,?,n, i≠k) mkk=1／akk. (5.14) 4消元 aij=aij+mikakj(i=1,2,?,n,且i≠k,j=k+1,?,n) bi=bi+mikbk (i=1,2,?,n,且i≠k) (5.15) 5主行计算 akj=akj×mkk (j=k,k+1,?,n) bk=bk×mkk (5.16) 当k=n时, 显然xi=bi,i=1,2,?,n 就是 Ax=b的解. 结束 * 高斯-若当消元法的消元过程比高斯消元法略复杂,但省去了回代过程,它的计算量约为n3/2,大于高斯消元法.也称为无回代的高斯消元法. 5.4.2 求方阵的逆 高斯-若当消元法解方程组并不比高斯消元法优越,但用于矩阵求逆是适宜的,实际上它是初等变换方法求逆的一种规范化算法: 例3 结束 * 结束 * 所以 将高斯-若当消元法算法略加改造便成了求逆的算法. AX=I 对k=1,2,?,n 1按列选主元,确定ik 2换行,交换第k行和ik行 (j=k,k+1,?,n) (j=1,2,?,n) 3计算乘数 mik=-aik／akk , (i=1,2,?, i≠k) mkk=1／akk, (5.16) 结束 * 结束 5主行计算 akj=akj×mkk (j=k,k+1,…,n) xkj=xkj×mkk (j=1,2,…,n) (5.19) A-1=(xij)n×n 应该注意到,在上述算法中都加进了选列主元的步骤,它可以保证在方阵A非奇异的情况下都能求得A-1.求逆算法中不宜使用全面选主元,否则会改变逆阵元素的排列,增加 程序设计的复杂性. 4消元 aij=aij+mik akj,(i=1,2,…,n,且i≠k,j=k+1,…,n) xij=xij+mik xkj,(i=1,2,…,n,且i≠k,j=1,2,…,n) (5.18) * Ax=b是线性方程组,A是n×n方阵,并设A的各阶顺序主 子式不为零。令 A(1)=A,当高斯消元法进行第一步后,相当于 用一个初等矩阵左乘A(1).不难看出,这个初等矩阵为 其中li1 (i=2,3,…,n) 由式(5.9)确定,即 结束 * §5.5 矩阵的 LU分解 同样第k步消元有 进行n-1步后,得到A(n) ,记U=A(n),显然U的下三角部分全化为零元素,它是一个上三角阵。整个消元过程可表达如下: 结束 已知U是上三角矩阵,下面讨论L的性态,首先指出 * 其次指出 结束 证明见习题。 * L相当是由各L-1k (k=1,2,…,n-1)的所有左下元素拼凑后加上对角元1而得. 当A进行LU分解后,Ax=b就容易解了. 即Ax=b等价于: 结束 这样,Ax=b分解为解两个三角形方程组,三角形方程组是极易求解的.(5.24)即是高斯消元法的矩阵表述. * 定理5.2 矩阵 An×n,只要A的各阶顺序主子式非零,则A可以分解为一个单位下三角阵L和一个上三角阵U的乘积,即A=LU,且这种分解是唯一的. L是下三角矩阵,且所有的对角元为1,故称为单位下三角阵.这样,A=LU称为A的LU分解,其中L是单位下三角阵,U是上三角阵. 5.5.2 直接LU分解 A的LU分解可以用高斯消元法完成,但也可以用矩阵乘法原理推出另一种方法,结果是完全一致的.设: 由矩阵乘法公式可推出直接LU分解算法,也称杜利特(Doolittl