数值分析实验(4).docVIP

实验四数值积分与数值微分 专业班级：信计131班 姓名：段雨博 学号：2013014907 、实验目的 1、 熟悉matlab编程。 2、 学习数值积分程序设计算法。 3、 通过上机进一步领悟用复合梯形、复合辛普森公式，以及用龙贝格求积方法计算积分的 原理。 二、实验题目 P137 1、用不同数值方法计算积分 「x In xdx - - ￡。 （1） 取不同的步长h ?分别用复合梯形及复合辛普森求积计算积分，给出误差中关于 h的函 数，并与积分精确值比较两个公式的精度， 是否存在一个最小的 h，使得精度不能再被改善？ （2） 用龙贝格求积计算完成问题（1 ）。 三、实验原理与理论基础 1.1复合梯形公式及其复合辛普森求解 h n 丄 h nd 「「了 〔f （xQ f（x「）匸；f a 2 f Xk f b 2 k z0 2 k j 误差关于h的函数：Rn f - _b a h2 f i I 12 h 一 nJ nJ 复合辛普森公式：S^ - f a 4^ f xk .1/2 f xk f b 6 - k=0 kT 误差关于h的函数: 误差关于h的函数: -Sn b _ a h 180 2 1.2龙贝格求积算法：龙贝格求积公式是梯形法的递推化, 和柯特斯公式之间的关系的基础上, 量的前提下提高误差的精度的特点。 计算过程如下： 1.2龙贝格求积算法： 龙贝格求积公式是梯形法的递