第四章数值积分介绍.pptVIP

第四章 数值积分 §4.1 梯形求积公式、Simpson求积公式和Newton-Cotes公式 §4.2 复化求积公式 §4.3 自动选取积分步长 §4.4 求积公式的误差 §4.5 龙贝(Romberg)方法 引 言 工程上经常用到积分运算，但用牛顿－莱布尼兹公式计算定积分时常遇到如下情况： ※常用插值多项式来构造数值求积公式 §4.1 梯形求积公式、Simpson求积公式和Newton-Cotes公式 二、辛卜生求积公式： §4.2 复化求积公式 二、复化梯形求积公式： 三、复化辛卜生公式 四、复化牛顿－柯特斯公式： ［例］计算： §4.3 自动选取积分步长 复化求积方法对提高精度是行之有效的，但是使用复化求积公式之前必须给出合适的步长，步长取的太大精度难以保证，步长太小则会导致计算量的增加。而事先给出一个恰当的步长是困难的。 实际计算时通常采用变步长的求积方案，即在步长逐次折半（或称步长二分）的过程中，反复利用复化的求积公式进行计算，直到二分前后两次积分近似值相当符合为止。 下面以梯形公式为例，探求变步长的计算规律。 ※ 算法设计与流程图： 算法： №1 输入初值； №2 按公式（4.9）求二分后的积分值，从第一个分点出发取h为步长，逐步向右跨，即可依次确定（4.9）中的各个分点，并将所得分点暂存于单元x中； №3 控制精度； №4 修改步长。 流程图：（略） 练习：分别用复化梯形法和复化辛甫生法计算下列积分： （1） ，4等分积分区间； （2） ，8等分积分区间 ； （3） ，6等分积分区间。 §4.4 求积公式的误差 一、梯形公式： 2、复化梯形求积公式 (4.5)截断误差： 二、辛卜生公式： 1、辛卜生公式(4.3)截断误差： 1、牛顿－柯特斯公式(4.4)截断误差： §4.5 龙贝(Romberg)方法 复化的柯特斯公式，结构简单，易于上机实现，但收敛速度慢，为节约计算工作量，产生龙贝方法。 Romberg算法是在积分区间逐次分半的过程中，对用复化梯形法产生的近似值进行加权平均，以获得准确程度较高的近似值的一种方法。 优点：公式简练，使用方便，结果可靠。 一、基本原理： ［例］计算： （§4.2 例题） 解法二：利用公式（ 4.19） 二、龙贝求积公式： ［例］利用龙贝算法计算： ——(4.21) ——(4.22) ——(4.23) ——(4.24) ※综上，可以在积分区间逐次分半的过程中，利用上述公式，将粗糙的近似值逐步精确，得到 这种算法称为龙贝算法。 3.141593 3.141592 3.141592 3.140941 4 3.141586 3.141594 3.141593 3.138988 3 3.142118 3.141569 3.131176 2 3.133333 3.1 1 3 0 K ［解］ * * 因此要求建立积分的近似计算方法。 这样建立的求积公式称作插值求积公式 ——(4.1) 并计算 作为积分值的近似值，得 ——(4.2) 即梯形求积公式 一、梯形公式： 从而求得求积公式 (4.3)称为辛卜生(Simpson)公式（或抛物线求积公式）。 ——(4.3) 三、牛顿－柯特斯公式： ——(4.4) 当积分区间[a，b]较大时，直接使用上述公式所得近似 值的精度很低，实际应用中往往采用复合求积的方法。 一、定义：将积分区间分成几个小区间，在每个小 区间上用上述公式计算积分的近似值，即为复化求积。 由梯形公式，在第i个小区间上 称为复化梯形公式 ——(4.5) ——(4.6) 则有复化牛顿－柯特斯公式： ——(4.7) 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 ※都需调用9次函数，但精度差别很大，故常用复化 辛卜生公式。 ※上机时，可以通过循环来实现公式4.6。 ——(4.8) ——(4.9) （1） ??????? （2） （3） 答案： 证明：由第二章余项公式，有： ——(4.10) 于是(4.10)成立。 ——(4.11) ——(4.12) 2、复化辛卜生公式(4.6)截断误差： ——(4.13) ——(4.14) ——(4.15) 三、牛顿－柯特斯公式： ——(4.16) 2、复化牛顿－柯特斯公式(4.7)截断误差： ——(4.17) 由上节的误差估计式 (4.11) 由公式 (4.11), 则由上两式求得： 前已说过（见§2.2），这种直接用计算结果估计 误差的方法称为事后估计法。 ——(4.18) 按估计式(4.18)， 可能是更好的结果。 ——(4.