数值分析实验报告(华科书本实验4.1-附C++程序).docVIP

华中科技大学 数值分析实验报告 考生姓名 考生学号 班 级 指导老师 路志宏 2013年4月15日实验4.1 实验目的：复化求积公式计算定积分 试验题目：数值计算下列各式右端定积分的近似值。 （1）； （2）； （3）； （4）； 实验要求： （1）若用复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式做计算，要求绝对误差限为，分别利用他们的余项对每种算法做出步长的事前估计。 （2）分别用复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式做计算。 （3）将计算结果与精确解做比较，并比较各种算法的计算量。 实验内容： 1.公式介绍 （1）复化梯形公式： =； 余项：； （2）复化Simpson公式： =； 余项：； （3）复化Gauss-Legendre I型公式： ； 余项：； 2.步长估计 （1）； 则可以得到：；； ； 估计步长：； 将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到： 复化梯形公式：；； 复化Simpson公式：；； 复化Gauss-Legendre I型公式：；； （2）； 则可以得到：；； ； 估计步长：； 将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到： 复化梯形公式：；； 复化Simpson公式：；； 复化Gauss-Legendre I型公式：；； （3）； 则可以得到：；； ； 估计步长：； 将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到： 复化梯形公式：；； 复化Simpson公式：；； 复化Gauss-Legendre I型公式：；； （4）； 则可以得到：；； ； 估计步长：； 将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到： 复化梯形公式：；； 复化Simpson公式：；； 复化Gauss-Legendre I型公式：；； 3.C++编程计算结果 （1）区间逐次分半求积法： 依据“事后误差法”，将区间逐次分半进行计算，并利用前后两次计算结果来判断误差的大小。在逐次二分进行计算时，可以用与来估计误差，这种直接用计算结果来估计误差的方法通常称作误差的事后估计法，若（为计算结果允许的误差），则停止计算，并取作为积分的近似值；否则将区间再次二分后算出，并检验不等式是否满足。 由于是区间分半，因此区间等分数必定是2的n次方。 结束输出结果选择题号N复化公式函数选择f(x)开始（2）流程图： 结束 输出结果 选择题号N 复化公式函数 选择f(x) 开始 （3）计算结果及误差： 分别对4题作复化Trapezoid、Simpson、Gauss_Legendre计算，并计算计算值与精确值之间的误差，结果如下表： （1）计算结果表 数据类型 求积类型 （1） 计算值 区间二分相邻误差 绝对误差 区间二分等分 复化Trapezoid公式 -02043509 4.18e-008 1.3935e-008 2^11=2048 复化Simpson公式 -09952908 2.76e-008 1.8447e-009 2^5=32 Gauss_Legendre I公式 -06878386 1.84e-008 1.2298e-009 2^5=32 （2）计算结果表 数据类型 求积类型 （2） 计算值 区间二分相邻误差 绝对误差 区间二分等分 复化Trapezoid公式 3.1415926 1385336 1.19e-007 3.9736e-008 2^11=2048 复化Simpson公式 3.1415926 5122482 1.49e-007 2.365e-009 2^3=8 Gauss_Legendre I公式 3.1415926 5527155 1.06e-007 1.6818e-009 2^3=8 （3）计算结果表 数据类型 求积类型 （3） 计算值 区间二分相邻误差 绝对误差 区间二分等分 复化Trapezoid公式 1.820478 4969086 1.31e-007 4.3655e-008 2^11=2048 复化Si