数值积分与数值微分ch04c.pptxVIP

1 第四章数值积分与数值微分 计算方法 —— Gauss 求积公式 2 本讲内容 一般理论: 公式, 余项, 收敛性, 稳定性 Gauss-Legendre 求积公式 Gauss-Chebyshev 求积公式 无限区间的 Gauss 求积公式 Gauss 求积公式 3 Gauss 型求积公式 考虑求积公式 含 2n+2 个参数 (节点与系数), 为了使该公式具有尽可能高的代数精度, 可将 f (x) = 1, x, x2, …, x2n+1 代入公式, 使其精确成立, 则可构造出代数精度至少为 2n+1 的求积公式! 4 举例 例：试确定节点 xi 和系数 Ai ，使得下面的求积公式具有尽可能高的代数精度，并求出此求积公式的代数精度。 易验证该公式对 f (x)＝x4 不精确成立， 所以此求积公式具有 3 次代数精度。 非线性方程组 求解较困难 5 Gauss 型求积公式 一般情形: 考虑机械带权求积公式 6 Gauss 点 如何计算Gauss点 xi 和 高斯系数 Ai 法一: 解非线性方程组 7 Gauss 点 定理：节点 xi (i = 0, 1, … , n) 是 Gauss点的充要条件是：多项式 与任意次数不超过 n 的多项式 p(x) 关于权函数 (x) 正交，即 且高斯系数 Ai 为 其中 li(x) 为以 xi 为节点的 Lagrange 基函数。 证明: 板书 8 Gauss 点 证明： x0 … xn 为 Gauss 点 设 p(x)Hn ，则 p(x)n+1(x) H2n+1 “” 设 要证 xi 为 Gauss 点，即公式对  p(x) H2n+1精确成立 “” p(x), r(x)Hn ① 正交性 ② 公式是插值型的 将 li(x) 为代入即可得 Ai 的表达式。 9 Gauss 公式 Gauss 点的计算 求出 n+1(x) 的表达式 计算其零点 Gauss 系数的计算 将 f (x) = 1, x, x2, …, xn 代入，解方程 或利用 Lagrange 基函数 10 举例 例：试确定节点 xi 和系数 Ai ，使得下面的求积公式具有尽可能高的代数精度。 11 举例 将 f (x)＝1, x 代入求积公式，使其精确成立，可得 求积公式为 12 余项 13 收敛性与稳定性 可以证明：当 a, b 为有限数，且 f (x) C[a, b] 时 Gauss 型公式是收敛的 14 G-L 公式 积分区间: [-1, 1] 权函数: (x) = 1 Gauss-Legendre 求积公式 15 简单 G-L 公式 n =0 时, G-L 求积公式: Gauss 点: 将 f (x)＝1 代入求出 A0 16 简单 G-L 公式 n =2 时, 三点 G-L 求积公式: Gauss 点: 17 更多 G-L 公式 当 n 3 时，可用数值方法计算 Pn+1(x) 的零点 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 018 余项 余项:   (-1, 1) 19 G-L 公式 一般区间上的 G-L 求积公式 20 举例 例：用四点G-L公式 (n=3) 计算定积分 21 G-C 公式 积分区间: [-1, 1] 权函数: Gauss-Chebyshev 求积公式 Gauss 点 = Chebyshev 多项式 Tn+1(x) 的零点 G-C 求积公式: 22 G-C 公式 Tn+1(x) 的零点 (i = 0, 1, … , n) Gauss 系数 (i = 0, 1, … , n) G-C 求积公式: 余项:   (-1, 1) 23 简单 G-C 公式 n = 0 n = 1 n = 2 两点 G-C 公式 三点 G-C 公式 24 举例 例：用五点G-C公式计算定积分 25 Gauss-Laguerre 公式