实验五:Romberg算法

上一次实验课我们练习使用复化求积公式的使用,先来看一下上一次实验课的内容:
实验内容:将[0,1]区间4等分,仅有等分点的函数值,分别用复化梯形公式、复化Simpson公式和复化Cotes公式计算定积分I=∫_01▒1/(1+x2 )dx的近似值.
实验结果:
复化梯形公式求得此定积分的近似值为:0.782794
复化Simpson公式求得此定积分的近似值为:0.785398
复化Cotes公式求得此定积分的近似值为:0.785529
本节课我们练习使用Romberg算法计算定积分I=∫_01▒1/(1+x2 )dx的近似值,要求误差小于10-6。
实验步骤:可以按照书上的学习步骤,先按不同的二分区间计算梯形公式的计算结果,可以先计算到T8,然后再得到S序列、C序列和R序列,误差的判断只需要判断两次的结果之间的误差就可以。

#include <stdio.h>
#include <math.h>
# define f(x) 1/(1+x*x)
void main()
{double t[5][5],a=0,b=1,temp,h,add;int n,k=0;t[0][0]=(b-a)*(f(a)+f(b))/2;for(k=1;k<=4;k++){n=int(pow(2,k));//n表示区间份数h=(b-a)/n;//h表示二分k次后的区间宽度t[k][0]=t[k-1][0]/2;add=pow(2,k-1);//add 表示每段区间新增加的节点for(int m=1;m<=add;m++){temp=(m-1)*(2*h)+h;t[k][0]+=(h*f(temp));}}//这一段是变步长的梯形公式/*for(int i=0;i<=4;i++)printf("t[%d][0]=%f\n",i,t[i][0]);*///接下来这一段是完成加速for(int j=1;j<=4;j++)//j表示列号{for(int i=j;i<=4;i++)//i表示行号{t[i][j]=(        );}}//这一段是按列输出结果for(j=0;j<=4;j++){for(int i=j;i<=4;i++){printf("t[%d][%d]=%f\n",i,j,t[i][j]);}}
}

实验五:Romberg算法相关推荐

  1. 数据结构实验五 查找算法的实现

    广州大学学生实验报告 开课实验室:计算机科学与工程实验(电子楼416B)     2019年6月11日 学院 计算机科学与教育软件学院 年级.专业.班 姓名 学号 实验课程名称 数据结构实验 成绩 实 ...

  2. 《机器学习与数据挖掘》实验五 编程实现误差逆传播算法(BP算法)

    前言: 摘要:本文对机器学习实验五 标准BP算法的代码进行实现,如果不了解的BP算法的话,可以自行上网搜索BP算法的详解. 实验题目:编程实现误差逆传播算法(BP算法) 实验目的:掌握误差逆传播算法( ...

  3. educoder算法设计与分析 实验五回溯法

    实验五 回溯法 第1关:排列 第2关:子集合 第3关:TSP问题 第4关:n皇后问题 第5关:0-1背包 第1关:排列 题目描述: 1.设计算法从前m个大写字母(m≤26)种取出n个字母的所有排列(组 ...

  4. 算法设计与分析 实验五 算法综合实验

    实验5.<算法综合实验> 一.实验目的 理解和复习所学各种算法的概念 掌握和复习所学各种算法的基本要素 掌握各种算法的优点和区别 通过应用范例掌握选择最佳算法的设计技巧与策略 二.实验内 ...

  5. SLR分析算法c语言实现,实验五 SLR语法分析器.doc

    实验五 SLR语法分析器 <编译原理>实验报告 实验序号: 05 实验项目名称: SLR语法分析器 学 号姓 名专业.班实验地点指导教师实验时间一.实验目的及要求 利用SLR文法的原理,掌 ...

  6. 20172329 2017-2018-2 《程序设计与数据结构》实验五报告

    这是这学期最后一次实验了,也是学到了很多东西,希望自己可以可以善始善终,加油! 让我们开始这一篇博客吧! 20172329 2017-2018-2 <程序设计与数据结构>实验五报告 课程: ...

  7. java 密钥工厂 desede_20145212 实验五《Java网络编程》

    20145212 实验五<Java网络编程> 一.实验内容1.运行下载的TCP代码,结对进行,一人服务器,一人客户端: 2.利用加解密代码包,编译运行代码,一人加密,一人解密: 3.集成代 ...

  8. 20175203 2018-2019 实验五《网络编程与安全》

    20175203 2018-2019 实验五<网络编程与安全> 知识重点(摘自实验资料) 栈 :(Stack)是一种只允许在表尾插入和删除的线性表,有先进后出(FILO),后进先出(LIF ...

  9. java 网络实验_Java实验五网络编程与安全

    实验五 网络编程与安全 实验准备 活动一 两人一组结对编程: 0. 参考http://www.cnblogs.com/rocedu/p/6766748.html#SECDSA 1. 结对实现中缀表达式 ...

最新文章

  1. c#如何将一个整数转换二进制,并进行位运算
  2. c#推箱子小游戏代码_推箱子小游戏V1.0制作
  3. plsql(轻量版)_异常处理机制
  4. excel 26进制 php,记录一次华为招聘的编程题-excel中的26进制
  5. 大数据技术周报第 005 期
  6. GIT学习----第六节:撤销修改
  7. NYOJ-数独(dfs)
  8. SPOJ LCMSUM - LCM Sum
  9. 软件漏洞-泛微OA-汇总
  10. yii2 init初始化脚本分析
  11. linux gif录制工具,教学?演示?在Linux下安装超好用的屏幕录像机来录制gif动画...
  12. 华为HCDA认证全套最新讲解视频教程(笔记)
  13. 常用颜色RGB、灰度值、取色值、透明度。
  14. web前端设计与开发大作业(一)----黄山风景区
  15. 万字详解车路协同、C-V2X通信协议
  16. FFmpeg将视频提取成帧
  17. Win7 系统下配置WinCE 5.0 模拟器网络环境
  18. IDEA 报错:Information:java: Errors occurred while compiling module ‘项目名‘
  19. 英语语言学论文选题 计算机,88个英语语言学、语言论文题目选题参考
  20. 栈应用 括弧匹配问题 题解(C++)

热门文章

  1. Messenger 麻烦了解一下
  2. 基于Java毕业设计在线图书超市源码+系统+mysql+lw文档+部署软件
  3. Python数据清洗-——(1)选择子集,列名重命名;(2)缺失值处理;(3)数据转换(重复值等);(4)异常值的处理
  4. 游戏玩家的图形显示设置指南(9)(The Gamer's Graphics Display Settings Guide)
  5. 小船过河 matlab,小船过河matlab实现.doc_蚂蚁文库
  6. java Synthetic
  7. 扫描识别系统能干些什么?
  8. mysql压力写入测试_mysql 压力测试
  9. C#上位机,信捷XD系列modbus485通信例子
  10. 渐变色怎么画?如何才能学好上色?