计算方法数值积分讲解.pptVIP

化工计算：数值积分  或写成: * * 一、数值积分的必要性 主要讨论如下形式的一元函数积分 在微积分里，按Newton-Leibniz公式求定积分 要求被积函数 ? 有解析表达式； ? 的原函数 为初等函数． 实际问题 1. 的原函数 不能用初等函数表示 例如函数: 2. 有些被积函数其原函数虽然可以用初等函数表示成有限 形式,但表达式相当复杂,计算极不方便. 例如函数: 并不复杂,但它的原函数却十分复杂: 3. 没有解析表达式，只有数表形式: 1 4 2 3 4 5 4.5 6 8 8.5 原来通过原函数来计算积分有它的局限性。那…… 怎么办呢？ 呵呵…这就需要积分的数值方法来帮忙啦。 二、数值积分的基本思想 1、定积分的几何意义 2、数值积分的理论依据 依据积分中值定理, 对于连续函数 , 在 内存在一点 ,使得 称 为区间 的平均高度. 3、求积公式的构造 ? 若简单选取区间端点或中点的函数值作为平均高度，则可得一点求积公式如下： 左矩形公式： 中矩形公式： 右矩形公式： 左矩形公式： 中矩形公式： 右矩形公式： ? 若取 两点，并令 ，则可得梯形公式（两点求积公式） 则可得Simpson公式(三点求积公式) ? 若取三点， 并令 ? 一般地 ，取区间 内 个点 处的高度 通过加权平均的方法近似地得出平均高度 这类求积方法称为机械求积: 数值积分公式 求积系数 求积节点 记 称为数值 求积公式 称为求积公式余项(误差). 三、几种常用的低阶求积公式 n = 1: 梯形公式 /* 令 x = a+th, h = b?a, 用中值定理 */ 代数精度 = 1 n = 2: Simpson 公式 代数精度 = 3 四、复化求积公式 高次插值有Runge 现象，怎么办？ 可采用分段低次插值来解决 高阶Newton-Cotes公式会出现 数值不稳定。 而低阶Newton-Cotes公式 有时又不能满足精度要求，怎么办？ 可将积分区间 分成若干小 区间，在每个小区间上用 低阶求积公式计算，然后求和。 ? 复化梯形公式： 在每个 上用梯形公式： = Tn /*中值定理*/ *