ch24高斯型求积公式.pptVIP

第二章 数值积分 高斯型求积公式 问题: 是否有比等距节点的Newton-Cotes型求积公式 更高代数精度的求积公式? 最高能达到多大? 度 5 求积公式为 (1) 问题 定义 如何构造高斯求积公式？ 令(1)精确成立， 可以由上式求出 试构造下列积分的高斯求积公式： 例 而从分析高斯点的特性来构造高斯求积公式. 高斯点的基本特性 尽管高斯点的确定原则上可以化为代数问题，但是由于所归结的方程组是非线性的，而它的求解存在实质性的困难，所以我们要从研究高斯点的基本特性着手解决高斯公式的构造问题。 高斯点与正交多项式的零点 （2） 11 即有 故(2)成立. 充分性. 由(2)可得 证明 必要性. （3） 12 即 再注意到 知 从而由(3)有 13 14 Gauss型求积公式的构造方法 (1)求出区间[a,b]上权函数为 正交多项式pn+1(x) . (2)求出pn+1(x)的n个零点x0 , x1 , … xn 即为Gauss点. (3)计算积分系数 。 常见的正交多项式及高斯求积公式 勒让德多项式(Legendre) 切比雪夫多项式(Chebyshev) 拉盖尔多项式(Laguerre) 埃尔米特多项式 (Hermite ) 高斯-勒让德求积公式 2. Legendre多项式的性质: 19 这样构造出的一点高斯-勒让德求积公式 是中矩形公式. 20 三点高斯-勒让德公式的形式是 列出了高斯-勒让德求积公式的节点和系数. 从而得到两点高斯-勒让德求积公式 n xk Ak n xk Ak 1 0 2 6 ±0.9324695142 ±0.6612093865 ±0.2386191861 0.1713244924 0.3607615730 0.4679139346 2 ±0.5773502692 1 3 ±0.7745966692 0 0.5555555556 0.8888888889 7 ±0.9491079123 ±0.7415311856 ±0.4058451514 0 0.1294849662 0.2797053915 0.3818300505 0.4179591837 4 ±0.8611363116 ±0.3399810436 0.3478548451 0.6521451549 8 ±0.9602898565 ±0.7966664774 ±0.5255324099 ±0.1834346425 0.1012285363 0.2223810345 0.3137066459 0.3626837834 5 ±0.9061798459 ±0.5384693101 0 0.2369268851 0.4786286705 0.5688888889 高斯-切比雪夫求积公式 2. Chebyshev多项式的性质: 一般积分区间[a,b]的处理 高斯积分公式的数值稳定性 27 Gauss求积公式的余项： /* 设P为f 的过x0 … xn的插值多项式 */ 插值多项式的余项 Q：什么样的插值多项式在 x0 … xn 上有 2n+1 阶余项？ 28 A：Hermite 多项式！