网站大量收购闲置独家精品文档,联系QQ:2885784924

数值积分(基于MATLAB)课件 chapter3 线性代数方程组的迭代法 .ppt

数值积分(基于MATLAB)课件 chapter3 线性代数方程组的迭代法 .ppt

  1. 1、本文档共29页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多

**第3章3.5迭代法迭代法基础?问题?在实际应用中遇到的系数矩阵多为大型稀疏矩阵,如用求解线性方程组的直接法求解,在计算机上会耗费大量的时间和存储单元。在许多应用问题中使用迭代法。思路将改写为等价形式,建立迭代。从初值出发,得到序列。研究内容:?如何建立迭代格式??收敛速度??向量序列的收敛条件??误差估计?一般迭代法迭代格式的构造把矩阵A分裂为则将上式写为迭代过程这种迭代过程称为逐次逼近法,B称为迭代矩阵。收敛性定义:若称逐次逼近法收敛,否则,称逐次逼近法不收敛或发散。给定初值就得到向量序列问题:定理1任意给定初始向量x0,如果由逐次逼近法产生的向量序列收敛于向量x*,那么,x*是方程组x=Bx+g的解。证明:是否为方程组Ax=b的解?迭代法的收敛条件定理当k??时,Bk?0??(B)1定理2设线性方程组x=Bx+g有惟一解,那么逐次逼近法对任意初始向量X0收敛的充分必要条件是迭代矩阵B的谱半径?(B)1。证明:因此,注:要检验一个矩阵的谱半径小于1比较困难,所以我们希望用别的办法判断收敛性。定理3若逐次逼近法的迭代矩阵满足‖B‖1,那么逐次逼近法收敛。Remark:因为矩阵范数都可以直接用矩阵的元素计算,因此,用定理3,很容易判别逐次逼近法的收敛性。第3章3.6迭代法的收敛性定理3.3(充分条件)若存在一个矩阵范数使得||B||1,则迭代收敛,且有下列误差估计:②①证明:②迭代法的误差估计误差表达式及收敛速度。停机准则。①(1)1.雅克比(Jacobi)迭代法设有n阶方程组几种常用的迭代格式若系数矩阵非奇异,且(i=1,2,…,n),将方程组(1)改写成然后写成迭代格式(2)(2)式也可以简单地写为(3)写成矩阵形式:A=LUDBJacobi迭代阵(4)Algorithm:JacobiIterativeMethodSolve.Givenaninitialapproximation.Input:thenumberofequationsandunknownsn;thematrixentriesa[][];theentriesb[];theinitialapproximationX0[];toleranceTOL;maximumnumberofiterationsMmax.Output:approximatesolutionX[]oramessageoffailure.Step1Setk=1;Step2While(k?Mmax)dosteps3-6 Step3Fori=1,…,nSet;/*computexk*/ Step4IfthenOutput(X[]);STOP;/*successful*/ Step5Fori=1,…,nSetX0[]=X[];/*updateX0*/ Step6Setk++;Step7Output(Maximumnumberofiterationsexceeded);STOP./*unsuccessful*/Whatifaii=0?迭代过程中,A的元素不改变,故可以事先调整好A使得aii?0,否则A不可逆。必须等

您可能关注的文档

文档评论(0)

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

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

1亿VIP精品文档

相关文档