高斯赛德尔迭代法解线性方程组

1、.数值分析实验五班级： 10信计二班 学号：59 姓名：王志桃 分数： 一实验名称高斯-赛德尔迭代法解线性方程组二实验目的1. 学会利用高斯赛德尔方法解线性方程组2. 明白迭代法的原理3. 对于大型稀疏矩阵方程组适用于迭代法比较简单三实验内容利用gauss-seidel迭代法求解下列方程组， 其中取。四、算法描述由jacobi迭代法中，每一次的迭代只用到前一次的迭代值，若每一次迭代充分利用当前最新的迭代值，即在计算第个分量时，用最新分量，代替旧分量，就得到所谓解方程组的gauss-seidel迭代法。其迭代格式为 (初始向量)， 或者写为 五、 编码#include#include精品.#in

2、clude#include#definemax_n 100#define precision0.0000001#define max_number1000void vectorinput(float x,int n) /输入初始向量int i;for(i=1;i=n;+i)printf(x%d=,i);scanf(%f,&xi);void matrixinput(float amax_n,int m,int n) /输入增广矩阵int i, j;printf(n=begin input matrix elements=n);for(i=1;i=m;+i)printf(input_line %d

3、 : ,i);for(j=1;j=n;+j)scanf(%f,&aij);void vectoroutput(float x,int n) /输出向量int i;for(i=1;i=n;+i)printf(nx%d=%f,i,xi);int issatisfypricision(float x1,float x2,int n) /判断是否在规定精度内int i;for(i=1;iprecision) return 1;精品.return 0;int jacobi_(float amax_n,float x,int n) /具体计算float x_formermax_n;int i,j,k;pr

4、intf(ninput vector x0:n);vectorinput(x,n);k=0;dofor(i=1;i=n;+i)printf(nx%d=%f,i,xi);x_formeri=xi;printf(n);for(i=1;i=n;+i)xi=ain+1;for(j=1;jprecision)xi/=aii;elsereturn 1;+k;while(issatisfypricision(x,x_former,n) & k=max_number)return 1;elseprintf(ng-s %d times!,k);return 0;int main() /主函数int n;float amax_nmax_n,xmax_n;精品.printf(ninput n=);scanf(%d,&n);if(n=max_n-1)printf(n007n must %d!,max_n);exit(0);matrixinput(a,n,n+1);if(jacobi_(a,x,n)printf(ng-s failed!);elseprintf(noutput solution:);vectoroutput(x,n);printf(nn007press any key to quit!n);getch();精品.精品.通过实验使我更加熟练的掌握和使用高斯赛贝尔迭代法来解线性方