复合梯形公式与复合辛普森公式对比.docx

PAGE \* MERGEFORMAT11 SHANGHAI JIAO TONG UNIVERSITY 题目名称：复合梯形公式与复合辛普森公式对比 学生姓名: 学生学号: 班 级: 学院(系): 目录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc469831364 1. 概述 PAGEREF _Toc469831364 \h 3 HYPERLINK \l _Toc469831365 2. 问题提出 PAGEREF _Toc469831365 \h 4 HYPERLINK \l _Toc469831366 3. 算法推导 PAGEREF _Toc469831366 \h 5 HYPERLINK \l _Toc469831367 4. 算法框图 PAGEREF _Toc469831367 \h 6 HYPERLINK \l _Toc469831368 4.1复合梯形公式算法流程图 PAGEREF _Toc469831368 \h 6 HYPERLINK \l _Toc469831369 4.2 复合辛普森公式算法流程图 PAGEREF _Toc469831369 \h 7 HYPERLINK \l _Toc469831370 5. MATLAB源程序 PAGEREF _Toc469831370 \h 8 HYPERLINK \l _Toc469831371 6. 结论与展望 PAGEREF _Toc469831371 \h 9 图表目录 TOC \h \z \c 图 HYPERLINK C:\\Users\\wxf\\Desktop\\115010910167_王学范_复合梯形公式与复合辛普森公式对比.docx \l _Toc469827360 图 41 复合梯形公式算法流程图 PAGEREF _Toc469827360 \h 6 HYPERLINK C:\\Users\\wxf\\Desktop\\115010910167_王学范_复合梯形公式与复合辛普森公式对比.docx \l _Toc469827361 图 42 复合辛普森公式算法流程图 PAGEREF _Toc469827361 \h 7 HYPERLINK \l _Toc469827362 图 61 MATLAB计算结果 PAGEREF _Toc469827362 \h 9 TOC \h \z \c 表 HYPERLINK \l _Toc469827370 表 21函数计算结果表 PAGEREF _Toc469827370 \h 4 概述 梯形求积公式和辛普森求积公式分别是牛顿-科斯特公式中n=1和n=2时的情形。其中梯形求积公式可表示为 其公式左端是以[a,b]区间上积分，右端为b-a为高、端点函数值为上下底的梯形的面积值，故通称为梯形公式，具有1次代数精确度。 类似的，辛普森求积公式可以表示为 该公式一般在立体几何中用来求拟柱体的体积，由于偶数n阶牛顿-科特斯求积公式至少具有n+1次代数精确度，所以辛普森公式实际上具有3次代数精确度。 由于牛顿-科斯特公式在n≥8时不具有稳定性，故不可能通过提高阶的方法来提高求积精度。为了提高精度通常可把积分区间分成若干子区间（通常是等分），再在每个子区间上用低阶求积公式。这种方法称为复合求积法。 本文主要讨论复合梯形公式和复合辛普森公式在同一数学问题中的应用。首先给出了复合梯形公式和复合辛普森公式的推导过程以及其余项的表达形式，然后用流程图的形式介绍算法思路，再运用MATLAB编写代码计算结果，最后对结果进行对比讨论。 希望通过两个算法在同一个算例中的应用对比，更好的理解和掌握复合梯形公式和复合辛普森公式的适用范围和适用条件。并且能够熟悉MATLAB编程求解问题的流程，掌握编程化的思想方法。同时对两种方法的计算结果对比分析，讨论两种求积方法的计算精度。 问题提出 对于函数fx=sinx 表 STYLEREF 1 \s 2 SEQ 表 \* ARABIC \s 1 1函数计算结果表 x f(x) 0 1 1/8 0.997397867081822 1/4 0.989615837018092 3/8 0.989615837018092 1/2 0.958851077208406 5/8 0.936155636704740 3/4 0.908851680031112 7/8 0.877192573984031 1 0.841470984807897 算法推导 3.1复合梯形公式 根据梯形公式， 将区间[a,b]划分为n