数值分析第讲矩阵的三角分解法课件.ppt

《数值分析第讲矩阵的三角分解法课件.ppt》由会员分享，可在线阅读，更多相关《数值分析第讲矩阵的三角分解法课件.ppt（30页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数值分析第讲矩阵的三角分解法第1页，此课件共30页哦第四讲矩阵的三角分解法第二章 线性方程组的解法第2页，此课件共30页哦In Scientific ComputingLarge Linear SystemsAx=bas sub-problems/as intermediate stepsGauss-Seidel methodJacobi methodSOR methodConjugate Gradient methodfor symmetric systemsGaussian eliminationLU factorizationCholesky factorizationGMRESGCRB

2、i-CGCGSBi-CGSTABBi-CGSTAB2GPBi-CGBi-CGSTAB(L)第3页，此课件共30页哦回顾高斯消去法回顾高斯消去法Gauss消元法的第k步：从矩阵理论来看，相当于左乘矩阵指标为k的初等下三角阵第4页，此课件共30页哦因此，整个Gauss消元法相当于左乘了一个单位下三角阵所以有L为单位下三角阵，U为上三角阵第5页，此课件共30页哦矩阵的直接三角分解(LU分解)对解线性方程组有什么帮助？三角方程组，易于求解第6页，此课件共30页哦1、Doolittle分解分解L为单位下三角，U为上三角比较第1行：比较第1列：第7页，此课件共30页哦比较第2行：比较第2列：第8页，此课件

3、共30页哦比较第k行：比较第k列：K-1次K-11次第9页，此课件共30页哦分解过程完毕，加上两次反代过程总运算量为：第10页，此课件共30页哦存储在矩阵的原来位置，且不影响计算第11页，此课件共30页哦第12页，此课件共30页哦L为下三角，U为单位上三角第13页，此课件共30页哦两次反代过程比较第k列：比较第k行：第14页，此课件共30页哦第15页，此课件共30页哦定理定理2.4 若矩阵若矩阵A非奇异，则存在置换矩阵非奇异，则存在置换矩阵Q，使得使得QA可以作可以作Doolittle分解分解 QA=LU其中其中L是单位下三角矩阵，是单位下三角矩阵，U是上三角矩阵。是上三角矩阵。第16页，此课

4、件共30页哦选主元选主元直接三角分解法直接三角分解法当需选主元时，PA=LU.设第r-1步已完成,就有第17页，此课件共30页哦第18页，此课件共30页哦选主元三角分解算法选主元三角分解算法：第19页，此课件共30页哦求解求解Ly=PbLy=Pb及及Ux=yUx=y的算法的算法：第20页，此课件共30页哦第21页，此课件共30页哦第22页，此课件共30页哦三对角阵的追赶法三对角阵的追赶法（A的前的前n-1个顺序主子式非零）个顺序主子式非零）在数值求解常微分方程边值问题、热传导方程和建立三次样条函数时，都会要解三对角方程组：AX=bAX=b第23页，此课件共30页哦所以，有计算过程如下：说明说明

5、:稳定性（对角占优）；运算量5n-4次乘除法；存贮（一维数组）.第24页，此课件共30页哦平方根法平方根法应用有限元法解结构力学问题时，最后归结为求解线性代数方程组，系数矩阵往往对称正定。平方根法是一种对称正定矩阵的三角分解法，广泛用于求解系数矩阵为对称正定的线性代数方程组。设A为对称矩阵,且顺序主子式不为零，则第25页，此课件共30页哦若A为对称正定矩阵,则第26页，此课件共30页哦第27页，此课件共30页哦解AX=b的平方根法平方根法:称为平方根法，因为带了开方运算，因此不常用第28页，此课件共30页哦Cramer ruleGauss eliminationLU factorization

6、Gauss-Jordan eliminationSquare root/improved square root追赶法(n+1)!n3/3n3/3n3/2n3/65n-4Row/column/complete PivotingOnly eliminating elements in a column below the diagonal oneNo pivotingDirectly factorization column pivotingEliminating elements in a row except for the diagonal elementA symmetric and p

7、ositive definiteNo pivotingA三对角，弱对角占优Summary第29页，此课件共30页哦直接法内容总结直接法内容总结1 1：用：用：用：用LULU分解（分解（分解（分解（DoolittleDoolittle）求解线性代数方程组等价于顺序）求解线性代数方程组等价于顺序）求解线性代数方程组等价于顺序）求解线性代数方程组等价于顺序Gauss Gauss 消消消消元法：元法：元法：元法：AxAxb b LUx LUxb b Ux Uxy y，LyLyb b；2 2：用选主元：用选主元：用选主元：用选主元LULU分解（分解（分解（分解（DoolittleDoolittle）求解

8、线性代数方程组等价于选列主元）求解线性代数方程组等价于选列主元）求解线性代数方程组等价于选列主元）求解线性代数方程组等价于选列主元GaussGauss消元法：消元法：消元法：消元法：AxAxb b LUx LUxQb Qb Ux Uxy y，LyLyQbQb；3 3：对上半带宽为：对上半带宽为：对上半带宽为：对上半带宽为s s、下半带宽为、下半带宽为、下半带宽为、下半带宽为r r的带状矩阵作的带状矩阵作的带状矩阵作的带状矩阵作LULU分解，那么分解，那么分解，那么分解，那么L L为下半带为下半带为下半带为下半带宽为宽为宽为宽为r r的下三角矩阵，的下三角矩阵，的下三角矩阵，的下三角矩阵，U U为上半带宽为为上半带宽为为上半带宽为为上半带宽为s s的上三角矩阵的上三角矩阵的上三角矩阵的上三角矩阵；第30页，此课件共30页哦