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

数值分析第三次实习报告--第1页

学生学号实验课成绩

武汉理工大学

学生实验报告书

实验课程名称：数值分析(第三次实习)

开课名称：计算机科学与技术学院

指导老师姓名：熊盛武

学生姓名：

学生专业班级：软件工程0803班

2009——2010学年第一学期

数值分析第三次实习报告--第1页

数值分析第三次实习报告--第2页

一、实习目的：

（1）通过编程计算实践，体会和理解复化梯形公式和复化Simpson公式。

（2）通过编程计算实践，搞清变步长梯形公式的计算流程。

（3）通过编程计算实践，掌握和提高Romberg算法流程的控制技术。

（4）通过各种方法对同一题目的求解，体会各种方法的精度差异。

二、实习步骤：

（1）分别画出复化梯形公式，复化Simpson公式，变步长梯形公式和Romberg

公式的算法流程图。

（2）分别用复化梯形公式和复化Simpson公式通过编程计算积分，并分析和比

较算法的效率差异和精度差异。

（3）用变步长梯形公式通过编程计算积分，并分析算法精度与二分法次数之间

的关系。

（4）用Romberg公式通过编程计算积分，并分析算法精度与二分法次数之间的

关系。

三、算法流程图：

（一）复化梯形求算法积流程图

（三）变步长梯形求积算法流程图

开始（二）复化Simpson求积算法流程图

开始

输入区间a,b.和区间等分数n

求出步长h输入区间a,b.和区间等分数n求出步长h

Fori=0ton-1Fori=0ton-1temp1=f(a+i*h);

计算节点函数值f(a+i*h);temp2=f(a+(i+0.5)*h);

S1=S1+2*temp1;

S2=S2+4*temp2;

计算积分值T=T+2*temp;

Tn=(f(a)+f(b)+T)*h/2;

计算积分值

Sn=(f(b)-f(a)+S1+S2)*h/6;

输出积分值Tn

结束输出积分值

结束

数值分析