下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验十五矩阵的分解【实验目的】1 了解矩阵LU,QR和Cholesky分解的基本概念。2 学习掌握MATLAB软件有关的命令。【实验准备】1 矩阵的LU分解LU分解是将一个方阵表示两个基本三角矩阵的乘积,其中一个三角矩阵为上三角阵,另一个为下三角阵,这里使用的是高斯变量消去法. 2矩阵的分解矩阵为实矩阵,且满足,其中为单位矩阵,则称为正交矩阵.QR分解能把任意长方阵分解为正交矩阵和上三角矩阵的初等变换形式的乘积,比如把矩阵进行如下分解,其中为初等变换矩阵.3矩阵的Cholesky分解Cholesky分解把矩阵分解为上三角矩阵和其转置的乘积,即,其中为上三角矩阵.如果复数矩阵满足Hermite正
2、定,也有Cholesky分解.3矩阵的MATLAB命令MATLAB中主要用lu,qr和chol分别完成矩阵的LU,QR和Cholesky分解。L,U=lu(A) 矩阵A的LU分解Q,R=qr(A) 矩阵A的QR分解R=chol(A) 矩阵A的Cholesky分解可以用help lu, help qr,help chol查阅有关这些命令的详细信息【实验方法与步骤】 练习1 求矩阵的LU分解相应的MATLAB代码和计算结果为:A=1 2 3;4 5 6;7 8 0; %创建矩阵L,U=lu(A) %进行LU分解L = 0.1429 1.0000 0 0.5714 0.5000 1.0000 1.0
3、000 0 0U = 7.0000 8.0000 0 0 0.8571 3.0000 0 0 4.5000练习求矩阵的QR分解.实际上是一个亏损矩阵,中间列是其他两列的平均,秩和亏损值由QR分解得出.相应的MATLAB代码和计算结果为: A=1 2 3;4 5 6;7 8 9;10 11 12A = 1 2 3 4 5 6 7 8 9 10 11 12Q,R=qr(A)Q = -0.0776 -0.8331 -0.2636 -0.4801 -0.3105 -0.4512 0.7093 0.4437 -0.5433 -0.0694 -0.6278 0.5530 -0.7762 0.3124 0.
4、1821 -0.5166R = -12.8841 -14.5916 -16.2992 0 -1.0413 -2.0826 0 0 -0.0000 0 0 0 本练习表明,当矩阵为长方形时,矩阵的末行元素为0,当矩阵的纵横方向上的元素个数相差较大时,中的零元素将被压缩速掉,中相应的列也被压缩掉.这种做法在大型复杂程序中很有用,不大可以节省磁盘空间,还能节省内存.练习3求5阶pascal矩阵的Cholesky分解.MATLAB代码和计算结果为:A=pascal(5) %5阶pascal矩阵A = 1 1 1 1 1 1 2 3 4 5 1 3 6 10 15 1 4 10 20 35 1 5 15 35 70R=chol(A) %Cholesky分解R = 1 1 1 1 1 0 1 2 3 4 0 0 1 3 6 0 0 0 1 4 0 0 0 0 1可以验证,有R*R=A.由于有了Cholesky分解,当解方程时,如果矩阵对称正定,Matlab会把方程代换为,让后按x=R(Rb)求解,计算速
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- (合作办学合同书)
- 民办学校竞业限制合同文本
- 2024年碱锰电池合作协议书
- 2024年住宿服务项目发展计划
- 2024年H-系列卷材涂料合作协议书
- 2024年石油产品添加剂:燃料油添加剂项目合作计划书
- 植物园装修拆除合同模板
- 2024版房产转让合同样本
- 污泥资源化运输服务合同
- 2024版购销合同与采购合同的合同履行义务
- 高铁隧道洞口陡坡危石清理专项方案(共17页)
- 综合布线 工程 11 综合布线系统的验收方案
- 信用证分析单(共享)
- 10KV配电调试报告
- 中医院新入职护士培训计划
- 水泵振动测试报告
- 小学高年级体育教案(全册)
- 涉密计算机审计报告
- 变电站工程质量评价报告范本
- 通用技术台灯作品方案
- 学校电子白板使用记录(样表)
评论
0/150
提交评论