高斯迭代法报告.docxVIP

高斯迭代法实验报告专业班级学号姓名实验内容和要求用jacobi迭代法解线性方程组Ax=b，其中A是非奇异矩阵。用编织出来的程序，求解线性方程组Ax=b方程如下：=;算法说明Jacobi迭代格式的引出是依据迭代法的基本思想：构造一个向量系列，使其收敛至某个极限，则就是要求的方程组的准确解。迭代将方程组：在假设，改写成如果引用系数矩阵，及向量，，，方程组（1）和（2）分别可写为：及，这样就得到了迭代格式用迭代解方程组时，就可任意取初值带入迭代可知式，然后求。但是，比较大的时候，写方程组和是很麻烦的，如果直接由，能直接得到，就是矩阵与向量的运算了，那么如何得到，呢实际上，如果引进非奇异对角矩阵将分解成：要求的解，实质上就有而是非奇异的，所以存在，从而有我们在这里不妨令就得到迭代格式：现在考虑迭代法的计算程序分别代表的系数和等号右边的常数项即 。由Jacobi迭代法中，每一次的迭代只用到前一次的迭代值，若每一次迭代充分利用当前最新的迭代值，即在计算第个分量时，用最新分量，代替旧分量，，就得到所谓解方程组的Gauss-Seidel迭代法。其迭代格式为 (初始向量)，或者写为源程序#includemalloc.h#includestdio.h#includeconio.h#includemath.h#define EPS 1e-6#define MAX 100float *Jacobi(float a[3][4],int n);float *Jacobi(float a[3][4],int n){ float *x,*y; float epsilon,s; int i,j,k=0; x=(float *)malloc(n*sizeof(float)); y=(float *)malloc(n*sizeof(float)); for(i=0;in;i++) x[i]=0; while(1) { epsilon=0; k++; for(i=0;in;i++) { s=0; for(j=0;jn;j++) { if(j==i) { continue; } s+=a[i][j]*x[j]; } y[i]=(a[i][n]-s)/a[i][i]; epsilon+=(float)fabs(y[i]-x[i]); } if(epsilonEPS) { printf(\n迭代次数为：%d\n,k); return x; } if(k=MAX) { printf(方程不收敛\n); return y; } for(i=0;in;i++) { x[i]=y[i]; } }}main(){ int i, k, l; float a[3][4]; float *x; printf (请输入增广矩阵：\n); for (k = 0; k 3; k++) for (l = 0; l 4; l++) scanf (%f, a[k][l]); x=(float *)malloc(3*sizeof(float)); x=Jacobi(a,3); for(i=0;i3;i++) { printf(x[%d]=%f\t,i,x[i]); } printf(\n); getch();}高斯迭代部分：float *Gauss-Seidel(float a[3][4],int n);float *Gauss-Seidel(float a[3][4],int n){ float *x,*y; float epsilon,s,t; int i,j,k=0; x=(float *)malloc(n*sizeof(float)); y=(float *)malloc(n*sizeof(float)); for(i=0;in;i++) { x[i]=0; y[i]=0; } while(1) { epsilon=0; k++; for(i=0;in;i++) { t=0; for(j=0;ji;j++) { t+=a[i][j]*y[j]; } s=0; for(j=i+1;jn;j++) { s+=a[i][j]*x[j]; } y[i]=(a[i][n]-s-t)/a[i][i];/*计算x^(k+1)*/ epsilon+=(float)fabs(y[i]-x[i]); } if(epsilonEPS) { printf(\n迭代次数为：%d\n,k); return x; } if(k=MAX) { printf(方程不收敛\n); return y;