数值分析实验报告三迭代法解线性方程组.docxVIP

数学与计算科学学院实 验 报 告实验项目名称 迭代法解线性方程组 所属课程名称 数值方法A 实 验 类 型 验证型 实 验 日 期 2014年12月3日 学 号 班 级 信计 成 绩 一、实验概述：【实验目的】通过实验，熟悉线性方程组的Gauss-Seidel迭代法和松弛迭代法求解【实验内容】用熟悉的计算机语言编程上机完成用Gauss-Seidel迭代法求解下列方程组，当满足时结束迭代2. 用松弛迭代法（取）求解下阿列方程组，当满足时结束迭代，说明迭代的收敛性。并讨论初值的变化对收敛的影响【实验环境】Windows 7Microsoft Virual C++ 6.0二、实验内容：【实验方案】Gauss-Seid迭代法在Jacobi方法迭代基础上对新计算出来的第k+1次近似的分量加以了利用，即（1）( k=0,1,2,…;i=1,2,…,n) （分量形式）（2）若D-L可逆k=0,1,2,… x(0)为初始向量（矩阵形式）2.松弛迭代法(1)在Gauss-Seid迭代法中若在修正项前加一个参数，便得到松弛法的迭代公式 i=1,2,…,n.k=0,1,2……(向量形式)（2）由索取=1.251,故为超松弛因子（矩阵形式）【实验过程】（实验步骤、记录、数据、分析）1.根据算法事先写出相应的程序；2.对Gauss-Sseid迭代法和松弛迭代法的程序分别进行编译和调试；4.调试通过，计算出正确结果；5.分析结果数据。【实验结果】1.Gauss_Seidel迭代次数为12次结果：X[1]= -3.090909 X[1]=1.094545X[1]=1.2654552.初值为（0,0,0,0） 初值为（10,10,10,10） 初值为（100,100,100,100）由以上结果可以看出初值的改变对结果影响不大。【实验小结】（收获体会） 这个实验中通过输入不同初值观测结果改变来查看初值对结果的影响，可以更直观看出初值大小对结果影响不大。在调试程序时出现错误但不知道怎么修改，感觉C++很酷炫 ，不容易，要多花工夫，多练习才行呢。三、指导教师评语及成绩：评 语评语等级优良中及格不及格1.实验报告按时完成,字迹清楚,文字叙述流畅,逻辑性强2.实验方案设计合理3.实验过程（实验步骤详细,记录完整,数据合理,分析透彻）4实验结论正确. 成 绩： 指导教师签名： 批阅日期：附录：源 程 序1.#includestdio.h#includestdlib.h#includeconio.h#includemath.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],intm,int n) //输入增广矩阵{ int i, j; printf(\n输入系数矩阵：\n); for(i=1;i=m;++i) { printf(增广矩阵行数%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]); } intIsSatisfyPricision(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]; inti,j,k; printf(\n初始向量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]