线性方程组求解的数值方法.pptxVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多

*;01高斯消元法02矩阵分解法0;理解各种线性方程组数值求解;掌;02直接法:计算过程没有舍入误;(2.1)43.1高斯消元法;解:(1)消元:整个计算过程可;对于一般情形:n阶线性方程组的;若记增广矩阵(2.2);(1)第1步(k=1),一般形;设已完成上述消元过程第1步,第;01其中:02设;(即用乘以2.;0102求解三角形方程组2.3;而采用(,1)作为主元素;错误在哪个地方呢?原因是我;01第一列中绝对值最大是10,;x3=6.2/6.2=1x2=;整个矩阵中绝对值最大是10,取;x1=0Ax2=-1Bx3=1;1完全主元素消除法与列主元素消;对矩阵A实行初等行变换相当于用;第k步的初等矩阵为01将上三角;030201将上三角矩阵;消元法:消元法与三角分解法间的;01直接三角分解法解线性方程组;2x的计算1再求解Ly=b,;故:解:(一)矩阵LU分解例;经计算:;01求解x:02从而;§3.2矩阵的Choles;uii0(i=1,··;01故:02由分解的唯一性知:;经过n步可直接求得01思路:分;Cholesky方法具体计算公;向量范数的定义C向量、矩阵与;容易证明向量范数满足如下性质:;如果向量范数取为2-范数,则矩;类似地,Matlab函数nor;求解线性代数方程组的方法*3.;即线性方程组为:*线性方程组的;如果对该方程组进行变形,将变量;所以(1)式可以表示为迭代形式;考虑线性方程组的一般形式:其中;上面式子就是一般形式的Jaco;在Jacobi迭代法中充分利用程组的精确解;kx1(k)x2(k)x3(k;Gauss-Seidel迭代法;为了进一步研究,从矩阵角度来讨;其中或写成x(k+1)=Bx(;其中01Gauss-Seide;3.5迭代法的分析01构造迭;kx1(k)x2(k)x3(k;1设是矩阵B任意一个特;01先引入两个引理(详细证明见;设迭代法产生的序列;由;那么当然;定理3.2如果A是严格对角;矩阵的范数为由A是;定理3.3系数矩阵为严格对角;03所以存在模大于1的特征值.;故,只有才能使得;收敛速度问题证明:定理3.4;而所以故,本页第一个式子可写为;得证01得证显然对于任意的正整;定理3.4的讨论给出的关系可以;3.6超松弛迭代(SOR);这时迭代矩阵:其中;新的迭代每步的计算代价与Gau;定理3.5SOR方法收敛的;分块矩阵以及分块迭代方法其中;求解方程组;分块的超松弛迭代法(BSOR);3.7线性方程组的条件???;例设有方程组已知系数测量数据;如果测量数据和系数矩阵的数据都;考虑由于右端的扰动所引起的解的个式子给出了;又Tylor展开,;对;从上面三种情况可见,问题扰动所;病态问题与条件数(针对问题本身;其中定理3.6对任意的非;证明:01另外,中的各条性质是

文档评论(0)

gongmengjiao + 关注
实名认证
文档贡献者

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

1亿VIP精品文档

相关文档