数值分析第三次实习报告.doc

学生学号 实验课成绩 武汉理工大学 学 生 实 验 报 告 书 实验课程名称： 数 值 分 析( 第三次实习) 开 课 名 称：计算机科学与技术学院 指导老师姓名： 熊盛武 学 生 姓 名： 学生专业班级：软件工程0803班 2009 —— 2010学年 第 一 学期 一、实习目的： （1）通过编程计算实践，体会和理解复化梯形公式和复化Simpson公式。 （2）通过编程计算实践，搞清变步长梯形公式的计算流程。 （3）通过编程计算实践，掌握和提高Romberg算法流程的控制技术。 （4）通过各种方法对同一题目的求解，体会各种方法的精度差异。 二、实习步骤： （1）分别画出复化梯形公式，复化Simpson公式，变步长梯形公式和Romberg公式的算法流程图。 （2）分别用复化梯形公式和复化Simpson公式通过编程计算积分，并分析和比较算法的效率差异和精度差异。 （3）用变步长梯形公式通过编程计算积分，并分析算法精度与二分法次数之间的关系。 （4）用Romberg公式通过编程计算积分，并分析算法精度与二分法次数之间的关系。 三、算法流程图： （一）复化梯形求算法积流程图 （三）变步长梯形求积算法流程图 （二）复化Simpson求积算法流程图  （四）Romberg公式求积流程图 四、编程计算积分 （一）复化梯形： （1）复化梯形算法描述 （1）输入求积区间a,b和区间等分数n （2）计算步长h h=(b-a)/n （3）对于i=0到i=n-1计算节点函数值f(a+i*h) （4）计算积分值T=T+2*temp; Tn=(f(a)+f(b)+T)*h/2; （5）输出积分值 （2）复化梯形求积程序源代码 #includeiostream #includecmath double f(double t) { return 1/(1+t*t); } using namespace std; int main() { int i,n; double h,a,b,T,temp,Tn; cout请输入求积区间:endl; cinab; cout请输入区间等分数:endl; cinn; h=(b-a)/n; temp=0; T=0; for(i=1;in;i++) { temp=f(a+i*h); T=T+2*temp; } Tn=(f(a)+f(b)+T)*h/2; cout经过复化梯形求积所得积分值Tn为: Tnendl; return 0; } （3）调试过程，实验结果及分析： 1/（x*x）在区间[0,1]上8等分后,程序运行结果如下: 对区间等分后求每个子区间上的积分值，然后将每个子区间的积分值相加，就得到整个积分区间上的积分值，梯形公式具有一次代数精度。 （二）复化Simpson： （1）复化Simpson算法描述 （1）输入求积区间a,b和区间等分数n （2）计算步长h h=(b-a)/n （3）对于i=0到n-1 计算节点函数值temp1=f(a+i*h); temp2=f(a+(i+0.5)*h); S1=S1+2*temp1; S2=S2+4*temp2; （4）计算积分值 Sn=(f(b)-f(a)+S1+S2)*h/6; （5）输出积分值. （2）复化Simpson算法程序源代码 #includeiostream #includecmath double f(double t) { return 1/(1+t*t); } using namespace std; int main() { int i,n; double h,a,b,S1,S2,temp1,temp2,Sn; cout请输入求积区间:endl; cinab; cout请输入区间等分数:endl; cinn; h=(b-a)/n; temp1=0; temp2=0; S1=0; S2=0; for(i=0;in;i++) { temp1=f(a+i*h); temp2=f(a+(i+0.5)*h); S1=S1+2*temp1; S2=S2+4*temp2; } Sn=(f(b)-f(a)+S1+S2)*h/6; cout经过复化Simpson求积所得积分值Sn为: Snendl; return 0; } （3）调试过程，实验结果及分析： 1/