多种数值积分的比较(Gauss抛物线龙贝格)要点.docxVIP

多种数值积分求积公式的分析比较 吴春晖 （中国海洋大学海洋环境学院 山东 青岛 266100） 摘要： 对于运用牛顿-莱布尼茨积分公式不能较好解决的定义在区间[a,b]上的可积函数，原函数并不能简单地用初等函数来表达，故需要构造定积分的近似计算公式。在本文中，主要构建了抛物线求积公式及其复化抛物线公式。在对抛物线类的求积公式进行应用检验后，再运用Gauss求积公式，构建Gauss-Laguerre求积公式，对相同的问题进行运用，并比较截断误差。之后再对求积过程进行优化，在限定误差范围的情况下，利用龙贝格算法，对求积加速收敛。 关键词： 抛物线求积 复化求积 Gauss-Laguerre 加速收敛 引言： 对于一些较为复杂的函数，在一定的误差要求下，需要通过构造的方式求给定函数的定积分。基本的替代法主要有梯形面积及抛物线近似代替曲边梯形。并通过划分更小的区间，减少截断误差从而提出了复化梯形及抛物线公式。为了提高运算效率，有加速收敛的Richardson外推法和Romberg求积公式。之后，针对节点数固定情况下，提出了Gauss公式，使其获得最大的精度。本文主要研究的是抛物线求积法与Gauss-Laguerre公式。 目录  TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc438471405 第一章 抛物线求积公式及应用  PAGEREF _Toc438471405 \h 3  HYPERLINK \l _Toc438471406 1.1 抛物线求积公式的算法  PAGEREF _Toc438471406 \h 3  HYPERLINK \l _Toc438471407 1.2抛物线求积公式的matlab程序  PAGEREF _Toc438471407 \h 3  HYPERLINK \l _Toc438471408 1.3复化抛物线求积公式的应用  PAGEREF _Toc438471408 \h 4  HYPERLINK \l _Toc438471409 第二章Gauss-Laguerre求积公式及应用  PAGEREF _Toc438471409 \h 5  HYPERLINK \l _Toc438471410 2.1 Gauss-Laguerre的算法  PAGEREF _Toc438471410 \h 5  HYPERLINK \l _Toc438471411 2.2Gauss-Laguerre公式的matlab程序  PAGEREF _Toc438471411 \h 5  HYPERLINK \l _Toc438471412 2.3 Gauss-Laguerre求积公式的应用  PAGEREF _Toc438471412 \h 6  HYPERLINK \l _Toc438471413 第三章 龙贝格算法与算法优化  PAGEREF _Toc438471413 \h 7  HYPERLINK \l _Toc438471414 3.1龙贝格算法及程序  PAGEREF _Toc438471414 \h 7  HYPERLINK \l _Toc438471415 3.2利用龙贝格算法优化求积  PAGEREF _Toc438471415 \h 7  HYPERLINK \l _Toc438471416 3.3 龙贝格算法的应用  PAGEREF _Toc438471416 \h 8  HYPERLINK \l _Toc438471417 第四章 数值积分的分析总结  PAGEREF _Toc438471417 \h 9  抛物线求积公式及应用 抛物线求积公式的算法 抛物线求积公式，是将区间二等分，以中点及两端点作为抛物线的三个点，并求出抛物线，在区间上对抛物线函数求积分。而复化抛物线公式，即在给定区间上，分成多个小份，并分别用抛物线公式求积。 所求得的数值积分余项如下： 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时，则为简化