《高斯积分法讲义》PPT课件.ppt

《《高斯积分法讲义》PPT课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、4.5高斯求积公式14.5.1一般理论求积公式含有个待定参数当为等距节点时得到的插值求积公式其代数精度至少为次.如果适当选取有可能使求积公式具有次代数精度，这类求积公式称为高斯(Gauss)求积公式.2为具有一般性，研究带权积分这里为权函数，类似(1.3)，求积公式为（5.1）为不依赖于的求积系数.使（5.1）具有次代数精度.为求积节点，可适当选取定义4如果求积公式(5.1)具有次代数精度，则称其节点为高斯点，相应公式(5.1)称为高斯求积公式.3根据定义要使(5.1)具有次代数精度，只要对（5.2）当给定权函数，求出右

2、端积分，则可由(5.2)解得令(5.1)精确成立，即4例5（5.3）解令公式(5.3)对于准确成立，试构造下列积分的高斯求积公式：得（5.4）5由于利用(5.4)的第1式，可将第2式化为同样地，利用第2式化第3式，利用第3式化第4式，分别得从上面三个式子消去有6进一步整理得由此解出从而7这样，形如(5.3)的高斯公式是由于非线性方程组(5.2)较复杂，通常就很难求解.故一般不通过解方程(5.2)求，而从分析高斯点的特性来构造高斯求积公式.8定理5是高斯点的充分必要条件是以这些节点为零点的多项式与任何次数不超过的多项式带权

3、正交，（5.5）证明即插值型求积公式(5.1)的节点必要性.设则9是高斯点，因此，如果精确成立，因即有故(5.5)成立.则求积公式(5.1)对于充分性.用除，记商为，余式为，即,其中.对于由(5.5)可得（5.6）10由于求积公式(5.1)是插值型的，它对于是精确的，即再注意到知从而由(5.6)有11可见求积公式(5.1)对一切次数不超过的多项式均精确成立.因此，为高斯点.定理表明在上带权的次正交多项式的零点就是求积公式(5.1)的高斯点.有了求积节点，再利用对成立，的线性方程.解此方程则得则得到一组关于求积系数12下面

4、讨论高斯求积公式(5.1)的余项.利用在节点的埃尔米特插值于是也可直接由的插值多项式求出求积系数即13两端乘，并由到积分，则得（5.7）其中右端第一项积分对次多项式精确成立，故由于（5.8）由积分中值定理得(5.1)的余项为关于高斯求积公式的稳定性与收敛性，有：14定理6证明它是次多项式，因而是次多项式，注意到高斯求积公式(5.1)的求积系数全是正的.考察故高斯求积公式(5.1)对于它能准确成立，即有上式右端实际上即等于从而有15由本定理及定理2，则得推论定理7定理得证.高斯求积公式(5.1)是稳定的.设即则高斯求积公

5、式(5.1)收敛，164.5.2高斯-勒让德求积公式在高斯求积公式(5.1)中，由于勒让德多项式是区间上的正交多项式，因此，勒让德多项式的零点就是求积公式(5.9)的高斯点.形如(5.9)的高斯公式称为高斯-勒让德求积公式.区间为则得公式若取权函数（5.9）17令它对准确成立，即可定出这样构造出的一点高斯-勒让德求积公式为是中矩形公式.若取的零点做节点构造求积公式再取的两个零点构造求积公式18令它对都准确成立，有由此解出三点高斯-勒让德公式的形式是表4-7列出高斯-勒让德求积公式(5.9)的节点和系数.从而得到两点高斯-

6、勒让德求积公式1920由(5.8)式，这里是最高项系数为1的勒让德多项式.由第3章(2.6)及(2.7)公式(5.9)的余项21得（5.10）当时，有它比区间上辛普森公式的余项还小，且比辛普森公式少算一个函数值.当积分区间不是，而是一般的区间时，只要做变换22可将化为,（5.10）对等式右端的积分即可使用高斯-勒让德求积公式.这时23例6用4点()的高斯-勒让德求积公式计算解先将区间化为，根据表4-7中的节点及系数值可求得由(5.11)有244.5.3高斯-切比雪夫求积公式若且取权函数则所建立的高斯公式为（5.12）称为

7、高斯-切比雪夫求积公式.25由于区间上关于权函数的正交多项式是切比雪夫多项式，因此求积公式(5.12)的高斯点是次切比雪夫多项式的零点，即为(5.12)的系数使用时将个节点公式改为个节点，（5.13）于是高斯-切比雪夫求积公式写成26由(5.9)，余项带权的高斯求积公式可用于计算奇异积分.（5.14）27例7用5点()的高斯-切比雪夫求积公式计算积分解当时由公式(5.13)由(5.14)式，误差这里可得284.6数值微分数值微分就是用函数值的线性组合近似函数在某点的导数值.294.6.1中点方法与误差分析按导数定义可以简

8、单地用差商近似导数，这样立即得到几种数值微分公式其中为一增量，称为步长.（6.1）30后一种数值微分方法称为中点方法，它其实是前两种方法的算术平均.但它的误差阶却由提高到较为常用的是中点公式.为利用中点公式计算导数的近似值，首先必须选取合适的步长，为此需要进行误差分析.分别将在处做泰勒展开有31代入中点公式得从截断误