第四章1936最全版.pptVIP

* Gauss 型求积公式 一般情形：考虑机械带权求积公式 定义：若存节点在 xi ?[a, b] 及系数 Ai ，使得上面的求积公式具有 2n+1 次代数精度，则称节点 xi 为高斯点，Ai 为高斯系数，求积公式为 高斯型求积公式 性质：上面的求积公式至多具有 2n+1 次代数精度 （将 代入验证即可） Gauss 求积公式在所有机械求积公式中代数精度最高 KUHKUK * Gauss 点 问题：如何计算 Gauss 点 xi 和 高斯系数 Ai 法一：解非线性方程组 太困难! ? 法二：分开计算 先确定 Gauss 点 再通过解线性方程组计算 Gauss 系数 KUHKUK * Gauss 点 定理：上面的插值型求积公式中的节点 xi (i = 0, 1, … , n) 是 Gauss点的充要条件是：多项式 与任意次数不超过 n 的多项式 p(x) 都关于权函数 ?(x) 正交，即 证明: 板书 推论：设 p0(x), p1(x), ?, pn(x) , ? 是 [a, b] 上带权 ?(x) 正交的多项式族，则 Gauss 点即为 pn+1(x) 的零点！ KUHKUK * 计算 Gauss 点的一般方法 求出 ?n+1(x) 的表达式 计算其零点 与 1, x, x2, ..., xn 带权正交 特殊情形： (1) [a, b]=[-1, 1], ?(x)=1， 则 Gauss 点即为 Legendre 多项式的零点 (2) [a, b]=[-1, 1], 则 Gauss 点即为 Chebyshev 多项式的零点 KUHKUK * Gauss 系数的计算 例：试确定节点 xi 和系数 Ai ，使得下面的求积公式具有尽可能高的代数精度。 将 f (x) = 1, x, x2, …, xn 代入，解方程 或利用 Lagrange 基函数 KUHKUK * 余项公式 设 p2n+1(x) 是 f(x) 在节点 x0, x1, ?, xn 上的 2n+1 次 Hermite 插值多项式, 即 KUHKUK * 收敛性与稳定性 可以证明：当 a, b 为有限数，且 f (x) ?C[a, b] 时 Gauss 型公式是收敛的 令 Gauss 型公式是稳定的 KUHKUK * Gauss 公式与正交多项式 积分区间: [-1, 1]，权函数： ?(x) = 1 利用正交多项式构造 Gauss 求积公式 积分区间: [-1, 1]，权函数： Gauss-Legendre 求积公式 Gauss-Chebyshev 求积公式 KUHKUK * Gauss-Legendre 求积公式 积分区间: [-1, 1]， 权函数: ?(x) = 1 Gauss 点 = Legendre 多项式 pn+1(x) 的零点 G-L 求积公式: KUHKUK * 低阶 G-L 公式 n =0 时, G-L 求积公式: Gauss 点: 将 f (x)＝1 代入求出 A0 n =1 时, 两点 G-L 求积公式: Gauss 点: 将 f (x)＝1, x 代入 求出 A0 , A1 KUHKUK * 低阶 G-L 公式 n =2 时, 三点 G-L 求积公式: Gauss 点: KUHKUK * 更多 G-L 公式 当 n 3时，可用数值方法计算 Pn+1(x) 的零点 (教材122页) n 节点个数 Gauss点 Gauss系数 0 1 0.0000000 2.0000000 1 2 ?0.5773503 1.0000000 2 3 ?0.7745967 0.0000000 0.5555556 0.8888889 3 4 ?0.8611363 ?0.3399810 0.3478548 0.6521452 4 5 ?0.9061798 ?0.5384693 0.0000000 0.2369269 0.4786287 0.5688889 5 6 ?0 ?0 ?0 0 0 0KUHKUK * G-L 公式余项 余项公式 ? ? (-1, 1) KUHKUK * 一般区间上的 G-L 公式 做变量代换 积分区间: [a, b]， 权函数: ?(x) = 1 KUHKUK * G-L公式举例 例：用四点G-L公式 (n=3) 计算定积分 解： 令 KUHKUK MKK