文档介绍：课程设计报告课程设计题目:求解的近似值课程名称:数值分析课程设计指导教师:XXX小组成员:XXXXXXXXX2013年12月31日目录目录 1题目 2一、摘要 2二、设计目的 2三、理论基础 31、复合矩形法求定积分的原理 32、复合梯形法求定积分的原理 33、复合辛普森法求定积分的原理 44、龙贝格求积公式原理 5四、程序代码及运算结果 51、复合矩形法求定积分:用函数 52、复合梯形法求定积分 6方法一 6方法二:用函数 73、复合辛普森法求定积分 7方法一 7方法二:用函数 74、龙贝格求定积分 85、Lobatto数值积分法 96、波尔文(Borwein)高阶公式 9五、结果分析 10六、设计心得 10七、参考文献 11题目:(1)已知:,所以。于是,我们可以通过计算上述定积分的近似值来得到π的近似值。(2)波尔文(Borwein)高阶公式在π值的高阶算法研究中,最好的结果来自两个都叫波尔文的数学家。他们在1984年发表了一个2阶收敛公式:,,,式中。试运用上述迭代算法,计算圆周率的近似值,并和前面传统方法进行比较。一、摘要借助matlab环境下的计算机编程语言,先用基本的积分函数对给出的题目进行求积分,然后基于给出的波尔文高阶收敛公式,在进行了连续迭代后,对运行结果做出分析,同时与之前的传统做法进行比较。二、设计目的用熟悉的计算机语言编程,上机完成用复合矩形法、复合梯形法、复合辛普森法、龙贝格法以及Lobatto数值积分方法,掌握各种方法的理论依据及求解思路,了解数值积分各种方法的异同与优缺点。三、理论基础1、复合矩形法求定积分的原理定积分的几何意义是计算曲边梯形的面积,若将积分区间[a,b]分为n段,每段都是一个小曲边梯形,用一个个小矩形代替这些小曲边梯形,然后把所有小矩形的面积加起来就近似得等于整个曲边梯形的面积,于是便求出了定积分的近似值,这就是复合矩形法求定积分的基本原理。设第i个小矩形的宽度,高度,根据可知其面积为,可得出复合矩形求积公式:如果积分区间被等分为段:,则。2、复合梯形法求定积分的原理把积分区间[a,b]分成若干小区间,在每个小区间上以梯形面积近似曲边梯形面积,即用梯形公式求小区间上积分的近似值定积分存在定理表明,只要被积函数连续,当小区间长度趋于零时,小区间面积之和趋于曲边梯形面积的准确值,即定积分的准确值。设第i个小梯形的宽度,两底高度分别为和,则定积分的近似值为设被积函数在上连续可导,把区间等分,令,于是有由于式中,,,代入上式得出复合梯形求积公式:表示区间分为等分时,用复合梯形求积法求出的定积分值。3、复合辛普森法求定积分的原理辛普森公式:因为辛普森公式用到了区间的中点,所以在构造复合辛普森公式时,把积分区间等分为偶数份。令,其中为正整数,节点为,在每两个小区间上用辛普森公式,则有记上式叫做复合辛普森公式,的下标表示将积分区间等分。4、龙贝格求积公式原理采用龙贝格求积公式,即逐次对分积分区间的方法,可以把前面计算的结果作为一个整体带入对分后的计算公式中,只需增加新的分点的函数值。龙贝格秋季攻势是一个很实用的公式。若已知与的关系:,其中,若记为新增分点函数值的和,则ξ。对而言,为下标为奇数位置的函数值的和。龙贝格求积公式:,用作为近似值,截断误差为。四、程序代码及运算结果1、复合矩形法求定积分:用函数编辑如下命令做出函数,并保存:functiony=qiupai(x);y=4./(1+x.^2);该函数的图像如右图所示,生成方法如下:x=0::1;plot(x,qiupai1(x),'linewidth',2),gridlegend('qiupai1(x)')在命令窗口输入并运行:h=1/100;x=0+h:h:1;s=h*sum(qiupai1(x));pi1=vpa(s,6);formatlongpi1运行结果:pi1=,值的精度相应提高,下表即为改变h的值得到的相应的pil的值:h1/10,0001/1,000,000/100,000,、复合梯形法求定积分方法一先编辑梯形函数并保存:functions=tixing(f,a,b,n)%f是被积函数;%a,b分别为积分的上下限;%n是子区间的个数;%s是梯形总面积;h=(b-a)/n;s=0;fork=1:(n-1)x=a+h*k;s=s+feval('f',x);endformatlongs=h*(feval('f',a)+feval('f',b))/2+h*s;然后在编辑窗口输入如下命令并运行:pi1=tixing('4/(1+x^2)',0,1,100)运行结果:pi1=,值的精度相应提