导数的数值计算方法文献综述.doc

文献综述 导数的数值计算方法 前言部分 导数概念的产生有着直觉的起源，．，．，，，，． ，，，，．，，， Rolle， Lagrange， Cauchy ， Taylor ， ，．，，，．，． ， ， 需要注意的是微分是非常敏感的问题，． 数学中研究导数，．．，，，，． (1820-1895)曾指出:“在一切理论成就中，17世纪下半叶微积分的发明那样被看作人类精神的最高胜利了”． ，，，．，16世纪才应运萌生，从15世纪初文艺复兴时期起，，，，．，，，，，，．，: (1) 求变速运动的瞬时速度； (2) ； (3) ． ，．(Fermat)为研究极值问题而引入的，:已知运动规律求速度和已知曲线求它的切线．(Newton)和德国数学家莱布尼茨(Leibniz)分别在研究力学和几何学过程中建立起来的． （一）利用差商表求导数 最简单的数值微分公式是用向前差商近似代替导数： ． （1．1） 类似地，也可用向后差商近似代替导数 ． （1．2） 或用中心差商近似代替导数 ． （1．3） T A B C 数值微分示意图 在几何图形上，这3种差商分别表示弦AB，AC和BC的斜率．将这3条弦同过A点的切线AT相比较，从上图可以看出，一般地说以BC的斜率更接近于切线AT的斜率，因此就精确度而言，以式（1．3）更为可取．称 （1．4） 为求的中点公式． 现在来考察用式（1．4）代替所产生的截断误差．首先分别将在处作Taylor展开，有 ． 然后代入中点公式（1．4），得 ． 所以截断误差 （1．5） 由此可以得到：从截断误差的角度来看，步长越小，计算结果越准确．但从计算角度看，越小，与越接近，直接相减会造成有效数字的严重损失，因此从舍入误差的角度看步长不宜去的太小．怎样选取合适的步长呢？可采用二分步长及误差事后估计法，即比较二分前后所得值与，若，则为所需的合适的步长且． （二）插值型求导公式 对于列表函数 … … 应用插值原理，．，， （2．1） 统称为插值型求导公式． ． 式中，且依赖于；． ． （2．2） 由于是的位置函数，因此求较麻烦，一般都限定求某个节点上的导数值，此时（2．2）右端的第2项由于而变为零，这时的截断误差为 ． （2．3） 由于以上的原因，以下仅考察节点处的导数值，为简化讨论，假定所给节点是等距的． ， ． 2． ． ． ． ， 我们对它不作深入讨论，但要指出的是，尽管与的值相差不多，，． ，．，；，．，，． 理查森外推法是科学计算领域提高算法精度的重要方法，，． ，． ． （3．1） 它由泰勒定理的两种情况导出： ， （3．2） ． （3．3） 我们现在就引入理查森外推的过程，并介绍如何利用它来巧妙地改进数值公式的精度．把（3．2）式和（3．3）式扩展到具有高阶项．假设用它的泰勒级数表示为 ， （3．4） ． （3．5） 如果第一个等式减去第二个等式，则消去了所有是偶数的项，得 ． 重新整理得 ． 这个等式具有形式 （3．6） 其中表示，表示数值微分公式（3．1）；即 ． 其中