数值计算方法第7章+数值微分与数值积分part1.ppt

第7章 数值微分与数值积分 差商型求导公式的精度分析 几点说明 插值型求导公式 常用的等距插值型求导公式 利用样条函数求数值微分 二. 牛顿-柯特斯求积公式 插值型求积公式 Newton-Cotes公式 几个特殊的NewtonCotes 积分公式 n=2 的情况 Cotes 系数表 (表7-2) 误差估计 求积公式的代数精度 (续) 梯形公式和Simpson公式的误差估计 定理7.3 定理7.3 证明 (续) 三. 复化求积公式 复化梯形公式的截断误差 复化Simpson公式 复化Simpson公式的截断误差 复化求积公式举例 例. (复化求积公式) (续) 逐次分半算法 梯形公式的逐次分半算法 (续) 算法7.1 梯形公式逐次分半算法的描述 Simpson公式的逐次分半算法 Simpson逐次分半算法的误差估计 算法7.2 Simpson逐次分半算法的描述 四. 龙贝格(Romberg)求积公式 构造加速收敛的积分序列（续） Romberg求积公式 就是Simpson积分 Richardson外推公式的一般形式 Romberg积分的实施步骤 算法7.3 Romberg积分 Romberg积分的例子 Romberg积分的例子（续） Romberg方法的收敛性 第7章 数值微分与数值积分 方法的基本思想 将Richardson外推法用于梯形积分公式序列 . 设 将区间[a,b]分成2k等分, 复化梯形公式近似: 其中, 为 积分的步长, (k=0,1,2,…) 由复化梯形公式的误差: 用Richardson外推法, 取α=1/2, 则: 即: 误差: 第7章 数值微分与数值积分 验证如下: 将 代入, 并记: 即, 2m等分的复化Simpson公式 第7章 数值微分与数值积分 即, Romburg 求积公式 仍取α=1/2, 重复使用Richardson外推法, 得 积分误差: 第7章 数值微分与数值积分 3 2 1 0 k m=0 m=1 m=2 m=3 m=4 逐次分半梯形公式 逐次分半的Simpson公式 Romberg方法 按上表逐行进行计算, 当 则满足精度要求, 作为积分的近似值. 等分区间 第7章 数值微分与数值积分 对书上的算法作了修改, 更易于编程时使用. 用2个数组T1 和T2 分别存放上面积分表中最近两行的结果. 输入积分区间[a,b], 函数f(x) 以及积分误差限ε; c?1.0; for j=1 to m step 1 do: 若 , 输出积分近似值T2(m), 停机; 否则 { T2(0:m)?T1(0:m); h/2?h; n+n?n; m+1?m; 转3; } Fortran程序… 置 计算 ; 计算: 逐次分半梯形公式 注. n是下一步逐次分半算法中区间的等分数, m是Romberg积分阶数. Romberg积分公式 end for (j). 置q?0; for j=1 to n-1 step 2 do: end for (j). 第7章 数值微分与数值积分 解. (1) 用梯形公式, h=1, 计算 计算积分 的近似值, 要求误差不超过 . 用逐次分半算法(h=0.5)计算: 用Romberg积分公式计算: 3 2 1 0 k (2) 用梯形逐次分半算法(h=0.25)计算: 用Romberg积分公式计算: 第7章 数值微分与数值积分 0.9460832 0.9460832 0.9460834 0.9456909 3 0.9460830 0.9460869 0.9445135 2 0.9461459 0.9397933 1 0.9207355 0 k m=0 m=1 m=2 m=3 用Romberg积分公式计算: (3) 用逐次分半算法(h=0.125)计算: 积分近似值: * 实际问题中, 很多函数仅给出一些离散点上的值, 因此无法用分析的方法求导与积分, 只能通过数值方法解决该类问题. 由导数的定义: 背景 一. 数值微分 理论依据 当h足够小, 作为近似的导数值. 差商型求导公式 向前差商 向后差商 中心差商 由Talor 展开公式: 向前差商公式的余项, 为O(h)精度. 第7章 数值微分与数值积分 向后差商公式的余项, 为O(h)精度. 进一步