第五章数值微积分一、内容分析与教学建议本意内容是数值微积分。
数值微分包括：用插值多项式求数值微分、用三次样条函数求数值微分和用Richardson外推法求数值微分。
数值积分包括：常见的Newton-Cotes求积公式，如：梯形公式、Simpson公式和Cotes公式；复化求积公式；Romberg求积公式和Gauss型求积公式等内容。
(一)数值微分1、利用Taylor展开式建立数值微分公式，实际上是利用导数的离散化，即用差商近似代替导数，在由Taylor公式的余项估计误差；由于当步长h很小时，回出现两个非常接近的数相减，因此，在实际运用中往往采用事后估计的方法来估计误差。
2、用插值多项式求数值微分，主要是求插值节点处的导数的近似值。
借助第二章的Lagrange插值公式及其余项公式，确定插值节点处的导数的近似值及其误差。
常用的有三点公式和五点公式。
3、阐明用三次样条函数s(x)求数值微分的优点：由第三章的三次样条函数s(x)的性质知：只要f(x)的4阶导数连续，则当步长h 0时，s(x)收敛到f (x) , s(x)收敛到f (x) , s (x) 收敛到f (x).因此，用三次样条函数s(x)求数值微分，效果是很好的。
指出其缺点是：需要解方程组，当h很小时，计算量较大。
4、讲解用Richardson外推法求数值微分时，首先阐明方法的理论基础是导数的离散化，即用差商近似代替导数；然后重点讲解外推法的思想和推导过程，因为这种方法和思路在后面的数值积分和微分方程数值解中还要用到。
（二）数值积分的一般概念1、由定积分的几何意义引入数值积分的思想，介绍求积公式、求积节点、求积系数、余项等基本概念。
2、重点介绍代数精度以及如何求一个判定积公式的代数精度，并举例说明。
3、介绍插值型求积公式以及插值型求积公式的代数精度的特点。
（三）等距节点的求积公式1、简单介绍一般的等距节点的插值型求积公式--- Newton-Cotes公式以及Cotes系数。
2、重点介绍几种常用的Newton-Cotes公式：梯形公式、Simpson公式和Cotes公式。
要求学生掌握上述三种求积公式的表达式，并了解三种求积公式各自的余项。
3、以Simpson公式为例，求出它的代数精度是3;并要求学生课后自己求出梯形公式和Cotes公式的代数精度。
（四）复化求积公式1、结合分段插值的思想阐明复化求积公式的思想。
2、重点介绍复化梯形公式、复化Simpson公式和复化Cotes公式以及它们各自的余项，并举一、两个例子加以说明。
3、简介事后估计和自适应Simpson方法。
（五）Romberg求积法1、Romberg求积法是一种逐步分半加速法，它是以复化梯形公式为基础构造高精度求积公式的方法，是一种快速、有效的求积法。
2、阐明Romberg公式的建立过程：利用事后估计的思想，从复化梯形公式建立一整套递推算法，进而得到Romberg公式，整个过程实际上是一个加速的过程。
3、可通过例子验证Romber评积法的加速效果。
（六）Gauss型求积公式1、Gauss型求积公式也是一种高精度的插值型求积公式，但它的节点不是等距的，因而Gauss型求积公式不届于Newton-Cotes公式的范畴。
2、阐明Gauss型求积公式的代数精度是插值型求积公式的最大值，介绍Gauss点的概念，并说明Gauss点实际上是某个正交多项式的零点。
3、讲活楚Gauss型求积公式的求积系数的特殊构造，并由此证明Gauss型求积公式是稳定的，以及Gauss型求积公式的收敛性。
4、介绍几种Gauss型求积公式：古典Gauss公式、Gauss-Tchebyshev公式、Gauss-Laguerre公式和Gauss-Hermite公式。
让学生了解上述四中Gauss型求积公式的表达式、表达式中的权函数、定积分的上、下限以及求积系数，并通过2— 3个例子具体阐述上述Gauss 型求积公式是如何求数值积分的，并和以前的方法比较它们的精度。
本意结束时，建议安排一次上机实习，让学生自己动手，根据书中的算法，编程计算各种数值积分的例子，加深和巩固学生对本意内容和方法的了解和掌握。
二、补充例题例1用三点公式求f (x)([ 1)2在x 1.0, 1.1, 1.2处的导数值，并估计误差，f(x)的函数值由下表给出：X j 1.0 1.1 1.2. f(x) 0.250000 0.2267570.206612解三点求导公式为取上表中X0 1.0, X1 1.1, X2 1.2,再分别将有关数值代入上式，即可得导数的近似值。
因为f ( i) max f (x) max —4—5弓0.75 ,所以可得误差估计及导数值如下表：'I 1.0 x 1.2 '， 1.0 x 1.2 (1 x)525例2从地面发射一枚火箭，在最初80秒内，记录其加速度如下表。
试求火箭在第80 秒时的速度。
分析：速度对时间t的导数等丁加速度，因此已知加速度求速度，只需把速度看作是加t 80、、速度的原函数即可。
右设速度为v(t)，则v(t) v(0) 0a(t)dt , 丁是v(80) v(0) 0 a(t)dt . N样就把问题转化为求积分的问题。
解应用复化Simpson求积公式计算。
此题中积分区间的长度是80,有9个节点，故n 4, h 80/4 20.由丁火箭从地面向上发射，因此v(0) 0. 丁是火箭在第80秒时的速度80 80为v(80) v(0) 0 a(t)dt 0 a(t)dt例3计算椭圆丈y2 1的周长，使结果具有5位有效数字。
4分析：这是一个求周长的问题，因此要用到线积分中的弧长公式。
在估计误差时，由丁弧长公式中含有根式，其高阶导数较复杂，故可用事后误差估计的方法来做；另外还必须把误差与有效数字结合起来使用。
解由丁在直角坐标系下求弧长表达式较复杂，因此采用极坐标来求解。
令x 2cos , y sin ,则椭圆弧长为l 4 : ~(y )2 d 4 j cos2 d 4 j ~~3sin2 d ,因为—I W 3sin2 d 2 —，所以I有一位整数。
故若要求结果有5位有效数字, 2 02则必须使截断误差 1 104.列表计算如下：2故可取I T8 2.4221可使I有5位有效数字，从而l 4 I 9.6884 .例4用反证法证明：不存在A k, X k (k 0,1,2,|||,n)，使得求积公式的代数精度超过2n 1次。
分析：只要能找到一个2n 2次的多项式，使求积公式两边不相等即可。
而具有2n 1次代数精度的求积公式的节点是 [a, b]上带权(x)的正交多项式的零点x k (k 0,1,2，川,n),可 n 考察2n 2次的多项式2 3i (x) (x X i )2.n解 构造多项式K(x) 2 i (x) (x x i )2，并令f (x) K(x),代入上述求积公式，贝U 左i 0b n n(x)K(x)dx 0;右端有 A k f(xQ AKg) 0; 即左端 右端。
这k 0度的求积公式。
1 ， ， ， (f2 8f i 8f i12h2 3 1⑵ f(x °) *1 2f 0 f 1) O(h2).证本题用Taylor 公式来证。
1(2h)4 f (4) (X 0) O(h 5), 4!上 h 4 f (4)(x °) O(h 5), 4!5f(x 0 h) 8f(x 0 h) f (X 0 2h) 12h f (X 0) O(h 5),即 f(x °) ^(f 2 8f 1 8f 1 f 2) O(h4).⑵ 利用⑴中f(X 0 h)的展开式，得3，,1 , , , - 2E ?(f 1 2f 0 f 1) O(h).b跚有 (x)f(x)dx说明：不存在具有 2n 2次代数精度的求积公式。
故Gauss 型求积公式是具有最高次代数精例5设f (x)5C 5[x ° 2h, x) 2h], h 0, x kx kh, f k f (x<), k 0, 1, 2 ,求证：⑴f (x °)f 2) O(h 4);(1)因为 f(x 。
2h) f (x)) 2h f (x 0)1 2.1 3 _(2h) f (X ))(2h) f (X 0)2!3!所以 f (X 0 2h)例6确定常数A, B, C, D (均用分数精确表示)，使求积公式@(f) I(f),其中I(f) ：(x a)f(x)dx, @(f) h 2[Af(a) Bf(b)] h 3[Cf (a) Df (b)]具有尽可能高的代数精确 度，并指出代数精确度是多少其中h b a .解 设该求积公式对f(x) 1, x a, (x a)2, (x a)3精确成立，得-(x a)2 b h 2[A 1 B 1] h 3[C 0 D 0],2a -(x a)3 b h 2[A 0 B h] h 3[C 1 D 1],3 a !(x a)4 b h 2[A 0 B h 2] h 3[C 0 D 2h],4 a l(x a)5 b h 2[A 0 B h 3] h 3[C 0 D 3h 2],5a化简得解得 A —, B —, C —, D —. 2020 30 20例7寻找合适的数值求积公式，计算出积分1x—x 2 —dx 的准确值。
3J 2x t 21；21 1「2 2x_4x x 2 3dx(t 2) .1 t 2dt — 、.1 t 2d(t 2)21-12 -1_ 2 _ 21t(t 2)21 1t(t 2)2 1 1 1—2) dt 1 L (： 2) dt - , 1 f (t)d t , 1-1.1 t 2 2 -1 ,1 t 2 2 -1 ,1 t 21(x)—，所以可取Gauss-Tchebyshev 求积公式-^=f(t)dt A<f (x k )，其中 A 一，k 1,2，川,n . 1 t 2k 1n解因为其中 f (t) t(t(t 2)2 ■ 1 t 22)2,权函数 11 -1乂因为f(t) t(t 2)6是3次多项式，且()具有2n 1次代数精度，所以取n 2,可计算出积分：x 4x X2 3dx 711 2 -1i t2f(t)dt的准确值。
此时Al — , A21 2 3 2—,x〔cos— --- , X2 cos— ---2 2 4 2 4 2f (X i)2松 1 2)2了 2 守，2 2 4f (X2) x2(x2 2)2M 号 2 2 *2 2 46 2 4 2 4 87 _ 3 2 4 3 .4 3 2。