数值微分的计算方法.docVIP

数值微分摘要 求解数值微分问题,就是通过测量函数在一些离散点上的值,求得函数的近似导数。实验结果表明了方法的有效性关键词 　　数值微分根据函数在一些离散点的函数值，推算它在某点的导数或高阶导数的近似值的方法。一个能够近似代替该函数的较简单的可微函数（如多项式或样条函数等）的相应导数作为能求导数的近似值当函数可微性不太好时，利用样条插值进行数值微分要比多项式插值更适宜。Taylor展开式方法 理论基础：Taylor展开式 我们借助Taylor展开式，可以构造函数在点的一阶导数和二阶导数的数值微分公式。取步长则 （1） 所以 （2） 同理 （3） （4） 式(2)和式(4)是计算的数值微分公式，其截断误差为，为提高精度，将Taylor展开式多写几项 两式相减得 （5） 上式为计算的微分公式，其截断误差为O(h2),比式（2）和(4)精度高。 两式相加，如果，则有 , (6) 式(6)是计算的数值微分公式,其截断误差为。 例1．设函数，试用数值微分公式计算的值。 解 由式(2)、式(4)和式(5)分别计算结果为 与真值相比，式(5)计算的结果精度较高。 2.数值微分的Lagrange插值方法 设函数具有个实验数据：，我们希望估计的值，特别时，估计的值。 基于插值方法的数值微分做法是，由已知建立Lagrange插值多项式或Newton插值多项式（这里以Lagrange插值方法为例），即 于是 当时，有 其中 略去误差项有 实际运用中，等距节点更为常见。设 ， 于是有 所以 3.多点数值微分公式 由于高阶插值的不稳定性，实际应用时多采用n=1,2,4的两点、三点和五点等多点插值型求导公式。 两点公式(n=1) (7) （2）三点公式(n=2) （8） (3)五点公式（n=4） （9） 例2．设，取，分别用三点公式和五点公式计算的近似值。 解:由式(8)有 由式(9)有 与真值相比，三点公式已有相当满意精度，而五点公式的结果是十分满意的。 4.数值微分的隐式格式 前述的Taylor展式法和数值微分格式均称为显式格式，即直接由已知的，经过适当的算术四则运算，立即可得的近似值。显式格式优点是计算方便，工作量小，缺点是数值不稳定。为克服后一缺点，隐式格式常常具有数值稳定性。 数值微分的隐式格式建立方法常用通过Taylor展开式方法或数值积分方法等不同途径。 首先我们用Taylor展开式方法来推导数值微分的隐式格式。 由式(5)和式(6)，我们用代替得： 所以也有 将最后表达式代入表达式可得 略去误差项，并且表示的近似值，则有 同样，用数值积分方法也可推导数值微分的隐式格式。 根据simpson公式有 上式是关于个未知量的个方程，如果已知 ，记，故有方程组： （10） 这就是求的线性方程组。由于系数矩阵是严格对角占优的三对角矩阵，因此非奇异，解存在且唯一，可由追赶法求解，且数值稳定。 例3．设函数在节点上的函数值，及如表1所示。试用数值微分隐式格式式(10)和式(11)求出相应节点上一阶导数的近似值。 表1 例3节点数据 x f(x) f (x) 1.5 0.405465108 0.666666667 1.6 0.470003629 1.7 0.530628251 1.8 0.587786664 1.9 0.641853886 2.0 0.693147182 0.500000000 解：由式(10)有 解方程组得： =0.62499828611483 =0.58823447854067 =0.55555484972249 =0.52631517256938 与在各相应节点上数值相比，上述结果约有五位有效数字，精度较高。h越小，精度越高。 附追赶法程序： %追赶法 %定义三对角矩阵A的各组成单元。方程为Ax=d % a为主对角线下面的次对角线，b为主对角线，c为主对角线上面的次对角线，d为右端向量% A=[4 1 0 0 % 1 4 1 0 % 0 1 4 1 % 0 0 1 4] function x=zhuiganfa(a,b,c,d) a=