ch06 数值积分.pptVIP

第6章 数值积分 引入 Newton-Cotes积分公式 变步长积分与Romberg求积公式 Gauss求积公式 数值微分 6.1 引入 问题提出 数值积分 常见的 几种求积公式 代数精度 1。问题的提出 若f(x)连续，则f(x)存在原函数F(x)，据Newton-Leibnits公式有： 但是：很多情况下用此法求定积分有困难： F(x)(原函数）存在但求不出： F(x)表达式很复杂 f(x）本身是一组离散的观测数据 需要用数值方法计算定积分的近似值 2。数值积分 定积分的极限定义： 设想：能否不求f(x)的原函数，利用(1)去掉极限，并将ξk取为具体的xk 数值积分公式In( f )： 定义：是离散点处函数值的线性组合，即： A k：求积系数，与f(x)的形式无关，只与积分区间、积分节点有关==避免求f(x)的原函数 Xk：求积节点 求积余项（截断误差）： 问题： 如何评价一个求积公式的精度？ 如何确定求积系数 A k，使 In( f )精度越高？ 若xk可任意选择，取什么点好？ 3。常见的几种求积公式 中距公式： 几何意义 梯形公式： 几何意义： Simpson公式： 如何评价上述求积公式的优劣？（例6-1: p 189) 4。代数精度 定义：若 对所有不超过m次的多项式均精确成立，对至少一个m+1次多项式不精确成立，则称此求积公式：In( f )（数值积分）具有m次代数精度。 即 In(f) 对不高于m次的多项式 f(x) 无误差 代数精度越高，求积公式越好 此定义验证麻烦 定理：求积公式 具有m次代数精度的充要条件是：公式对 f(x)=xi (i=0,1,2,…,m)精确成立，对f(x)=xi+1 不精确成立 故有以下结论： 此方程组为范德蒙行列式，故方程组有且仅有一解向量，解之即可确定求积公式 例1：考查： 的代数精度 例2：确定求积公式 的系数Ai使之具有最高的代数精度 思考：试确定一个至少具有2次代数精度的求积公式： 积分公式的代数精度与积分节点个数的关系 n+1个互异的积分节点xk，k=0,1,…,n，则存在求积系数Ak，使得求积公式 具有至少 n 次代数精度 证明过程同于：前面充要条件的证明 6.2 Newton-Cotes公式 思想 插值型求积公式的代数精度与截断误差 Newton-Cotes公式（等距节点的求积公式） 截断误差 1。思想 已知一组节点：a≤x0 ≤x1 ≤… ≤ xn ≤b,及函数值f(x0),f(x1)…f(xn) 构造求积公式的思想：用被积函数f(x)的Lagrange插值函数Ln(x)?逼近f(x),作定积分?近似f(x)的定积分。 因： 故： 求积系数： 则数值积分公式为： 2。 插值型求积公式的代数精度与截断误差 1）截断误差： 2）代数精度： ∵：f (x)为任意次数小于等于n的多项式时，f(n+1)(x)=0 ∴：R(f)＝0，即In(f)=I(f)，求积公式精确成立 ∴：插值型求积公式至少具有n次代数精度 若 至少具有n次代数精度，则上式对任意n次多项式精确成立 而 是n次多项式 故： 所以： 此即插值型求积公式 定理：插值求积公式的代数精度： 求积公式： 至少具有n次代数精度的充要条件是该公式是插值型的，即： 例题：p194: 6-4~6-6 例2：在区间 [-h, h] 上取节点 -α,0,α ，确定 α 及求积系数，构造代数精度尽可能高的求积公式，并确定其代数精度。 3。等距节点的求积公式(Newton-Cotes公式） 1）构造：当a=x0 ≤x1 ≤… ≤ xn =b是区间[a,b]上的n等分点，插值求积公式：