矩阵的三角分解_第1页
矩阵的三角分解_第2页
矩阵的三角分解_第3页
矩阵的三角分解_第4页
矩阵的三角分解_第5页
已阅读5页,还剩43页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

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

2、tle分解9Doolittle分解10Doolittle分解11Doolittle分解12例题13例题14例题15例题16例题17Doolittle分解18Crout分解若矩阵A有分解:A=LU,其中L为下三角阵,U为单位上三角阵,则称该分解为Crout分解,若矩阵A的Doolitlle分解为A=LU ,则矩阵AT的Crout分解为UTLT。所以得到计算Crout分解的计算方法如下:19Cruou分解20Crout分解213.2.3 对称正定矩阵的Cholesky分解在应用数学中,线性方程组大多数的系数矩阵为对称正定这一性质,因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法,且

3、分解过程无需选主元,有良好的数值稳定性。 22对称正定矩阵的Cholesky分解 A对称:AT=A A正定:A的各阶顺序主子式均大于零。即 23对称正定矩阵的Cholesky分解24对称矩阵的Cholesky分解定理3.2.4 设A为对称正定矩阵,则存在唯一分解A=LDLT,其中L为单位下三角阵,D=diag(d1,d2,dn)且di0(i=1,n)25对称矩阵的Cholesky分解证明: 26对称矩阵的Cholesky分解27对称矩阵的Cholesky分解28对称矩阵的Cholesky分解推论:设A为对称正定矩阵,则存在唯一分解 其中L为具有主对角元素为正数的下三角矩阵。29对称矩阵的Cholesky分解证明: 30Cholesky分解的求法31Cholesky分解的求法32Cholesky分解的求法33Cholesky分解法Cholesky分解法缺点及优点 优点:可以减少存储单元。 缺点:存在开方运算,比较耗时。34改进Cholesky分解法改进的cholesky分解A=LDLT35改进的cholesky分解36改进的cholesky分解37改进的cholesky分解算法38改进的cholesky分解算法39例题40例题41例题42例题43A=LDLT分解,既适合于解对称正定方程组,也适合求解A为对称,而各阶顺序主子式不为零的方程组而对A=LLT只适合于对称正定方程组44

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论