矩阵论课件:矩阵分解_第1页
矩阵论课件:矩阵分解_第2页
矩阵论课件:矩阵分解_第3页
矩阵论课件:矩阵分解_第4页
矩阵论课件:矩阵分解_第5页
已阅读5页,还剩51页未读 继续免费阅读

下载本文档

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

文档简介

矩阵分解高斯消去法的矩阵表示LU分解法(LUDecomposition)平方根法

(CholeskyDecomposition)3.1

直接三角分解法矩阵三角分解

将一个矩阵分解成结构简单的三角形矩阵的乘积称为矩阵的三角分解。

顺序高斯消去法其实就是一个矩阵的三角分解过程。

LU分解先求解y再求解x则A(k)

与A(k+1)之间的关系式可以表示为:其中:(i=k+1,…,n),将Gauss消去过程中第k-1步消元后的系数矩阵记为:(k=1,…,n-1)LU分解记:,则其中:L---单位下三角矩阵,U---上三角矩阵LU分解于是有:容易验证:(k=1,…,n-1)(杜利脱尔Doolittle分解)LU

分解的存在唯一性LU分解存在顺序高斯消去法不被中断?定理

顺序高斯消去法求解方程组Ax=b时,所有主元

的充要条件是:A的所有顺序主子式不为零。定理若A的所有顺序主子式不等于0,则A存在唯一的LU分解。(LU分解的唯一性

)LU分解紧凑方式直接利用矩阵乘法来计算LU分解(待定系数法)

比较等式两边的第一行得:u1j=a1j比较等式两边的第一列得:

比较等式两边的第二行得:比较等式两边的第二列得:(j=1,…,n)(i=2,…,n)(j=2,…,n)(i=3,…,n)U的第一行L的第一列U的第二行L的第二列待定系数法LU分解计算顺序

LU分解紧凑算法(续)第k

步:此时U的前k-1行和

L

的前k-1列已经求出比较等式两边的第k行得:比较等式两边的第k列得:直到第n

步,便可求出矩阵L和U的所有元素。(j=k,…,n

)(i=k+1,…,n

)LU分解紧凑算法(续)LU分解的算法:Fork=1,2,...,nEndForj=k,…,ni=k+1,…,n为了节省存储空间,通常用A

的绝对下三角部分来存放L(对角线元素无需存储),用

A

的上三角部分来存放U。矩阵分解3.2矩阵的QR分解

定理

任何实的非奇异n阶矩阵A可以分解成正交矩阵Q和上三角矩阵R的乘积,且除去相差一个对角线元素之绝对值全等于1的对角矩阵因子D外,分解式是唯一的。

矩阵分解3.3矩阵的满秩分解

矩阵分解3.4矩阵的奇异值分解

温馨提示

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

最新文档

评论

0/150

提交评论