数值积分_数值积分原理__matlab实现.docVIP

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