雅克比迭代法与高斯-赛德尔.docVIP

实验三 线性方程组的迭代法 班级：**计本 ** 班 姓名：** 座号： ** 时间：2010/6/2 实验目的 熟悉VC++开发平台和开发语言。 掌握雅可比及高斯-塞德尔迭代法解方程组的迭代法，并能根据给定的精度要求计算（包括迭代过程）；比较两种方法的优劣。 培养编程和上机调试能力。 实验设备 一台PC机，XP操作系统，VC++软件 实验内容 用雅可比迭代法和高斯-塞德尔迭代法解方程组，其精度为10-5。 算法描述 雅克比迭代法公式： X[0][k+1] = ( 1 / a[0][0] )*[ b[0] –a[0][1]*X[1][k] – a[0][2]*X[2][k] ]; X[1][k+1] = ( 1 / a[1][1] )*[ b[1] –a[0][0]*X[0][k] – a[0][2]*X[2][k] ]; X[2][k+1] = ( 1 / a[2][2] )*[ b[2] –a[0][0]*X[0][k] – a[0][1]*X[1][k] ]; 高斯-赛德尔迭代法公式： X[0][k+1] = ( 1 / a[0][0] )*[ b[0] –a[0][1]*X[1][k] – a[0][2]*X[2][k] ]; X[1][k+1] = ( 1 / a[1][1] )*[ b[1] –a[0][0]*X[0][k+1] – a[0][2]*X[2][k] ]; X[2][k+1] = ( 1 / a[2][2] )*[ b[2] –a[0][0]*X[0][k+1] – a[0][1]*X[1][k+1] ]; 流程图 雅克比流程图 高斯-赛德尔流程图 程序代码 /*雅克比方法类Jacobi.h*/ #ifndef A #define A #define MAX_SIZE 50 class Jacobi{ public: void Solve(int n,double a[MAX_SIZE][MAX_SIZE],double b[MAX_SIZE],int MAX_XunHuan,double e,int pre); }; #endif /*雅克比方法类方法实现Jacobi.cpp*/ #ifndef B #define B #include Jacobi.h #include iostream #include math.h #include iomanip using namespace std; double X[MAX_SIZE][MAX_SIZE]={0}; double E[MAX_SIZE]; void Jacobi::Solve(int n,double a[MAX_SIZE][MAX_SIZE],double b[MAX_SIZE],int MAX_XunHuan,double e,int pre) { int i,j=1,z; double s=0; int flag=1; cout请输入初始值X[i][0]:; for (i=0;in;i++) { cinX[i][0]; } for (j=1;jMAX_XunHuanflag;j++) { flag=0; for (i=0;in;i++) { s=0; for (z=0;zn;z++) { if (i!=z) { s+=a[i][z]*X[z][j-1]; } } X[i][j]=1.0/a[i][i]*(b[i]-s); } E[j]=fabs(X[0][j]-X[0][j-1]); for (i=1;in;i++) { if(fabs(X[i][j]-X[i][j-1])E[j]) E[j]=fabs(X[i][j]-X[i][j-1]); } if (E[j]e) { flag=1; } } if(flag==0) { for (i=0;in;i++) { cout Xi(k) ; } cout 误差endl; for (i=0;ij;i++) { coutsetiosflags(ios::left)setiosflags(ios::fixed)setprecision(pre)X[0][i] X[1][i] X[2][i] ; if(i0) { if(i!=j-1) coutE[i]endl; else coutE[i] eendl; } else cout ∞endl; } } else cout没有结果,最大循环次数不够\n; cout迭代次数为jendl; } #