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

数值分析实验五 班级： 10信计二班 学号：59 姓名：王志桃 分数： 一．实验名称 高斯-赛德尔迭代法解线性方程组 二．实验目的 学会利用高斯赛德尔方法解线性方程组 明白迭代法的原理 对于大型稀疏矩阵方程组适用于迭代法比较简单 三．实验内容 利用Gauss-Seidel迭代法求解下列方程组 ， 其中取。 四、算法描述 由Jacobi迭代法中，每一次的迭代只用到前一次的迭代值，若每一次迭代充分利用当前必威体育精装版的迭代值，即在计算第个分量时，用必威体育精装版分量，代替旧分量，，就得到所谓解方程组的Gauss-Seidel迭代法。 其迭代格式为 (初始向量)， 或者写为 五、 编码 #include stdio.h #include stdlib.h #include conio.h #include math.h #define MAX_n 100 #define PRECISION 0.0000001 #define MAX_Number 1000 void VectorInput(float x[],int n) //输入初始向量 { int i; for(i=1;i=n;++i) { printf(x[%d]=,i); scanf(%f,x[i]); } } void MatrixInput(float A[][MAX_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 : ,i); for(j=1;j=n;++j) scanf(%f,A[i][j]); } } void VectorOutput(float x[],int n) //输出向量 { int i; for(i=1;i=n;++i) printf(\nx[%d]=%f,i,x[i]); } int IsSatisfyPricision(float x1[],float x2[],int n) //判断是否在规定精度内 { int i; for(i=1;i=n;++i) if(fabs(x1[i]-x2[i])PRECISION) return 1; return 0; } int Jacobi_(float A[][MAX_n],float x[],int n) //具体计算 { float x_former[MAX_n]; int i,j,k; printf(\nInput vector x0:\n); VectorInput(x,n); k=0; do{ for(i=1;i=n;++i) { printf(\nx[%d]=%f,i,x[i]); x_former[i]=x[i]; } printf(\n); for(i=1;i=n;++i) { x[i]=A[i][n+1]; for(j=1;j=n;++j) if(j!=i) x[i]-=A[i][j]*x[j]; if(fabs(A[i][i])PRECISION) x[i]/=A[i][i]; else return 1; } ++k; }while(IsSatisfyPricision(x,x_former,n) kMAX_Number); if(k=MAX_Number) return 1; else { printf(\nG-S %d times!,k); return 0; } } int main() //主函数 { int n; float A[MAX_n][MAX_n],x[MAX_n]; printf(\nInput n=); scanf(%d,n); if(n=MAX_n-1) { printf(\n\007n must %d!,MAX_n); exit(0); } MatrixInput(A,n,n+1); if(Jacobi_(A,x,n)) printf(\nG-S Failed!); else { printf(\nOutput Solution:); VectorOutput(x,n); } printf(\n\n\007Press any key to quit!\n); getch(); } 通过实验使我更加熟练的掌握和使用高斯赛贝尔迭代法来解线性方程组，而且这个方法可以用来解大型的稀疏