柯特斯公式复化梯形公式积分法2牛顿.ppt

定义 称某求积公式具有m次代数精度,如果它满 足如下两个条件: (1)对所有次数≤ m次的多项式 ,有 (2)存在m+1次多项式 ,使 定理 n为偶数时, Newton-Cotes公式至少具有n+1次代数精度 误差估计 在每个子区间上 例 对于函数 , 用复化Simpson法得 例 设 带余项的差分公式为 重复使用外推法，得到误差为 的外推公式 当j=2时 当K=2时,有 这是n=2的复化Simpson公式的 . 计算精确度的要求来确定区间逐次分半的次数，常用不等式 作为达到精度要求的判断准则。?是给定的一个小的正数. 例2 用Romberg求积方法计算 的近似值,给定?=0.001 解 首先令区间长度 h1=1 ,用梯形求积公式计算 其误差为 . * 数值积分是否必要？ 本章主要问题 Newton-Leibniz公式 f(x)应满足 有解析表达式； f(x)的原函数F(x)为初等函数． 的数值解法 例 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的. 若要求波纹瓦长4英尺,每个波纹的高度(从中心线)为1英寸,且每个波纹以近似2π英寸为一个周期. 求制做一块波纹瓦所需铝板的长度L. 解 f(x)=sin x从x=0到x=48英寸间的弧长即为L 该积分称为第二类椭圆积分,不能用普通方法来计算 积分值 的几何意义　　　　　　　　　 1-1 数值积分的基本思想 积分中值定理 ξ是[a,b]内一点 提供一种平均高度 f(ξ) 算法,相应获得一种数值求积方法 梯形公式 I(f)≈(b-a)[f(a)+f(b)]/2， f(ξ) =[f(a)+f(b)]/2 中矩形公式 I(f)≈(b-a)f[(a+b)/2]， f(ξ) = f[(a+b)/2] 取区间[a,b]内n+1个点 {xi},(i=0,1,2,…n)处的高度{f(xi)} (i=0,1,…,n)的加权平均作为平均高度 f(ξ） 机械求积方法 作为 求积节点 数值积分公式 求积系数 称为求积公式余项 1-1 代数精度的概念 插值型求积公式 在[a,b] 上取n+1个节点xi，i=0,1,2,…,n,作f(x)的n次插值多项式 则有 插值余项为 于是 Ak 插值型 误差 定理 形如 的求积公式至少有 n 次代数精度 的充分必要条件是 该公式为插值型 即 取等距节点 令 Cotes系数 注：Cotes 系数仅取决于 n 和 k，与 f (x) 及区间[a, b]均无关。 构造的插值型求积公式称为Newton-Cotes公式, 求积系数为 2-1 柯特斯系数 证明 只需验证当n为偶数时,Newton-Cotes公式对 f(x)=xn+1 的余项为零. 由 f(x)=xn+1 得 f(n+1)(x)=(n+1)! 误差为 引进变换t=u+n/2,因为n为偶数,故n/2为整数, 于是有 R(f)=0 2-2 偶阶求积公式的代数精度 2-3 几种低阶求积公式的余项 n = 1: 梯形公式 n = 2: 代数精度= 1 辛甫生公式 代数精度= 3 复化梯形公式 在每个 上用梯形公式 = Tn （中值定理） 2-4 复化求积法及其收敛性 则 复化梯形公式积分法 收敛性 由误差估计式可知，当 f(x) ?C2[a,b] 时,只要 h?0，则 Tn(f) ? I(f)，且收敛速度为O(h2). 当f(x) ?C2[a,b] 时条件相对苛刻 现假定f(x)在[a,b]上Riemann可积 复化 Simpson 公式 复化Simpson公式积分法 相加得误差 由介值性质可知 当f(x)有四阶导数时,复化Simpson公式具有4阶收敛 试用数据表计算积分 x f (x) 0 1 1/8 0.9973978 2/8 0.9896158 3/8 0.9767267 4/8 0.9588510 5/8 0.9361556 6/8 0.9088516 7/8 0.8771925 1 0.8414709 解 将[0,1]分为8等分，用复化梯形法得 准确值 I(