Chap__NewtonCotes公式.ppt

第6章 数值积分（Numerical Integration） 数理学院 内容提纲 数值积分的必要性 求积公式及其代数精度 插值型求积公式 Newton-Cotes公式及数值稳定性 复化求积公式及误差估计  数值积分的必要性 本章主要讨论如下形式的一元函数积分 在微积分里，按Newton-Leibniz公式求定积分 要求被积函数f(x) ? 有解析表达式； ? f(x)的原函数F(x)为初等函数． 实际问题 1.2 f(x)的原函数F(x)不能用初等函数表示 例如函数: 假若要求波纹瓦长4英尺,每个波纹的 高度(从中心线)为1英寸,且每个波纹以近 似2π英寸为一个周期. 求制做一块波纹 瓦所需铝板的长度L. 这个问题就是要求由函数f(x)=sin x给定 的曲线,从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 2. 有些被积函数其原函数虽然可以用初等函数表示成有限形式,但表达式相当复杂,计算极不方便.例如函数  3.f(x)没有解析表达式，只有数表形式: 求积公式及其代数精度 求积公式的概念 积分值　　　　　　　　　 　 在几何上可解释为由x=a, x=b, y=0和 y=f(x) 所围成的曲边梯形的面积.积分计算之所以有 困难，就是因为这个曲边梯形有一条边y=f(x) 是曲的. 　　 依据积分中值定理,对于连续函数f(x) ,在 [a,b]内存在一点ξ,使得 如果简单地选取区间[a,b]的一个端点或 区间中点的高度作为平均高度,这样建立的 求积公式分别是: 左矩形公式: I(f)≈(b-a)f(a) 右矩形公式: I(f)≈(b-a)f(b) 中矩形公式: I(f)≈(b-a)f[(a+b)/2] 此外,众所周知的梯形公式: I(f)≈(b-a)[f(a)+f(b)]/2 和 Simpson公式: I(f)≈(b-a)[f(a)+4f((a+b)/2)+f(b)]/6 则分别可以看作用 a, b, c=(a+b)/2, 三点 高度的加权平均值 [f(a)+f(b)]/2 和 [f(a)+4f(c)+f(b)]/6 作为平均高度f(ξ)的近似值. 更一般地,取区间[a,b]内n+1个点 {xi},(i=0,1, 2,…n)处的高度{f(xi)} (i=0,1,…,n)通过加权平 均的方法近似地得出平均高度f(ξ),这类求积 方法称为机械求积: 或写成: 记 求积公式的代数精度 定义1 称求积公式(2)具有m次代数精度,如果它满 足如下两个条件: (i)对所有次数≤ m次的多项式 ,有 (ii)存在m+1次多项式 ,使得 插值型求积公式 在积分区间[a,b] 上取n+1个节点xi，i=0,1, 2,…,n,作f(x)的n次代数插值多项式（拉格朗 日插值公式）: 则有 为插值余项 于是有 取 称(4)式为插值型求积公式,其中求积系Ak由(5) 式确定. 推论1 求积系数满足: 记 定理2 当阶数n为偶数时, Newton-Cotes 公式(8)至少具有n+1次代数精度. Newton-Cotes公式的数值稳定性 现在讨论舍入误差对计算结果产生的影响.设 用公式 近似计算积分 时,其中计算函数值f(xj)有误差εj (j=0,1,2,…,n).设 计算Cj(n)没有误差,中间计算过程中的舍入误差也 不考虑,则在式(10 )的计算中,由εj引起的误差为 收敛性 由上述的误差估计式可知，当 f(x) ?C2[a,b] 时,只要 h?0时 数列Tn(f) ? I(f),且收敛速度为二阶O(h2). 误差估计 每个子区间上的误差估计式为 例 对于函数 , 应用复化Simpson法计算,得 复化Simpson公式积分法 将n个子区间的误差相加得 由闭区间上连续函数的介值性质可知在[a,b]上至少