雅可比迭代法和高斯-塞德尔迭代法求解线性方程组.doc

实验报告内容 一 实验目的与要求（实验题目）1分别利用雅可比迭代法和高斯-塞德尔迭代法求解以下线性方程组 使得误差不超过 2.用不动点迭代法求方程的实根： 二 模型建立（相关主要计算公式）1. 雅可比迭代法 其中为初始向量.2.高斯-塞德尔迭代法3.不动点迭代法 三、实验过程、步骤（程序）1. 雅可比迭代法#include stdio.h#include math.h#include string.hmain()int i,j,k;float m1=0.0,m2=0.0;float a34=8,-3,2,20,4,11,-1,33,6,3,12,36;float x3=0.0,0.0,0.0;for(k=1;k=10;) for(i=0;i=2;i+) for(j=0;ji;j+) m1=m1+aij*xj; for(j=i+1;j=2;j+) m2=m2+aij*xj; xi=(ai3-m1-m2)/aii; m1=0,m2=0; k+; printf(雅可比迭代法计算结果为：n);for(i=0;i=2;i+)printf(x%2d=%8.9fn,i+1,xi);2高斯-塞德尔迭代法#include#include# define n 3void main()int i,j,k=1;float xn=0,0,0,mn=0,0,0,s=1;float ann=8,-3,2,4,11,-1,6,3,12,dn=20,33,36; printf(高斯-塞德尔迭代法运算结果为:n);for(k=0;fabs(s-x0)1e-6;k+)s=x0;for(i=0;in;i+) mi=0; for(j=0;jn;j+) mi=mi-aij*xj; mi=mi+di+aii*xi; xi=mi/aii; printf(Y1=%f Y2=%f Y3=%fn,x0,x1,x2); getchar() ;3.#include #include double f( double x ) return x * x * x + 2 * x * x + 10 * x - 20;double fdx( double x ) return 3 * x * x + 18.4 * x + 16.7;int main( ) int t1 = 0, t2 = 1; double x 2 , ep = 1e-8; x 0 = 0; do t1 = 1 - t1; t2 = 1 - t2; x t1 = x t2 - f( x t2 ) / fdx( x t2 ); while( fabs( x t1 - x t2 ) ep ); printf(解得x=%lfn, x t1 ); return 0;四实验结果：1.雅可比迭代法：2.高斯-塞德尔迭代法：. 3.不动点迭代法：五实验小结 通过这次上机，学会了用Jacobis迭代法，高斯-塞德尔迭代法求解线性方程组，算法程序比较复杂，特别是要多次使用数组条件及