c语言编程求解线性方程组论文.docVIP

数值计算小报告 题目：线性方程组求解方法比较 姓 名 和 丽 专 业 软件工程 班 级 11级软件（2）班 完成日期：2013 年 5月 18日 摘 要 目前在许多实际应用领域,诸如航空、造船以及其它结构工程中,常遇到求解大型线性代数方程的问题。本文根据方程的。在第中方程 目 录 第一章 绪论 …………….………………………………………………………………………………………………………4 第二章 求解线性方程组的基本理论 2.1 迭代法 …………………………………………………………………………………………………………….5 2.2 直接三角分解法 …………………………………………………………………………………………….6 2.3 高斯消去法 …………………………………………………………………………………………………….7 第三章 三种算法求解方程组实例 3.1 迭代法 ……………………………………………………………………………………………………………8 3.2 直接三角分解法 …………………………………………………………………………………………..10 3.3 高斯列主元消去法 ……………………………………………………………………………………….14 3.4 三种方法的优缺点比较 ……………………………………………………………………..…….16 参考文献 ………………………………………………………………………………………………………..17  第一章 绪 论 计算数学是数学学科的一大分支，它研究如何借助于计算机求解各类数值问题。应用计算机求解各类数值问题需要经历以下几个主要过程：1、实际问题2、数学模型3、计算方法4、算法设计5、计算求解 目前已有的数学软件可以帮助我们实现上机计算，基本上已经将数值分析的主要内容设计成简单的函数，只要调用这些函数进行运算便可得到数值结果。 数值分析的内通包括线性代数方程组求解、非线性代数方程（组）求解、矩阵的特征值与特征值向量的计算、函数插值、函数逼近、数值积分与数值微分以及微分方程数值解法。 线性方程组的求解从理论上可分为两类：直接法和迭代法。直接法是不考虑计算过程中的舍入误差，经过有限次的运算得到方程组精确解的方法，常见的方法是高斯顺序消去法、高斯列主元消去法和矩阵的LU分解法。迭代法是采用某种极限过程，用线性代数方程组的近似解逐步逼近精确解的方法。迭代法中常见的方法有简单迭代法、J-迭代法、GS-迭代法和SOR-迭代法。 本文主要是分析高斯列主元消去法、矩阵的LU分解法和简单迭代法理论上的异同，并用C语言程序通过具体实例进行了分析比较。 本文将线性方程组的求解过程用计算机实现，本文的编写由以下几个特点： 1、对于难点问题从具体模型引入，淡化抽象的概念与定理，通俗易通； 2、对于具体模型本文给出了多种解题的思想及方法； 3、对问题进行简洁易懂的理论证明，突出了线性代数的理论和基本思想，使数学方法更加利于理解掌握。 4、简要分析了算法的计算效果、稳定性、收敛效果、计算精度以及优劣性。 第1章 求解线性方程组的基本理论 迭代法的基本思想：是将线性方程组转化为便于迭代的等价方程组，对任选一组初始值（i=1,2…n），按某种计算规则，不断地对所得到的值进行修正，最终获得满足精度要求的方程组的近似解。 对于线性方程组Ax=b 其中，A为非奇异矩阵。 将A分裂为A=M-N，其中，M为非奇异矩阵，且要求线性代数方程组Mx=d容易求解，一般选择为A的某一部分元素构成的矩阵，称M为A的分裂矩阵。于是，求解Ax=b转化为求解Mx=Nx+b,由此可构造一个迭代法： x(0)(初始向量) ， x(k+1)=Bx(k)+f (k=0,1,2…) 其中，f=b/M,B=I-A/M为迭代法的迭代矩阵。 选取M为A的对角元素组成的矩阵，即选取M=D，可得到解Ax=b的迭代法： x(0)(初始向量),x(k+1)=Bx(k)+f (k=0,1,2…) BJ为求解Ax=b的雅克比迭代法的迭代矩阵。 解迭代法的计算公式为： （k=0,1,2，……:i=1,2,3,……..n） 迭代法方法是求对称矩阵的全部特征值以及相应的特征向量的一种方法，，它是基于以下两个结论： 1)任何实对称矩阵A可以通过正交相似变换成对角型，即存在正交矩阵Q，使得 AQ=diag() 其中i（i=1,2,…,n）是A的特征值，Q