合肥工业大学宣城校区计算方法实验二

1实验目的和要求

实验目的：
(1)熟悉复化梯形方法、复化Simpson方法、梯形递推算法、龙贝格算法；
(2)能编程实现复化梯形方法、复化Simpson方法、梯形递推算法、龙贝格算法；
(3)理解并掌握自适应算法和收敛加速算法的基本思想；
(4)分析实验结果体会各种方法的精确度，建立计算机求解定积分问题的感性认识

程序任务：
（1）设计复化梯形公式求积算法，编制并调试相应的函数子程序
（2）设计复化辛卜生求积算法，编制并调试相应的函数子程序
（3）用龙贝格算法计算
输入：积分区间，误差限
输出：序列Tn，Sn，Cn，Rn及积分结果（参考书本P55的表2-5）

取n=2，4，8，16，精确解为0.9460831

2实验环境和工具

用matlab或C语言实现

3实验结果

3.1算法流程图

复化辛卜生公式算法流程：

3.2程序核心代码

#include<stdio.h>
#include<math.h>
double ctl(double a,double b,int n);
double smp(double a,double b,int n);
void lbg(double a,double b,double e);
int main(void)
{int m,n;double a,b,s,e;printf("1.复化梯形公式\n2.复化辛卜生\n3.龙贝格算法\n\n");scanf("%d",&m);switch(m){case 1:{printf("输入下限、上限、分段数。\n\n");scanf("%lf %lf %d",&a,&b,&n);s = ctl(a,b,n);printf("%lf\n",s);break;}case 2:{printf("输入下限、上限、分段数。\n\n");scanf("%lf %lf %d",&a,&b,&n);s = smp(a,b,n);printf("%lf\n",s);break;}case 3:{printf("输入下限、上限、误差限。\n\n");scanf("%lf %lf %lf",&a,&b,&e);lbg(a,b,e);}}return 0;
}double f(double x)//原函数f(x)
{if(x==0){return 1;}else{return sin(x)/x;}
}double ctl(double a,double b,int n)//复化梯形公式
{double h = (b-a)/n;double t = 0;int i;for(i=1;i<n;i++){t = t+f(a+i*h);}return h/2*(f(a)+2*t+f(b));
}double smp(double a,double b,int n)//复化辛卜生
{double h = (b-a)/n;double x1 = f(a+h/2);double x2 = 0;int i;for(i=1;i<n;i++){x1 = x1+f(a+i*h+h/2);x2 = x2+f(a+i*h);}return h/6*(f(a)+4*x1+2*x2+f(b));
}void lbg(double a,double b,double e)//龙贝格算法
{double h = b-a;double t1,t2=0,s1=0,s2=0,c1=0,c2=0,r1=0,r2=1;double x,s;int k = 1;t1 = h/2*(f(a)+f(b));printf("\n%.7lf\n",t1);while(fabs(r2-r1)>=e){s = 0;x = a+h/2;while(x<b){s = s+f(x);x = x+h;}t2 = t1/2+h/2*s;s2 = t2+(t2-t1)/3;if(k==1){k = k+1;h = h/2;t1 = t2;s1 = s2;printf("%.7lf\t%.7lf\n",t1,s1);continue;}c2 = s2+(s2-s1)/15;if(k==2){k = k+1;h = h/2;t1 = t2;s1 = s2;c1 = c2;printf("%.7lf\t%.7lf\t%.7lf\n",t1,s1,c1);continue;}r2 = c2+(c2-c1)/63;k = k+1;h = h/2;t1 = t2;s1 = s2;c1 = c2;r1 = r2;printf("%.7lf\t%.7lf\t%.7lf\t%.7lf\n",t1,s1,c1,r1);}
}

3.3运行结果

3.4运行结果分析
根据运行结果可以得出结论：
复化梯形公式、复化辛卜生公式、龙贝格算法的精确度是逐渐提高的，复化梯形公式的精确度较低，复化辛卜生公式的精度较前者有所提高，龙贝格算法的精度是三者之中最高的。

4实验心得

在本次实验之中，我又一次回顾了复化梯形公式、复化辛卜生公式和龙贝格算法。这三种算法不仅仅在数学上有很大的作用，在计算机编程方面也有极大的帮助。
我认为在计算机编程的过程中，算法是极为重要的一环，甚至可以说是最为重要的一环，没有好的算法，就难以构造好的函数，就难以得到符合要求的程序。而计算方法课程教会了我们很多好的算法，这些算法对于一个合格的程序员是必不可少的。将这三种算法牢记于心，对我们编程有莫大的好处。
这次实验开始的时候，我对于三种算法的印象已经有些模糊了，而这次实验让我重新回忆起了这三种算法，这也给我敲响了警钟，计算方法课程上的知识对于我们是必备的，我们应该经常回顾这些知识，温故而知新，只有在不断固化对过去知识的记忆，并去学习新的知识，才能在未来做到不后悔。
每一次实验，都是对过去知识的回顾，同样也是在应用过去的知识，我应该珍惜每一次的实验机会，利用好每一次实验，充分提高自己，巩固根基，才能在未来的学习生活中取得更好的成绩。