合肥工业大学计算方法复化梯形公式实验
更新时间:2024-05-14 00:39:01 阅读量: 综合文库 文档下载
- 合肥工业大学计算方法试卷推荐度:
- 相关推荐
《计算方法》实验报告
学号 实验项目名称 一、实验名称 姓名 班级 实验二 数值积分 实验二 数值积分 二、实验目的: (1) 熟悉复化梯形方法、复化Simpson方法、梯形递推算法、龙贝格算法; (2) 能编程实现复化梯形方法、复化Simpson方法、梯形递推算法、龙贝格算法; (3) 理解并掌握自适应算法和收敛加速算法的基本思想; (4) 分析实验结果体会各种方法的精确度,建立计算机求解定积分问题的感性认识 三、实验内容及要求 (1)设计复化梯形公式求积算法,编制并调试相应的函数子程序 (2)设计复化辛浦生求积算法,编制并调试相应的函数子程序 (3)用龙贝格算法计算?10sinxdx x输入:积分区间,误差限 输出:序列Tn,Sn,Cn,Rn及积分结果(参考书本P81的表2-5) 取n=2,4,8,16,精确解为0.9460831 四、实验原理及算法描述 在许多实际问题中,常常需要计算定积分?af(x)dx的值。根据微积分学基本定理,若被积函b数f(x)在区间[a,b]上连续,只要能找到f(x)的一个原函数F(x),便可利用牛顿-莱布尼兹公式?af(x)?F(b)?F(a)求得积分值。 但是在实际使用中,往往遇到如下困难,而不能使用牛顿-莱布尼兹公式。 (1) 找不到用初等函数表示的原函数 (2) 虽然找到了原函数,但因表达式过于复杂而不便计算 (3) f(x)是由测量或计算得到的表格函数 由于以上种种困难,有必要研究积分的数值计算问题。 b计算机科学与工程学院
b利用插值多项式 Pn(x)?f(x)则积分?af(x)dx转化为?P(x)dx,显然易算。?P(x)dx称为ananbb插值型求积公式。最简单的插值型求积公式是梯形公式和Simpson公式,。当求积结点提供较多,可以分段使用少结点的梯形公式和Simpson公式,并称为复化梯形公式、复化Simpson公式。如步长未知,可以通过误差限的控制用区间逐次分半的策略自动选取步长的方法称自适应算法。梯形递推公式给出了区间分半前后的递推关系。由梯形递推公式求得梯形序列,相邻序列值作线性组合得Simpson序列, Simpson序列作线性组合得柯特斯序列, 柯特斯序列作线性组合的龙贝格序列。若|R2-R1|,则输出R2;否则…依此类推。如此加工数据的过程叫龙贝格算法,如下图所示: n?11复化梯形公式Tn?h[f(x0)?2?f(xk)?f(xn)] 2k?1n?1n?1h复化Simpson公式Sn?[f(x0)?2?f(xk)?f(xn)?4?f(x1)] k?6k?1k?02Tnhn?1梯形递推公式T2n???f(x1) k?22k?024T2n?Tn42S2n?Sn43C2n?Cn?Sn ?Cn ?Rn 加权平均公式:4?142?143?1龙贝格算法大大加快了误差收敛的速度,由梯形序列O(h2) 提高到龙贝格序列的O(h8) 五、程序代码及实验结果 1. 主程序 int main() { cout << \ << \ << \ << \ << \ << \ << \ << \ << cout << \ << \ < 2 计算机科学与工程学院 0.00000001) << endl;; cout << \<<4<<\ << longbeige(0, 1, 0.00000001); return 0; } 2. 复化梯形公式子程序: double tixing(double a, double b, int n) { } } double res = (h / 2.0)*(fa + 2.0*fxk + fb); return res; double fa = f(a); double fb = f(b); double h = (b-a) / n; double fxk = 0.0; for (int k = 1; k <= n - 1; k++) { double xk = a + k*h; fxk = f(xk) + fxk; 3. 复化辛浦生公式子程序: double xingbusheng(double a, double b, int n) { double h = (b - a) /n; double fa = f(a); double fb = f(b); //double s = fb - fa; double x = a; double fxk12 = 0.0; double fxk = 0.0; /*for (int k = 1; k <=n; k++) { } s = (h / 6.0)*s; return s;*/ 3 x = x + h / (2.0); double fx = f(x); s = s + 4.0*f(x); x = x + h / (2.0); s = s + 2.0*f(x); 计算机科学与工程学院 } double s = (h / 6)*(fa + 4 * fxk12 + 2 * fxk + fb); return s; for (int k = 1; k <= n - 1; k++) { } double xk = a + k*h; fxk = fxk + f(xk); for (int k = 0; k <= n - 1; k++) { } double xk1 = a + k*h; x = xk1 + h / 2.0; fxk12 = fxk12 + f(x); 4. 龙贝格公式子程序: double longbeige(double a, double b, double wuchaxian) { S2 = T2 + (1 / 3)*(T2 - T1); if (k == 1) { Tck[k] = T2; k++; h = h / 2; T1 = T2; 4 double h = b - a; double T1 = (h / 2.0)*(f(a) + f(b)); int k = 1; double S,x,T2,S2,S1=0,C1=0,C2,R1=0,R2; double Tck[100],Sck[100],Cck[100],Rck[100]; S = 0; x = a + h / 2.0; while (x } } R2 = C2 + (1 / 63)*(C2 - C1); if (k == 3) { } if (fabs(R2 - R1) >= wuchaxian) { } else { } return R2; R1 = R2; C1 = C2; k++; h = h / 2; T1 = T2; S1 = S2; goto loop1; Tck [k]= T2; R1 = R2; C1 = C2; k++; h = h / 2; T1 = T2; S1 = S2; goto loop1; h = h / 2; T1 = T2; S1 = S2; goto loop1; } C2 = S2 + (1 / 15)*(S2 - S1); if (k == 2) { C1 = C2; Tck[k]=T2; k++; S1 = S2; goto loop1; 5 计算机科学与工程学院 实验结果: 图1 六、实验总结 1. 梯形公式的收敛速度太慢,所以我们才会选择之后的几种公式加快收敛速度。 2. 通过实验可以看出,龙贝格算法大大加快了误差收敛的速度,由梯形序列O(h2) 提高到龙贝格序列的O(h8) 3. 在编程的过程中可以切身体验到龙贝格算法是基于前几种算法的实现才能实现的。从而可以知道其实现可以运用前几种算法。 五、教师评语(或成绩) 教师签字 : 6 计算机科学与工程学院 实验结果: 图1 六、实验总结 1. 梯形公式的收敛速度太慢,所以我们才会选择之后的几种公式加快收敛速度。 2. 通过实验可以看出,龙贝格算法大大加快了误差收敛的速度,由梯形序列O(h2) 提高到龙贝格序列的O(h8) 3. 在编程的过程中可以切身体验到龙贝格算法是基于前几种算法的实现才能实现的。从而可以知道其实现可以运用前几种算法。 五、教师评语(或成绩) 教师签字 : 6
正在阅读:
合肥工业大学计算方法复化梯形公式实验05-14
高考化学重要知识点总结!04-05
民用建筑设计通则07-19
概率论与数理统计及其应用课后答案(浙江大学 - 盛骤版)12-14
马哲论文 当代大学生素质调查报告06-13
洞口县醪田镇中心小学为民办实事汇报材料01-13
绩效制定原则制度05-30
研究注册商标和社会主义道德风尚的关系12-11
微信小程序开发合同协议书09-20
- 清真菜谱
- 我国国民经济和社会发展十二五规划纲要(全文)
- 高三物理机械振动和机械波复习2
- 浙江省公路山岭隧道机械化装备应用指导手册 doc - 图文
- 2018届高三数学文科二轮复习:专题检测(九) 导数的简单应用
- 2015年上海市公务员录用考试《行政职业能力测验》试卷(B类)
- 七年级道德与法制下册
- 大班户外游戏教案
- 病虫害预警 - 图文
- 某养鱼场为了提高经营管理水平
- 汉中市勉县尧柏余热汽机规程 10
- 烹饪试卷
- 事业单位考试公共基础知识专项分类题库训练
- 语文:第2课 走一步,再走一步 课堂导学案(人教版 七上)
- 天汉使用手册
- 人教版小学三年级数学下册教学计划
- 房地产销售管理完全操作手册122页
- 2009年评审通过具有中学高级教师专业技术资格人员名单...
- 《15秋公共关系学》作业1
- 2017最新版监理公司三标一体管理手册
- 合肥工业大学
- 梯形
- 公式
- 实验
- 计算
- 方法