函授 数值 课件 chap3：数值积分2(复化求积公式).pptVIP

华长生制作 * 第四章 数值积分与数值微分 4.2 复合求积法 4.2 复合求积法 直接使用Newton-Cotes公式的余项将会较大 公式的舍入误差又很难得到控制 为了提高公式的精度,又使算法简单易行,往往使用复合方法 然后在每个小区间上使用低阶Newton-Cotes公式 最后将每个小区间上的积分的近似值相加 一、复合求积公式 各节点为 记为 由积分的区间可加性,可得 复合求积公式 复合梯形公式 复合Simpson公式 复合抛物线公式 复合Cotes公式 复合梯形公式分解 复合Simpson公式分解 例1. 解: 为简单起见,依次使用8阶复合梯形公式、4阶 复合Simpson公式和2阶复合Cotes公式 可得各节点的值如右表 0 1 0.125 00.25 00.375 00.5 00.625 00.75 00.875 0 1 0 复合求积 公式的程序 newtoncotes.m 函数程序 func.m 分别由复合Trapz、Simpson、Cotes公式有 原积分的精确值为 精度最高 精度次高 精度最低 比较三个 公式的结果 那么哪个复合求积公式的收敛最快呢？ 二、复合求积公式的余项和收敛的阶 我们知道,三个求积公式的余项分别为 单纯的求积公式 复合求积公式的每个小区间 则复合梯形公式的余项为 由于 即有 又由 比较三种复合公式的的余项 为此介绍收敛阶的概念 定义1. 不难知道,复合梯形、Simpson、Cotes公式的收敛阶分别为 2阶、4阶和6阶 通常情况下,定积分的结果只要满足所要求的精度即可 三、复合求积公式步长的自动选取 复合梯形公式的余项为 因此有 即 依此类推 步长自动选取的步骤: 依此类推 以上这种方法称为自适应求积法 有时也去掉 精度会更高 以复合Simpson求积公式的特点为例 具有以下特点: 旧节点 新节点 步长折半 四、复合自适应求积法的算法设计 规定步长的求积法的算法比较简单,这里只以自动选取步 长的复合Simpson求积公式为例介绍自适应法的算法设计 (一) 算法名称 求定积分 (二) 存储方式 (三) 自然语言 华长生制作 (四) 程序实例 程序名： Autosimpson.m 命令格式：autosimpson(fun,a,b,eps) 例2. 用自适应Simpson公式计算下列定积分,并比较 Li42.m 解: m = n=8 trapz n=20 trapz n=50 trapz n=100 trapz n=4 simpson n=20 simpson n=4 cotes autosimp 1e-4 autosimp 1e-6 autosimp 1e-10 z = 0.745865614845695 0.746670836939873 0.746799607189351 0.74681800146797 0.746826120527467 0.746824136005348 0.746824133229615 0.746826120527467 0.746824140606985 0.74682413281433 t = 0 0 0 0.06 0 0 0 0.11 0.16 0.22 d = -9.6e-04 -1.5e-04 -2.5e-05 -6.1e-06 2.0e-06 3.2e-09 4.2e-10 2.0e-06 7.8e-09 1.9e-12 h =0.25 4 h =0.0625 16 h =0.0078125 128 z1= 0.746824132812427 区间数及方法 积分结果 运行时间 误差 * * * * *