原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
线性方程组AX=B的数值计算方法实验
实验描述：
随着科学技术的发展，线性代数作为高等数学的一个重要组成部分，在科学实践中得到广泛的应用。本实验的通过C语言的算法设计以及编程，来实现高斯消元法、三角分解法和解线性方程组的迭代法（雅可比迭代法和高斯-赛德尔迭代法），对指定方程组进行求解。
实验原理：
1、高斯消去法：
运用高斯消去法解方程组，通常会用到初等变换，以此来得到与原系数矩阵等价的系数矩阵，达到消元的目的。初等变换有三种：(a)、(交换变换)对调方程组两行；(b)、用非零常数乘以方程组的某一行；(c)、将方程组的某一行乘以一个非零常数，再加到另一行。
通常利用(c)，即用一个方程乘以一个常数，再减去另一个方程来置换另一个方程。在方程组的增广矩阵中用类似的变换，可以化简系数矩阵，求出其中一个解，然后利用回代法，就可以解出所有的解。
2、选主元：
若在解方程组过程中，系数矩阵上的对角元素为零的话，会导致解出的结果不正确。所以在解方程组过程中要避免此种情况的出现，这就需要选择行的判定条件。经过行变换，使矩阵对角元素均不为零。这个过程称为选主元。选主元分平凡选主元和偏序选主元两种。平凡选主元：如果，不交换行；如果，寻找第行下满足的第一行，设行数为，然后交换第行和第行。这样新主元就是非零主元。偏序选主元：为了减小误差的传播，偏序选主元策略首先检查位于主对角线或主对角线下方第列的所有元素，确定行，它的元素绝对值最大。然后如果，则交换第行和第行。通常用偏序选主元，可以减小计算误差。
3、三角分解法：
由于求解上三角或下三角线性方程组很容易所以在解线性方程组时，可将系数矩阵分解为下三角矩阵和上三角矩阵。其中下三角矩阵的主对角线为1，上三角矩阵的对角线元素非零。有如下定理：
如果非奇异矩阵可表示为下三角矩阵和上三角矩阵的乘积：
(1)
则存在一个三角分解。而且，的对角线元素为1，的对角线元素非零。得到和后，可通过以下步骤得到：
(1)、利用前向替换法对方程组求解。
(2)、利用回代法对方程组求解。
4、雅可比迭代：
考察一般形式的线性方程组：
(2)
设从(2)中分离变出变量将它改写成：
(3)
据此建立雅可比迭代公式：
(4)
5、高斯-赛德尔迭代：
通过对雅可比迭代的观察，由于通常被认为是比更好的的近似值，所以在计算时用来替换是合理的。所以对雅可比迭代进行了改进，得到了高斯-赛德尔迭代，其收敛速度更快。它的迭代公式为：
(5)
实验内容：
许多科学应用包含的矩阵带有很多零。在实际情况中很重要的三角形线性方程组有如下形式:
(6)
构造一个程序求解三角形线性方程组。可假定不需要行变换，而且可用第k行消去第k+1行的。
使用C语言编写一个程序求解线性方程组，其中：
使用2中的程序求解线性方程组，其中= ；而且，当时，。对的情况分别求解。精确解为。对得到的结果与精确解的差异进行解释。
修改2中的程序，使得它可以通过重复求解个线性方程组
其中 (7)
来得到 ，
则
(8)
而且
(9)
设有如下三角线性方程组，而且系数矩阵具有严格的对角优势： (10)
（1）、根据解线性方程组的迭代法原理，设计一个算法来求解上述方程组，算法必须有效利用系数矩阵的稀疏性。
、用（1）中的程序解下列的三角线性方程组：
(a) (b)
利用高斯-赛德尔迭代法求解下列带状方程：
四、计算结果及分析：
此题由于没有提供具体的实例计算，所以,以《数值方法(MATLAB版)》(第四版)第107页的11、12题为例，进行代码的计算验证。此题以高斯消去法作为算法的核心，并依此来进行程序的编写。
实例：
第11题:求解下列线性方程组。
第12题: 求解下列线性方程组。
第11题的正确解为:
第12题的正确解为:
经程序的运算：
题号
X值
(11)、此方程组的解为：
(12)、此方程组的解为：
X[0]
1.0000
2.0000
X[1]
3.0000
3.0000
X[2]
2.0000
-2.0000
X[3]
-2.0000
1.0000
经过与正确答案比较，可知计算结果在保留四位小数的时候，得出比较精确的结果。
此题以带选主元的三角分解法作为算
文档评论（0）