数值计算方法龙贝格积分报告.doc

. . 姓名 班级 学号 【实验目的】： 1通过本实验加深对“龙贝格公式”的认识 2通过c++编写程序求解下列积分（E=1E-6） 3预期达到题中结果 4注意实验给出的程序有两个子程序为：f f1使用时只需该程序名即可。 【实验题目】： 复合求积公式计算定积分 【实验原理及理论基础】： 用梯形公式得到的积分近似值的误差是，因此，人们希望用这个误差作为对的一种补偿，则得到求积公式 的代数精度会有所提高。 通过直接验证可知 梯形公式二分前后的两个积分值与按如上公式做线性组合其结果正好是抛物线公式得到的积分值 同理可得到科茨公式积分值 和龙贝格公式 计算 【实验内容】： 龙贝格求积的计算步奏如下： 计算，，算出； 把[a,b]2等分，计算，算出与； 把[a,b]4等分，计算，，算出、与； 把[a,b]8等分，计算，算出、、与； 把[a,b]16等分，计算，算出、、与，继续重复进行，直到 （允许误差） 时停止运算，就是所求积分值。 输入 输入a,b 按公式计算 按公式 按公式(4－2)计算，，…， 输出近似值 输出近似值 非 非 输出失败信息 输出失败信息 【实验结果】： 【实验结果分析】： 实验结果比较理想，程序输出结果与预期结果相近。 本实验有助于我们学习龙贝格公式，掌握c++的使用方法进一步掌握结构化设计程序的好处（只需在外围改变子程序即可）。 下附程序： /*romberg method*/ #include stdio.h #include math.h #define n 1 float f(float x); float tt(int m,float t1,float t2); main() { float t[20][20]; int i,m,k; float p,a,b; printf(\nplease input zone [a,b]); scanf(%f,%f,a,b); p=f(a); for(i=1;in;i++) p=p+2*f(a+i*(b-a)/n); p=p+f(b); t[0][0]=(b-a)*p/(2*n); p=0; k=1; m=1; do { for(i=0;ipow(2,k-1);i++) p=p+f(a+(b-a)/(n*pow(2,k))+i*2*(b-a)/(n*pow(2,k))); t[0][k]=0.5*t[0][k-1]+(b-a)*p/(pow(2,k)*n); p=0; for(i=1;i=m;i++) t[i][k-i]=tt(i,t[i-1][k-i],t[i-1][k-i+1]); k=k+1; m=m+1; if(m=20) printf(m20,stop acclerate,I=%f,t[m-1][0]); } while((t[m-2][1]-t[m-2][0])1e-6); printf(\nthe integral value of sin(x)/x in zone [%f,%f] is I=%f,a,b,t[m-1][0]); /*printf(\nthe integral value of exp(-x^2) in zone [%f,%f] is I=%f,a,b,t[m-1][0]);*/ getch(); } float f(float x) { float y; if(x1e-10) y=1; else y=sin(x)/x; return y; } float tt(int m,float t1,float t2) { float tm; tm=(pow(4,m)*t2-t1)/(pow(4,m)-1); return tm; } float f1(float x) { float y; y=exp(-x*x); return y; }  程序由太原理工大学老师编写，我只做部分修改。