第8講(2014-11-19).docVIP

数值积分与数值微分 上次课要点： §4 龙贝格求积公式 固定复化求积公式中的步长会导致： （1）步长取得太大，精度难以保证； （2）步长太小，则计算量太大，有时不必要。 解决思路：采用变步长方法。 一、梯形法的递推化 将区间进行等分，设节点为，，则梯形法求积公式为 再将每一个子区间二等分，记相应的中点为，则此时梯形法求积公式为 于是有递推公式 。 二、龙贝格算法的思路 上述算法过于依赖梯形公式，效率较低。 问题：能否利用和的某种组合，得到精度更高的求积公式？ 理查德外推法：用精度较低的近似公式，组合成精度较高的近似公式。 龙贝格将它应用到数值积分上。 等于辛普森法积分值，即 . 用同样方法，可得 三、龙贝格算法 步1 初始化：计算，置（记录区间的二分次数），； 步2（二分） 计算积分值： ； 步3（加速） 求加速值：计算 ，； 步4 精度检验：对指定的精度，若 ，则终止计算，并取作为所求的结果；否则置，，转步2。 §5 自适应积分方法 问题： 如果被积函数在积分区间上，有的部分变化剧烈，有的部分变化平缓，则将整个积分区间一起细分是非理性的。 解决方法： 根据各个子区间满足精度的情况，采用不同的步长——这就是自适应积分方法。 §8 数值微分 研究内容： 用数值方法求函数在离散点处的导数值。 解决方法： 用某些点处函数值的线性组合，近似指定点处的导数值。 问题分为两类： （1）求一点处的导数值； （2）求等距分布的多个点处的导数值。 一、泰勒展开法 1、显式求导公式（用于求一个点处的一阶、二阶导数） （1）一阶向前差商公式（用于求区间左端点的近似导数）： （2）一阶向后差商公式（用于求区间右端点的近似导数）： （3）一阶中心差商公式（用于求区间内点的导数值）： （4）二阶中心差商公式（用于求区间内点的二阶导数值）： 注：不能太小，也不能太大。 2、隐式求导公式 节点： 函数值： 边界条件：和 求：。 （1）用表示的近似值，得 这是由个未知数、个方程组成的线性方程组。 上述方程组主对角占优，解存在且唯一。用追赶法求解是数值稳定的。 本次课继续。。。。。。 再看二阶导数。由在处的泰勒展开式，可得 而 ， 于是 用表示的近似值，得近似等式 即 再加上两个边值条件 ， 并令 便得（大约）具有精度的二阶导数隐格式方程组 上述方程组同样是主对角占优的，解存在且唯一，用追赶法求解数值稳定。 二、插值型求导公式 过已知点，，建立插值多项式，用的导数作为的导数的近似值，即 。 （一）用Lagrange插值多项式求近似导数 常用的有三点求导公式和五点求导公式。 1、三点公式 已知三个点 ， 则 的导函数 因此，各点处的一阶导数近似公式： 2、五点公式 已知五个点 用三点类似的方法可得 （二）用三次样条插值求导（求多个点处的导数值） （略） 三、利用数值积分求导 设，， ， 则 上式右边的积分可以采用不同的求积公式，从而得到不同的求导公式。 1、中矩形法（用于求一点处的导数） 对 用中矩形公式 则得 即 。 2、辛普森法（用于求多点处的导数） 对 用辛普森公式 则得 用表示的近似值，得 这是由个未知数、个方程组成的线性方程组。 若已知，则可得到前面相同的一阶导数隐格式。 四、数值微分的外推算法 外推法是一种逐步提高公式精度的方法，使截断误差逐步减小。由泰勒展开式 （1） （2） 得 （3） 其中只与有关，与无关。 令，则 （4） （5） 由得 （6） 可见，用精度更高，余项为。 令，可以继续外推。一般地，令 ， 可得 ， 误差： 。 当较大时，截断误差很小。考虑到舍入误差，一般不能太大。 五、复数法求近似导数 William Squire , George Trapp. Using Complex Variables to Estimate Derivatives of Real Functions. SIAM Review, Vol.40, No.1, PP110-112, March 1998. 算法要点： 将按复数展开，表示虚数单位。 则虚部 令 。 该式计算的截断误差为。 没有两个相近数相减产生的有效数字位数的损失。 上述算法与中心差商法进行数值比较 （1） （2） 例1 ， 公式（1） 公式（2） 0.1D-01 0.1D-02 0.1D-03 0.1D-04 0.1D-05 0.1D-06 0.1D