《第五章数值微积分.doc》由会员分享，可在线阅读，更多相关《第五章数值微积分.doc（7页珍藏版）》请在金锄头文库上搜索。
1、第五章 数值微积分一、内容分析与教学建议本章内容是数值微积分。数值微分包括：用插值多项式求数值微分、用三次样条函数求数值微分和用Richardson外推法求数值微分。数值积分包括：常见的Newton-Cotes求积公式，如：梯形公式、Simpson公式和Cotes公式；复化求积公式；Romberg求积公式和Gauss型求积公式等内容。（一） 数值微分1、利用Taylor展开式建立数值微分公式，实际上是利用导数的离散化，即用差商近似代替导数，在由Taylor公式的余项估计误差；由于当步长很小时，回出现两个非常接近的数相减，因此，在实际运用中往往采用事后估计的方法来估计误差。2、用插值多项式求数值
2、微分，主要是求插值节点处的导数的近似值。借助第二章的Lagrange插值公式及其余项公式，确定插值节点处的导数的近似值及其误差。常用的有三点公式和五点公式。3、阐明用三次样条函数求数值微分的优点：由第三章的三次样条函数的性质知：只要的4阶导数连续，则当步长时，收敛到，收敛到，收敛到. 因此，用三次样条函数求数值微分，效果是很好的。指出其缺点是：需要解方程组，当很小时，计算量较大。4、讲解用Richardson外推法求数值微分时，首先阐明方法的理论基础是导数的离散化，即用差商近似代替导数；然后重点讲解外推法的思想和推导过程，因为这种方法和思路在后面的数值积分和微分方程数值解中还要用到。（二） 数
3、值积分的一般概念1、由定积分的几何意义引入数值积分的思想，介绍求积公式、求积节点、求积系数、余项等基本概念。2、重点介绍代数精度以及如何求一个判定积公式的代数精度，并举例说明。3、介绍插值型求积公式以及插值型求积公式的代数精度的特点。（三） 等距节点的求积公式1、 简单介绍一般的等距节点的插值型求积公式Newton-Cotes公式以及Cotes系数。 2、重点介绍几种常用的Newton-Cotes公式：梯形公式、Simpson公式和Cotes公式。要求学生掌握上述三种求积公式的表达式，并了解三种求积公式各自的余项。3、以Simpson公式为例，求出它的代数精度是3；并要求学生课后自己求出梯形公
4、式和Cotes公式的代数精度。（四） 复化求积公式1、结合分段插值的思想阐明复化求积公式的思想。2、重点介绍复化梯形公式、复化Simpson公式和复化Cotes公式以及它们各自的余项，并举一、两个例子加以说明。3、简介事后估计和自适应Simpson方法。（五） Romberg求积法1、Romberg求积法是一种逐步分半加速法，它是以复化梯形公式为基础构造高精度求积公式的方法，是一种快速、有效的求积法。2、阐明Romberg公式的建立过程：利用事后估计的思想，从复化梯形公式建立一整套递推算法，进而得到Romberg公式，整个过程实际上是一个加速的过程。3、可通过例子验证Romberg求积法的加速
5、效果。（六） Gauss型求积公式1、Gauss型求积公式也是一种高精度的插值型求积公式，但它的节点不是等距的，因而Gauss型求积公式不属于Newton-Cotes公式的范畴。2、阐明Gauss型求积公式的代数精度是插值型求积公式的最大值，介绍Gauss点的概念，并说明Gauss点实际上是某个正交多项式的零点。3、讲清楚Gauss型求积公式的求积系数的特殊构造，并由此证明Gauss型求积公式是稳定的，以及Gauss型求积公式的收敛性。4、 介绍几种Gauss型求积公式：古典Gauss公式、Gauss-Tchebyshev公式、Gauss-Laguerre公式和Gauss-Hermite公式。
6、让学生了解上述四中Gauss型求积公式的表达式、表达式中的权函数、定积分的上、下限以及求积系数，并通过23个例子具体阐述上述Gauss型求积公式是如何求数值积分的，并和以前的方法比较它们的精度。本章结束时，建议安排一次上机实习，让学生自己动手，根据书中的算法，编程计算各种数值积分的例子，加深和巩固学生对本章内容和方法的了解和掌握。二、补充例题例1 用三点公式求在处的导数值，并估计误差，的函数值由下表给出： .解 三点求导公式为取上表中，再分别将有关数值代入上式，即可得导数的近似值。因为，所以可得误差估计及导数值如下表：例2 从地面发射一枚火箭，在最初80秒内，记录其加速度如下表。试求火箭在第8
7、0秒时的速度。分析：速度对时间的导数等于加速度，因此已知加速度求速度，只需把速度看作是加速度的原函数即可。若设速度为，则，于是.这样就把问题转化为求积分的问题。解 应用复化Simpson求积公式计算。此题中积分区间的长度是80，有9个节点，故.由于火箭从地面向上发射，因此. 于是火箭在第80秒时的速度为 例3 计算椭圆的周长，使结果具有5位有效数字。分析：这是一个求周长的问题，因此要用到线积分中的弧长公式。在估计误差时，由于弧长公式中含有根式，其高阶导数较复杂，故可用事后误差估计的方法来做；另外还必须把误差与有效数字结合起来使用。解 由于在直角坐标系下求弧长表达式较复杂，因此采用极坐标来求解。
8、令，则椭圆弧长为，因为，所以有一位整数。故若要求结果有5位有效数字，则必须使截断误差. 列表计算如下：等分012312482.35619452.419920782.422103102.422112060.02124210.000727440.000002986故可取可使有5位有效数字，从而.例4 用反证法证明：不存在，使得求积公式的代数精度超过次。分析：只要能找到一个次的多项式，使求积公式两边不相等即可。而具有次代数精度的求积公式的节点是上带权的正交多项式的零点,可考察次的多项式.解 构造多项式，并令，代入上述求积公式，则左端有；右端有； 即左端右端。这说明：不存在具有次代数精度的求积公式。故Gauss型求积公式是具有最高次代数精度的求积公式。例5 设，求证：(1) ；(2) .证 本题用Taylor公式来证。(1) 因为 ，所以，即 .(2) 利用(1)中的展开式，得.例6 确定常数（均用分数精确表示），使求积公式，其中具有尽可能高的代数精确度，并指出代数精确度是多少？其中.解 设该求积公式对精确成立，得，化简得解得 例7 寻找合适的数值求积公式，计算出积分的准确值。解 因为，其中，权函数，所以可取Gauss-Tchebyshev求积公式 ，其中. 又因为是3次多项式，且具有次代数精度，所以取，可计算出积分的准确值。此时， .
