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

数值分析数值分析第讲矩阵的三角分解法第1页，共30页，星期日，2025年，2月5日第四讲

矩阵的三角分解法第二章线性方程组的解法第2页，共30页，星期日，2025年，2月5日InScientificComputing↓LargeLinearSystemsAx=bassub-problems/asintermediatestepsGauss-SeidelmethodJacobimethodSORmethodConjugateGradientmethodforsymmetricsystemsGaussianeliminationLUfactorizationCholeskyfactorizationGMRESGCRBi-CGCGSBi-CGSTABBi-CGSTAB2GPBi-CGBi-CGSTAB(L)第3页，共30页，星期日，2025年，2月5日回顾高斯消去法Gauss消元法的第k步：从矩阵理论来看，相当于左乘矩阵指标为k的初等下三角阵第4页，共30页，星期日，2025年，2月5日因此，整个Gauss消元法相当于左乘了一个单位下三角阵所以有L为单位下三角阵，U为上三角阵第5页，共30页，星期日，2025年，2月5日矩阵的直接三角分解(LU分解)对解线性方程组有什么帮助？三角方程组，易于求解第6页，共30页，星期日，2025年，2月5日1、Doolittle分解L为单位下三角，U为上三角比较第1行：比较第1列：第7页，共30页，星期日，2025年，2月5日比较第2行：比较第2列：第8页，共30页，星期日，2025年，2月5日比较第k行：比较第k列：K-1次K-1＋1次第9页，共30页，星期日，2025年，2月5日分解过程完毕，加上两次反代过程总运算量为：第10页，共30页，星期日，2025年，2月5日存储在矩阵的原来位置，且不影响计算第11页，共30页，星期日，2025年，2月5日第12页，共30页，星期日，2025年，2月5日L为下三角，U为单位上三角第13页，共30页，星期日，2025年，2月5日两次反代过程比较第k列：比较第k行：第14页，共30页，星期日，2025年，2月5日第15页，共30页，星期日，2025年，2月5日定理2.4若矩阵A非奇异，则存在置换矩阵Q，使得QA可以作Doolittle分解QA=LU其中L是单位下三角矩阵，U是上三角矩阵。第16页，共30页，星期日，2025年，2月5日选主元直接三角分解法当需选主元时，PA=LU.设第r-1步已完成,就有第17页，共30页，星期日，2025年，2月5日第18页，共30页，星期日，2025年，2月5日选主元三角分解算法：第19页，共30页，星期日，2025年，2月5日求解Ly=Pb及Ux=y的算法：第20页，共30页，星期日，2025年，2月5日第21页，共30页，星期日，2025年，2月5日数值分析