第7章数值微积分1.ppt

【定理2】 若求积公式（7.3）的系数 对任意 ，若取 ，且 则有 证毕. 则它是收敛的、稳定的. 证明 只证明公式（7.3）的稳定性. 【注】 若牛顿-柯特斯求积公式中系数 则它是收敛的、稳定的. 的牛顿-柯特斯求积公式是不稳定的. 小结： 牛顿-柯特斯公式 其中 称为柯特斯系数， 当 相应的求积公式称为梯形公式 当 时相应的求积公式称为辛普森公式 时称为柯特斯公式 当 其中， 截断误差 牛顿-柯特斯公式的截断误差（即余项）为 当 时，梯形公式的截断误差为 如果 连续， 辛普森公式的截断误差为 柯特斯公式的截断误差为 * * 第7章 数值微积分 基础教学部数学教研室 彭 晓 华 立体化教学资源系列——数值分析 波形屋顶平材的长度：一个波形屋 顶是通过将一张平的铝材料压成横 断面具有正弦波形式的材料而构造 出来的.现在需要一个48英寸长的 波形屋顶，每个波的高度均离开中 心线1英寸，每个波的周期大约为 英寸.求原来平材的长度. （从 英寸到 从而这个问题归结为求数值积分问题. 7.1数值积分与数值微分 确定此曲线的长度.根据微积分理论，此长度为 引言： 英寸）， 问题为:给定 年份 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 人口 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 人口相对增长率：已知20世纪美国人口统计数据如下表， 为了计算表中这些年份的人口相对增长率，记时刻 的 ，则人口相对增长率为 它表示每年人口增长的比例.从而这个问题归结为求数值 微分问题. 表7-1 20世纪美国人口统计数据（ ） 人口为 在许多实际工程中，直接或间接地涉及计算导数和计算定积分.在这些微分积分的计算过程中存在如下一些问题： 1、牛顿—莱布尼兹公式 大量的被积函数找不到用初等函数表示的原函数； 例如，积分 2、当 微积分理论无法精确求积分，也无法精确求导数. 是由测量或数值计算给出的一张数据表时， * 对于计算导数的问题，从导数定义想到用差商近似的算 法，即 虽然这种近似计算的精度较差，但是它启示我们可以用两个 点上的函数值来近似求导，如果用有限个点上的函数值，能 否建立一个求导公式并估计误差呢？这是数值微分研究的问 题. 需要研究的问题： 从几何方面看，公式（7.1）表示以区间 的长度为底而高为 恰等于 在 边梯形的面积，见图7-2.问题在于 一般是不知道的，因而难以准确算出 的值. 区间 只要对平均高度 数值求积方法.如果近似地取 . 上的积分值，即曲 的具体位置 称为 提供一种算法，相应地便获得一种 对于计算定积分的问题，根据积分中值定理：存在点 使 . （ 7.1） 图7-2 矩形公式几何意义 的矩形面积， 上函数的平均高度， 如果近似地取 则由公式（7.1）得梯形公式，几何 意义见图7-3. 图7-3 梯形公式几何意义 则公式（7.1）称为中矩形公式 如果取 上有限个节点 的函数值 的加权平 ，则公式（7.1）称为机械求积公式 （7.2） 称为求积节点， 称为求积系数. 仅仅与节点 的选取有关，而不 的具体形式. 均值为 其中 依赖于被积函数 问题: 的位置并确定求积系数 才能使求得的积分值具有预先给定的任意的精确度？ (2)怎样估计数值计算的误差？ (1)如何安排求积节点 寻找便于数值计算，又能满足精度要求的微积分公式和方法. 数值积分与数值微分的基本内容： 复习：拉格朗日插值多项式 一、满足插值条件 Pn(xi)=f(xi), ( i=0,1,2,…,n) n次插值多项式Pn(x)=a0+a1x+a2x2+……+anxn 存在而且惟一。 二、Lagrange插值多项式： 称为Lagrange插值基函数。 三、插值余项: Rn (x)= f (x) - Pn (x)= 7.2 牛顿-柯特斯求积公式 上取定 个点 ，经过这些点作插值多项式，用插值多项式 ，进而确定 如果利用Lagrange插值多项式 型求积公式 其中 插值型求积公式的基本思想：在 7.2.1 牛顿-柯特斯求积公式 代替被积函数 ，则有插值 (7.3) 等距节点情形:将区间 划分为 等分，步长 ，选取等距节点 ，此时， 即记 将其代入公式（7.3），得到牛顿-柯特斯公式 其中 称为柯特斯系数， . （7.5） ， （7.4） 对不同的 ,柯特斯系数可按公式（7.5）计算. 当 时（有两个等距求积节点）， 相应的求积公式就是梯形公式 当 时（有三个等距求积节点