Nton-ots公式.pptVIP

(1)复化梯形公式： 在每个 上用梯形公式： = Tn /*中值定理*/ 复化梯形公式积分法 (2)收敛性 由上述的误差估计式可知,当f(x) ?C2[a,b] 时,只 要h?0时 数列Tn(f) ? I(f),且收敛速度为二阶 O(h2). 但是f(x)?C2[a,b]条件相对苛刻, 现假定f(x)在[a,b]上Riemann可积,讨论复化求积公式的收敛性 (3)复化 Simpson 公式 复化Simpson公式积分法 每个子区间上的误差估计式为 将n个子区间的误差相加得 由闭区间上连续函数的介值性质可知在[a,b]上至少存在一点?，使 可见,当f(x)有四阶导数时,复化Simpson公式具有4阶收敛. (4)误差估计 将区间[0,1]划分为8等分,应用复化梯形法求得 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.9456909 例1 试用数据表计算积分 对于函数 解 应用复化Simpson法计算,得 比较上面两个结果T8和S4,它们都需要提供9个点上的函数值工作量基本相同,然而精度却差别很 大. 同积分的准确值I(f)=0.9460831比较,复化梯形法的结果T8=0.9456909只有两位有效数字, 而复化 Simpson法的结果S4=0.9460832却有六位有效数字. =0.9456909 第四章 数值积分(Numerical Integration)  ? f(x)的原函数F(x)为初等函数． 本章主要讨论如下形式的一元函数积分 在微积分里，按Newton-Leibniz公式求定积分 要求被积函数f(x) ? 有解析表达式； 一、数值积分的必要性 §1 引言 /*Introduction*/ 例如函数: 考虑一个实际问题: 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的. 1、f(x)的原函数F(x)不能用初等函数表示 假若要求波纹瓦长4英尺,每个波纹的高度(从中心 线)为1英寸,且每个波纹以近似2π英寸为一个周 期. 求制做一块波纹瓦所需铝板的长度L. 这个问题就是要求由函数f(x)=sin x给定的曲 线,从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 上述积分称为第二类椭圆积分,它不能用普通方法来计算. 例如函数 并不复杂,但它的原函数却十分复杂: 2. 有些被积函数其原函数虽然可以用初等函数表示成有限形式,但表达式相当复杂,计算极不方便.  3. f(x)没有解析表达式，只有数表形式: x 1 2 3 4 5 f(x) 4 4.5 6 8 8.5 这些都说明,通过原函数来计算积分有它的局限性,因而,研究关于积分的数值方法具有很重要的实际意义. 在几何上可解释为由x=a, x=b, y=0和 y=f(x) 所围 成的曲边梯形的面积.积分计算之所以有困难， 就是因为这个曲边梯形有一条边y=f(x)是曲的. (1)求积公式的概念 积分值 二、求积公式及其代数精度 依据积分中值定理,对于连续函数f(x) ,在[a,b]内 存在一点ξ,使得 称f(ξ)为区间[a,b]的平均高度. 问题在于点ξ的具体位置一般是不知道的.这样,只要对平均高度f(ξ)提供一种算法,相应地便获得一种数值求积 方法. 左矩形公式: I(f)≈(b-a)f(a) 右矩形公式: I(f)≈(b-a)f(b) 中矩形公式: I(f)≈(b-a)f[(a+b)/2] 如果简单地选取区间[a,b]的一个端点或区间中 点的高度作为平均高度,这样建立的求积公式分 别是: 作为平均高度f(ξ)的近似值. 此外,众所周知的梯形公式: 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 更一般地,取区间[a,b]内n+1个点 {xi},(i=0,1,2,…n) 处的高度{f(xi)} (i=0,1,…,n)通过加权平均的方法 近似地得出平均高度f(ξ),这类求积方法称为 机械求积: 或写成: 数值积分公式 求积系数 求积节点 记 称为求积公式余项(误差). 构造或确