数值分析实验二 数值积分3.docxVIP

PAGE

PAGE1

数值分析实验二 数值积分

组号 班级 学号 姓名 分数

一：实验目的

学会用复化梯形公式，复化辛普森公式，Romberg求积公式，Guass求积公式，高斯—勒让德求积公式，Newton—Cots求积公式等利用函数在若干点的函数值，近似的计算该函数在某个区间上的积分值

二：实验内容及基本知识介绍

数值积分的基本思想

实际问题当中常常需要积分，有些数值方法，如微分方程和积分方程的求解，也都和积分计

算相联系。依据人们所熟知的微积分基本定理，对于积分I=?b

a

f(x)dx,理论上可以用牛顿—

莱布尼茨公式?b

a

f(x)dx=F(b)-F(a),其中F(x)是被积函数f(x)某个原函数，但对很多实际问

题，上述公式却无能为力。

这是因为：

被积函数f(x)原函数理论上存在，但无法知道它的可用于计算的表达式，如ex2，sinx

x

等初等函数。

被积函数f(x)的本身没有可用于计算的表达式，而仅仅是一种数表函数，即只知道该函数的部分特点的函数值。

因此，借助插值理论是解决数值计算定积分的有效途径之一。

,复化梯形公式

将区间［a,b］划分为几等份，分点x

=a+kh, h=

k

b?an

， k=0.1.2……n，

在每个子区间［x,x ］ (k=0.1.2…..n－1)

k k+1

上采用梯形公式：

bb?a

b

2? f(x)dx≈ ［f(a)+f(b)］，

2

a

则得

I=?b

f(x)dx= ?k?1

f(x)d(x)

?n?1 xa k

?n?1 x

=h?n?1

［f(x)+f(x )］+R(f)

2

k?0

k k+1 n

记 T=h?n?1

［f(x

)+f(x

)］=h［f(a)+2?n?1

f(x

)?f(b)］

n 2 k k+1 2 k

k?0 k?1

称为复化梯形公式。

b?a

其余项 R(f)=－

n 12

h2f＂(η) ?η?（a,b）

，复化辛普森求积公式

将区间［a,b］分为几等份，在每个子区间［x,x ］上采用辛普森公式即

k k+1

b?a a?b

S= 6

1

［f(a)+4f( 2

b

)+f(b)］

?n?1 x

若记x

=x+ h，则得 I=? f(x)dx=

?k?1

f(x)d(x)

1k? k 2

1

2

a x

?k0 k

?

=h?n?1

［f(x

)+4f(x

)+f(x

)］+R(f),

16 k k?

1

k+1 n

1k?0 2

1

记 S=h?n?1

［f(x

)+4f(x

)+f(x

)］=

n 6 k

k+1

k?

k?0 2

h［f(a)+4?n?1

6

k?0

f(x

1k?

1

2

)?2?n?1

k?1

f(x

k

)?f(b)］

( )

( )

其余项R

h

(f)=I－S=－

h4?n?1

f(? ) η

?(x x )

n n

，Romberg求积公式

180 2

k

k?0

k k,

k?1

设以T

0

(k)表示二分k次后求得的梯形值，且以T

m

(k)表示序列｛T

0

(k)｝的m次加速值，则

4m h 1

依递推公式： T

(h)=

T ( )?

T (h)

m 4m?1

m?1

2 4m?1

m?1

4m 1

可得 T

m

(k)=

4m

4

T

?1 m?1

(k?1)?

4m

T

?1 m?1

(k), k=1.2

称为Romberg求积算法。

三：实验问题及方法、步骤

例1，利用复化梯形公式或Simpson公式计算积分I=?1sinxdx时，若欲使结果的误差为10-4，

0 x

则应将［0,1］几等分并计算结果。算法与程序：定步长梯形法function t=natrapz(fname,a,b,n)

%用途：定步长梯形求函数的积分

%格式：t=natrapz(fname,a,b,n)其中fname是被积函数a,b分别为上下限，n为等分数h=(b-a)/n;

fa=feval(fname,a);fb=feval(fname,b);f=feval(fname,a+h：h：b-h+0.001*h);t=h*(0.5*(fa+fb)+sum(f));

用于解例题1，在MATLAB命令窗口执行：

formatlong;natrapz(i