高斯积分学习.ppt

高斯积分;在进行有限元分析时，一般需要大量的数值积分，通常都是通过坐标变换把被积函数全部化为局部坐标的函数。经过这样的变换，计算将会变得简单。 若被积函数不是多项式函数，比如采用等参单元就是这样，此时就需要求助于数值积分，在有限元分析中，经常采用的就是高斯数值积分。;对于积分区间[a,b]，通过变换 可化为区间[-1,1] 所求积分则化为： 不失一般性，取a=-1，b=1，讨论区间[-1,1]上的积分即可，则有 其中 为权函数 ; 求积公式 的代数精度最高不超2n-1次。 证明：分别取 f(x)=1, x，x2，...xn 时代入公式，并让其成为等式得 A1 + A2 + …… + An =∫ab1dx.= b-a x1 A1 + x2 A2+ …… +xn An =∫abxdx.= （b2-a 2)/2 ...... x1 rA1 + x2 rA2+ …… +xn rAn =∫abxr dxr =(br+1-a r+1)/ (r+1) 上式共有 r 个 等式，2n个待定系数(变元),要想如上方程组有唯一解，应有方程组中方程的个数等于变元的个数,即 r=2n,这样求出的解答应的求积公式的代数精度至少是2n-1,下面证明代数精度只能是2n-1. [ 如果事先已选定[a ，b]中求积节点xk如下a?x1 ?…x n?b,上式成为n个未知数 A1、...An的n元线性方程组，此时要r=n 时方程组有唯一解] ;;;　Gauss - Legendre 求积公式 　 　　　　　　　　　　　　　　 　其中高斯点为Legendre多项式的零点 Ln(x)= 对于一般有限区间[a,b]，用线性变换x=(a+b)/2+(b-a)t/2使它变成为[-1,1]。;;例题：利用高斯求积公式计算;;