第3章 数值积分与数值微分教案.pptVIP

图 3.7 例 4 计算积分 精确到10-4。 解 于是 由于 实际上 §4 数值微分 4.1 用插值法求数值微分 用插值多项式Pn(x)近似地表示函数f(x),即 f(x)≈Pn(x) 于是有 f(k)(x)≈P(k)n(x) 其余项相应地为R(k)n(x)。 设插值基点为等距分布,由牛顿前差插值多项式 其中 由于 于是 即 因为 而当x=xi时,s=i,此时 (3―32) (3―33) 特别当x=x0时,s=0,则 (3―34) 1. 两点公式(n=1) 于是在区间［x0,x2］上有 (3―35) 2. 三点公式(n=2) 于是在区间［x0,x2］上有 9.2 用三次样条函数求数值微分 设s(x)是f(x)在各区间［xk,xk+1］上的三次样条插值函数,则在区间［xk, xk+1 ］上可通过三次样条函数来求f(x)的数值微分。 1.一阶数值微分公式 (3―36) 若只求基点xk(k=0,1,…,n-1)上的一阶导数值,则 (3―37) 2. 二阶数值微分公式 特别,若只求基点上的二阶导数值,则 (3―38) (3―39) 例2 若用复合梯形公式计算积分 问积分区间要等分多少才能保证有五位有效数字 解 由余项(3―21)式 则当0＜x＜1时,有 因为 又 故 由于原积分的准确值具有一位整数,因此要使近似积分值有五位有效数字,只需取n满足 两边取对数得 整理后得到 2.2 复合抛物线公式 类似复合梯形公式的做法,把区间［a,b］分成n个相等的子区间［x2i,x2i+2］(i=0,1,…,n-1),设每个子区间上的中点为x2i+1(i=0,1,…,n-1),且 在每一个子区间［ x2i,x2i+2 ］上利用抛物线公式得 (3―22) 相加后得 (3―23) 图 5.4 图 3.4 若f(4)(x)在［a, b］上连续,则 从而得到复合抛物线公式 (3―24) 其余项为 (3―25) 复合抛物线公式的计算框图见3.4。 例3 根据给出的函数 的数据表3―2,分别用复合梯形公式和复合抛物线 公式计算 表 3―2 解 用复合梯形公式,这里 用复合抛物线公式可得 而I的准确值为0.946 083 1…,可见用复合抛物线 公式比用复合梯形公式精确。 2.3 变步长公式 前面介绍的复合梯形公式和复合抛物线公式的步长都是预先确定的。它的主要缺点是事先很难估计出n的大小(或步长h的大小),使结果达到预先给定的精度。在实际计算中,我们常常借助于计算机来完成积分步长h的自动选择,即采用变步长求积公式。具体地讲,就是将步长逐次折半,反复利用复合求积公式,直到满足精度要求为止。 下面介绍变步长复合抛物线公式(变步长复合梯形公式留给读者作为练习)。 逐次将区间［a,b］分成21,22,…,2m等分,并按复合抛物线公式逐次计算积分得到S1,S2,…,Sm,而 (3―26) 其中 再把每个子区间分成两半,用 图 3.5 图 3.5 作步长,按复合抛物线公式计算出积分的近似值S2m。对于相邻两次的积分近似值Sm、S2m,考察 当｜S2m｜＜1 当｜S2m｜≥1 (3―27) 设预先给定的精度为ε,若 ｜d｜＜ε 则以S2m作为所要求的积分近似值,否则继续将区间 分半,利用复合抛物线公式求积分,直到满足预给的精度 为止。 §3 龙贝格(Romberg)积分方法 我们已经知道,当被积函数f(x)在区间［a,b］上连续时,要使得复合梯形公式或复合抛物线公式比较精确地代替定积分