第二章.数值积分.pptVIP

2. Simpson公式 例8 依次用n=8的复化梯形公式、n=4的复化 Simpson公式计算定积分 解:首先计算出所需各节点的函数值,n=8时， 由复化梯形公式得 由复化Simpson公式得 （积分准确值I=0.9460831） 例9 用复化梯形公式、Simpson公式计算定积分 已知 2 2.26549 2.56 2.87640 3.2 3.50685 3.76471 3.93846 4 1 0.875 0.75 0.625 0.5 0.375 0.25 0.125 0 解：复化梯形公式n＝8， 复化Simpson公式n=4, 二、变步长的梯形法 -------------变步长梯形公式 思考题：导出二分前后的 与 的关系 变步长梯形法的流程图 例10 用变步长梯形求积法计算定积分 ， ． 解: 由 得 由于 ,故有 新分点上的函数值 二分一次 二分二次 第五节 Romberg算法 一、梯形法的加速 则有 是否存在？ 复化梯形公式 故 二、 法的再加速 其中 则有 三、Cotes法的进一步加速 其中 －－－－－－ Romberg公式 只有7阶的精度 四、 Romberg算法的实现 用梯形公式计算积分近似值 按变步长梯形公式计算积分近似值 将区间逐次分半,令区间长度 计算 ③ 按加速公式求加速值 梯形加速公式： 辛卜生加速公式： 龙贝格求积公式： ④ 精度控制 （其中ε为允许的误差） 表示二分前后的差 例11 用Romberg算法计算定积分 要求误差不超过 ． 解:由题意 3.141593 3.141593 3.141593 3.140942 3.141586 3.141595 3.141593 3.138988 3.142118 3.141569 3.131176 3.133333 3.1 3 解：构造Romberg算法表 ?(1) 被积函数f(x)并不一定能够找到用初等函数表示的原函数F(x)，例如： (2) 还有被积函数f(x)的原函数能用初等函数表示，但表达式太复杂，例如函数 其原函数F(x)为： 第二章 数值积分 (3) 被积函数f(x)没有具体的解析表达式, 其函数 关系由表格或图形表示． 一、求积方法的历史变迁 第一节 机械求积 二、机械求积的概念 如 －－－－－－ 左矩形公式 －－－－－－ 右矩形公式 平均高度 三个求积分公式 ① y=f(x) y x a b a b y x (a+b)/2 ② y=f(x) a b a b ----------------梯形公式 －－－－－－中矩形公式 y=f(x) y a b ③ Simpson公式： a b －－－－－－－－ 机械求积公式 求积系数 权系数 求积节点 1．求积系数及节点如何确定？ 2．公式的计算精度如何判断？提高计算精度？ 3．此公式与Lagrange插值多项式有何关系？ 例1 设积分区间[a, b]为[0, 2]，分别用梯形和Simpson公式 解:梯形公式和Simpson公式的计算结果与准确值比较如下表所示 f(x) 1 x x2 x3 x4 ex 准确值 2 2 2.67 4 6.40 6.389 梯形公式 2 2 4 8 16 8.389 辛卜生公式 2 2 2.67 4 6.67 6.421 计算 时积分结果并与准确值进行比较. 定义 设求积公式 对于一切次数小于等于m的多项式是准确的，而对于次数为m+1的多项式是不准确的，则称该求积公式具有m阶代数精度（简称精度） 三、代数精度的概念 等价定义 设求积公式 对于 是准确的，而对于