已阅读5页,还剩26页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第四节Gauss消去法的变形,7.4.1直接三角分解法,我们知道对矩阵进行一次初等变换,就相当于用相应的初等矩阵去左乘原来的矩阵。因此我们这个观点来考察Gauss消元法用矩阵乘法来表示,即可得到求解线性方程组的另一种直接法:矩阵的三角分解。,1.不选主元的三角分解法,矩阵三角分解原理,应用高斯消去法解n阶线性方程组Ax=b,经过n步消元之后,得出一个等价的上三角型方程组A(n)x=b(n),对上三角形方程组用逐步回代就可以求出解来。上述过程可通过矩阵分解来实现。将非奇异阵A分解成一个下三角阵L和一个上三角阵U的乘积A=LU称为对矩阵A的三角分解,又称LU分解。,L为单位下三角阵而U为一般上三角阵的分解称为Doolittle分解;L为一般下三角阵而U为单位上三角阵的分解称为Crout分解。,比较第1行:,比较第1列:,Doolittle分解,比较第2行:,比较第2列:,比较第k行:,比较第k列:,由Ly=b可得:,Doolittle法在计算机上实现是比较容易的,但如果按上述流程运算仍需要较大的存储空间:,因此可按下列方法存储数据:,直接三角分解的Doolittle法可以用以下过程表示:,存储单元(位置),紧凑格式的Doolittle法,杜丽特尔分解的紧凑格式,计算原则:,上表中U(y)的第一行等于A的第一行,其余元素等于A(b)的对应元素减去同一行左边L的元素与同一列上边U的元素对应乘积之和;L的元素求法与U的类似,只是最后要除以同一列U的对角元。,由Ux=y可得:,由此,求方程组的解就只剩下回代过程了。,即求解方程组Ux=y,例:用Doolittle分解求解方程组,解:由紧凑格式可得,2,1,2,2,-2,6,6,3,0,3,-1,-1,Ux=y,即,2.列主元的三角分解法,在Doolittle法(包括紧凑格式)中,反复用到公式,仍有可能为小主元做除数,为此,我们也要考虑在算法中加入选取列主元,我们下面介绍Doolittle列主元法,符号因换行只
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 餐饮业节能环保改造项目方案
- 呼叫中心员工激励方案设计
- 护理流程中电解质紊乱的识别与干预
- 大气污染防治工程技术方案
- 安全保工知识题库及答案解析
- 证劵从业考试吧及答案解析
- 装卸管理员从业资格考试及答案解析
- 环境监测数据采集与分析方法
- 企业内部审计流程与规范标准
- 护理专业急救知识题库大全及答案解析
- GB/T 32260.1-2015金属材料焊缝的破坏性试验焊件的冷裂纹试验弧焊方法第1部分:总则
- GB/T 15601-2013管法兰用金属包覆垫片
- 【原创】课题专题讲座-《抓好朗读训练播下语感种子》PPT
- FZ/T 01002-2010印染企业综合能耗计算办法及基本定额
- CJJ28-2014城镇供热管网工程施工及验收规范
- 美术技法理论试卷A
- 防爆电器基础知识
- 减员减支方案一阶段
- 苏教版四年级上册数学第三单元观察物体试卷【含答案】
- 银行笔试真题
- 安全检查记录 受限空间作业
评论
0/150
提交评论