第3章积分数值方法.pptVIP

* 第3章积分的数值方法 计算定积分有微积分基本公式： 但很多函数找不到原函数，如 等。而实际上，有很多函数只知一些离散点的函数值，并无表达式，这就需要利用已知条件求出近似值。 3.1概述 3.2 梯形积分法 （简介） 3.2 抛物线积分法 （简介） 的被积函数?(x)在节点a?x0x1…xn?b上的函数值 yk= ?(xk) ,k=0,1,2,…,n 所以有 设?(x)在[a,b]具有n+1阶连续导数,则有 则可以构造n次Lagrange插值多项式Ln(x),因为?(x)=Ln(x)+Rn(x) 3.4 龙贝格积分法 若已知定积分 3.4.1 牛顿－科特斯积分 求积公式(3－10)称为求积公式的一般形式. 称之为插值型求积公式. 为了简化计算,取等距节点xk=a+kh,(k=0,1,2,…,n, 令 其中 则有 称式(3-13)为Newton-Cotes公式.Ck(n)称为Cotes系数. 例1 设?(x)?C2[a,b],求n=1时的Newton-Cotes公式并估计误差. 解 计算Cotes系数 于是有 从几何上看: 所以公式 ,则有误差估计 若记 o x y a b 也称为梯形公式,记为T. =T 例2 设?(x)?C4[a,b],求n=2时的Newton-Cotes公式并估计误差. 解 计算Cotes系数 y=?(x) 称之为Simpson公式或抛物线公式，记为S. 于是有 =S. 若记 则有 由于构造Newton-Cotes公式需要Cotes系数,将其列表如下: 1 2 3 4 5 6 7 8 8 7 6 5 4 3 2 1 0 k n 例3 求n=4的Newton-Cotes公式及误差. 解 查表可得 于是有 称之为Cotes公式，记为C。其误差为 其中, xk=a+kh , k=0,1,2,3,4 , h=(b-a)/4 . 例4 用梯形公式、Simpson公式和Cotes公式求积分 的近似值。 解 I?T=1/2(4+2)=3 I?S=1/6(4+12.8+2)=3.13333 I?C=1/90(28+…+14)=3.14212 在区间[a，b]上,取等距节点 若在每个小区间[xk,xk+1]用梯形公式，则有 称为复化梯形公式，记为 3.4.2 复化求积公式 xk=a+kh , k=0,1,2,…,n , h=(b-a)/n 由定积分的区间可加性,有 可见,复化梯形公式是收敛的。而且,要使|I-Tn|? ,只要 如果记M2= 复化梯形公式的误差为 若将New-Cotes公式中,每个小区间上的积分采用Simpson公式,则可得到复化Simpson公式: ，则有 其中， ，而且误差为 如果记M4= ，则有 复化Simpson公式也是收敛的,而且,要使|I-Sn|? ,只要 其中, 误差为 类似地可得复化Cotes公式: 复化Cotes公式也是收敛的,而且,要使|I-Cn|? ,只要 已知函数 分别用复化梯形公式、复化Simpson公式计算积分 例5 的数据表 0.8414710 1 0.9361556 5/8 0.9896158 1/4 0.8771926 7/8 0.9588511 1/2 0.9973978 1/8 0.9088517 3/4 0.9767267 3/8 1 0 ?(xk) xk ?(xk) xk ?(xk) xk 的近似值。 定积分I精确到小数点后7位的值是0.9460831。 解 由此得 由于 一方面，若|T2n-Tn|3? ,则有近似误差|I-T2n|? . 所以有 另一方面，(4T2n-Tn)/3应比Tn和T2n的近似程度更好.事实上,有 3.4 .3 龙贝格积分法 其中,xk=a+kh,k=0,1,2,…,n,h=(b-a)/n, 而且有 于是有 因此有逐次分半的复化梯形公式的递推公式：