数值积分(论文)

《数值积分(论文)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、目录第一章数值积分计算的重述11.1弓丨言11.2问题重述2第二章复化梯形公式22」复化梯形公式的算法描述32.2复化梯形公式在C语言中的实现32.3测试结果4笫三章复化simpson公式53」复化simpson公式的算法描述53.2复化simpson公式在C语言中的实现63.3测试结果6第四章复化cotes公式84」复彳匕cotes公式的算V去描述84.2复化cotes公式在C语言屮的实现84.3测试结果9第五章Romberg积分法105」Romberg积分法的算法描述105.2Romberg积分法在C中的实现115.3测试结果12第六章结果对比分析和体会134参考

2、文献14附录14数值积分(—)第一章数值积分计算的重述1.1引言数值积分是积分计算的重要方法，是数值逼近的重要内容，是函数插值的最直接应用，也是工程技术计算中常常遇到的一个问题。在应用上，人们常要求算出具体数值，因此数值积分就成了数值分析的一个重要内容。在更为复杂的计算问题中，数值积分也常常是一个基本组成部分。在微积分理论屮，我们知道了牛顿■莱布尼茨(Newton-Leibn⑵公式^f(x)dx=F(b)-F(a)其中F")是被积函数/(X)的某个原函数。但是随着学习的深入，我们发现一个问题：对很多实际问题，上述公式却无能为力。这主耍是因为:它们或是被积函数没有解析形

3、式的原函数，或是只知道被积函数在一些点上的值，而不知道函数的形式，对此，牛顿一莱布尼茨(Newton-Leibn⑵公式就无能为力了。此外,即使被积函数存在原函数,但因找原函数很复杂，人们也不愿花费太多的吋间在求原函数上，这些都促使人们寻找定积分近似计算方法的研究,特别是有了计算机后，人们希望这种定积分近似计算方法能在计算机上实现，并保证计算结果的精度,具冇这种特性的定积分近似计算方法称为数值积分。由定积分知识，定积分只与被积函数和积分区间冇关，而在对被积函数做插值逼近时，多项式的次数越高，对被积函数的光滑程度要求也越高，月.会出现Runge现彖。如〃〉7时,Newto

4、n-Cotes公式就是不稳定的。因而，人们把目标转向积分区间，类似分段插值,把积分区间分割成若干小区间，在每个小区间上使用次数较低的Newton-Cotes公式，然后把每个小区间上的结果加起来作为函数在整个区间上积分的近似，这就是复化的基木思想。木文主要研究的公式有：复化梯形公式、复化Simpson公式、复化Cotes公式、Romberg积分法。1.2问题重述本文主耍介绍微积分方程的复化解法。通过运用复化梯形公式、复化Simpose公式、复化cotes公式和Romberg积分法这四种积分法方法，解出微分方程的近似解。并进行误差分析和结果比较。当积分区间［a,b］的长度

5、较大，而节点个数n+1固定时，直接使用Newton-Cotes公式的余项将会较大，而如呆增加节点个数即n+1增加时，公式的舍入误差又很难得到控制，为提高公式的精度，又使算法简单易行，往往使用复化方法。即将积分区间［a,b］分成若干个子区间，然后在每个小区间上使用低阶Newton-Cotes公式，最后将每个小区间上的积分的近似值相加。将定积分的积分区间［ab］分割为n等份h—Z7n・l）上各节点为忑=a+kh,k=0丄…nh=^—^在子区间［耳，忑+J(k=O丄1n使用NewtonCotes公式将子区间分割为1等份，节点为Xk,xkh2hIh+7,心+了=林+】在子区间

6、上作f(x)的1^Newton-Cotes求积公式打⑴d⑴胡—(s-b心0牛#0巧由积分的区间可加性，可得an-1畑J7(x)d(x)二工f(x)dxb"()I"一1"一1I复化求积公式u⑴.£(“+•)=inR=()Jt=()/=()第二章复化梯形公式2.1复化梯形公式的算法描述复化求积公式£/$)“空士cPfgj=Ink=0BO/=O当L=1时可得复化梯形公式：

7、7⑴d(QS=hYYC^f(xk+i)/k=()i=O/i-11"工刃/g)+/(畑)]k=O厶h—n"一1复化梯形公式二〒lm)+2£/(九)+,(b)]2/7“02.2复化梯形公式在C语言中的实现复

8、化梯形公式运用的程序如下：T0=(a-b)*(f_x(a)+f_x(b))/2;//n=1时的cotes公式即梯形公式for(i=1;i<=100;i++){〃计算sum_num、xishu>s_point(startpoint)>d_pointsum_num=pow(23-1)；xishu=double(a-b)/sum_num;s_point=double(b)+double(a-b)/pow(24);d_point=clouble(a-b)/pow(2,i-1);for(j=l;j<=sum_num;j++){add_T=add_T+f_x(s