3.数值积分-2说课.ppt

4.5.2 高斯-勒让德求积公式 勒让德多项式 勒让德多项式 N=0~5时的勒让德多项式 例 用两点(n=1)高斯-勒让德求积公式计算 积分系数： A0=1，A1=1 高斯点(积分结点)： t0=0.5773502692， t1=-0.5773502692 作积分变换： b=1, a=0 积分结点： x0=0.7886751346，x1=0.2113248654 高斯-勒让德求积公式结果 4.6 数 值 微 分 4.6.2 插值型的求导公式 4.6.3 利用数值积分求导 4.6.4 三次样条求导 本章小结 数值积分的基本思想 插值型求积公式及其余项 牛顿-科特斯公式 复合求积公式 龙贝格求积公式 高斯求积公式 数值微分 对于列表函数 运用插值原理，可以建立插值多项式 作为它的近似. 由于多项式的求导比较容易，我们取 的值作为 的近似值，这样建立的数值公式 （6.3） 统称插值型的求导公式. 即使 与 的值相差不多， 与导数的真值 仍然可能差别很大. 导数的近似值 因而在使用求导公式(6.3)时应特别注意误差的分析. 依据插值余项定理，求导公式(6.3)的余项为 式中 但如果限定求某个节点 上的导数值，那么第二项中 由于 是 的未知函数，所以对随意给出的点 ， 误差是无法预估的. 因式 变为零，这时余项公式为 （6.4） 下面仅考察节点处的导数值并假定所给节点是等距的. 1. 两点公式 设已给出两个节点 上的函数值 对上式两端求导，记 ，有 做线性插值 于是有下列求导公式： 利用余项公式(6.4)知，带余项的两点公式是 2. 三点公式 设已给出三个节点 上的函数值， 做二次插值 令 上式可表示为 两端对 求导，有 （6.5） 式中撇号（′）表示对变量 求导数. 分别取 得到三种三点公式： 带余项的三点求导公式为 （6.6） 其中的公式(6.6)是中点公式. 它比其余两个三点公式少用 了一个函数值. 用插值多项式 作为 的近似函数，还可以建立 高阶数值微分公式： 例如，将式(6.5)再对 求导一次，有 于是有 而带余项的二阶三点公式如下： （6.7） 微分是积分的逆运算，因此可利用数值积分的方法来 计算数值微分. 设 是一个充分光滑的函数， （6.8） 对上式右边积分采用不同的求积公式就可得到不同的数值 微分公式. 则有 例如，用中矩形公式(1.2)，则得 从而得到中点微分公式 若对(6.8)右端积分用辛普森求积公式，则有 略去上式余项，并记 的近似值为 则得到 辛普森数值微分公式 这是关于 的 个方程组， 已知， （6.9） 若 则可得 这是关于 的三对角方程组，且系数矩阵为严格 对角占优的，可用追赶法求解(见第5章5.4节). 如果端点导数值不知道，那么对(6.9)中第1个和第 个方程可分别用 及 的中点微分公式近似， 然后求 即为 的近似值. 即取 例8 给定 的一张数据表(表4-9左部)， 并给定 及 的值(见表4-9). 解 解之得 利用辛普森数值微分公式求 在 上的一阶导数. 结果见表4-9. 根据(6.9)有 三次样条函数 与 ，不但函数值很接近，而且 导数值也很接近，并有 （6.10） 因此利用三次样条函数 接得到 根据第2章(7.8)，(7.9)可求得 这里 为一阶均差. 其误差由(6.10)可得 4.5.1 一般理论 求积公式 含有 个待定参数 当 为等距节点时得到的插值求积公式其代数精度至少为 次. 如果适当选取 有可能使求积公式 具有 次代数精度，这类求积公式称为高斯(Gauss) 求积公式. 可以出现：x0 ≠ a