- 1、本文档共8页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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)