数值分析.南京电大 36讲 szfx04.doc

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
4、LU分解的紧凑格式 我们不必在高斯消去法过程中产生L和U,而是直接用矩阵A来进行LU分解。下面导出计算公式,因A=LU则 利用矩阵乘法,得 由此可以得到和的计算公式: 先计算U的第一行。在U的第行和L的第列计算出来后再计算L的第j列: 在U的第行和L的第列均计算出来后再计算U的第i行:。计算的过程是:U第1行,L第1列,U第2行,L第2列,……顺序计算。 1.计算流程 2.存储问题和全部流程 U可以存储在A的上三角部分,L可以存储在A的下三角部分,对角元存U的对角元,A的对角元为1而不必存储,如下图: 用LU分解来解线性方程组的全部流程如下: 例1:利用LU分解求解线性方程组 解:第一步,A的LU分解 第二步,求解 解得: ,, ,。 第二步,求解 解得: , , , 。 二、追赶法 1. 三对角矩阵:在实际问题中,经常会遇到如下形式的方程组, 其系数矩阵为 没有写的部分均为零(下同),对角元为,后次对角线上的元素为,前次对角线上的元素为。 我们称矩阵A为三对角矩阵,相应的方程组称三对角方程组。如果对A进行LU分解,则将有如下形式: 并称它们为二对角矩阵。 2.定理3:设上述三对角矩阵A满足 且 则对矩阵A的LU分解能进行,且分解是唯一的。 3.追赶法的计算公式 利用矩阵乘法可得: 从而可以得到: 解得: 解得: , 4.追赶法的计算流程 , 第一个循环称之为追的过程,相当于消元过程; 第二个循环称之为赶的过程,相当于回代过程。 例2:用追赶法求解方程组 解:追的过程: ,; , ,; ,,;,,。 赶的过程: , , , 。 14

文档评论(0)

0520 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档