c语言程序编程线性方程,C语言编程求解线性方程.doc

C语言编程求解线性方程

本 科 专 业 学 年 论 文

题目：线性方程组求解方法比较

姓 名 郭 凤

专 业 计算机科学与技术专业

班 级 08级本科(2)班

指导教师 刘 晓 娜

完成日期：2010 年 1月 8日

题目：线性方程组求解方法比较

摘 要

目前在许多实际应用领域,诸如航空、造船以及其它结构工程中,常遇到求解大型线性代数方程的问题。本文根据方程的。在第中方程

目 录

第一章 绪论 …………….………………………………………………………………………………………………………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、简要分析了算法的计算效果、稳定性、收敛效果、计算精度以及优劣性。

第二章 求解线性方程组的基本理论

迭代法的基本思想：是将线性方程组转化为便于迭代的等价方程组，对任选一组初始值(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的雅克比迭代法的迭代矩阵。

解雅克比迭代法的计算公