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

课程设计报告课程设计题目：求解 的近似值

课程名称：数值分析课程设计

指导教师： 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、复合矩形法求定积分：用sum函数 (5)

2、复合梯形法求定积分 (6)

方法一 (6)

方法二：用trapz函数 (7)

3、复合辛普森法求定积分 (7)

方法一 (7)

方法二：用quad函数 (7)

4、龙贝格求定积分 (8)

5、Lobatto数值积分法 (9)

6、波尔文(Borwein)高阶公式 (9)

五、结果分析 (10)

六、设计心得 (10)

七、参考文献 (11)

题 目：

（1）已知：411

02π=+⎰

x dx ，所以 ⎰+=10214

dx x π 。于是，我们可以通过计算上述定积分的近似值来得到π的近似值。 （2）波尔文(Borwein )高阶公式

在π值的高阶算法研究中，最好的结果来自两个都叫波尔文的数学家。他们在1984年发表了一个2阶收敛公式：

20=a ，00=b ，220+=p ，

⎪⎪⎪

⎩⎪

⎪

⎪⎨⎧++=++=+=++++++1

111

11

1)1()1(2)

1(k k k k k k

k k k k k k k b a b p p b a b a b a a a

式中π→k p 。试运用上述迭代算法，计算圆周率的近似值，并和前面传统方法进行比较。

一、摘要

借助matlab 环境下的计算机编程语言，先用基本的积分函数对给出的题目进行求积分，然后基于给出的波尔文高阶收敛公式，在进行了连续迭代后，对运行结果做出分析，同时与之前的传统做法进行比较。

二、设计目的

用熟悉的计算机语言编程，上机完成用复合矩形法、复合梯形法、复合辛普森法、龙贝格法以及Lobatto 数值积分方法，掌握各种方法的理论依据及求解思路，了解数值积分各种方法的异同与优缺点。

三、理论基础

1、复合矩形法求定积分的原理

定积分的几何意义是计算曲边梯形的面积，若将积分区间[a,b]分为n 段，每段都是一个小曲边梯形，用一个个小矩形代替这些小曲边梯形，然后把所有小矩形的面积加起来就近似得等于整个曲边梯形的面积，于是便求出了定积分的近似值，这就是复合矩形法求定积分的基本原理。

设第i 个小矩形的宽度),,2,1(1n i x x h i i i =-=-，高度)(i i x f y =，根据

)()()()()()(00)

0(00x f a b x f c a b dx x P dx x f b

a

b

a

-=-==⎰⎰

可知其面积为i i i i i y h x f x x =--)()(1，可得出复合矩形求积公式：

i n

i i n n b

a

h y h y h y h y f I dx x f ∑⎰

==+++≈=1

2211)()(

如果积分区间被等分为段：h n a b x x b x x x a i i n =-=-=<<<=-/)()(,110 ，则

i n

i i n

i i b

a

h y y n a b f I dx x f ∑∑⎰

===-≈=1

1)()(。 2、复合梯形法求定积分的原理

把积分区间[a,b]分成若干小区间，在每个小区间上以梯形面积近似曲边梯形面积，即用梯形公式求小区间上积分的近似值定积分存在定理表明，只要被积函数连续，当小区间长度趋于零时，小区间面积之和趋于曲边梯形面积的准确值，即定积分的准确值。

设第i 个小梯形的宽度),,2,1(1n i x x h i i i =-=-，两底高度分别为

)(11--=i i x f y 和),,2,1)((n i x f y i i ==，则定积分的近似值为

=+≈=∑∑⎰

⎰

=-=-i n

i i i n

i x x b

a

h y y x f dx x f i i

)(2

1

)()(111

1

])()()()[(2

1

1112221110n n n n n n h y y h y y h y y h y y ++++++++---- 设被积函数)(x f 在],[b a 上连续可导，把],[b a 区间n 等分，令n a b h /)(-=，于是

有

⎪⎪⎭

⎫ ⎝⎛

++=+++++≈∑⎰

-=-110121022)222(2)(n j j n n n b

a

y y y h h y y y y y h dx x f 由于式中)(0a f y =，)(b f y n =，)()(jh a f x f y j j +==，代入上式得出复合梯形求积公式：

n n j b

a

T jh a f b f a f h dx x f =⎪⎪⎭

⎫ ⎝⎛

+++=∑⎰

-=11)(2)()(2)( n T 表示区间分为n 等分时，用复合梯形求积法求出的定积分值。

3、复合辛普森法求定积分的原理 辛普森公式：)]()2

(4)([6)(b f b

a f a f a

b f I +++-=

因为辛普森公式用到了区间的中点，所以在构造复合辛普森公式时，把积分区间],[b a 等分为偶数份。令m n 2=，其中m 为正整数，节点为

)12,,2,1,0(-=+=m k kh a x k ，m

a

b h 2-=

在每两个小区间],[122+k k x x 上用辛普森公式，则有

==

∑⎰

⎰

-=+n m k x

x b

a

k k

x f dx x f 10

21

22)()(

=+++∑∑-=+-=+1

221

0122][)]()(4([3m k k k m k k k f R x f x f x f h

∑∑∑-=-=-=++⎥⎦⎤⎢⎣⎡+++1

101

0122][)(4)(2)()(3m k k m k m k k k f R x f x f b f a f h

记

⎥⎦

⎤⎢⎣⎡+++=∑∑-=-=+101

01222)(4)(2)()(3m k m k k k m

x f x f b f a f h S 上式叫做复合辛普森公式，m S 2的下标m 2表示将积分区间],[b a m 2等分。