jacobi迭代法求解线性方程组.doc

//======================================================== //此程序用gauss-seidel迭代法对方程组进行近似求解 //======================================================== #includestdio.h #includestdlib.h #includemath.h //======================================================================== //全局变量 #define N 20 double a[N][N]; //声明系数矩阵 double b[N]; //声明方程组右边的值（数组） double eps; //声明相对误差变量 double x[3]={0,0,0}; int iteration; //函数声明 void GetData(); //此函数用以进行对系数矩阵赋值 void jacobi(); /*声明jacobi()函数， 用以对方程的增广矩阵进行简单迭代*/ void show(); //显示最后的结果 //===================================================================== //以下是主函数 void main() { printf(============================================================\n\n); printf( 此程序运用Jacobi迭代法求解方程组的解 \n\n); printf(============================================================\n\n); GetData(); jacobi(); show(); printf(============================================================\n\n); } //====================================================================== //一下定义GetData()函数 void GetData() { FILE *filename;//定义文件的指针 int i,j; if((filename=fopen(D:\\gauss_seidel.dat,r+))==NULL) { printf(文件打开失败，退出程序!); exit(0); } for(i=0;i3;i++) fscanf(filename,%lf,%lf,%lf,%lf,a[i][0],a[i][1],a[i][2],b[i]); printf(请输入误差允许量：); scanf(%lf,eps); printf(\n\n); printf(对应的系数和方程右边的结果是：\n\n); for(i=0;i3;i++) { for(j=0;j3;j++) { printf(%10lf ,a[i][j]); } printf(%10lf\n\n,b[i]); } printf(求解方程的准许误差是：%8lf\n\n,eps); } //================================================================================== //以下是对gauss-seidel函数的定义 void jacobi() { int i,j,flag; double y[3],sum,temp[3],k; do { flag=1; iteration++; for(i=0;i3;i++) { sum=0.0; for(j=0;j3;j++) { if(j!=i) sum+=a[i][j]*x[j]; } y[i]=(b[i]-sum)/a[i][i]; } {//此段代码求y[i]-x[i]绝对值得最大值 for(i=0;i3;i++) temp[i]=fabs(y[i]-x[i]); for(i=0;i2;i++) for(j=1;j3;j++)