第五章直接法_数值积分.pptVIP

第五章. 解线性代数方程组的直接方法 §1.引言; ;（一般有限步内得不到精确解）;§2. Gauss (高斯)消去法;相当于第i个方程-第一个方程×数→新的第i方程—同解！第一方程不动！ ; 在［A(1):b(1)］中,红方框中的元素是要化为0的部分；［A(2):b(2)］中,红方框中的元素全部已发生变化,故上标由(1)改(2).; (i=2,3,?,n); 本次消元的目的是对框内部分作类似第一次消元的处理,消掉第k+1个方程到第n个方程中的xk项,即把 到  化为零.计算公式如下:; 它的方阵部分A(n)是一个上三角形矩阵,它对应的方程组是一个上三角形方程组,只要 ≠0,就可以回代求解公式为:;(i=k+1,k+2,?,n);;系数矩阵与常数项：;返回;(即为P169 定理6);;消去第一列的 n-1 个系数要计算n*(n-1） 个乘法。;§3. 高斯主元素消去法; 若 =0,则必有方框内的元素全为零,此时易证; 全主元消元法 在列主元消去法中,若每次选主元不局限所在列中,而在整个主子矩阵; 值得注意的是,在全主元的消元法中,由于进行了列交换, x各分量的顺序已被打乱.因此必须在每次列交换的同时,让 机器“记住”作了一次怎样的交换,在回代解出后将x各分 量换回原来相应的位置,这样增加了程序设计的复杂性. 此外作①步比较大小时,全主元消元法将耗用更多的机时. 但全主元消元法比列主元消元法数值稳定性更好一些.实际 应用中,这两种选主元技术都在使用.; 高斯-若当(Gauss-Jordan)消元法 很多问题都需要计算方阵的逆阵.线性代数中有多种求逆方法,本节讨论求A-1的数值方法.A非奇异,求A-1;在第k步计算时,考虑将第k行上下的第k列元素都化为零,且akk化为1.对k=1,2,?,n 1按列选主元,确定ik使;3.计算乘数 mik= - aik／akk… mkk=1／akk. ; 高斯-若当消元法的消元过程比高斯消元法略复杂,但省去了回代过程,它的计算量约为n3/2,大于高斯消元法.也称为无回代的高斯消元法.;; 所以; 5主行计算 akj = akj × mkk ( j=k,k+1,…,n) xkj = xkj × mkk ( j=1,2,…,n ) A-1=(xij)n×n;; Ax=b是线性方程组,A是n×n方阵,并设A的各阶顺序主子式不为零。 令A(1)=A,当高斯消元法进行第一步后,相当于用一个初等矩阵左乘A(1).不难看出,这个初等矩阵为;同样第 k 步消元有;其次指出;L相当是由各L-1k (k=1,2,…,n-1)的所有左下元素拼凑 后加上对角元1而得.;例1;选主元素的矩阵表示; 直接LU分解法 A的LU分解可以用高斯消元法完成,但也可以用矩阵乘法原理推出另一种方法,结果是完全一致的．;例如：;2.解 L y = b;杜利特算法可用文字描述如下: 将A的元素划分为形如“┏”的n框, 如A的第一行和第一列元素为第一框,紧靠第一框内部的第二行和第二列元素为第二框,依此类推,ann一个元素为第n框.为便于描述我们称原A矩阵中元素为a元素,计算后每框中的每行元素为u元素(包含对角元素),每列元素为l元素(不包含对角元素).;也可以用增广矩阵[A:b]进行以上操作,这时每框要多算一个 u元素,结果矩阵中的最后一列就是 Ly=b 的解,回代时只要解 Ux=y 就行了.;;方阵行列式求法 在实际问题中,有时会遇到求方阵的行列式.在线性代数中讲到的行列式定义算法和展开算法均不适用于阶数较高的行列式.而LU分解是计算行列式的十分方便和适用的算法. A=LU;克劳特(Crout)分解 将LU分解换一个提法:要求L为一般下三角阵,U为单位上 三角阵,就是克劳特分解.; 平方根法 实际问题中Ax=b,A若是对称正定矩阵,则高斯消元法简化为 平方根法或改进的平方根法.;定理１１ 设 A 对称正定,则存在三角分解 A=LLT,L是非 奇异下三角矩阵,且当限定 L 的对角元为正时, 这种分解是唯一的.;形如 A=LLT 的分解称为对称正定矩阵的乔累斯基分解。;对称正定矩阵的A=LLT分解对应于解对称正定方程组Ax=b 的平方根法.