m教材《数值分析》陈晓江主编武汉理工大学出版社Y参考书目《数值分析》李庆扬、王能超、易大义编清华大学出版社施普林格出版社《科学和工程计算基础》施妙根、顾丽珍编清华大学出版社上页下页返圆
上页 下页 返回 教材 《数值分析》 陈晓江主编 武汉理工大学出版社 参考书目 《数值分析》 李庆扬、王能超、易大义编 清华大学出版社 施普林格出版社 《科学和工程计算基础》 施妙根、顾丽珍编 清华大学出版社
第四章数值积分与数值微分第一节问题的提出第二节机械求积法和代数精度第三节 牛顿一柯特斯求积公式第四节复化求积公式第五节龙贝格求积公式第六节#高斯求积公式第七节数值微分上页下页返圆
上页 下页 返回 第四章 数值积分与数值微分 第一节 问题的提出 第三节 牛顿—柯特斯求积公式 第四节 复化求积公式 第五节 龙贝格求积公式 第六节 高斯求积公式 第七节 数值微分 第二节 机械求积法和代数精度
86高斯求积公式一、一般理论考察两个点的求积公式：L f(x)dx ~ A(xo)+ Af(x) f(x)dx~ f(-1)+f(I) 具有1次代数精度,它可以是：也可以是：(x)dx~(-)+(具有3次代数精度它表明，有2个节点的求积公式，最高可具有3次代数精度。ZA(x)一般的，有n+1个节点的求积公求f(x)dx~)k=0最高可以具有2n+1次代数精度。上页这类求积公式称为高斯（Gauss）求积公式下页返圆
上页 下页 返回 §6 高斯求积公式 考察两个点的求积公式： 1 1 0 0 1 1 f (x)dx A f (x ) A f (x ) 它可以是： 1 1 f (x)dx f ( 1) f (1) 具有1次代数精度， 也可以是： 1 1 ) 3 1 ) ( 3 1 f (x)dx f ( f 具有3次代数精度。 它表明,有2个节点的求积公式,最高可具有3次代数精度。 ( ) ( ) 0 k b a n k f x dx Ak f x 最高可以具有2n+1次代数精度。 这类求积公式称为高斯(Gauss)求积公式。 一般的，有n+1个节点的求积公式 一、一般理论
f'f(x)dx~ZArf(x)因为求积公式k=0含有2n+2个待定参数xk、Ak(k=0, 1, ..., n)当x为给定节点时，插值型求积公式至少有n次代数精度，如果适当选取xk(k=0，1，…，n)，有可能使求积公式具有2n+1次代数精度。为使问题更具一般性，我们研究带权积分I = f" f(x)p(x)dx这里p(x)为权函数，类似机械求积公式，它的求积公式为J' f(x)p(x)dx ~ZAxf(x)k=0Ak(k=0，1，...，n)为不依赖于f(x)的求积系数，xr(k=0，1，.….，n)为求积节点,上页可适当选取x及Ak（k=0，1，，n)使求积公式下页具有2n+1次代数精度返圆
上页 下页 返回 Ak (k＝0，1，.，n)为不依赖于f (x)的求积系数， b a n k k xk f x x A f 0 ( )d ( ) b a I f (x)(x)dx n k k k b a f x x x A f x 0 ( )( )d ( ) 因为求积公式 含有2n+2个待定参数xk、Ak (k＝0，1，.，n)． 当 xk 为给定节点时，插值型求积公式至少有n次代数精度， 如果适当选取 xk (k＝0，1，.，n)，有可能使求积公式 具有 2n+1 次代数精度。 为使问题更具一般性，我们研究带权积分 这里 (x)为权函数，类似机械求积公式，它的求积公式为 xk (k＝0，1，.，n)为求积节点， 可适当选取 xk 及 Ak (k＝0，1，.，n)使求积公式 具有2n+1次代数精度．
如果求积公式具有2n+1次代数精度，则称其节点xk(k=0，1，，n)为高斯点，相应的公式称为高斯求积公式·要使求积公式具有2n+1次代数精度，只要取f(x）=xm对m=0，1，，2n+1，求积公式精确成立，则得ZAxx" = I"' x" p(x)dxm=0,1,..,2n+1.当给定权函数P(x)，求出右端积分，则可由上面非线性方程组解得Ak及x(k=0，1，，n)从教材p140看到，求解上面非线性方程组较复杂，通常n>2就很难求解故一般不通过解方程组求xr及 Ak(k=0, 1, ..., n),上页下页而从分析高斯点的特性来构造高斯求积公式返圆
上页 下页 返回 b a m n k m Ak xk x (x)dx m 0, 1, , 2n 1. 0 而从分析高斯点的特性来构造高斯求积公式． 如果求积公式具有2n+1次代数精度， 则称其节点 xk (k＝0，1，.，n)为高斯点， 相应的公式称为高斯求积公式 . 要使求积公式具有2n+1次代数精度， 只要取f(x)＝x m ， 对m＝0，1，.，2n+1，求积公式精确成立，则得 则可由上面非线性方程组解得 Ak 及 xk (k＝0，1，.，n)． 当给定权函数 (x) ，求出右端积分， 从教材p140看到,求解上面非线性方程组较复杂， 通常n≥2就很难求解． 故一般不通过解方程组求 xk 及 Ak (k＝0，1，.，n)