第五章数值微积分

第五章数值微积分一、内容分析与教学建议本章内容是数值微积分。数值微分包括：用插值多项式求数值微分、用三次样条函数求数值微分和用Richardson外推法求数值微分。数值积分包括：常见的Newton-Cotes求积公式，如：梯形公式、Simpson公式和Cotes公式；复化求积公式；Romberg求积公式和Gauss型求积公式等内容。（一）数值微分1、利用Taylor展开式建立数值微分公式，实际上是利用导数的离散化，即用差商近似代替导数，在由Taylor公式的余项估计误差；由于当步长h很小时，回出现两个非常接近的数相减，因此，在实际运用中往往采用事后估计的方法来估计误差。2、用插值多项式求数值微分，主要是求插值节点处的导数的近似值。借助第二章的Lagrange插值公式及其余项公式，确定插值节点处的导数的近似值及其误差。常用的有三点公式和五点公式。3、阐明用三次样条函数()sx求数值微分的优点：由第三章的三次样条函数()sx的性质知：只要()fx的4阶导数连续，则当步长0h时，()sx收敛到()fx，()sx收敛到()fx，()sx收敛到()fx.因此，用三次样条函数()sx求数值微分，效果是很好的。指出其缺点是：需要解方程组，当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、可通过例子验证Romberg求积法的加速效果。（六）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用三点公式求21()(1)fxx在1.0,1.1,1.2x处的导数值，并估计误差，()fx的函数值由下表给出：1.01.11.2()0.2500000.2267570.206612iixfx.解三点求导公式为200120210212201221()[3()4()()](),231()[()()](),261()[()4()3()]();23hfxfxfxfxfhhfxfxfxfhhfxfxfxfxfh取上表中0121.0,1.1,1.2xxx，再分别将有关数值代入上式，即可得导数的近似值。因为551.01.21.01.24!4!()max()max0.75(1)2ixxffxx，所以可得误差估计及导数值如下表：()1.01.11.20.247920.216940.18596()0.250000.215960.187830.002500.001250.002500.002080.000980.00187fxxfx用三点公式求出的的近似值准确值理论误差限实际误差限例2从地面发射一枚火箭，在最初80秒内，记录其加速度如下表。试求火箭在第80秒时的速度。2()01020304050607080(/)30.0031.6333.4435.4737.7540.3343.2946.6950.67ta秒米秒分析：速度对时间t的导数等于加速度，因此已知加速度求速度，只需把速度看作是加速度的原函数即可。若设速度为()vt，则0()(0)()tvtvatdt，于是800(80)(0)()vvatdt.这样就把问题转化为求积分的问题。解应用复化Simpson求积公式计算。此题中积分区间的长度是80，有9个节点，故4,80420nh.由于火箭从地面向上发射，因此(0)0v.于是火箭在第80秒时的速度为808000(80)(0)()()vvatdtatdt120[30.004(31.6335.4740.3346.69)2(33.4437.7543.29)50.67]63087.03333(/).米秒例3计算椭圆2214xy的周长，使结果具有5位有效数字。分析：这是一个求周长的问题，因此要用到线积分中的弧长公式。在估计误差时，由于弧长公式中含有根式，其高阶导数较复杂，故可用事后误差估计的方法来做；另外还必须把误差与有效数字结合起来使用。解由于在直角坐标系下求弧长表达式较复杂，因此采用极坐标来求解。令2cos,sinxy，则椭圆弧长为222222220004()()44sincos413sinlxyddd，因为22013sin222Id，所以I有一位整数。故若要求结果有5位有效数字，则必须使截断误差41102.列表计算如下：k等分2k2kT12213kkTT012312482.35619452.419920782.422103102.422112060.02124210.000727440.00000298641210故可取82.4221IT可使I有5位有效数字，从而49.6884lI.例4用反证法证明：不存在,(0,1,2,,)kkAxkn，使得求积公式0()()()nbkkakxfxdxAfx的代数精度超过21n次。分析：只要能找到一个22n次的多项式，使求积公式两边不相等即可。而具有21n次代数精度的求积公式的节点是[,]ab上带权()x的正交多项式的零点(0,1,2,,)kxkn,可考察22n次的多项式2210()()nniixxx.解构造多项式2210()()()nniiKxxxx，并令()()fxKx，代入上述求积公式，则左端有()()()()0bbaaxfxdxxKxdx；右端有00()()0nnkkkkkkAfxAKx；即左端右端。这说明：不存在具有22n次代数精度的求积公式。故Gauss型求积公式是具有最高次代数精度的求积公式。例5设5000()[2,2],0,,(),0,1,2kkkfxCxhxhhxxkhffxk，求证：(1)4021121()(88)()12fxffffOhh；(2)2010121()(2)()fxfffOhh.证本题用Taylor公式来证。(1)因为230000011(2)()2()(2)()(2)()2!3!fxhfxhfxhfxhfx4(4)501(2)()()4!hfxOh，230000011()()()()()2!3!fxhfxhfxhfxhfx4(4)501()()4!hfxOh，所以500000(2)()8()(2)12()()fxhfxhfxhfxhhfxOh，即4021121()(88)()12fxffffOhh.(2)利用(1)中0()fxh的展开式，得2010121()(2)()fxfffOhh.例6确定常数,,,ABCD（均用分数精确表示），使求积公式()()IfIf，其中23()()()d,()[()()][()()]baIfxafxxIfhAfaBfbhCfaDfb具有尽可能高的代数精确度，并指出代数精确度是多少？其中hba.解设该求积公式对23()1,,(),()fxxaxaxa精确成立，得2231()[11][00]2baxahABhCD，3231()[0][11]3baxahABhhCD，42231()[0][02]4baxahABhhCDh，523321()[0][03]5baxahABhhCDh，化简得1,21,312,413,5ABBCDBDBD解得3711,,,.20203020ABCD例7寻找合适的数值求积公式，计算出积分32143dxxxx的准确值。解因为31122221-1-12143d(2)1d1d(2)2xtxxxxttttt令22111122-1-1-122211(2)1(2)11(2)1dd()d222111ttttttttfttttt，其中2()(2)fttt，权函数21()1xt，所以可取Gauss-Tchebyshev求积公式1-1211()d()1nkkkfttAfxt，其中,1,2,,kAknn.()又因为2()(2)fttt是3次多项式，且()具有21n次代数精度，所以取2n，可计算出积分3121-121143d()d21xxxxfttt的准确值。此时12121232,,cos,cos224242AAxx，221112222222324()(2)2,22422324()(2)2,224fxxxfxxx31211221-1211143d()d[()()]221xxxxfttAfxAfxt132432432224248.