数值计算方法（第3章） ppt课件.ppt

《数值计算方法（第3章） ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、3.2矩阵的三角分解法我们知道对矩阵进行一次初等变换，就相当于用相应的初等矩阵去左乘原来的矩阵。因此我们这个观点来考察Gauss消元法用矩阵乘法来表示，即可得到求解线性方程组的另一种直接法：矩阵的三角分解。3.2.1Gauss消元法的矩阵形式3.2.2Doolittle分解Doolittle分解定理3.2.1.若矩阵A有分解：A=LU，其中L为单位下三角阵，U为上三角阵，则称该分解为Doolittle分解，可以证明，当A的各阶顺序主子式均不为零时，Doolittle分解可以实现并且唯一。Doolittle分解定理3.2.1.的证明分两步:1)构造L,U;2)唯一性的证明。A的各阶顺序

2、主子式均不为零，即Doolittle分解Doolittle分解Doolittle分解Doolittle分解Doolittle分解Doolittle分解例题例题例题例题例题Doolittle分解Crout分解定理3.2.2.若矩阵A是非奇异的,则存在置换矩阵P,使得:PA=LU，其中L为单位下三角阵，U为上三角阵(对角元为1)，则称该分解为Crout分解。证明(略)3.2.3对称矩阵的Cholesky分解在应用数学中，线性方程组大多数的系数矩阵为对称正定这一性质，因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法，且分解过程无需选主元，有良好的数值稳定性。对称矩阵的Cho

3、lesky分解A对称：AT=AA正定：A的各阶顺序主子式均大于零。即对称矩阵的Cholesky分解由Doolittle分解，A有唯一分解对称矩阵的Cholesky分解定理3.2.4设A为对称正定矩阵，则存在唯一分解A=LDLT，其中L为单位下三角阵，D=diag(d1,d2,…,dn)且di>0(i=1,…,n)对称矩阵的Cholesky分解证明：对称矩阵的Cholesky分解对称矩阵的Cholesky分解对称矩阵的Cholesky分解推论：设A为对称正定矩阵，则存在唯一分解其中L为具有主对角元素为正数的下三角矩阵。对称矩阵的Cholesky分解证明：Cholesky分解的求法Cho

4、lesky分解的求法Cholesky分解的求法Cholesky分解法Cholesky分解法缺点及优点优点：可以减少存储单元。缺点：存在开方运算，可能会出现根号下负数。例题例题改进Cholesky分解法改进的cholesky分解A=LDLT改进的cholesky分解改进的cholesky分解改进的cholesky分解算法改进的cholesky分解算法例题例题例题例题A=LDLT分解，既适合于解对称正定方程组，也适合求解A为对称，而各阶顺序主子式不为零的方程组而对A=LLT只适合于对称正定方程组3.2.4三对角方程组求解的追赶法三对角方程组求解的追赶法三对角方程组求解的追赶法三对角方程组

5、求解的追赶法三对角方程组求解的追赶法其计算工作量为5n-4次乘除法。工作量小，其实现的条件为qi不为零。有以下定理可得证三对角矩阵求解的充分性条件。解三对角矩阵线性方程组的追赶法程序框图