数值积分与数值微分 高斯型求积公式.pptVIP

第2章 数值积分与数值微分 高精度的求积公式 举例（一） 高斯点与高斯公式 高斯点的确定 高斯点的确定 高斯公式代数精度 Gauss-Legendre 公式 G-L 公式的余项 几个简单的 G-L 公式 G-L公式的Gauss系数 更多 G-L 公式 一般区间上的 G-L 公式 G-L 公式的优缺点 例题1 例题2 2.5 数值微分 插值型的求导公式 两点公式 三点公式（等距） 例题2 例题3 * * 高斯型求积公式 考虑积分 将节点 x0 , … , xn 以及系数 A0 , … , An 都作为待定系数。令 f (x) = 1, x, x2, …, x2n+1 代入可求解，得到的公式具有2n+1 次代数精度。 能否利用 n+1 个节点 x0 ,… , xn 构造出具有 2n+1 次代数精度的求积公式 ？ 例：试确定 x0 , x1 以及系数 A0, A1，使得下面的求积公式具有尽可能高的代数精度。 解： 将 f (x) = 1, x, x2, x3 代入，使其精确成立得 解得 不是线性方程组，不易求解 定义 若存在节点 xi ?[a, b]及求积系数 Ai ，使得下面的求积公式具有 2n+1 次代数精度，则称节点 xi 为高斯点，Ai 为高斯系数，求积公式为高斯(Gauss)求积公式。 注：(1)Gauss求积公式仍然是插值型求积公式； (2)Gauss系数可通过Gauss点和Lagrange基函数得到； 定理 节点 xi (i = 0, 1, … , n) 是求积公式(2-30)的Gauss点的充要条件是：多项式 与任意次数不超过 n 的多项式 p(x) 正交，即 且高斯系数 Ai 为 其中 li 为以节点 xi 为节点的 Lagrange 基函数。 证明： “?” x0 … xn 为 Gauss 点, 则公式 至少有 2n+1 次代数精度。 对任意次数不大于n 的多项式 p (x)， p (x) w(x) 的次数不大于 2n+1，则代入公式应精确成立： 0 = 0 ? “?” 要证明 x0 … xn 为 Gauss 点，即要证公式对任意次数不大于2n+1 的多项式 p(x) 精确成立，即证明： 设 0 ? 定理 用 n+1 个点 x0 , … , xn 构造的插值型求积公式 的代数精度不超过 2n+1。 即Gauss公式是插值型求积公式中代数精度最高的。 设 f (x) ?C[-1, 1] ，考虑 Gauss型 求积公式 在 [-1, 1] 上的正交多项式为Legendre多项式 取其 n+1 个零点作为 Gauss 点，即可得 Gauss-Legendre 求积公式。 定理 设 f (x) ?C 2n+2[-1, 1] ，则 G-L求积公式的余项为 n = 0: Pn+1(x) = x, x0 = 0, A0 = 2 n = 1: Pn+1(x) = (3x2 - 1)/2, A0=A1=1 n = 2: Pn+1(x) = (5x3 - 3x)/2, 两点G-L公式 三点G-L公式 定理 G-L求积公式中的 Gauss点为 Pn+1(x) 的 n+1 个零点，Gauss系数为 (i = 0, 1, … , n) 当 n 3 时，可用数值方法计算 Pn+1(x) 的零点 (mygl.m) 0 0 0 ?0 ?0 ?06 5 0.2369269 0.4786287 0.5688889 ?0.9061798 ?0.5384693 0.0000000 5 4 0.3478548 0.6521452 ?0.8611363 ?0.3399810 4 3 0.5555556 0.8888889 ?0.7745967 0.0000000 3 2 1.0000000 ?0.5773503 2 1 2.0000000 0.0000000 1 0 Gauss系数 Gauss点 节点个数 n 设 f (x) ?C[a, b] 作变量替换 x = (b- a) t/2+(b + a)/2，则 t ? [-1, 1] 其中 xi 和 Ai 分别为 Gauss点 和 Gauss系数。 Ai 与前面求积方法的比较 复合梯形公式：用了 210+1 个节点达到 7 位有效数字 Romberg公式：用