jacobi.c

平时完成作业做的几个小程序。有数值计算方面的梯形公式求积、simpson算法求积、Jacobi算法求解线性方程组、Gauss-Seidel法求解线性方程组

/*
*C语言实现Jacobi迭代法解方程组;
*实现特定方程组的求解并比较他们的收敛的快慢，eps=0.0000001;
*author:黄翔 date:2008-6-5
*/

#include <math.h>
#include <stdio.h>

#define eps 0.000001
#define true 1


void main()
{	
	double A[10][10]={{4,-1,0,0,0,0,0,0,0,0},  //系数矩阵
	                  {-1,4,-1,0,0,0,0,0,0,0},
					  {0,-1,4,-1,0,0,0,0,0,0},
                      {0,0,-1,4,-1,0,0,0,0,0},
					  {0,0,0,-1,4,-1,0,0,0,0},
					  {0,0,0,0,-1,4,-1,0,0,0},
					  {0,0,0,0,0,-1,4,-1,0,0},
					  {0,0,0,0,0,0,-1,4,-1,0},
					  {0,0,0,0,0,0,0,-1,4,-1},
					  {0,0,0,0,0,0,0,0,-1,4,}};
	double B[10]={4.0,11.0,15.0,18.0,21.0,24.0,27.0,30.0,37.0,45.5};//常数项
	
	double x[10]={0}; //初始解
	double X[10];     //迭代解
	double E[10];     //迭代误差
	double e;         //误差无穷范数，用于停机条件判断
	int i,j;          //用于循环条件中数组下标
	unsigned count=0; //迭代次数计数器
	double M,N,G;	  //迭代系数

//Jacobi迭代;

	while(true){
		for(i=0;i<10;i++){
			for(j=0,M=0;j<=i-1;j++){
				M=M-A[i][j]*x[j];
			}
			for(j=i+1,N=0;j<10;j++){
				N=N-A[i][j]*x[j];
			}
			G=1/A[i][i];

			X[i]=G*(M+N+B[i]);
		
			E[i]=fabs(X[i]-x[i]);
		}
	
		for(i=0,e=E[0];i<10;i++){//
			if(E[i]>e){
				e=E[i];
			}
		}

	    count++;  //计数器加一

		if(e<eps) break;  //判断范数误差是否达到精度要求，是则跳出循环停机
			
		for(i=0;i<10;i++){//否则以迭代求得的解赋给初始解作为出发点，继续进行迭代求解
			x[i]=X[i];
		}
	}

	printf("迭代次数\t%d\n",count);//输出求解结果
	printf("求得方程组的解为:\n");
	for(i=0;i<10;i++)
		printf("%.10f\n",X[i]);	
}