数值分析课程数值积分实验报告.docVIP

数值分析课程—数值积分实验报告 实验名称 数值积分 班级 学号 姓名 实验地点 实验时间 任课教师 评分 实验目的 体会数值积分的基本概念； 掌握低阶的插值型数值积分公式； 掌握区间逐次分半的复化求积方法； 掌握龙贝格算法的基本思路和迭代步骤； 二、用文字或图表记录实验过程和结果 计算积分： 精度要求。 编写C++程序，分别用梯形公式、Simpson公式和Cotes公式计算上述积分的近似值。并对计算结果作一比较。 答： 梯形公式 Simpson公式 Cotes公式 0.68394 0.74718 0.746834 0.920735 0.946146 0.946083 分析：随着结点数目的增加，积分公式的代数精度逐渐增高，截断误差减小，得到的结果也越加接近精确值。 编写C++程序，分别用区间逐次分半的复化梯形公式和区间逐次分半的复化Simpson公式计算上述积分的近似值，比较它们的迭代次数。 答: 精度为1e-6 积分函数 复化梯形公式 复化Simpson公式 0.746823 （迭代次数）8 0.746825 （迭代次数）22 0.946082 7 0.946084 22 分析：在相同精度条件下，由于迭代次数越少，收敛速度越快，所以复化Simpson公式的收敛速度比复化梯形公式的收敛速度慢。 编写龙贝格算法的C++程序，并用以计算上述积分，与（2）比较迭代次数。 注：迭代次数是收敛快慢的指标之一 积分函数 龙贝格 0.746824 4 0.946083 4 分析：在相同精度要求下，龙贝格公式的迭代次数比复化梯形公式的迭代次数少，比复化Simpson公式的迭代次数少，因此龙贝格公式的收敛速度比复化梯形公式的收敛速度快，比复化Simpson公式的收敛速度快。 三、练习与思考题分析解答 练习思考 为什么要引进复化求积的方法？ 答：应用高阶的newton-cotes型求积公式计算积分会出现数值不稳定，而低阶公式又因为积分区间步长过大使得离散误差很大，然而，若积分区间愈小，则离散误差愈小，因此，为提高求积公式的精确度，可以将积分区间分成若干个子区间，在每个子区间上使用低阶求积公式，然后将计算的结果加起来，据成了复化求积公式。 简述为什么引进区间逐次分半的思想，什么是区间逐次分半法（即自适应求积步长法）？ 答：在应用复合求积公式时，为了保证计算结果的精度，往往需要确定步长取多大，我们希望计算机自动选择积分步长，计算出满足精度要求的积分近似值。我们将积分区间逐次分半，使用复合求积公式计算出积分近似值后随时比较相邻两次结果，若二者之差小于所允许的误差界限，则最后计算结果作为积分近似值。即为区间逐次分半法。 四、本次实验的重点难点分析 1、本次实验的重点是理解如何将积分公式转化为数值计算，并且熟练掌握基本的梯形，simpson和龙贝格公式及他们的复化形式，分别区分开来他们各自的作用和优缺点，便于以后更方便的使用。 2、难点在于对复化梯形、复化simpson、龙贝格基本思路的熟悉掌握以及熟悉它们迭代过程，能够通过算法，将其转化为C++格式，然后在计算机上实现计算过程。 3、在编写过程中发现，精度控制是很重要的，对于不同精度下的计算，结果会出现偏差，因此要选择合适的精度，才能得到和精确解比较接近的结果。另外，在比较收敛速度时也要注意所比较的对象是不是在相同的精度下，才能达到比较的效果，换句话说，在学习上，也要有“细节决定成败”的思维。 4、在输出迭代次数时，要注意迭代次数与那个变量对应，在变量多的情况喜爱更要分清楚。 4