5数值积分与数值微分.docVIP

第5章数值积分与数值微分方法 1 基本概念 梯形公式 b f(x)dx a 中矩形公式 n A i f (x」 i 二 0 则上式为一个数值求积公式. Ai称为求积系数，X称为求积节点；而称 n 瓦 Af(x) i=0 为求积余项或求积公式的截断误差. 从定义可以看到，数值求积公式依赖于求积节点 个数n、求积节点｛人｝和求积系数｛A｝，这三个量有一 个发生变化，则产生不同的求积公式. 定义1若求积公式对于次数不超过 m 的多项式准确成立，而对于m 1次多项式 不准确成立，则称该求积公式具有 m次代 数精度为. 一般，一个求积公式的代数精度越大， 则该求积公式越好. 确定代数精度的方法 依次取f(x)=xk (k= 0,1,111)代入公式  解 本题要先求出具体的求积公式，然后 再判断所求公式的代数精度. 公式有3个待定参数，故利用3个条件 得到的3个等式关系就可以解决求出具体求 积公式的问题. 依次取f (x)二 z 1,x,x2代入求积公式并取等号，有 A+B+C=4h A-2C = 0 、A+4C = 16h/3 解之得 八 16“ 4h 8h A = ——h,B =——,C =—— 9 3 9 故所求的求积公式为 2h 2hf(x)dx —2h 詈 f(-h) +号 f(0) +詈 f(2h) 为确定其代数精度，再取f(x^ x3代入求出的公式继 16 续计算，有R(x3)「§h4 = 0，故所求的求积公式具有 二次代数精度• 插值型求积公式 考虑f(x)关于n 1个节点 a兰沧 x Ilk x/ b的Lagrange插值多项式 Ln(x)与f(x)的余项，有  b A 厂 l i ( x)dX a 则有 b n b f x dx = 2 Aj f （X） a i=0 a （n 1）! b f（时）（匕） 若舍去七眉莎」飞皿，得求积公式 n f ( x) dx Ai f ( xi ) i= 0 Aj= J：li(x)dx )该公式是 插值型 求积公式. 插值型求积公式的求积余项 b f (竹) R f a (n + 1)! 当f （X）为次数不超过n次的多项式时， 有 f n+1（X）= 0，对应的 R（f） = 0.因此 n 1 个节点的插值型求积公式的代数精度至少 为n.若求积公式 b n f ( x) dx 、 A , f (x,) a i =0 的代数精度至少是n,则该公式是插值型求 积公式• 2. Newto n-Cotes 求积公式 n 1 点的 Newton-Cotes 公式 将求积节点K取为［a,b］上的等距节点 xk二 a+kh,h= — ,k=0,1, |||nn n Z k=0  称ckn)为Cotes系数.求积公式n (b a) Ckn)f (Xk) k=0 称为Newton-Cotes求积公式. 易验证 n x cknj 1 k=0 2 点的 Newton-Cotes 公式 f x dx f a f b a 2 这正是我们熟悉的梯形公式. 3点的Newton-Cotes公式为 b b - a f x dx f a 4f a 6 称它为Simpson公式. 例 算 “ 1沁dx 0 x 用梯形公式计算，有 1-0 1 2 F【K°)+f(l)l = 7(l + sinl)出 0.920 73S49 L L 用Simpson公式计算，有 I ―[f(c)+二 £ (1 + 8 sin；+血1)尢 0.946145 88 梯形公式与Simpson公式的余项 梯形公式余项为 b f气巴) Rr f i I i x - a x - b dx a 2! 利用积分中值定理可有 fr)b (b- a) R f x- a x- b dx=- 2! a 12 梯形公式余项 (a,b) 3 b - a Rt f 12 “…(a,b) Simpson公式的余项 (b- a)5 (a,b) n G（n） 1 1 t 2 2 2 3 1 1 4 — — 5 € 6 6 1 3 1 7 3 8 8 8 8 1 16 2 16 7 90 45 15 45 90 部分Cotes系数 19 2S 25 25 25 19 289 96 144 144 9G 283 41 9 9 34 9 9 41 m 35 289 m 280 35 540 7S1 3577 1323 2989 2989 1323 3577 7S1 172801728017280172801728017280172801728D 989 588B 928 10496 4S40 10496 92B 58B8 989 2B3502035fi 2835028350 2835028350 2B3502835028350 当n较大时Cote