- 1、本文档共18页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
解线性方程组的迭代法讲解
PAGE \* MERGEFORMAT1
MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 1 \h \* MERGEFORMAT 解线性方程组的迭代法
Haha
送给需要的学弟学妹
摘要:因为理论的分析表明,求解病态的线性方程组是困难的,但是实际情况是否如此,需要我们来具体检验。系数矩阵H为Hilbert矩阵,是著名的病态问题。因而决定求解此线性方程组来验证上述问题。
详细过程是通过用Gauss消去法、J迭代法、GS迭代法和SOR迭代法四种方法求解线性方程组。
关键词:病态方程组、Gauss消去法、J迭代法、GS迭代法、SOR迭代法
目录:
一、问题背景介绍
二、建立正确额数学模型
三、求解模型的数学原理
1、Gauss消去法求解原理
2、Jacobi迭代法求解原理
3、G-S迭代法求解原理
4、SOR迭代法求解原理
5、Jacobi和G-S两种迭代法收敛的充要条件
四、计算过程
(一)Hilbert矩阵维数n=6时
1、Gauss消去法求解
2、Jacobi迭代法求解
3、G-S迭代法求解
4、SOR迭代法求解
(二)Hilbert矩阵维数n=20、50和100时
1、G-S迭代法求解图形
2、SOR迭代法求解图形
五、编写计算程序
六、解释计算结果
1、Gauss消去法误差分析
2、G-S迭代法误差分析
3、SOR迭代法误差分析
G-S迭代法与SOR迭代法的误差比较
七、心得体会
正文:
一、问题背景介绍。
理论的分析表明,求解病态的线性方程组是困难的。实际情况是否如此,会出现怎样的现象呢?
二、建立正确的数学模型。
考虑方程组的求解,其中系数矩阵H为Hilbert矩阵,
这是一个著名的病态问题。通过首先给定解(为方便计算,笔者取x的各个分量等于1),再计算出右端这样的解就明确了,再用Gauss消去法、J迭代法、GS迭代法和SOR迭代法四种方法分别求解将求解结果与给定解比较,而后求出上述四种方法的误差,得出哪种方法比较好。
三、求解模型的数学原理。
1、Gauss消去法求解原理
对于(A非奇异)求解时,可以先将A分解成一个下三角矩阵L和一个上三角矩阵U的乘积,即,就可以通过
MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1)
求解出的值。
接下来就具体讲讲如何将A分解成L和U,也就是Gauss消去法。
欲把一个给定的矩阵A分解为一个下三角阵L与一个上三角阵U的乘积,最自然的做法便是通过一系列的初等变换,逐步将A约化为一个上三角阵,而又能保证这些变换的乘积是一个下三角阵。这可归结为:对于一个任意给定的向量找一个尽可能简单的下三角阵,使经这一矩阵作用之后的第至第个分量均为零。能够完成这一任务的最简单的下三角阵便是如下形式的初等下三角阵:
其中
即
这种类型的初等下三角阵称作Gauss变换,而称向量为Gauss向量。
对于一个给定的向量我们有
由此立即可知,只要取
便有
当然,这里我们要求
而后经过多次变换可以得到
MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2)
从而求出上三角阵U,而后通过
求得下三角阵
MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3)
将(1.2)和(1.3)带入到(1.1)式中求出的值即可。
2、J迭代法求解原理
考虑非奇异线性代数方程组
令
MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGE
文档评论(0)