龙贝格数值分析作业c语言,数值分析龙贝格实验报告.doc

数值分析龙贝格实验报告

实验三龙贝格方法

【实验类型】验证性

【实验学时】 2学时

【实验内容】

1.理解龙贝格方法的基本思路

2.用龙贝格方法设计算法，编程求解一个数值积分的问题。

【实验前的预备知识】

1．计算机基础知识2．熟悉编程基本思想3．熟悉常见数学函数；

【实验方法或步骤】

龙贝格方法的基本思路龙贝格方法是在积分区间逐次二分的过程中，通过对梯形之值进行加速处理，从而获得高精度的积分值。

龙贝格方法的算法

步骤1 准备初值和，用梯形计算公式计算出积分近似值

步骤2 按区间逐次分半计算梯形公式的积分近似值令

，计算，

步骤3 按下面的公式积分梯形公式：

辛普生公式：龙贝格公式：

步骤4 精度控制

当，(为精度)时，终止计算，并取为近似值否则将步长折半，转步骤2。

[实验程序]

#include

# define Precision 0.00001//积分精度要求

# define e 2#define MAXRepeat 10 //最大允许重复

double function(double x)//被积函数

{

double s;

s=2*pow(e,-x)/sqrt(3.1415926);

return s;

}

double Romberg(double a,double b,double f(double x))

{

int m,n,k;

double y[MAXRepeat],h,ep,p,xk,s,q;

h=b-a;

y[0]=h*(f(a)+f(b))/2.0;//计算T`1`(h)=1/2(b-a)(f(a)+f(b));

m=1;

n=1;

ep=Precision+1;

while((ep>=Precision)&&(m

{

p=0.0;

for(k=0;k

{

xk=a+(k+0.5)*h; // n-1

p=p+f(xk); //计算∑f(xk+h/2),T

} // k=0

p=(y[0]+h*p)/2.0; //T`m`(h/2),变步长梯形求积公式

s=1.0;

for(k=1;k<=m;k++)

{

s=4.0*s;// pow(4,m)

q=(s*p-y[k-1])/(s-1.0);//[pow(4,m)T`m`(h/2)-T`m`(h)]/[pow(4,m)-1],2m阶牛顿柯斯特公式，即龙贝格公式

y[k-1]=p;

p=q;

}

ep=fabs(q-y[m-1]);//前后两步计算结果比较求精度

m=m+1;

y[m-1]=q;

n=n+n; // 2 4 8 16

h=h/2.0;//二倍分割区间

}

return q;

}

main()

{

double a,b,Result;

cout<

cin>>a;

cout<

cin>>b;

Result=Romberg( a, b, function);

cout<

return 0;

}

龙贝格数值分析作业c语言,数值分析龙贝格实验报告.doc相关推荐

c语言课设航班信息管理软件管理系统,c语言航班信息查询系统实验报告.doc
c语言航班信息查询系统实验报告.doc 下载提示(请认真阅读)1.请仔细阅读文档,确保文档完整性,对于不预览.不比对内容而直接下载带来的问题本站不予受理. 2.下载的文档,不会出现我们的网址水印. 3 ...
c语言图形时钟编程,C语言图形时钟课程设计实验报告.doc
C语言图形时钟课程设计实验报告.doc 目录 1.系统功能要求 . 2. 数据结构设计及说明 . 3.程序结构(画流程图) . 4.各模块的功能 . 5.试验结果(包括输入数据和输出结果) . 6. ...
c语言应用报告,C语言：数组的应用实验报告.doc
C语言:数组的应用实验报告.doc -1- 课程名称课程名称结构化程序设计与结构化程序设计与 C C 语言语言开课实验室开课实验室年年月月日日年级年级.专业.班.专业.班学号学号姓名姓名 ...
c语言程序设计航空查询系统,c语言航班信息查询系统实验报告.doc
c语言航班信息查询系统实验报告 PAGE 13 - 软件学院大作业任务书题目:航班信息查询系统专业: 班级: 姓名: 学号: 完成人数: 起讫日期: 任课教师: 职称: 讲师部分管主任 ...
数值分析matlab实验报告,数值分析第一次作业matlab实验报告.doc
数值分析第一次作业matlab实验报告.doc 几种线性方程组迭代算法的MATLAB实现和性能比较用有限差分方法(五点差分格式)求解正方形域上的Poisson方程边值问题用MATLAB语言编写算法程序 ...
北京电大c语言实验作业二,大学大一c语言程序设计实验室上机题全部代码答案(实验报告).doc...
大学大一c语言程序设计实验室上机题全部代码答案(实验报告).doc C语言实验报告实验1-1:helloworld程序:源代码:#includemain(){printf("hellowor ...
c语言学生档案管理课设作业,2019-2020年c语言课程设计学生档案管理系统实验报告.doc...
2019-2020年c语言课程设计学生档案管理系统实验报告.doc 还剩 10页未读, 继续阅读下载文档到电脑,马上远离加班熬夜! 亲,喜欢就下载吧,价低环保! 内容要点: *********C 语 ...
c语言飞机票实验报告格式,c语言航班信息查询系统实验报告
<c语言航班信息查询系统实验报告>由会员分享,可在线阅读,更多相关<c语言航班信息查询系统实验报告(21页珍藏版)>请在人人文库网上搜索. 1.希昌丸會软件学院大作业任务书题目 ...
c程序语言第四版实验报告,C语言程序设计第四次实验报告
C语言程序设计第四次实验报告姓名:熊毅实验地点:家实验时间:2020.04.09 实验项目: 5.3.1练习2 求数列的前n项和 5.3.2练习2 求水仙花数 5.3.4 十进制转换 5.3.5 ...

龙贝格数值分析作业c语言,数值分析龙贝格实验报告.doc

龙贝格数值分析作业c语言,数值分析龙贝格实验报告.doc相关推荐

最新文章

热门文章