Chapter4_3_Gauss积分讲述.ppt

§3 Gauss积分/*Gauss Integration*/ 一、引言 二、Gauss型求积公式 三、求积公式的构造与应用 求积公式 当求积系数Ak(k=1,2,…,n)、求积节点xk(k=1,2,…,n) 都可以自由选取时,其代数精确度最高可以达到多 少次? 当求积系数Ak(k=1,2,…,n)和求积节点 xk(k=1,2,…,n)都可以自由选取时,n点的求积公式 (1)的代数精确度最高可以达到2n-1次. (1) 一、引言 引理1 假设求积公式(1)具有m次代数精确度,即对任意 的m次代数多项式 求积公式(1)精确成立. 即 (2) 证明 于是 记 由于系数 的任意性 则(2)式成为 (3) 故使(3)式成为恒等式的充要条件是 (4) (4)式的待定系数有2n个,所以确定待定系数的 独立条件至多给出2n个,从而可知m至多为2n-1 n点的求积公式(1)具有2n-1次代数精确度(或称为具有最高的代数精确度)时,称为Gauss型求积公式. 定义1 Gauss型求积公式的求积节点 , 称为Gauss点,它们可以通过求区间[a,b]上带权?(x)的n次正交多项式gn(x)的n个根获得.所以先介绍正交多项式及其性质.然后讨论Gauss型求积公式的构造,等等. 二、Gauss型求积公式 n点的求积公式若具有最高的代数精确度,或具有 2n-1次的代数精确度成为Gauss型求积公式. 求积节点 和求积系数 如何选取,才能 使之成为Gauss型求积公式? 定理1 求积公式(1)中的n个求积节点 ,取在区 间[a,b]上带权函数?(x)的n次正交多项式 gn(x)的n 个根成为Gauss型求积公式. [a,b]上带权函数?(x)的n次正交多项式gn(x)的n个 根记为 在两边同乘 ?(x),并从a到b积分. 由正交多项式的 性质可知,含 项的积分为零. ,首项系数为dn ,由定义有 证明 设 因此, 所以 至少具有n-1次代数精确度, 又可知 即 注意到当 作为插值节点时建立的n点插值 求积公式 而 所以 综合可知当 时,求积公式 成立. 证明略. 之值近似积分值, In(f) 有下面的误差估计 用n点Gauss求积公式 定理2 则Gauss型求积公式(1)的误差 估计R(?,f)为 其中 在稍后讨论Gauss积分值数列的收敛性等问题时, 需要用到Gauss型求积公式的求积系数大于零的 结论. 这里用下面的定理给出 . , 定理3 Gauss型求积公式的求积系数 大于零 证明 这里 为区间 [a,b]上带权函数?(x) 的n次正交多项式 gn(x)的n个根 显然 由于 ,所以对f(x)求积公式(1)精确成 立, 即 因为 所以 存在 ,使得 Gauss型求积公式的收敛性 若f(x)?C[a,b],则Gauss型求积公式所求积分值 序列{In(f)}收敛于积分值I(f). 因为f(x)?[a,b], 由Weierstrass定理对任意的?10, 定理4 即 证明 对任意的x?[a,b]成立. 由于公式(1)为Gauss型求积公式时具有2n-1次代 数精确度,取 N(m+1)/2 故当 nN时,即m2N-12n-1时, 有 成立. 于是 因为 记 故 即 三、Gauss型求积公式的构造与应用 定理实际上给出了构造Gauss型求积公式的一种 方法. 先求出区间[a,b]上带权函数?(x)的n次正交多项式 然后求出多项式的n个根,从而获得了求积节点。 最后代入下式求求积系数。 Gauss-Lagendre求积公式的积分区间为[-1,1],权 函数?(x)=1. Gauss-Legendre求积公式的积分区间为[0,+?),权 函数 Gauss-Hermitte求积公式的积分区间为(-?,+?), 权函数 表4-5为Gauss-Legendre求积公式的求积系数和 求积节点 注意： 需要指出的是,不同求积公式的求积系数与求积 节点,积分区间和权函数是不同的. (1) 例如,n点的Gauss-Chebyshev求积公式,它的n个 求积系数: 而n个求积节点则为: . 正是因为Gauss-Chebyshev求积公式的求积系数 相同,所以在实际计算时,乘