《zfeAAAGauss-Seidel迭代法》由会员分享，可在线阅读，更多相关《zfeAAAGauss-Seidel迭代法（9页珍藏版）》请在金锄头文库上搜索。
1、 数值分析课程论文 姓名: 学号: Gauss-Seidel迭代法求解线性方程组摘要 线性方程组的求解在许多的工程技术中是一个极为常见的问题，对于线性方程组的求解无论从理论上还是实践应用上都已经成熟.对于一般线性方程组的求解有Gauss消元法为基础的直接法，也有迭代法.其中Gauss-Seidel是一个重要的组成部分.鉴于此，本论文细致地研究了用Gauss-Seidel迭代法求解线性方程组. 论文的第一部分先介绍了迭代法求解线性方程组的一般模式，并给出这种迭代法的收敛性条件，Gauss-Seidel迭代法求解线性方程组的基本原理.这一部分是Gauss-Seidel迭代法的理论基础.论文的第二部
2、分给出了Gauss-Seidel迭代法的具体操作步骤，以伪代码的形式细致的描绘如何使用Gauss-Seidel迭代法的求解方程组.同时，为了验证算法的有效性，在这一部分，还引入一个简单的算例，用于MATLAB编程发现计算结果完全正确.论文的第三部分给出了关于Gauss-Seidel迭代法的MATLAB程序，用于计算线性方程组.关键词：Gauss-Seidel迭代法，基本原理，算例，MATLAB程序目 录1 Gauss-Seidel迭代法的基本理论.11.1线性方程组的迭代法求解.11.2Gauss-Seidel迭代法的原理.22.具体的算例和操作步骤.32.1. Gauss-Seidel迭代法
3、的伪代码.32.2.具体的算例验证算法的有效性.33MATLAB程序.4参考文献.6 Gauss-Seidel迭代法求解线性方程组一. Gauss-Seidel迭代法的基本理论1.1线性方程组的迭代法求解在考虑求解线性方程组Ax=b时，其中A为非奇异矩阵.尽管Guass消元法通过有限次运算可以求解此问题，其对应的计算复杂度为.但是对于工程技术中和某些偏微分方程过程中出现的大型稀疏型矩阵利用迭代法可以更快的收敛，找到解.另外一方面，由于迭代法占用的计算机内存少，且便于计算.这两方面的优势促成了迭代法求解线性方程组的研究.关于迭代法的收敛的几个判定条件1（迭代法基本原理）设有方程组，对于任意初始向
4、量及任意，解此方程组的迭代法（即）收敛的充要条件是.2（迭代法收敛的充分条件）如果方程组的迭代公式为（为任意初始向量），且迭代矩阵的某一种范数，则：迭代法收敛；.定理3如果为严格对角占优阵或为不可约弱对角占优阵，则对于任意的，解方程组的Jacobi迭代法，Gauss-Seidel迭代法均收敛.定理4如果为对称正定矩阵，且，则解式的SOR方法收敛.1.2Gauss-Seidel迭代法的原理由Jacobi方法迭代公式可知，迭代的每一步计算过程，都是用的全部分量来计算的所有分量，显然在计算第个分量时，已经计算出的最新分量，没有被利用.从直观上看，最新计算出的分量可能比旧的分量要好些.因此，对这些最新
5、计算出来的第次近似的分量加以利用，就得到所谓解方程组的Gauss-Seidel迭代法（简称G-S方法）： （初始向量）， 或写为上面第2个式子利用了最新计算出的分量，第个式子利用了计算出的最新分量.还可写成矩阵形式 ，若设存在，则 ，于是Gauss-Seidel迭代公式的矩阵形式为 ，其中 ，.由此可以看出，应用Gauss-Seidel迭代法解式，就是对方程组应用迭代法.称为解式的Gauss-Seidel迭代法的迭代矩阵.Gauss-Seidel迭代法的一个明显优点是，在用计算机计算时，只需一组工作单元，以便存放近似解.由式可以看出，每迭代一步只需计算一次矩阵与向量的乘法.二具体的算例和操作步
6、骤2.1. Gauss-Seidel迭代法的伪代码1.输入问题的参数A,b2.分解A为D,L,U.3.计算迭代方程G，f.4.开始迭代，随机设定一个初值.5.以迭代方程更新x的值.6.如果到达迭代次数，则进入步骤7；否则，回到步骤5.7.输出x，结束.2.2.具体的算例验证算法的有效性求解如下的线性方程组这个方程的真实解为（3,2,1）.程序运行结果：情况1：输入GS（A,b）GS(A,b)xhis = 0 0 0 2.5000 2.0909 1.2273 2.9773 2.0289 1.0041 3.0098 1.9968 0.9959 2.9998 1.9997 1.0002 2.9998
7、 2.0001 1.0001 3.0000 2.0000 1.0000 3.0000 2.0000 1.0000 3.0000 2.0000 1.0000 3.0000 2.0000 1.0000 3.0000 2.0000 1.0000ans = 3.0000 2.00001.0000 图一。解的迭代过程程序说明：xhis为迭代过程，ans为最终迭代结果.本次迭代次数为10.图一中红点标记每次迭代的解的轨迹。从图中可以看出，初始值选择在（0,0,0），经过几次就应经到达很高的精度，表明GaussSeidel的迭代法收敛很快。三MATLAB程序function y=GS(A,b)if(narg
8、in=0) disp(你没有输入任何参数); A=ceil(10*rand(3,3);b=ceil(10*rand(3,1); A,b GS(A,b)elseMaxInteration=10;n=length(A(:,1);D=zeros(n,n);L=zeros(n,n);U=zeros(n,n);for i=1:n for j=1:n if(i=j) D(i,j)=A(i,j); else if(ij) L(i,j)=-A(i,j); else U(i,j)=-A(i,j); end end endendG=inv(D-L)*U;f=inv(D-L)*b;x1=zeros(n,1);xhis=x1;for i=1:MaxInteration x2=G*x1+f; x1=x2; xhis=xhis,x1;endxhis=xhis;y=x1;if(n=3) plot3(xhis(:,1),xhis(:,2),xhis(:,3),or);title(解的迭代情况);endend参考文献1刘卫国，matlab程序设计教程，中国水利水电出版社，2010年2 李庆扬，王能超，易大义，数值分析，华中科技大学出版社，2010年3卓金武 魏永生 秦健 李必文，matlab在数学建模中的应用 北京航天航空出版社 2011年
