多种数值积分的分析比较(Gauss抛物线龙贝格)..docx

多种数值积分求积公式的分析比较吴春晖（中国海洋大学海洋环境学院 山东 青岛 266100）摘要：对于运用牛顿-莱布尼茨积分公式不能较好解决的定义在区间[a,b]上的可积函数，原函数并不能简单地用初等函数来表达，故需要构造定积分的近似计算公式。在本文中，主要构建了抛物线求积公式及其复化抛物线公式。在对抛物线类的求积公式进行应用检验后，再运用Gauss求积公式，构建Gauss-Laguerre求积公式，对相同的问题进行运用，并比较截断误差。之后再对求积过程进行优化，在限定误差范围的情况下，利用龙贝格算法，对求积加速收敛。关键词： 抛物线求积 复化求积Gauss-Laguerre加速收敛引言：对于一些较为复杂的函数，在一定的误差要求下，需要通过构造的方式求给定函数的定积分。基本的替代法主要有梯形面积及抛物线近似代替曲边梯形。并通过划分更小的区间，减少截断误差从而提出了复化梯形及抛物线公式。为了提高运算效率，有加速收敛的Richardson外推法和Romberg求积公式。之后，针对节点数固定情况下，提出了Gauss公式，使其获得最大的精度。本文主要研究的是抛物线求积法与Gauss-Laguerre公式。抛物线求积公式及应用抛物线求积公式的算法抛物线求积公式，是将区间二等分，以中点及两端点作为抛物线的三个点，并求出抛物线，在区间上对抛物线函数求积分。而复化抛物线公式，即在给定区间上，分成多个小份，并分别用抛物线公式求积。所求得的数值积分余项如下：1.2抛物线求积公式的matlab程序具体的程序代码如下：function?s=simpr1(f,a,b,M)?h=(b-a)/(2*M);?s1=0;?s2=0;?for?k=1:M?????x=a+h*(2*k-1);?????s1=s1+feval(f,x);?end?for?k=1:(M-1)?????x=a+h*2*k;?????s2=s2+feval(f,x);?end?s=h*(feval(f,a)+feval(f,b)+4*s1+2*s2)/3;?说明：f代表的是原函数，该程序为抛物线复化积分公式。当M值为1时，则为简化的抛物线求积公式。最后s即为用matlab语言所表述的抛物线求积公式。1.3复化抛物线求积公式的应用在应用部分，所计算的积分为1/(1+x^2)。将函数与节点，积分区间输入函数后，得出的结果如下图表所示。图1.3.1 复化抛物线求积（所取节点为5-50）Fig.1.3.1Compound parabolic quadrature (from the node for 5-50)图1.3.2 积分误差（节点数为0-50）Fig.1.3.2quadrature error (from the node for 0-50)节点数05101520253035404550误差0.0020659.91E-091.55E-101.36E-112.42E-126.35E-132.13E-138.44E-143.77E-141.87E-141.01E-14表1.3.1 每隔5节点的误差Table 1.3.1 every 5 node error分析总结：复化抛物线积分对给定的函数有较好的适用性，在节点数为3时，误差就处于E-3数量级。这主要与函数的特性及选择区间有关，所求得积分为1/4圆，故抛物线较能吻合特征。从图像我们可以观察到：随着节点数的增加，积分的准确性得到增强，且增加速度较快。在50个节点的情况下，误差达到了E-14的数量级。从误差观察得出，误差减少到下一数量级的节点间隔数增大。第二章Gauss-Laguerre求积公式及应用2.1 Gauss-Laguerre的算法Gauss-Laguerre求积公式是Gauss求积公式的一种建立在无穷区间上的特殊求积公式。Laguerre多项式：在关于权函数为的正交多项式（1）。故在求积分时，我们主要使用的求积公式为：（2）其中,（k=1,…,n）是的n个零点，求积系数k=1,2,…,n （3）2.2Gauss-Laguerre公式的matlab程序程序代码如下：syms xf;syms l;n=20;l(n)=exp(x)*diff(x^(n)*exp(-x),x,n);[xs]=solve([char(l(n)),'=0'],x);ll=diff(l(n),1);for i=1:n a(i)=(factorial(n))^2/(double(xs(i))*(subs(ll,x,double(xs(i)))^2));ends=double(a)*feval(f1,double(xs));程序分析：n代表Laguerre的项数，f为函数。该程序先求出相应的Laguerre多项式后，通过solve解出节点，再用节点算出系数，代入