[高等教育]东北大学数值分析实验报告.docVIP

数值分析实验报告 课题一 迭代格式的比较 问题提出 设方程f(x)=x- 3x –1=0 有三个实根 x=1.8793 , x=-0.34727 ,x=-1.53209现采用下面三种不同计算格式，求 f(x)=0的根 x 或x x = x = x = 二、要求 1、编制一个程序进行运算，最后打印出每种迭代格式的敛散情况； 2、用事后误差估计来控制迭代次数，并且打印出迭代的次数； 3、初始值的选取对迭代收敛有何影响； 4、分析迭代收敛和发散的原因。 三、目的和意义 1、通过实验进一步了解方程求根的算法； 2、认识选择计算格式的重要性； 3、掌握迭代算法和精度控制； 4、明确迭代收敛性与初值选取的关系。 程序代码： #includeiostream #includemath.h int k=1,k1=1,k2=1,k3=1; float x1,x2,x3; float one(float x0) { x1=(3*x0+1)/(x0*x0); return(x1); } float two(float x0) { x2=(pow(x0,3)-1)/3; return(x2); } float three(float x0) { x3=pow(3*x0+1,0.33333); return(x3); } main() { float x,x0; printf(输入x0=); scanf(%f,x); x0=x; x1=one(x0); printf(第一个公式迭代结果: \n); while(fabs(x0-x1)1e-5) { printf(x1=%6.5f\n,x1); x0=x1; x1=one(x0); k1++; } printf(x1=%6.5f \n,x1); printf(k1=%i\n,k1); x0=x; x2=two(x0); printf(第二个公式迭代结果: \n); while(fabs(x0-x2)1e-5) { printf(x2=%6.5f\n,x2); x0=x2; x2=two(x0); k2++; } printf(k2=%i\n,k2); x0=x; x3=three(x0); printf(第三个公式迭代结果: \n); while(fabs(x0-x3)1e-5) { printf(x3=%6.5f\n,x3); x0=x3; x3=three(x0); k3++; } printf(x3=%6.5f\n,x3); printf(k3=%i\n,k3); scanf(%); } 实验结果： 四、程序运行结果讨论和分析： 对于第一种迭代格式，收敛区间[-8.2 -0.4]，在该收敛区间内迭代收敛于-1.53209，只能求得方程的一个根； 对于第二种迭代格式，收敛区间[-1.5 1.8]，在该收敛区间内迭代收敛于-0.34730，同样只能求得方程的一个根； 对于第三种迭代格式，收敛区间[-0.3 +∞)，在该收敛区间内迭代收敛于1.87937，只能求得方程的一个根； 由以上结果很容易发现，初值的选取对迭代敛散性有很大影响。以第一种迭代格式为例，当初值大于等于-0.3时，迭代格式发散；当初值小于等于-8.3时，迭代格式也发散；只有初值在-0.3和-8.3之间时，迭代格式才收敛于—1.53209。其他迭代格式也有这样的性质，即收敛于某个数值区间，超出这个区间迭代格式就是发散的，这就是所谓迭代格式的收敛性。 对于不同迭代格式在不同区间具有不同的敛散性的原因，我认为可以从一下两方面理解：1、迭代法是一种逐次逼近法，其基本思想是将隐式方程归结为一组显式的计算公式，就是说，迭代过程实质上是个逐步显式化的过程。2、我们可以用几何图像来更好地理解迭代过程。由图可知，在某些区间选取的初始值随着迭代次数的增加会越来越逼近精确值，即收敛于精确值，而在另外一些区间选取的初始值随着迭代次数的增加却离精确值越来越远，即不会收敛于一个确定值。 课题二 线性方程组的直接算法 一、问题提出 给出下列几个不同类型的线性方程组，请用适当算法计算其解。 1、设线性方程组 = x= ( 1, -1, 0, 1, 2, 0, 3, 1, -1, 2 ) 2、设对称正定阵系数阵线方程组 = x = ( 1, -1, 0, 2, 1, -1, 0, 2 ) 3、三对角形线性方程组 = x= ( 2, 1, -3, 0, 1, -2, 3, 0, 1, -1 ) 二、要求 对上述三个方程组分别利用Gauss顺序消去法与Gauss列主元消去法；平方根法