数值逼近实习报告八.docVIP

数 值 逼 近 实 习 报 告 班级 姓名 指导教师 实习题目 龙贝格积分算法 计算公式： 目的意义 掌握龙贝格积分公式的计算方法 数学模型 如果将区间等分，则复化梯形积分公式为 现变步长，将区间等分，则有 令 算法描述 对于给定的和初始 Step1：用复化梯形公式 求得外推出 Step2: 对于已经计算出来的 计算出 Step3: 用外推法一次计算 Step4： 输出 （流程图）程序 #include <iostream.h> #include <math.h> const int k=4; //次数 const double a=0; //起点 const double b=1; //终点 double h=b-a; //步长 double (*data)[k]=new double[k][k]; double f(const double &x); //求积函数 double sum(const int &n,const double &h); //多项式求和 int main() { cout<<" ///////////////////////////////////////"<<endl; cout<<" // Ronberg积分公式 //"<<endl; cout<<" ///////////////////////////////////////"<<endl; cout<<"利用Ronberg积分公式计算圆周率π"<<endl; cout<<"步长为(b-a)的时候积分结果为："<<endl; data[0][0]=h/2*(f(a)+f(b)); cout<<"T0="; cout<<data[0][0]<<endl; cout<<"变步长之后的计算结果为："<<endl; for (int i=1;i<k;++i) { data[i][0]=data[i-1][0]/2+h/2*sum(i,h); //复化梯形积分计算公式 h=h/2; } for(i=1;i<k;i++) for (int j=i;j<k;j++) data[j][i]=(pow(4,i)*data[j][i-1]-data[j-1][i-1])/(pow(4,i)-1); for (i=0;i<k;i++) { for (int j=0;j<=i;j++) cout<<"T"<<j+1<<"="<<data[i][j]<<' '; cout<<endl; } } double f(const double &x) { return 4/(1+x*x); } double sum(const int &n,const double &h) { int t=pow(2,n-1); double s=0; for (int i=1;i<=t;i++) { s=s+f(a+(2*i-1)*(h/2)); } return s; } 数值算例 对计算结果进行分析评价 参考文献 [1] ***，数值逼近，***：****大学***学院，20**.