三点Gauss公式

三点gauss-chebyshev 求积公式
三点Gauss-Chebyshev求积公式是一种数值积分方法，用于计算函数在[-1, 1]区间上的积分。

它基于Chebyshev多项式的零点，并使用三个等距节点进行插值和积分计算。

具体的求积公式如下：
积分近似值 = (b - a) / 3 * [f(a) + 4f((a + b) / 2) + f(b)]
其中，a和b是积分区间的起始和终止点，f(x)是需要计算的函数。

请注意，这个求积公式特别适用于在[-1, 1]区间上具有较光滑性质的函数。

当函数在该区间上不是很光滑或具有较大的变化时，可能需要采用更精确的高阶Gauss-Chebyshev求积公式或其他数值积分方法。

对于其他节点数，可以使用更高阶的Gauss-Chebyshev求积公式，如五点、七点或更多，以提高积分的精度。

这些公式的推导和使用方法是类似的，只需根据具体的节点数和相应的节点位置进行调整。

需要注意的是，当使用数值积分方法时，应根据具体应用需求和被积函数的特点选择合适的公式和节点，并对数值误差进行评估和控制，以获得准确的积分结果。

一、 引言介绍高斯型求积公式，并使用其求积分⎰=1sin I xdx 。

要求：数值实验结果要体现出随高斯点的增加误差的变化。

我们知道,求积公式⎰∑=≈bani i ix f Adx x f 0)()( （1.1）含有22+n 个待定常数i x 及),,2,1,0(n i A i =,如果它具有n 次代数精确度,则它应使1+m 个方程mk dx x x A bakni ki i ,,2,1,0,==⎰∑= (1.2)精确成立。

作为插值型求积公式(1.1)它至少具有n 次代数精确度；另一方面,令)())(()(101n n x x x x x x x ---=+ ω,则对22+n 次多项式)()(21x x f n +=ω而言，(7.5.1)右端为零,而左端严格大于零,即(7.5.1)式对22+n 次多项式)(21x n +ω不准确成立。

但要确定方程组(7.5.2)中的22+n 个待定常数i x 与i A ,最多需要给出22+n 个独立条件,所以m最大取12+n 。

因此,插值型求积公式(1.1)的代数精确度最小是n ,最大是12+n ．由此可见，高斯公式的代数精度比牛顿-科特斯公式高，求解高斯求积公式的关键就是解出上述2n+2个待定常数。

为解决上述问题，首先要先给出三个定理：定理一：以n x x x ,,,10 为节点的插值型求积公式(7.5.1)具有12+n 次代数精确度的充要条件是以这些节点为零点的多项式)())(()(101n n x x x x x x x ---=+ ω与任意次数不超过n 的多项式)(x P 均在区间],[b a 上正交,即⎰=+ban dx x x P 0)()(1ω (1.3)定理二：高斯公式(1.1)的求积系数k A 全为正,且nk dx x l dx x l A bak bak k ,1,0,)()(2===⎰⎰(1.4)定理三：对于高斯公式(1.1),其余项为dxx fn f R ban n ⎰+++=)()()!22(1)(21)22(ωη (1.5)其中).())(()(],,[101n n x x x x x x x b a ---=∈+ ωη证明 以n x x x ,,,10 为节点构造)(x f 的埃尔米特插值多项式)(x H),()(i i x f x H = ni x f x H i i ,1,0),()(='='因为)(x H 是12+n 次多项式,而它的余项是)()()!22(1)()(21)22(x fn x H x f n n +++=-ωξ所以高斯公式(7.5.1)对)(x H 能准确成立,即∑∑⎰====ni i in i iibax f Ax H A dx x H 0)()()(从而dxx fn dxx H dx x f x f A dx x f f R n ban babani i i ba)()()!22(1)()()()()(21)22(0++=⎰⎰⎰∑⎰+=-=-=ωξ若)()22(x fn +在区间],[b a 上连续,由于)(21x n +ω在],[b a 上不变号,故应用积分中值定理可得],[,)()()!22(1)(21)22(b a dx x fn f R ban n ∈+=⎰++ηωη上述定理说明,与牛顿—科兹公式进行比较,高斯公式不但具有高精度,而且它还是数值稳定的,但是节点和求积系数的计算比较麻烦。