4-5-1高斯求积公式 数值微分.ppt

4.7.4 Gauss-Chebyshelv quadrature formula According to the error estimation of the Gauss- Type formula,we have: 4.1 Numerical Differentiation Balance between round-off truncated error 4.2 Richardson’s Extrapolation(1927) Numerical Differentiation Revisit-------Using Extrapolation Method 构造高斯公式的一般方法： 1、构造正交多项式，继而求其零点，再按插值求积公式获得高斯公式； 2、待定系数法 此外，还可涉及到无穷区间上的广义积分等。例如： ---拉盖尔-高斯积分 举例 要构造下列形式的高斯公式 解 则其代数精度应为 即 求解…?! 定理（稳定性） 高斯求积公式的求积系数Ak0. 证明：事实上 这表明高斯求积法是稳定的。 关于积分余项和收敛性有： 积分余项： 收敛性： 设f(x) ∈C[a,b],则有： However , (i) There is no error estimation; (ii) Are there any other numerical methods for ND? How to construct them what about error? To answer these questions, we observe first: Error Bound Called forward difference central difference formula. There are also backward difference formulas. Five-point formula below can be obtained similarly: It then be called compact form. For higher order derivatives, it can also be obtained by interpolation like to the 1st order derivative using more points. Alternately, we can obtain the formulas which are algebraically tedious by Taylor’s expansion such as: Cf. the results obtained by the two methods. Richardson’s Extrapolation is used to generate high-accuracy results while using low-accuracy formulas. Then combined with the formula of N2(h) to eliminate the h2 term, we obtain: Which posses higher order truncated error! 一、高斯点 定义：高斯公式 机械求积公式 含有2n+2个待定参数 若适当选择这些参数使求积公式具有尽量高次（2n+1次?!）代数精度，则这类公式称为高斯公式。 (4.1) 定义：高斯公式的求积节点称为高斯点。 请回顾: 以前学过的梯形公式、辛甫生公式、柯 特斯公式、中矩形公式是高斯公式吗？ 除中矩形公式外都不是！ 注：机械型高斯求积公式一定是插值求积公式。 举例 求 [a,b]上的两点高斯公式。 解 设两点高斯公式为 这是关于四个未知数的非线性方程组，是否有解？一般难于求解… 要求其代数精度最高，四个未知数，可列出4个方程: 高斯点具有以下性质： 定理 插值型求积公式(4.1)成为Gauss求积公式的充要条件： 求积节点 为n+1次正交多项式的零点。 如何求高斯公式? 正交多项式概述： 首先证明对于任给节点 x0， x1，… ， xn，均存在某个次数为2n+2的多项式f(x),机械型求积公式不能精确成立，即其最高代数精度不能达到2n+2 。如取： 证明 则有： 设求积节点 为n+1次正交多项式ωn+1(x) 的零点。 现证充分性。即 求积公式是高斯