数学实验题目5 相对Gauss列主元消去法.doc

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

数学实验题目5 相对Gauss列主元消去法 摘要 由一般线性方程组在使用Gauss消去法求解时,从求解过程中可以清楚地看到,若,必须施以行交换的手续,才能使消去过程继续下去。有时既使,但其绝对值很小,由于舍入误差的影响,消去过程也会出现不稳定现象。因此,为使这种不稳定现象发生的可能性减至最小,在施行消去过程时每一步都要选主元素,即要寻找行,使 并将第行与第行交换,以使的当前值(即的数值)远大于0。 这种列主元消去法的主要步骤如下: 1.消元过程 对,做 1o 选主元,记 若,说明方程组系数矩阵奇异,则停止计算,否则进行2o。 2o 交换(增广矩阵)的两行元素 3o 计算 2.回代过程 对,计算 前言 利用Gauss列主元消去法、显式相对Gauss列主元消去法、隐式相对Gauss列主元消去法求解线性方程组 程序设计流程 问题1 程序运行如下: x = GaussSysSolve(Mat1_1,b1_1) x = 1.0000 1.0000 1.0000 1.0000 x = GaussExpSysSolve(Mat1_1,b1_1) x = 1.0000 1.0000 1.0000 1.0000 x = GaussIneSysSolve(Mat1_1,b1_1) x = 1.0000 1.0000 1.0000 1.0000 程序运行如下: x = GaussSysSolve(Mat1_2,b1_2) x = 1.0000 1.0000 1.0000 1.0000 x = GaussExpSysSolve(Mat1_2,b1_2) x = 1.0000 1.0000 1.0000 1.0000 x = GaussIneSysSolve(Mat1_2,b1_2) x = 1.0000 1.0000 1.0000 1.0000 程序运行如下: x = GaussSysSolve(Mat1_3,b1_3) x = 1.0000 1.0000 1.0000 1.0000 x = GaussExpSysSolve(Mat1_3,b1_3) x = 1.0000 1.0000 1.0000 1.0000 x = GaussIneSysSolve(Mat1_3,b1_3) x = 1.0000 1.0000 1.0000 1.0000 程序运行如下: x = GaussSysSolve(Mat1_4,b1_4) x = 1.0000 1.0000 1.0000 1.0000 x = GaussExpSysSolve(Mat1_4,b1_4) x = 1.0000 1.0000 1.0000 1.0000 x = GaussIneSysSolve(Mat1_4,b1_4) x = 1.0000 1.0000 1.0000 1.0000 问题2 程序运行如下: = GaussSysSolve(Mat2_1,b2_1) x = 1.0915 0.2832 1.1463 -0.1008 x = GaussExpSysSolve(Mat2_1,b2_1) x = 1.0915 0.2832 1.1463 -0.1008 x = GaussIneSysSolve(Mat2_1,b2_1) x = 1.0915 0.2832 1.1463 -0.1008 程序运行如下: x = GaussSysSolve(Mat2_2,b2_2) x = 0.5162 0.4152 0.1100 1.0365 x = GaussExpSysSolve(Mat2_2,b2_2) x = 0.5162 0.4152 0.1100 1.0365 x = GaussIneSysSolve(Mat2_2,b2_2) x = 0.5162 0.4152 0.1100 1.0365 程序运行如下: x = GaussSysSolve(Mat2_3,b2_3) x = 1.0000 1.0000 1.0000 x = GaussExpSysSolve(Mat2_3,b2_3) x = 1 1 1 x = GaussIneSysSolve(Mat2_3,b2_3) x = 1.0000 1.0000 1.0000 程序运行如下: x = GaussSysSolve(Mat2_4,b2_

文档评论(0)

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

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

1亿VIP精品文档

相关文档