第6章数值微积分精选.docVIP

第章 数值微积分 .1 引言 的导数（或微分）和求函数在闭区间上的定积分等问题. 由高等数学知道：若函数为初等函数，则可由基本求导公式和运算法则求出的导函数；若函数在区间上连续且原函数为，则由Newton-Leibniz公式 可求得在区间上的定积分.但是，在实际问题中遇到的求导数、求定积分的计算中，经常会有这样的情况： (1) 函数的原函数无法用初等函数给出，例如定积分 等，从而无法用Newton-Leibniz公式计算定积分。 (2) 只给出了函数的若干数据，函数表达式未知，因而无法用导数公式或定积分公式。 (3) 函数的导函数或原函数虽然能够求出，但形式过于复杂，不便使用. 由此可见，利用求导公式及求导法则求导数或利用原函数求定积分有它的局限性。因此，研究如何利用函数在有限个点处的函数值计算其导数的近似值，以及定积分的近似值的数值计算方法——数值微分(numerical differentiation)和数值积分(numerical integration)，既有理论意义又有实际应用价值。 这些方法也是微分方程和积分方程数值解法的基础。 6.2 数值微分 .2.1 插值型的数值微分公式 已知数据. 记 是包含插值节点的区间, ，则的n次插值多项式为 , 其中, ，误差为 , 故 . (6.2.1) 因为是的未知函数，所以无法对上式右端的第二项作出判断，因而对于任意给定的点，误差是无法预估的. 如果只是求某个节点上的导数值，这时有 . (6.2.2) 以为例，为简便起见，设，已知 . 设满足插值条件的2次插值多项式为，即 则 . (6.2.3) 于是得三点数值公式 (6.2.4) 类似上述推导，在等距节点的情形，即时，可得如下常用的数值公式. (1) (two-point formulas) (6.2.5) (2) 三点公式(three-point formulas) 一阶求导公式： (6.2.6) 二阶求导公式： (6.2.7) (3) 五点公式(five-point formulas) 一阶求导公式： (6.2.8) 二阶求导公式： (6.2.9) 注 上述数值求导公式不仅是计算导数值的有效方法，在微分方程数值解及样条插值的数据扩充方面都很有用处. 例6.2.1 已知函数的数值表6.2.1 0.6 0.7 0.9 0.826446 0.591716 0.510204 0.444444 0.390265 0.346021 0.277008 分别取，试用两点、三点和五点公式计算处的一、二阶导数值。 时 , 当时 , , , , .公式和和6.2.2中。 表6.2.2 两点公式 三点公式 五点公式 精确值 0.2 -0.736360 -0.614238 -0.590051 -0 1.221225 1.181288 1 0.1 -0.657600 -0.599695 -0.594848 -0 1.158100 1.137058 1 注意和的值是. 上面的计算表明五点公式比三点公式与两点公式准确，步长越小结果越准确. 一般情况下，这个结论也是对的. 但由误差表示式可见，如果高阶导数无界，或者舍入误差超过截断误差时，这个结论就不对了.2.2 利用三次样条求数值微分 前述数值微分公式不便于计算节点之间的导数值；而且使用多点公式时由于截断误差包含高阶导数，当高阶导数数值较大时难以保证截断误差很小，此时最好用样条函数求数值微分。由于三次样条是最常用的样条函数，本节只讨论利用三次样条求数值微分。 由.6节知，对于剖分 ， 三次样条函数在子区间上的表达式为 , (.2.10) 其中 对(6.2.10)求导得 . (.2.11) 于是得 . 若只要求节点上的值，则由上式可直接得到 对(6.2.11)求导得 , (.2.12) 于是得 , 由定理4.6.3用样条函数建立的数值公式. 其缺点是需要解方程组，当很小时，计算量较大.2.3 变步长的中点方法 因为，所以若精度要求不高，我们可以简单地取差商作为导数的近似值：分别可用差商向后差商或中心差商近似代替导数，即 . (6.2.15) , (6.2.16) (