第数值微积分

第五章数值微积分

一、内容分析与教学建议

本意内容是数值微积分。数值微分包括：用插值多项式求数值微分、用三次样条函数求数值微分和用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.226757

0.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

分析：这是一个求周长的问题，因此要用到线积分中的弧长公式。在估计误差时，由

丁弧长公式中含有根式，其高阶导数较复杂，故可用事后误差估计的方法来做；另外还必

须把误差与有效数字结合起来使用。