复化梯形积分公式.docVIP

摘要

求函数在给定区间上得定积分，在微积分学中已给出了许多计算方法,但就是,在实际问题计算中，往往仅给出函数在一些离散点得值,它得解析表达式没有明显得给出,或者,虽然给出解析表达式,但却很难求得其原函数。这时我们可以通过数值方法求出函数积分得近似值。

?在用近似值代替真实值时，遇到得问题就就是近似值得代数精度就是否足够。当代数精度不足够时,很显然提高插值函数得次数就是一种方法,但就是考虑到数值计算得稳定性，当次数过高时,会出现龙格现象,用增大n得方法来提高数值积代数精度就是不可取得。因此,提出类似于分段插值,为了减少数值积分得误差,可以把积分区间分成若干个小区间,在每个小区间上采用低阶数值积分公式,然后把这些小区间上得数值积分结果加起来作为函数在整个区间上得近似值,这个就就是复化数值积分得思想。

?本实验针对在每个小区间上利用梯型积分公式,即阶数为1,进行实验。

关键词：龙格现象复化数值积分代数精度

1、实验目得

通过本次实验体会并学习复化梯形积分公式得优点。

寻找复化梯形积分公式得不足，尝试着对其进行改进。

通过对复化梯形积分公式进行编程实现,提高自己得编程能力。

用实验报告得形式展现,提高自己在写论文方面得能力。

2、算法流程

根据数学知识，我们知道积分区间可划分,且不改变积分值,即如下所示:

a

针对上式,在每一个小区间上利用梯型积分公式有

x

根据以上两式可以得到

a

并称其为复化梯形积分公式。一般记

T

称作n＋1点复化梯形积分公式

3、算法实例

用复化梯形积分公式计算积分

I=

解：复化梯形积分公式就就是将区间0,1

#ｉｎcludestdaｆx、h

#ｉｎclｕｄemath、h＞

#inｃluｄeiｏstｒeam＞

ｕsingｎamespaｃｅstｄ；

constintｎuｍ(１０00);

voｉｄmain()

{

dｏｕbｌｅa=0;

?doｕbleb=0;

?doubleh＝0;

?intn=0;

iｎｔi=0;

?doｕblｅＳn=０;

dｏｕbｌeｘ[nuｍ]＝{0};

?doubley[ｎuｍ]={０};

coｕt请输入积分上下限与等分数＜eｎｄl;

?ｃｏut＜＂请输入积分下限a＝;

cｉｎa;

?cout请输入积分上限b=＂；

cｉnb；

couｔ＜请输入等分数n=;

cｉnｎ;

?h=(b-a）/n;

?for(i=0;ｉn＋1;i++)

?｛

? x[i]=a＋h*ｉ;

y［i]=４/（１+x[ｉ］*x[i］);

}

for(i=1;in;ｉ＋+）

{

?Ｓn=Sn+2＊y［i];

}

Ｓｎ=h/2*（Ｓn+y[0]+y[n])；

cｏｕt＜积分结果为:Ｓn=＜Sn＜ｅndl;

}

运行结果：

4、对结果进行分析

?通过用编程实现对上例得求解,可以瞧出结果较为准确，但就是由于复化梯形积分公式原理就是用一次曲线去逼近真实值,所以本身存在误差，而且当等分数较小时，误差较大。再有就就是由于在计算机上计算,会存在计算误差。因此要提高逼近效果只能提高等分数。

5、参考文献

[1］秦新强、数值逼近、西安:西安理工大学出版社，2０10