Gauss 求积公式.ppt

6.5 Gauss 求积公式 1. Gauss 求积公式 设插值型的数值积分公式： 前面讲述的方法（lagrange插值型数值积分法）是事先给定积分节点xk 。例如Newton-Cotes公式把区间[a , b] 的等分点作为求积节点，这样所求积分公式的代数精度至多为n+1。 1. Gauss 求积公式 但是，考虑2n+2次多项式： 它只在节点xi（i = 0,1,…,n）处为零，在其它点处均大于零，所以 1. Gauss 求积公式 例 6.5.1 确定下列求积公式中的待定参数 A0 , A1 , x0 和 x1： 关于 Gauss 点的求法 , 有如下定理。 2. Gauss-Legendre 求积公式 这是因为: 常见Gauss-Legendre求积公式 常见Gauss-Legendre求积公式 常见Gauss-Legendre求积公式 1~5 个节点的Gauss-Legendre求积系数 1~5 个节点的Gauss-Legendre求积系数（真值） 续Gauss-Legendre求积系数（真值） 对于一般的区间[a, b]，可作坐标变换 例 6.5.3 利用三点Gauss-Legendre求积公式计算积分 3. 带权的Gauss 求积公式 从而有 Gauss-Chebyshev 求积公式 4. Gauss 求积公式的余项, 收敛性与稳定性 Gauss 求积公式求积系数有如下性质: 对于 Gauss 求积公式的收敛性, 有如下的定理 6.6 数值微分 本节讨论数值微分。即对于定义在区间[a, b] 上，由列表给出的函数 y = f (x)： 1. 插值型数值微分公式 由上述列表函数,我们可以建立f (x)的n次Lagrange插值多项式Ln(x), 根据 f (x) ? Ln (x), 可以得到插值型数值微分公式 由于Lagrange插值多项式Ln(x) 的余项为 但当计算插值节点处的导数时, 因ωn+1(xk) = 0, 这时数值导数的余项公式可以表示为 2. 两点数值微分公式 已知在两个插值节点 x0, x1上的函数值 f (x0), f (x1), 要计算 f (x) 在x0, x1处的导数. 由线性插值公式 由余项公式 (6.6.5) 3. 三点数值微分公式 从而由 若略去余项, 有三点数值微分公式 还可以建立高阶导数的数值微分公式. 由余项公式 (6.6.3) 4. 利用三次样条插值函数求数值导数 基本思想就是在区间[a , b] 上, 根据互异的节点 ?a = x0 x1 x2 . . . xn = b 及函数值yk = f (xk) (k = 0, 1, … , n) , 构造三次样条函数S(x), 于是 f (x) ? S (x). 以三弯矩公式为例, 有 这样 从而, 在插值节点处的导数有 用三次样条插值函数计算数值导数, 在一定条件下,（如f (x)有连续四阶导数）有如下误差估计式： 因此，用三次样条插值函数求数值导数比用插值多项式可靠性大． 具体解法如下： 首先构造一个三次样条插值函数， 然后对其求导，用计算出的导函数作为f(x)的近似值. 5. 数值微分的外推算法 由三点数值微分公式 ( 6.6.11 ) 由Taylor公式有 上两式相减整理后, 得 由外推算法, 选取 本章小结 数值积分的基本概念、思想与理论 数值积分公式、插值型数值积分公式、余项、代数精度、收敛阶 插值型数值积分及其数值稳定性 Newton-Cotes 公式，复化求积公式、变步长的求积公式、 Richardson外推算法与Romberg求积公式、Gauss 求积公式 数值微分的思想与各种方法 两点、三点数值微分公式、样条插值函数求数值导数、外推算法 本章结束 准备复习考试了！ 能否举个例子说明这样做,不总是合理的？ 其中 , 且依赖于x. 对式 (6.6.3) 两边求导得 (6.6.4) 误差估计式中含有不确定函数 ξ(x) 的导数 可记为 从而可以对节点处误差作出适当的估计. 下面, 我们在等距节点情况下讨论函数 f (x) 在插值节点处导数的求法. (6.6.4) 令 h = x1－x0 , 将上式求导得 可得 带余项的两点数值微分公式 若略去余项, 可得两点数值微分公式 截断误差为O(h) . 注：h = x1－x0 类似可得 已知 f (x) 在三点 处的函数值分