《数值分析gauss积分》由会员分享，可在线阅读，更多相关《数值分析gauss积分（38页珍藏版）》请在金锄头文库上搜索。
1、数值分析数值分析 前面介绍的前面介绍的 n 1个节点的个节点的 newton -cotes求积公式，求积公式，其特征是节点是等距的。这种特点使得求积公式便于其特征是节点是等距的。这种特点使得求积公式便于构造，复化求积公式易于形成。但同时也限制了公式构造，复化求积公式易于形成。但同时也限制了公式的精度。的精度。 n是偶数时，是偶数时，代数精度为代数精度为n 1， n是奇数时，是奇数时，代数精度为代数精度为n 。 我们知道我们知道 n 1个节点的插值型求积公式的代数精个节点的插值型求积公式的代数精确度不低于确度不低于n 。设想：设想：能不能在区间能不能在区间a,b上适当选择上适当选择n 1个节点个
2、节点 x 0x1,x2,xn ,使插值求积公式使插值求积公式的代数精的代数精度高于度高于n？ 答案是肯定的，适当选择节点，可使公式的精度答案是肯定的，适当选择节点，可使公式的精度最高达到最高达到2n 1，这就是本节所要介绍的高斯求积公式。，这就是本节所要介绍的高斯求积公式。第四节第四节 高斯高斯(gauss)(gauss)求积公求积公式式数值分析数值分析1数值分析数值分析考虑更一般形式的数值积分问题考虑更一般形式的数值积分问题定义：定义：若求积公式若求积公式 对一切对一切不高于不高于m次的多项式次的多项式p(x)都等号成立，即都等号成立，即r(p)=0;=0;而对而对于某个于某个m 1 1次多
3、项式等号不成立，则称此求积公式的次多项式等号不成立，则称此求积公式的代数精度为代数精度为m. .一、构造高斯型求积公式的基本原理和方法一、构造高斯型求积公式的基本原理和方法数值分析数值分析2数值分析数值分析定理定理1：设节点设节点x0, x1,xna,b，则求积公式，则求积公式 的的代数精度最高为代数精度最高为2n 1次。次。 分别取分别取 f(x)=1, x，x2，.xr 代入公式，并让其成为代入公式，并让其成为等式，得：等式，得： a0 a1 an =ab1dx.= b-ax0 a0 x1 a1 xn an =abxdx.= （b2-a 2)/2 .x0 ra0 x
4、1 ra1 xn ran =abxr dxr =(br 1-a r 1)/(r 1)数值分析数值分析3数值分析数值分析 事实上事实上,取取 2n 2次多项式次多项式g(x)=(x-x0)2(x-x1)2.(x-xn)2 代入求积公式代入求积公式,这里这里 x0, x1,xn是节点，是节点，有有左左 右右,故等式不成立故等式不成立,求积公式求积公式的的代数精度最高为代数精度最高为2n 1次。次。 证毕证毕. 上式共有上式共有 r 1个个 等式，等式，2n 2个待定系数个待定系数(变元变元),要想要想如如上方程组有唯一解，应有方程的个数等于变元的个数上方程组有唯一解，应有方程的个数等于变元的
5、个数,即即 r 1=2n 2, 这样导出求积公式的代数精度至少是这样导出求积公式的代数精度至少是2 n 1,下面证明代数精度只能是下面证明代数精度只能是2n 1. 数值分析数值分析4数值分析数值分析定义定义: 使求积公式使求积公式达到最高代数精度达到最高代数精度2n 1的求积公式称为的求积公式称为guass求积公式。求积公式。guass求积公式的节点求积公式的节点xk称为称为guass点点,系数系数ak称为称为guass系数系数.因为因为guass求积公式也是插值型求积公式求积公式也是插值型求积公式,故有故有结论结论: n 1个节点的个节点的插值型求积公式的代数精度插值型求积公式的代数精度 d
6、 满足满足： n d 2n 1。数值分析数值分析5数值分析数值分析例：例：选择系数与节点，使求积公式（选择系数与节点，使求积公式（1） 成为成为gauss公式。公式。解：解：n=1, 由定义，若求积公式具有由定义，若求积公式具有3次次代数精度，则代数精度，则 其是其是gauss公式。公式。 为此，分别取为此，分别取 f(x)=1, x，x2，x3 代入公式，并让代入公式，并让 其成为等式，得其成为等式，得c1 c2=2c1 x1 c2 x2=0c1 x12 c2 x22 =2/3c1 x13 c2 x23 =0求解得：求解得：所求所求gauss公式为：公式为：(1) 用待定系数法构造
7、高斯求积公用待定系数法构造高斯求积公式式数值分析数值分析6数值分析数值分析 设设pn(x)，n=0,1,2,为正交多项式序列，为正交多项式序列， pn(x)具有如下性质：具有如下性质：1）对每一个）对每一个n ,pn(x)是是 n 次多项式。次多项式。 n=0,1,2）(正交性正交性)3）对任意一个次数）对任意一个次数n-1的多项式的多项式p(x)，有，有4）pn(x)在在(a,b)内有内有n个互异零点。个互异零点。（2）利用正交多项式构造高斯求积公式）利用正交多项式构造高斯求积公式数值分析数值分析7数值分析数值分析定理定理2 设设x0,x1, ,xn 是是n 1次次正交多项式正交多项式pn
8、1(x)的的n 1 个零点个零点,则插值型求积公式则插值型求积公式是是guass型型求积公式。求积公式。证明：证明：只要证明只要证明求积公式的代数精确度为求积公式的代数精确度为2n 1,即即对对任意一个次数任意一个次数2n 1的多项式的多项式求积公式求积公式都精确成立。都精确成立。设设 f(x)为任意一个次数为任意一个次数2n 1的多项式，则有的多项式，则有 f(x)=q(x)pn 1(x) r(x)，满足，满足 f(xk)=r(xk)这里，这里， pn 1(x)是是 n 1次次正交多项式，正交多项式， q(x)、r(x)均是均是次数次数n的多项式。的多项式。数值分析数值分析8数值分析数值分析
9、由性质由性质3）及）及(4)式，有式，有由于由于n 1个节点的插值型求积公式的代数精确度不低个节点的插值型求积公式的代数精确度不低于于n，故有，故有即即对对 f(x)为任意一个次数为任意一个次数2n 1的多项式的多项式求积公式求积公式都都精确成立精确成立。 证毕证毕数值分析数值分析9数值分析数值分析利用正交多项式构造高斯求积公式利用正交多项式构造高斯求积公式的基本步骤：的基本步骤：代入积分式代入积分式因此，求积系数为因此，求积系数为数值分析数值分析10数值分析数值分析数值分析数值分析11数值分析数值分析数值分析数值分析12数值分析数值分析常用的高斯求积公式常用的高斯求积公式1.gauss -
10、legendre 求积公式求积公式 (1)其中高斯点为其中高斯点为legendre多项式的零点多项式的零点 guass点点xk, guass系数系数ak都有表可以查询都有表可以查询.数值分析数值分析13数值分析数值分析数值分析数值分析7 714数值分析数值分析数值分析数值分析15数值分析数值分析数值分析数值分析16数值分析数值分析一般区间的一般区间的gauss - gauss - legendre legendre 求积公式求积公式 如果积分区间是如果积分区间是a,b，用线性变换，用线性变换这样就可以用这样就可以用gauss - gauss - legendrelegendre求积公式计算一求
11、积公式计算一般区间的积分般区间的积分.将积分区间从将积分区间从a,b变成变成-1,1,由定积分的换元积由定积分的换元积分法有分法有数值分析数值分析17数值分析数值分析数值分析数值分析18数值分析数值分析数值分析数值分析19数值分析数值分析数值分析数值分析20数值分析数值分析数值分析数值分析21数值分析数值分析例例 利用高斯求积公式计算利用高斯求积公式计算解解: 令令x=1/2 (1 t), 则则用用高斯高斯-legendre求积公式计算求积公式计算.取取n=4 积分精确值为积分精确值为i=ln2=0.69314718由此可见，高斯公式精确度是很高的由此可见，高斯公式精确度是很高的.数值分析数值
12、分析22数值分析数值分析例例:分别用不同方法计算如下积分分别用不同方法计算如下积分,并做比较并做比较各种做法比较如下：各种做法比较如下：1、用、用newton-cotes公式公式当当n=1时，即用梯形公式，时，即用梯形公式，i0.9270354当当n=2时时, 即用即用simpson公式公式, i 0.9461359当当n=3时时, i 0.9461090当当n=4时时, i 0.9460830当当n=5时时, i 0.9460830i准准=0.9460831=0.9460831数值分析数值分析23数值分析数值分析2:用复化梯形公式用复化梯形公式 令令h=1/8=0.1253：用复化辛卜生公式
13、：用复化辛卜生公式 令令h=1/8=0.125i准准=0.9460831数值分析数值分析24数值分析数值分析4、用用romberg公式公式k tn sn cn rn0 0.9207355 1 0.9397933 0.94614592 0.9445135 0.9460869 0.94008303 0.9456906 0.9460833 0.9460831 0.9460831 i准准=0.9460831数值分析数值分析25数值分析数值分析5、用、用gauss公式公式解：解：令令x=(t 1)/2, i准准=0.9460831（2）用）用3个节点的个节点的gauss公式公式（1）用）用2个节点的个节
14、点的gauss公式公式数值分析数值分析26数值分析数值分析算法比较算法比较n此例题的精确值为此例题的精确值为0.9460831.n由例题的各种算法可知：由例题的各种算法可知：n对对newton-cotes公式，当公式，当n=1时只有时只有1位有效位有效数字，当数字，当n=2时有时有3位有效数字，当位有效数字，当n=5时有时有7位有效数字。位有效数字。n对复化梯形公式有对复化梯形公式有2位有效数字，对复化辛卜位有效数字，对复化辛卜生公式有生公式有6位有效数字。位有效数字。n用复合梯形公式，对积分区间用复合梯形公式，对积分区间0，1二分了二分了11次用次用2049个函数值，才可得到个函数值，才可得
15、到7位准确数字。位准确数字。n用用romberg公式对区间二分公式对区间二分3次，用了次，用了9个函个函数值，得到同样的结果。数值，得到同样的结果。n用用gauss公式仅用了公式仅用了3个函数值，就得到结果。个函数值，就得到结果。数值分析数值分析27数值分析数值分析2.gauss-chebyshev2.gauss-chebyshev公式公式常用的高斯求积公式常用的高斯求积公式数值分析数值分析28数值分析数值分析3.gauss-laguerre3.gauss-laguerre公式公式数值分析数值分析7 729数值分析数值分析4.gauss-hermite4.gauss-hermite公式公式数值
16、分析数值分析7 730数值分析数值分析二、高斯型求积公式的截断误差和稳定性分析二、高斯型求积公式的截断误差和稳定性分析数值分析数值分析31数值分析数值分析已知已知hermitehermite插值误差是插值误差是因为对因为对2n 12n 1次多项式求积公式准确成立，即次多项式求积公式准确成立，即代入上式代入上式即有即有数值分析数值分析32数值分析数值分析以下将证明高斯形求积公式的求积系数恒正以下将证明高斯形求积公式的求积系数恒正数值分析数值分析33数值分析数值分析数值分析数值分析由此说明高斯求积公式的计算误差可以控制，由此说明高斯求积公式的计算误差可以控制，高斯型求积公式是数值稳定的方法。高斯型
17、求积公式是数值稳定的方法。34数值分析数值分析 将积分区间将积分区间a , b n等分，在每个小子区间上等分，在每个小子区间上使用一个节点数较少的使用一个节点数较少的gaussgauss型求积公式型求积公式, ,然后把它然后把它们加起来，就得到整个区间上们加起来，就得到整个区间上gaussgauss型求积公式的复型求积公式的复化形式。化形式。 复化复化gauss求积公式的基本思想：求积公式的基本思想： 下面用下面用gauss-legender求积公式推导求积公式推导复化复化gauss型求积公式型求积公式.将积分区间将积分区间a , b n等分，等分，三、三、复化复化gauss求积公式求积公式数值分析数值分析35数值分析数值分析数值分析数值分析36数值分析数值分析 例如例如,用用2点的点的gauss-legender求积公式复合求积公式复合,由表由表9-4,取取n=1,得得aj =1,xj=0.5773502692代入到上式代入到上式中中,得得2点的复化点的复化gauss-legender求积公式求积公式 再将上式应用再将上式应用gauss-legender求积公式就得求积公式就得到了到了复化复化gauss型求积公式型求积公式.数值分析数值分析37数值分析数值分析数值分析数值分析二版习题 p276-11(1)，16，三版习题三版习题 p250-11(1)，16，38