数值微分和数值积分总结复习+习题课(陈)祥解.ppt

* o a=x0 x1 x2 x3 xn-2 xn-1 xn=b y=f(x) y x …… * * …… …… y=f(x) o a=x0 x1 x2 x3 x4 x2i x2i+1 x2i+2 x2m-2 x2m-1 x2m=b y x * 6. 解：利用给定的精度和两种方法的余项公式进行计算 0≤x≤1 课堂练习E：若用复化梯形公式与复化辛卜生公式计算积分 ，问区间[0，1]应多少等分才能使截断误差不超过0.5×10-5？若改用复化辛普生公式，要得到同样精度区间[0，1]应多少等分？ 由复化梯形截断误差公式，得 解得 所以积分区间应等分213份，节点为n+1=214个 由复化辛卜生截断误差公式，得 所以积分区间应等分8份，节点为2m+1=9个 解得 * y=f(x) …… …… y x o a=x0 x1 x2 xi xi+1 xn-1 xn=b x0+1/2 x1+1/2 xi+1/2 xn-1+1/2 * 课堂练习F： 利用复化梯形积分的自适应算法计算 下面的积分，使误差不超过0.5?10-6。 等分数n Tn( f ) 等分数n Tn( f ) 1 0.9207355 32 0.9460596 2 0.9397933 64 0.9460769 4 0.9445135 128 0.9460815 8 0.9456909 256 0.9460827 16 0.9459850 512 0.9460830 * 课堂练习G：用龙贝格方法计算下列积分 要求精确到10-2 解：（1） （2）将区间[0，1]二等分，x=0. 5是新 分点，得 在区间[0，1]上使用梯形公式，得 （3）再将区间等分，增加的分点为 （4）将区间再等分，增加的分点为 （5）将区间再等分，增加的分点为 注意:什么时候停止计算? k 区间等分数 T S C R 0 1 0.5 1 2 0.6036 0.6380 2 4 0.6433 0.6565 0.6578 3 8 0.6582 0.6632 0.6536 0.6536 4 16 0.6636 0.6654 0.6656 0.6658 5 32 0.6656 0.6663 0.6663 0.6663 * 数值微分 差商求导 插值型求导 两点公式 三点公式 数值积分 插值型积分（代数精度） 梯形公式 辛普森公式 复化梯形公式 复化辛普森公式 自适应积分公式 龙贝格积分公式 上节课我们介绍了计算方法这门课程的基本信息，应用背景，发展前景以及误差理论。今天开始，我们来介绍数值逼近这一个大方向。这一部分内容最接近计算数学所讲述的实质性的东西，也是应用范围最为广泛的一个方向。包含以下几个内容：1.函数插值；2.数值积分和数值微分；3.曲线拟合。 特别是函数插值，这一部分让大家掌握的是数值逼近的思想； 数值积分和微分，使大家对于连续和离散，极限与逼近这些数学上的概念有更深刻的了解。我们可以学到如何在计算机上实现微积分中学过的很多东西。对于一个连续的复杂的问题，如何用一些简单的模型在满足一定的精度范围内来近似的替代，然后呢，使得它在计算机上的操作成为可行。数值逼近学完了以后，下一部分，也就是我们教材的第8章，微分方程数值解就是这一部分和数值代数的应用。 第三部分，曲线拟合，相当于函数插值已知点本身误差较大的一种数值逼近处理方式。 大部分数值计算的教材只讲这么两部分，微分方程数值解呢只讲常微。我们只要干净利落的干掉将这一部分的内容，大家在学微分方程数值解时就不会有太大的困难。 * * * * * * * * * * * * * * * * * 数值逼近 * 1.函数插值 2.数值积分和数值微分 3.曲线拟合 * 数值微分复习课 o P Q x0 x0+h f(x) y x A类：差商求导 B类：插值型求导 * 数值微分复习课 A类：差商求导 * 1、向前差商截断误差 2、中心差商截断误差 1、向后差商截断误差 * 课堂练习A ：函数f(x)由下表给出，用差商公式计算x=0.00、0.20、0.40处的二阶导数值 x 0.00 0.10 0.20 0.30 0.40 f(x) 1.70 1.50 1.60 2.00 1.90 解 为求二阶导数值，须