原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《数值分析07解线性代数方程组的直接解法
第三章 线性代数方程组的数值解法 第一节 求解线性代数方程组的基本定理 第二节 高斯消元法及其计算机实现 一、三角形方程组的解法 二、高斯消元法 三、选主元高斯消元法 再交换第k,ik两个方程和第k,jk列,使 成为主元. 称这个过程为完全选主元. 不论是哪种方式选出主元,而后再按上面介绍的计 算步 骤进行消元的计算,一般都称为选主元的高斯消元法.在 实际计算中,常用按列选主元的高斯消元法. 例4:讨论下面方程组的解法 0.0001x1+x2=1 x1+x2=2 假设求解是在四位浮点十进制数的计算机上进行 0.1000?10-3 x1 + 0.1000 ?101 x2 = 0.1000 ?101 0.1000 ?101 x1 +0.1000 ?101 x2 = 0.2000 ?101 解:本题用机器数系表示为 a11 =0.0001, m21=a21/a11=1/0.0001= 104, 消元得 回代解得 x2=1 , x1=0 严重失真! (本题的准确解为 x1= 10000/9999, x2=9998/9999 ） a22(2)= 0.1000 ?101 - 104 ? 0.1000 ?101 = 0.00001 ?105 - 0.1000 ?105 (对阶计算) = - 0.1000 ?105 0.1000?10-3 x1 + 0.1000 ?101 x2 = 0.1000 ?101 -0.1000 ?105 x2 = -0.1000 ?105 主元a11过小 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 选主元基本思想 用高斯消元法求解线性方程组时,为避免小的主元.在进行第k步消元前,应该在第k列元素 (i=k,…,n)中找出第一个出现的绝对值最大者,例如 ， 再把第ik个方程与第k个方程组进行交换,使 成为主元.我们称这个过程为选主元.由于只在第k列元素中选主元,通常也称为按列选主元. 如果在第k步消元前,在第k个方程到第n个方程所有的xk到xn的系数 (i=k,…,n;j=k,…,n)中,找出绝对值最大者,例如 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. * 数值分析 数值分析 * 第二节 高斯消元法及其计算机实现 第三节 用矩阵分解法求解线性方程组 第四节 误差分析和解的精度改进 第五节 大型稀疏方程组的迭代法 第一节 求解线性代数方程组的基本定理 第六节 极小化方法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 线性代数方程组的一般形式 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. MATLAB实现: x=A\b Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 数
文档评论（0）