第三章 第六节 Gauss型求积公式.ppt

* * 第六节 Gauss型求积公式 一、问题的提出 三、 常用的Gauss型求积公式 二、 Gauss型求积公式的构造 四 Gauss型求积公式的余项 一、问题的提出 为了一般性，考虑积分 （6.1) 何普通积分 ，都可写成 其中， 为权函数，当 时，即为普通积分。对任 （6.2） 从而都可化成（3.1）形式的积分。 对于（3.1）中的 ，用 个不等距节点 作插值，则可得到 其中 它不依赖于函数 ,但可以依赖于权函数 。由于代数精度的次数是积分公式近似程度的一种量度，所以对于包含同样节点数的积分公式，自然希望采用一种准确度次数比较高的积分公式。因为这样可以用同等的代价，获得较高近似度的计算结果。于是提出如下问题： 能否选择 个节点 和 个系数 ，使得求积公式（6.2）具有尽可能高的代数精度？ 为了回答这个问题，先固定 值，并设（6.2）的代数精度为m（待定），即设（6.2）对所有的 次多项式 是准确的。于是有 即 （6.3） 令 则（6.3）成为 （6.4） 由于系数 是任意的，故使（6.4）成立的充要条件是 （6.5） 这是一个具有 个未知数， 个方程的方程组 ， 当 时，方程组（6.5）是可解的。 因此确实可以 找到一组 和 使求积公式（6.2）的代数 精度从 次提高到 次。 定义2 若对于次数 的多项式 ，关系式 （6.6） 为恒等式，则称（6.6）为Gauss型求积公式，相应节点 称为Gauss型点。 显然，由方程组（6.5）就可以把Gauss型求积公式的节点和系数求出来。但在一般情况下，要解这个方程组是相当困难的，下面我们给出一种行之有效的方法。 二、 Gauss型求积公式的构造 先从一个简单例子入手 对任何三次多项式 精确成立。 解： 以 除 得 于是， （6.7） 例4???????? 取 ，积分区间 为，求 和 使 由于 而 从而有 故要使 必须有 （6.8） 也就是说，当节点满足条件（6.8）时，能使（6.7）式 精 确成立。 由于 是任意三次多项式 ， 所以 是任意一次多项式，据正交多项式的性质知，任意一个一次多项式都与二 次Legendre多项式在 上带权正交，又因为 是最高次项系数为1的多项式，故（3.8）中的 应取为 就是二次Legendre多项式的根，即 求出节点后，可利用求积公式（6.7）确定 。为此， 取 则有 解得 故求积公式（3.7）变成 此式对于任何不高于3次的多项式 都准确成立。这 就是 时的Gauss---Legendre求积公式。 下 面 讨 论 一 般 情 形 。即 求 节 点 和 系数 ，使 求 积 公 式 对任何不高于 次的多项式 都准确。 令 取 ，用 除 得 于是 由于 而 从而有 故要使 必须有 即任意 次多项式 与 次多项式 在 上关于权函数 正交。因此有 定理3 Gauss求积公式(6.6)的n个节点