《合肥工业大学计算方法复化梯形公式实验》由会员分享，可在线阅读，更多相关《合肥工业大学计算方法复化梯形公式实验（5页珍藏版）》请在装配图网上搜索。
1、计算方法实验报告学号姓名班级实验项目名称实验二 数值积分一、实验名称实验二数值积分二、 实验目的:(1) 熟悉复化梯形方法、复化impon方法、梯形递推算法、龙贝格算法;(2) 能编程实现复化梯形方法、复化Simpon方法、梯形递推算法、龙贝格算法；(3) 理解并掌握自适应算法和收敛加速算法的基本思想；(4) 分析实验结果体会各种方法的精确度,建立计算机求解定积分问题的感性认识三、 实验内容及要求(1）设计复化梯形公式求积算法，编制并调试相应的函数子程序()设计复化辛浦生求积算法，编制并调试相应的函数子程序（）用龙贝格算法计算输入：积分区间,误差限输出：序列Tn,Sn,Cn,Rn及积分结果(参
2、考书本P81的表25) 取n=2,4，8，1,精确解为0946081四、 实验原理及算法描述在许多实际问题中，常常需要计算定积分的值。根据微积分学基本定理,若被积函数(x)在区间a，b上连续，只要能找到f（x）的一个原函数F()，便可利用牛顿莱布尼兹公式求得积分值。但是在实际使用中，往往遇到如下困难,而不能使用牛顿莱布尼兹公式。（1） 找不到用初等函数表示的原函数（2） 虽然找到了原函数，但因表达式过于复杂而不便计算（3） f(x）是由测量或计算得到的表格函数由于以上种种困难,有必要研究积分的数值计算问题.利用插值多项式 则积分转化为，显然易算。称为插值型求积公式。最简单的插值型求积公式是梯形
3、公式和Simpso公式,。当求积结点提供较多,可以分段使用少结点的梯形公式和ison公式，并称为复化梯形公式、复化isn公式。如步长未知，可以通过误差限的控制用区间逐次分半的策略自动选取步长的方法称自适应算法。梯形递推公式给出了区间分半前后的递推关系。由梯形递推公式求得梯形序列，相邻序列值作线性组合得Smpson序列， Simpson序列作线性组合得柯特斯序列， 柯特斯序列作线性组合的龙贝格序列.若R2R1|e，则输出R；否则依此类推。如此加工数据的过程叫龙贝格算法,如下图所示：复化梯形公式复化Smpson公式梯形递推公式加权平均公式: 龙贝格算法大大加快了误差收敛的速度，由梯形序列O（h2）
4、 提高到龙贝格序列的O(h8）五、 程序代码及实验结果1 主程序t mai()/cout longbeige(0， 1, 0。0000)；cout ”请输入 你的区间 和 误差限度: y z；cut ”根据龙贝格算法 求出的精确值为： lonbi（x, y,z) enl;out ”K ” ”T” ”S ”C” ” nd；co 0” ” tiing(0， 1, 1） endl；out 1” txing(0, 1, 2) ” xinbuhen(， 1, 1) enl;cou ”2” ” tixin(0， 1， 4) ” xgbuheng（, , ） (（doble（16/5） xigbusn（0,
5、1, 2)） (oule （1/1)）*xingbheg（0， 1,1) enl;ut 3 tixing（0，1,） xiushen（0,1,4）” (（double(1 /15)* xgbushe（0， ,） (double(1 / 1)）*inbuhen（，1， 2）” ”longege(0，1， 0。00000001) endl;out ”4” tixn(0， 1, 16） ” xngbusheng（0， 1,8) ” (（doble（16 15）) xingbusheng（0, 1， 4）) （oule(1 / 15)）xigbushn（0， 1，2） longbege(， 1， 0。
6、0000001) ndl;；ot ”h resul ”4” i longbeie(0， ，0.00001)；rrn 0；2 复化梯形公式子程序:doble iing（dolea, ublb，int）doublfa =f(a)；doul fb = f（b）；db (-） / n;dub fxk 0.0；f(ntk= ； k 1； k+）doublexk = a+ ;fk = f(k) + fxk;double res=（h / 。0)(fa + 。*fxk +fb)；return res；3 复化辛浦生公式子程序：ouble ingbusheg(doule a, doble b， it n）do
7、be h (b ） n；dbe fa (a）；doubl fb f（b）；/dobe s= b f;doue x= ;doube fx12 = 0；doul fxk = 0。;/or （itk =； k=；k+)x x + h （2.0);ube fx= f（x）；s = s 4。0*f（x)；x= x +h/(0)； = s + 2(x）； (h 60)s;retur s；*/for (ink = 0；k = uci）1= R2；1 = C2;k+;h= 2；T1 = T2；S1=S2；goto lop;elsrturn R;实验结果：图1 六、 实验总结1. 梯形公式的收敛速度太慢，所以我们才会选择之后的几种公式加快收敛速度。2. 通过实验可以看出，龙贝格算法大大加快了误差收敛的速度,由梯形序列(h2) 提高到龙贝格序列的O（8)3. 在编程的过程中可以切身体验到龙贝格算法是基于前几种算法的实现才能实现的。从而可以知道其实现可以运用前几种算法。五、教师评语（或成绩） 教师签字 ： 文中如有不足，请您见谅！5 / 5
合肥工业大学计算方法复化梯形公式实验
