《雅克比迭代法与高斯-赛德尔》由会员分享，可在线阅读，更多相关《雅克比迭代法与高斯-赛德尔（14页珍藏版）》请在金锄头文库上搜索。
1、实验三实验三 线性方程组的迭代法线性方程组的迭代法班级：*计本 * 班 姓名：* 座号： * 时间：2010/6/2一、一、实验目的实验目的（1）熟悉 VC+开发平台和开发语言。（2）掌握雅可比及高斯-塞德尔迭代法解方程组的迭代法，并能根据给定的精度要求计 算（包括迭代过程） ；比较两种方法的优劣。（3）培养编程和上机调试能力。二、二、实验设备实验设备一台 PC 机，XP 操作系统，VC+软件三、三、实验内容实验内容用雅可比迭代法和高斯-塞德尔迭代法解方程组，其精度为 10-1231231238322041133631236xxxx axxxxx 5。四、四、算法描述算法描述1.雅克比迭代法公
2、式：X0k+1 = ( 1 / a00 )* b0 a01*X1k a02*X2k ;X1k+1 = ( 1 / a11 )* b1 a00*X0k a02*X2k ;X2k+1 = ( 1 / a22 )* b2 a00*X0k a01*X1k ;2.高斯-赛德尔迭代法公式：X0k+1 = ( 1 / a00 )* b0 a01*X1k a02*X2k ;X1k+1 = ( 1 / a11 )* b1 a00*X0k+1 a02*X2k ;X2k+1 = ( 1 / a22 )* b2 a00*X0k+1 a01*X1k+1 ;3.流程图雅克比流程图开始输入n,系数矩阵A，矩阵B，精度e，循
3、 环次数Max_XuanHuan，输出精度控制 prefor (z=0;zis+=aiz*XGzj-1;是s+=aiz*XGzj;误差e是否输出结果高斯-赛德尔流程图输入 n,a,b,MAX_ XunHuan,e ，keyKey=1，2，0？雅克比迭代法高斯-赛德尔迭代 法是1是2其他Key=0,.是0非0雅克比迭代公式高斯-赛德尔迭代公式输出结果结束开始五、五、程序代码程序代码/*雅克比方法类雅克比方法类 Jacobi.h*/#ifndef A#define A#define MAX_SIZE 50class Jacobipublic:void Solve(int n,double aMAX
4、_SIZEMAX_SIZE,double bMAX_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 XMAX_SIZEMAX_SIZE=0;double EMAX_SIZE;void Jacobi:Solve(int n,double aMAX_SIZ
5、EMAX_SIZE,double bMAX_SIZE,int MAX_XunHuan,double e,int pre)int i,j=1,z;double s=0;int flag=1;coutXi0;for (j=1;jEj)Ej=fabs(Xij-Xij-1);if (Eje)flag=1;if(flag=0)for (i=0;i0)if(i!=j-1)coutXGi0;for (j=1;ji)s+=aiz*XGzj-1;elses+=aiz*XGzj;XGij=1.0/aii*(bi-s);EGj=fabs(XG0j-XG0j-1);for (i=1;iEGj)EGj=fabs(XGij-XGij-1);if (EGje)flag=1;if(flag=0)for (i=0;i0)if(i!=j-1)coutn;coutaij;coutbi;coutMAX_XunHuan;coute;coutpre;coutkey;if (key=1)J.Solve(n,a,b,MAX_XunHuan,e,pre);if (key=2)G.Solve(n,a,b,MAX_XunHuan,e,pre); while (key);return 0;六、六、实验结果实验结果七、七、实验小结实验小结1.这次试验很简单，主要是公式问题，做得很顺利。
