内容发布更新时间 : 2024/11/7 23:48:19星期一 下面是文章的全部内容请认真阅读。
SHANGHAI JIAO TONG UNIVERSITY
题目名称:复合梯形公式与复合辛普森公式对比
学生姓名: 学生学号: 班 级: 学院(系):
1
目录
1. 2. 3. 4.
概述 .............................................................................................................. 3 问题提出 ...................................................................................................... 4 算法推导 ...................................................................................................... 5 算法框图 ...................................................................................................... 6 4.1复合梯形公式算法流程图 ...................................................................... 6 4.2 复合辛普森公式算法流程图 ................................................................. 7 5. 6.
MATLAB源程序 ......................................................................................... 8 结论与展望 .................................................................................................. 9
图表目录
图 4-1 复合梯形公式算法流程图 ...................................................................... 6 图 4-2 复合辛普森公式算法流程图 .................................................................. 7 图 6-1 MATLAB计算结果 .................................................................................. 9
表 2-1函数计算结果表 ....................................................................................... 4
2
1. 概述
梯形求积公式和辛普森求积公式分别是牛顿-科斯特公式中n=1和n=2时的情形。其中梯形求积公式可表示为
?baf(x)dx?b?a(f(a)?f(b))2 其公式左端是以[a,b]区间上积分,右端为b-a为高、端点函数值为上下底的梯形的面积值,故通称为梯形公式,具有1次代数精确度。 类似的,辛普森求积公式可以表示为 S?b?aa?b[f(a)?4f()?f(b)]62 该公式一般在立体几何中用来求拟柱体的体积,由于偶数n阶牛顿-科特斯求积公式至少具有n+1次代数精确度,所以辛普森公式实际上具有3次代数精确度。 由于牛顿-科斯特公式在n≥8时不具有稳定性,故不可能通过提高阶的方法来提高求积精度。为了提高精度通常可把积分区间分成若干子区间(通常是等分),再在每个子区间上用低阶求积公式。这种方法称为复合求积法。
本文主要讨论复合梯形公式和复合辛普森公式在同一数学问题中的应用。首先给出了复合梯形公式和复合辛普森公式的推导过程以及其余项的表达形式,然后用流程图的形式介绍算法思路,再运用MATLAB编写代码计算结果,最后对结果进行对比讨论。
希望通过两个算法在同一个算例中的应用对比,更好的理解和掌握复合梯形公式和复合辛普森公式的适用范围和适用条件。并且能够熟悉MATLAB编程求解问题的流程,掌握编程化的思想方法。同时对两种方法的计算结果对比分析,讨论两种求积方法的计算精度。
3