数值分析实验报告.doc

华中科技大学 数值分析实验报告 考生姓名 考生学号 班 级 指导老师 路志宏 2013年4月15日实验 实验目的：复化求积公式计算定积分 试验题目：数值计算下列各式右端定积分的近似值。 （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型公式：；； ++编程计算结果 （1）区间逐次分半求积法： 依据“事后误差法”，将区间逐次分半进行计算，并利用前后两次计算结果来判断误差的大小。在逐次二分进行计算时，可以用与来估计误差，这种直接用计算结果来估计误差的方法通常称作误差的事后估计法，若（为计算结果允许的误差），则停止计算，并取作为积分的近似值；否则将区间再次二分后算出，并检验不等式是否满足。 由于是区间分半，因此区间等分数必定是2的n次方。 结束输出结果选择题号N复化公式函数选择f(x)开始（2）流程图： 结束 输出结果 选择题号N 复化公式函数 选择f(x) 开始 （3）计算结果及误差： 分别对4题作复化Trapezoid、Simpson、Gauss_Legendre计算，并计算计算值与精确值之间的误差，结果如下表： （1）计算结果表 数据类型 求积类型 （1） 计算值 区间二分相邻误差 绝对误差 区间二分等分 复化Trapezoid公式 -0. 2043509 2^11=2048 复化Simpson公式 -0. 9952908 2^5=32 Gauss_Legendre I公式 -0. 6878386 2^5=32 （2）计算结果表 数据类型 求积类型 （2） 计算值 区间二分相邻误差 绝对误差 区间二分等分 复化Trapezoid公式 1385336 2^11=2048 复化Simpson公式 5122482 2^3=8 Gauss_Legendre I公式 5527155 2^3=8 （3）计算结果表 数据类型 求积类型 （3） 计算值 区间二分相邻误差 绝对误差 区间二分等分 复化Trapezoid公式 4969086 2^11=2048 复化Simpson公式 4541318 2^5=32 Gauss_Legendre I公式 4438881 2^4=16 （4）计算结果表 数据类型 求积类型 （4） 计算值 区间二分相邻误差 绝对误差 区间二分等分 复化Trapezoid公式 1970612 2^13=8192 复化Simpson公式 0756376 2^5=32 Gauss_Legendre I公式 9317526 2^5=32 由上表中的误差分析可知，利用题目所要求的