4.1数值积分概论教材.pptVIP

第四章 数值积分 与数值微分 §4.1 数值积分概论 一、数值积分的必要性 本章主要讨论如下形式的一元函数积分 在微积分里，按Newton-Leibniz公式求定积分 要求被积函数 ☞ 有解析表达式； ☞ 的原函数 为初等函数．  实际问题 1. 的原函数 不能用初等函数表示 例如函数: 考虑一个实际问题: 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的. 假若要求波纹瓦长4英尺, 每个波纹的高度(从中心线)为1英寸, 且每个波纹以近似 英寸为一个周期. 求制做一块波纹瓦所需 铝板的长度L. 从 到 英寸间的弧长L. 这个问题就是要求由函数 给定的曲线, 由微积分学我们知道,所求的弧长可表示为: 上述积分称为第二类椭圆积分。 What’s the Original function?! It’s so complex that we can not get it. 2. 有些被积函数其原函数虽然可以用初等函数表示成有限 形式,但表达式相当复杂,计算极不方便. 例如函数: 并不复杂,但它的原函数却十分复杂: 3. 没有解析表达式，只有数表形式: 原来通过原函数来计算积分有它的局限性。那…… 怎么办呢？ 呵呵…这就需要积分的数值方法来帮忙啦。 数值积分----是计算定积分的具有一定精度的近似值的各种计算方法。 二、数值积分的基本思想 1、定积分的几何意义 2、数值积分的理论依据 依据积分中值定理, 对于连续函数 , 在 内存在一点 ,使得 称 为区间 的平均高度. 3、求积公式的构造  若简单选取区间端点或中点的函数值作为平均高度，则可得一点求积公式如下： 左矩形公式： 中矩形公式： 右矩形公式： 左矩形公式： 中矩形公式： 右矩形公式： 则可得Simpson公式(三点求积公式)  一般地 ，取区间 内 个点 处的高度 通过加权平均的方法近似地得出平均高度 这类求积方法称为机械求积: 或写成: 数值积分公式 求积系数 求积节点 记 上式可作为一般的求积公式，其特点是将积分问题归结为函数值的计算，从而避开了使用牛顿一莱布尼兹公式需要求原函数的困难，适合于函数给出时计算积分，也非常便于设计算法。便于上机计算。 三、求积公式的代数精度 数值积分是一种近似方法，但其中有的公式能对较多 的函数准确成立，而有的公式只对较少的函数准确成立。 为了反映数值积分公式在这方面的差别，引入代数精度的概念。 1、代数精度 定义：如果对于所有次数不超过 m 的多项式 f (x) ，公式 精确成立，但对某个次数为 m +1 的多项式不精确成立，则称该求积公式具有 m 次代数精度 例：试确定 Ai ，使得下面的求积公式具有尽可能高的代数精度 具有至少 n 阶代数精度 2、举 例 例：试确定系数 Ai ，使得下面的求积公式具有尽可能高的代数精度，并求出此求积公式的代数精度。 易验证该公式对 f (x)＝x3 也精确成立，但对 f (x)＝x4 不精确成立，所以此求积公式具有 3 次代数精度。 例：(教材100页) 试确定下面求积公式中的系数，使其具有尽可能高的代数精度。 将 f (x)＝x3 代入，等号成立，故公式具有 2 次代数精度。 容易验证： 左矩形公式 和 右矩形公式 具有 零次 代数精度 中矩形公式 和 梯形公式 具有 一次 代数精度 特别地，任意具有 m ( 0 ) 次代数精度的求积公式一定满足： 我们也可以从下面的角度考虑构造数值积分公式： 四、插值型求积公式 插值型求积公式 设求积节点为：a  x0 x1 ··· xn  b 若 f (xi) 已知，则可做 n 次多项式插值： 插值型求积公式 插值型求积公式 当 f (x)＝ 1, x, x2, … , xn 时，有 即公式精确成立 证:必要性 设n+1个节点的求积公式 为插值型求积公式,求积系数为 又 当f(x)为不高于n次的多项式时, f(x)=P(x),其余项R(f)=0。因而这时求积公式至少 具有n次代数精度。 充分性： 若求积公式至少具有n次代数精度, 则对n次多项式 精确成立,即 而 取 时 所以有 ,即求积公式为插值型求积公式 (1) 在积分区间[a,b]上选取节点xk (2) 求出f(xk)及利用