梯形多步法和辛普森积分
分别利用矩形法梯形法辛普森法对定积分进行近似计算并比较计算效果
分别利用矩形法梯形法辛普森法对定积分进行近似计算并比较计算效果定积分是微积分中重要的概念之一，表示在一个区间上函数的面积。
在计算定积分时，有时候我们无法通过解析方法求得精确的结果，这时候可以利用数值方法来进行近似计算。
常见的数值方法包括矩形法、梯形法和辛普森法。
本文将分别对这三种方法进行介绍并进行比较。
1.矩形法（矩形近似法）：矩形法是最简单的数值方法之一，它的基本思想是将函数曲线上每个小区间的面积近似为一个矩形的面积，然后将这些矩形的面积相加，即可得到函数曲线下的面积。
根据矩形法的计算公式可以得到：∫f(x)dx ≈ Δx·(f(x₁)+f(x₂)+...+f(xₙ))其中，Δx为区间的长度，f(x)为函数在区间上的值。
2.梯形法（梯形近似法）：梯形法同样是利用近似的思想，将函数曲线上每个小区间的面积近似为一个梯形的面积，然后将这些梯形的面积相加，即可得到函数曲线下的面积。
梯形法的计算公式为：∫f(x)dx ≈ （Δx/2）·[f(x₀)+2f(x₁)+2f(x₂)+...+2f(xₙ-1)+f(xₙ)]其中，Δx为区间的长度，f(x)为函数在区间上的值。
3.辛普森法（抛物线近似法）：辛普森法是一种基于三次多项式插值的数值积分方法，它通过将函数曲线上每个小区间的面积近似为一个抛物线的面积，然后将这些抛物线的面积相加，即可得到函数曲线下的面积。
辛普森法的计算公式为：∫f(x)dx ≈ （Δx/3）·[f(x₀)+4f(x₁)+f(x₂)+4f(x₃)+...+4f(xₙ-1)+f(xₙ)]其中，Δx为区间的长度，f(x)为函数在区间上的值。
例：计算函数f(x)=√(1+x²)在区间[0,1]上的定积分。
接下来，我们分别利用矩形法、梯形法和辛普森法对这个定积分进行近似计算，并比较计算结果。
1)矩形法：将区间[0,1]平均分为n个小区间，取xᵢ=i/n，其中i=0,1,2,...,n。
梯形公式和辛普森求解定积分
梯形公式和辛普森求解定积分
数学中，定积分是较为常见的运算，既可以通过梯形公式、辛普森公式等方式求解。
梯形公式是用来计算定积分的一种常用方法，主要就是把封闭的积分区间[a,b]分
成若干等分，每一等分长度相等，每一等分的两端点函数值分别用各积分一次积分而得小梯形面积来等价近似。
然后，将所得的小梯形面积加起来，即可求到积分的近似值。
辛普森求解定积分是将积分区间[a,b]表示为一组有限个点，然后利用辛普森公式
来近似计算函数在该组点上的值，最后加起来就可以得到整个积分区间的值。
一般而言，当积分区间越小而越窄，辛普森公式所得的积分结果的接近的越精确，且求解的速度最快。
定积分的求解方法有多种，梯形公式和辛普森求解定积分就是其中的两种求解方式，一般情况下，梯形公式用于在积分区间中间定点数多，但是积分段数相对较少的情况下，而辛普森求解定积分则用于积分区间窄，但是积分段数稍微多的情况下。
由于梯形公式和辛普森求解定积分有其各自的优点，在实际应用中可以根据不同的情况，灵活选用二者的优点，以达到最优的结果。
(完整版)复合梯形公式与复合辛普森公式对比
SHANGHAI JIAO TONG UNIVERSITY题目名称：复合梯形公式与复合辛普森公式对比学生姓名:学生学号:班级:学院(系):目录1.概述 (3)2.问题提出 (4)3.算法推导 (5)4.算法框图 (6)4.1复合梯形公式算法流程图 (6)4.2 复合辛普森公式算法流程图 (7)5.MATLAB源程序 (8)6.结论与展望 (9)图表目录图4-1 复合梯形公式算法流程图 (6)图4-2 复合辛普森公式算法流程图 (7)图6-1 MATLAB计算结果 (9)表2-1函数计算结果表 (4)1.概述梯形求积公式和辛普森求积公式分别是牛顿-科斯特公式中n=1和n=2时的情形。
其中梯形求积公式可表示为由于牛顿-科斯特公式在n≥8时不具有稳定性，故不可能通过提高阶的方法来提高求积精度。
为了提高精度通常可把积分区间分成若干子区间（通常是等分），再在每个子区间上用低阶求积公式。
这种方法称为复合求积法。
本文主要讨论复合梯形公式和复合辛普森公式在同一数学问题中的应用。
首先给出了复合梯形公式和复合辛普森公式的推导过程以及其余项的表达形式，然后用流程图的形式介绍算法思路，再运用MATLAB编写代码计算结果，最后对结果进行对比讨论。
希望通过两个算法在同一个算例中的应用对比，更好的理解和掌握复合梯形公式和复合辛普森公式的适用范围和适用条件。
并且能够熟悉MATLAB编程求解问题的流程，掌握编程化的思想方法。
同时对两种方法的计算结果对比分析，讨论两种求积方法的计算精度。
2.问题提出对于函数f(x)=sinxx给出的函数表如下，试用复合梯形公式和复合辛普森公式计算积分I=∫sinxx dx1。
表 2-1函数计算结果表3. 算法推导3.1复合梯形公式根据梯形公式，将区间[a,b]划分为n 等份，分点x k =a +kℎ,h =b−a n,k =0,1,…,n ，在每个子区间[x k ,x k+1](k =0,1,…,n −1)上采用梯形公式，则得：记则T n 为复合梯形公式。
复化梯形公式和复化Simpson公式
一、计算定积分的近似值：221x e xe dx =⎰ 要求：（1）若用复化梯形公式和复化Simpson 公式计算，要求误差限71021-⨯=ε，分别利用他们的余项估计对每种算法做出步长的事前估计；（2）分别利用复化梯形公式和复化Simpson 公式计算定积分；（3）将计算结果与精确解比较，并比较两种算法的计算量。
1.复化梯形公式程序：程序1（求f （x ）的n 阶导数：syms xf=x*exp(x) %定义函数f （x ）n=input('输入所求导数阶数:')f2=diff(f,x,n) %求f(x)的n 阶导数结果1输入n=2f2 =2*exp(x) + x*exp(x)程序2：clcclearsyms x %定义自变量xf=inline('x*exp(x)','x') %定义函数f(x)=x*exp(x)，换函数时只需换该函数表达式即可f2=inline('(2*exp(x) + x*exp(x))','x') %定义f(x)的二阶导数,输入程序1里求出的f2即可。
f3='-(2*exp(x) + x*exp(x))'%因fminbnd（）函数求的是表达式的最小值，且要求表达式带引号，故取负号，以便求最大值e=5*10^(-8) %精度要求值a=1 %积分下限b=2 %积分上限x1=fminbnd(f3,1,2) %求负的二阶导数的最小值点，也就是求二阶导数的最大值点对应的x值for n=2:1000000 %求等分数nRn=-(b-a)/12*((b-a)/n)^2*f2(x1) %计算余项if abs(Rn)<e %用余项进行判断break% 符合要求时结束endendh=(b-a)/n %求hTn1=0for k=1:n-1 %求连加和xk=a+k*hTn1=Tn1+f(xk)endTn=h/2*((f(a)+2*Tn1+f(b)))z=exp(2)R=Tn-z %求已知值与计算值的差fprintf('用复化梯形算法计算的结果 Tn=')disp(Tn)fprintf('等分数 n=')disp(n) %输出等分数fprintf('已知值与计算值的误差 R=')disp(R)输出结果显示：用复化梯形算法计算的结果 Tn=等分数 n=7019已知值与计算值的误差 R=2. Simpson公式程序：程序1：（求f（x）的n阶导数）：syms xf=x*exp(x) %定义函数f（x）n=input('输入所求导数阶数:')f2=diff(f,x,n) %求f(x)的n阶导数结果1输入n=4f2 =4*exp(x) + x*exp(x)程序2：clcclearsyms x%定义自变量xf=inline('x*exp(x)','x') %定义函数f(x)=x*exp(x)，换函数时只需换该函数表达式即可f2=inline('(4*exp(x) + x*exp(x))','x') %定义f(x)的四阶导数，输入程序1里求出的f2即可f3='-(4*exp(x) + x*exp(x))'%因fminbnd（）函数求的是表达式的最小值，且要求表达式带引号，故取负号，一边求最大值e=5*10^(-8) %精度要求值a=1 %积分下限b=2 %积分上限x1=fminbnd(f3,1,2) %求负的四阶导数的最小值点，也就是求四阶导数的最大值点对应的x值for n=2:1000000 %求等分数nRn=-(b-a)/180*((b-a)/(2*n))^4*f2(x1) %计算余项if abs(Rn)<e %用余项进行判断break% 符合要求时结束endendh=(b-a)/n %求hSn1=0Sn2=0for k=0:n-1 %求两组连加和xk=a+k*hxk1=xk+h/2Sn1=Sn1+f(xk1)Sn2=Sn2+f(xk)endSn=h/6*(f(a)+4*Sn1+2*(Sn2-f(a))+f(b)) %因Sn2多加了k=0时的值，故减去f(a)z=exp(2)R=Sn-z %求已知值与计算值的差fprintf('用Simpson公式计算的结果 Sn=')disp(Sn)fprintf('等分数 n=')disp(n)fprintf('已知值与计算值的误差 R=')disp(R)输出结果显示：用Simpson公式计算的结果 Sn=等分数 n=24已知值与计算值的误差 R=用复化梯形公式计算的结果为：，与精确解的误差为：。
梯形法则和辛普森法则
梯形法则与辛普森法则1. 概述梯形法则和辛普森法则是数值积分中常用的近似方法。
它们通过将曲线或曲面分割成若干个由直线或弧线组成的小区间，并在每个小区间内估计函数值来求解定积分。
这两种方法具有较高的精度和广泛的应用领域，本文将详细介绍这两种方法的原理、应用和优缺点。
2. 梯形法则梯形法则是一种利用梯形来近似曲线下面积的方法。
假设我们要计算函数f(x)在区间[a, b]上的定积分，梯形法则将该区间分成n个小区间，每个小区间的宽度为Δx=(b-a)/n。
然后，我们可以将每个小区间内的函数值连接起来，形成若干个梯形，计算每个梯形的面积，并将它们相加，即可得到函数在整个区间上的近似积分值。
梯形法则的近似积分公式如下：∫[a, b] f(x) dx ≈ Δx/2 * [f(a) + 2f(x1) + 2f(x2) + ... + f(b)]其中，x1, x2, …, xn-1 是分割点。
梯形法则的优点是简单易懂，容易推广到多维积分，适用于一般的函数。
然而，它的缺点是存在一定的误差，特别是在曲线弯曲的区域，误差较大。
3. 辛普森法则辛普森法则是一种利用拟合曲线来近似曲线下面积的方法。
与梯形法则类似，辛普森法则也将函数f(x)在区间[a, b]上分成若干个小区间，每个小区间的宽度为Δx=(b-a)/n，但不同的是，辛普森法则采用二次多项式来拟合每个小区间内的曲线。
辛普森法则的近似积分公式如下：∫[a, b] f(x) dx ≈ Δx/3 * [f(x0) + 4f(x1) + 2f(x2) + 4f(x3) + ... + f(xn)]其中，x0, x1, x2, …, xn 是分割点，且n为偶数。
辛普森法则的优点是相比于梯形法则，它的精度更高，尤其适用于曲线较为平滑的情况。
然而，辛普森法则的缺点是计算量较大，对于需要较高精度的计算而言，分割的区间数需要相对较多。
4. 应用场景梯形法则和辛普森法则在数值计算中有广泛的应用，特别是在求解无法用解析方法求得精确解的积分问题时，这两种方法成为重要的工具。
复化梯形公式和复化辛普生公式
1)复化梯形
#include<iostream.h>
#include<math.h>
class tixing
{
private:
int n,k;//n是等分数，
double h,f;//h步长
double a,b,integral,sum;//a上限，b下限，integral积分值
public:
《数值分析》实验报告
姓名
学号
日期
2012．11.20
实验室
设备编号
实验题目
用复化梯形公式和复化辛普生公式求I=∫01sinx/xdx
一实验目的
1．了解复化梯形公式和复化辛普生公式。
2.用复化梯形公式和复化辛普生公式求I=∫01sinx/xdx。
二实验内容
算法：复化梯形公式是Tn=∑h/2[f(xi)+ f(xi+1)]=(b-a)/2n[f(a)+2∑f(xi)+f(b)]记子段[xi,xi+1]的中点为xi+1/2,则复化Simpson公式为Sn=∑h/6[f(xi)+4f(xi+1/2)+ f(xi+1)]=b-a/6n[f(a)+4∑f(xi+1/2)+2f(xi)+f(b)]
double a,b,f,integral,interval;//a上限，b下限，integral积分值，interval每次应用法则的区间长度（2h）
double sum_odd_terms,x,sum_even_terms;//sum_even_terms偶数结点处的函数值之和，sum_odd_terms奇数结点处的函数值之和
复合梯形公式和复合辛普森公式
复合梯形公式和复合辛普森公式1.复合梯形公式步骤1：将积分区间[a,b]均匀地分成n个小区间，每个小区间的长度为h=(b-a)/n，其中n为正整数。
步骤2：定义一个函数f(x)，并在每个小区间上求出函数在小区间两个端点的函数值，记作f(x0), f(x1), f(x2), ..., f(xn)。
步骤3：根据梯形公式，每个小区间上的定积分近似值为(h/2) * (f(x0) + 2f(x1) + 2f(x2) + ... + 2f(xn-1) + f(xn))。
步骤4：将所有小区间上的定积分近似值相加，得到最终的近似值。
复合辛普森公式是通过将积分区间划分成若干个小区间，然后在每个小区间上应用辛普森公式来逼近定积分的值。
具体的步骤如下：步骤1：将积分区间[a,b]均匀地分成n个小区间，每个小区间的长度为h=(b-a)/n，其中n为正偶数。
步骤2：定义一个函数f(x)，并在每个小区间上求出函数在小区间三个点的函数值，记作f(x0),f(x1),f(x2)。
步骤3：根据辛普森公式，每个小区间上的定积分近似值为(h/3)*(f(x0)+4f(x1)+f(x2))。
步骤4：将所有小区间上的定积分近似值相加，得到最终的近似值。
复合辛普森公式的误差主要取决于小区间的数量和函数在每个小区间上的变化情况。
与复合梯形公式相比，复合辛普森公式的精度更高，但计算复杂度也更高。
综上所述，复合梯形公式和复合辛普森公式是数值积分中常用的近似计算方法。
它们通过将积分区间划分成小区间，并在每个小区间上应用相应的公式来逼近定积分的值。
这两种方法都可以通过增加小区间的数量来提高近似的精度，但相应地也会增加计算的复杂度。
根据实际情况，我们可以选择合适的方法来计算需要求解的定积分。
梯形公式和辛普森公式计算积分
梯形公式和辛普森公式是用于计算定积分的数值逼近方法。
在实际应用中，我们经常会遇到需要计算函数在某个区间上的定积分值的情况，而对于一些复杂的函数，直接进行积分计算可能会十分困难，甚至是不可能的。
我们需要借助数值逼近方法来得到积分的近似值。
梯形公式和辛普森公式都是数值积分的基本方法，它们的原理都是通过将被积函数在积分区间上进行分割，然后利用分割后的小区间上的函数值，以及各个小区间的长度来进行计算，从而得到积分的近似值。
梯形公式是一种线性插值法，它的原理是将积分区间等分成n个小区间，然后用每个小区间的两个端点处的函数值进行线性插值，将每个小区间的面积近似为一个梯形，再将所有梯形的面积相加就得到了整个积分的近似值。
具体地，对于被积函数f(x)在区间[a, b]上的积分，我们可以利用梯形公式进行近似计算：1. 将区间[a, b]等分成n个小区间，记每个小区间的长度为h，即h=(b-a)/n。
2. 根据梯形面积的计算公式，我们可以得到每个小区间上梯形的面积为(h/2)*(f(x[i])+f(x[i+1]))，其中x[i]和x[i+1]分别为第i个小区间的两个端点。
3. 将所有小区间上梯形的面积相加得到整个积分的近似值，即I ≈(h/2)*(f(a)+2*f(x[1])+2*f(x[2])+...+2*f(x[n-1])+f(b))。
梯形公式的优点在于其较为简单易懂，且可以很容易地通过计算机进行程序实现。
但是需要注意的是，当被积函数在积分区间上变化较大时，梯形公式可能会产生较大的误差。
与梯形公式类似，辛普森公式也是一种数值积分的方法，它是一种二次插值法，其原理是将积分区间等分成n个小区间，然后利用每个小区间的三个节点处的函数值进行二次插值，将每个小区间的面积近似为一个二次多项式曲线下的面积，再将所有小区间的面积相加就得到了整个积分的近似值。
具体地，对于被积函数f(x)在区间[a, b]上的积分，我们可以利用辛普森公式进行近似计算：1. 将区间[a, b]等分成n个小区间，记每个小区间的长度为h，即h=(b-a)/n。
梯形多步法和辛普森积分
2020/4/23
华南师范大学数学科学学院 谢骊玲
7
闭型牛顿－科特斯面积公式
定理7.1 设xk=x0+kh为等距节点，且fk=f(xk)。前4个闭型N－C 面积公式为
x1
(
f0
f1 )
（梯形公式）
x2 x0
f
(x)dx
h 3
(
f
0
4 f1
f2)
（辛普森公式）
x3 x0
2020/4/23
华南师范大学数学科学学院 谢骊玲
6
基于多项式插值的面积公式
通过M+1个等距点{(xk , f (xk ))}kM0存在唯一的 次数小于等于M的多项式PM(x)。当用该多 项式来近似[a,b]上的f (x)时，PM(x)的积分 就近似等于f (x)的积分，这类公式称为牛 顿－科特斯公式。当使用采样点x0=a和 xM=b时，称为闭型牛顿－科特斯公式
3
1
3.5
1.5
2
-0.5[x0,x3]上y=P3(x)的辛普森3/8积分公式
-0.5 [x0,x4]上y=P4(x)的布尔积分公式
-1 2020/4/23
华南师范大学数学科学-1学院 谢骊玲
10
N－C公式的精度
推论7.1 设f(x)充分可微，则N－C面积公式的E[f]包含一个高 阶的导数项。
梯形公式的精度为n=1，如果f∈C2[a,b]，则
意i次多项式Pi(x)=aixi+ai-1xi-1+…+a1x+a0，如果i≤n，则对所有x，
有Pi(n+1)(x)≡0，并且对所有的x，式
成立 P(n1) n1
(x)
(n
复化梯形公式和复化辛普森公式
复化梯形公式和复化辛普森公式1. 引言嘿，大家好！今天我们来聊聊数学里那些看似高深莫测的公式，尤其是复化梯形公式和复化辛普森公式。
这些名字听起来就像是从某部科幻片里蹦出来的角色，但其实它们是我们在数值积分中不可或缺的好帮手。
你知道吗？它们就像是数学世界里的“超能英雄”，让我们轻松搞定积分，简直是妙不可言。
2. 复化梯形公式2.1 你知道什么是梯形吗？首先，咱们得聊聊复化梯形公式。
说白了，就是把一个复杂的积分任务，分解成几个小的梯形来求解。
想象一下，你在河边钓鱼，河水流得可欢了。
为了找一个合适的钓鱼点，你可能得把河分成几段，然后每一段的宽度就是你的小梯形。
你看，这就是复化梯形的魅力所在！2.2 如何运用复化梯形公式？用这个公式的时候，你只需把整个区间分成N个小区间，每个区间的宽度都是一样的。
然后，把每个小区间的函数值拿来加一加，再乘上宽度的一半，最后再把头尾的函数值加上。
这听起来是不是很简单？比如，你想算从0到1的某个函数的积分，只要把这个区间分成若干段，像切蛋糕一样，每一片都求个函数值，然后把结果合起来就行了。
简单得就像吃个冰淇淋，大家都喜欢。
3. 复化辛普森公式3.1 辛普森是谁？接下来，让我们来看看复化辛普森公式。
辛普森这个名字，大家可能都听过，或者说过“这是辛普森家的事儿”。
其实，他是一位牛逼的数学家，专门研究如何让积分变得更加简单。
辛普森公式就像是对梯形公式的一次升级，像换了个新款手机，功能更强大，效果更好。
3.2 如何运用复化辛普森公式？用复化辛普森公式的时候，我们也是把整个区间分成N个小区间，不过这里的N必须是偶数哦！每个小区间的宽度仍然是一样的。
然后，用函数值的加权平均法来计算。
换句话说，你把每个小区间的头尾和中间的函数值都考虑进来，像是为你的冰淇淋加上各种口味的配料。
最后，你的结果就会比单纯用梯形公式得来的要精准多了，仿佛一口下去，味蕾都在舞蹈。
4. 比较与应用4.1 谁更强？说到这儿，很多人就会问，复化梯形公式和复化辛普森公式，谁更厉害呢？其实，这就像问“苹果和橘子，哪个更好吃”。
梯形公式和辛普森公式例题
梯形公式和辛普森公式例题1.概述本文将介绍梯形公式和辛普森公式两个数值积分方法，并通过例题演示其具体应用。
2.梯形公式梯形公式是一种数值积分方法，用来估计定积分的近似值。
它的思想是将曲线下的面积近似为一系列梯形的面积之和。
2.1公式推导考虑函数f(x)在区间[a,b]上的定积分，可以将[a,b]等分成n个小区间，每个小区间的宽度为h=(b-a)/n。
梯形公式的计算公式如下：∫[a,b]f(x)dx≈h/2*[f(a)+2f(a+h)+2f(a+2h)+...+f(b)]2.2例题假设要计算函数f(x)=x^2在区间[0,2]上的定积分，可以使用梯形公式进行近似计算。
首先将区间[0,2]等分成n个小区间，选择n=4。
则每个小区间的宽度为h=(2-0)/4=0.5。
根据梯形公式的计算公式，可以得到近似值为：∫[0,2]x^2d x≈0.5/2*[0^2+2(0.5)^2+2(1)^2+2(1.5)^2+2(2)^ 2]≈2.5因此，函数f(x)=x^2在区间[0,2]上的定积分的近似值为2.5。
3.辛普森公式辛普森公式是一种数值积分方法，用来估计定积分的近似值。
它的思想是将曲线下的面积近似为一系列抛物线的面积之和。
3.1公式推导考虑函数f(x)在区间[a,b]上的定积分，同样将[a,b]等分成n个小区间，每个小区间的宽度为h=(b-a)/n。
辛普森公式的计算公式如下：∫[a,b]f(x)dx≈h/3*[f(a)+4f(a+h)+2f(a+2h)+4f(a+3h)+...+ f(b)]3.2例题继续以函数f(x)=x^2在区间[0,2]上的定积分为例，使用辛普森公式进行近似计算。
同样选择n=4，计算每个小区间的宽度为h=(2-0)/4=0.5。
根据辛普森公式的计算公式，可以得到近似值为：∫[a,b]x^2d x≈0.5/3*[0^2+4(0.5)^2+2(1)^2+4(1.5)^2+2(2)^ 2]≈2.6667因此，函数f(x)=x^2在区间[0,2]上的定积分的近似值为2.6667。
复化梯形公式和复化Simpson公式
数值计算方法上机题目3一、计算定积分的近似值：221x e xe dx =⎰ 要求：（1）若用复化梯形公式和复化Simpson 公式计算，要求误差限71021-⨯=ε，分别利用他们的余项估计对每种算法做出步长的事前估计；（2）分别利用复化梯形公式和复化Simpson 公式计算定积分；（3）将计算结果与精确解比较，并比较两种算法的计算量。
1.复化梯形公式程序：程序1（求f （x ）的n 阶导数：syms xf=x*exp(x) %定义函数f （x ）n=input('输入所求导数阶数:')f2=diff(f,x,n) %求f(x)的n 阶导数结果1输入n=2f2 =2*exp(x) + x*exp(x)程序2：clcclearsyms x%定义自变量xf=inline('x*exp(x)','x') %定义函数f(x)=x*exp(x)，换函数时只需换该函数表达式即可f2=inline('(2*exp(x) + x*exp(x))','x') %定义f(x)的二阶导数,输入程序1里求出的f2即可。
f3='-(2*exp(x) + x*exp(x))'%因fminbnd（）函数求的是表达式的最小值，且要求表达式带引号，故取负号，以便求最大值e=5*10^(-8) %精度要求值a=1 %积分下限b=2 %积分上限x1=fminbnd(f3,1,2) %求负的二阶导数的最小值点，也就是求二阶导数的最大值点对应的x值for n=2:1000000 %求等分数nRn=-(b-a)/12*((b-a)/n)^2*f2(x1) %计算余项if abs(Rn)<e %用余项进行判断break% 符合要求时结束endendh=(b-a)/n %求hTn1=0for k=1:n-1 %求连加和xk=a+k*hTn1=Tn1+f(xk)endTn=h/2*((f(a)+2*Tn1+f(b)))z=exp(2)R=Tn-z %求已知值与计算值的差fprintf('用复化梯形算法计算的结果 Tn=')disp(Tn)fprintf('等分数 n=')disp(n) %输出等分数fprintf('已知值与计算值的误差 R=')disp(R)输出结果显示：用复化梯形算法计算的结果 Tn= 7.3891等分数 n=7019已知值与计算值的误差 R= 2.8300e-0082. Simpson公式程序：程序1：（求f（x）的n阶导数）：syms xf=x*exp(x) %定义函数f（x）n=input('输入所求导数阶数:')f2=diff(f,x,n) %求f(x)的n阶导数结果1输入n=4f2 =4*exp(x) + x*exp(x)程序2：clcclearsyms x%定义自变量xf=inline('x*exp(x)','x') %定义函数f(x)=x*exp(x)，换函数时只需换该函数表达式即可f2=inline('(4*exp(x) + x*exp(x))','x') %定义f(x)的四阶导数，输入程序1里求出的f2即可f3='-(4*exp(x) + x*exp(x))'%因fminbnd（）函数求的是表达式的最小值，且要求表达式带引号，故取负号，一边求最大值e=5*10^(-8) %精度要求值a=1 %积分下限b=2 %积分上限x1=fminbnd(f3,1,2) %求负的四阶导数的最小值点，也就是求四阶导数的最大值点对应的x值for n=2:1000000 %求等分数nRn=-(b-a)/180*((b-a)/(2*n))^4*f2(x1) %计算余项if abs(Rn)<e %用余项进行判断break% 符合要求时结束endendh=(b-a)/n %求hSn1=0Sn2=0for k=0:n-1 %求两组连加和xk=a+k*hxk1=xk+h/2Sn1=Sn1+f(xk1)Sn2=Sn2+f(xk)endSn=h/6*(f(a)+4*Sn1+2*(Sn2-f(a))+f(b)) %因Sn2多加了k=0时的值，故减去f(a)z=exp(2)R=Sn-z %求已知值与计算值的差fprintf('用Simpson公式计算的结果 Sn=')disp(Sn)fprintf('等分数 n=')disp(n)fprintf('已知值与计算值的误差 R=')disp(R)输出结果显示：用Simpson公式计算的结果 Sn= 7.3891等分数 n=24已知值与计算值的误差 R= 2.7284e-008用复化梯形公式计算的结果为：7.3891，与精确解的误差为：2.8300e-008。
几种常用数值积分方法的比较汇总
几种常用数值积分方法的比较汇总数值积分是一种用计算机逼近求解定积分的方法，它通过将区间划分为多个小区间，并在每个小区间上进行数值计算，最后将结果相加以得到整个区间上的定积分近似值。
在实际应用中，常用的数值积分方法有梯形法则、辛普森法则和复化求积法。
下面将详细介绍这几种方法，并对它们进行比较汇总。
1.梯形法则是一种基本的数值积分方法。
它的原理是将每个小区间视为一条梯形，并用该梯形的面积来近似表示该小区间的积分值。
具体而言，对于求解区间[a,b]上的定积分，梯形法则的计算公式为：∫[a,b]f(x)dx≈ (b-a)[f(a) + f(b)]/2梯形法则的优点是简单易懂、计算速度较快，但它的缺点是精度较低，特别是当被积函数曲线较为陡峭时。
2.辛普森法则是一种比梯形法则更精确的数值积分方法。
它的原理是将每个小区间视为一个二次曲线，并用该曲线下的面积来近似表示该小区间的积分值。
具体而言，对于求解区间[a,b]上的定积分，辛普森法则的计算公式为：∫[a,b]f(x)dx ≈ (b-a)[f(a) + 4f((a+b)/2) + f(b)]/6辛普森法则的优点是精度较高，特别是对于曲线比较平滑的函数，它能给出较为准确的积分近似值。
然而，辛普森法则的计算量较大，因为它需要在每个小区间上计算3个点的函数值。
3.复化求积法是一种综合性的数值积分方法，它基于划分区间的思想，将整个求积区间划分为多个小区间，并在每个小区间上采用其中一种数值积分方法来进行计算。
具体而言，复化求积法可以采用梯形法则或辛普森法则来进行计算。
它的计算公式如下：∫[a,b]f(x)dx ≈ ∑[i=0,n-1] (b-a)/n * [f(a + i(b-a)/n) +f(a + (i+1)(b-a)/n)]/2复化求积法的优点是能够灵活地根据被积函数的特点选择合适的数值积分方法，从而提高求积的准确性。
但它的计算量较大，尤其在需要高精度的情况下，需要划分较多的小区间。
梯形公式和辛普森公式计算积分
梯形公式和辛普森公式计算积分梯形公式和辛普森公式都是数值积分的常用方法。
数值积分是一种近似计算函数积分值的方法，适用于无法求得解析解的复杂函数。
梯形公式是数值积分中最简单的一种方法。
它的原理是将要积分的函数曲线划分为一系列梯形，然后计算每个梯形的面积并相加得到积分值。
具体步骤如下：1. 将积分区间[a, b]等分为n个小区间，每个小区间的宽度为h = (b - a) / n。
2. 计算函数在每个小区间的端点处的函数值，得到f(a), f(a+h), f(a+2h), ..., f(b)。
3. 使用梯形公式计算每个小区间的面积：Area = (f(i) + f(i+1)) * h / 2，其中i为当前小区间的序号。
4. 将每个小区间的面积相加得到最终的积分值：Integral = Area1 + Area2 + ... + Arean。
梯形公式的优点是简单易懂、易于实现，但是在求解复杂函数时可能会出现较大误差。
为了提高积分精度，可以增加划分的小区间数n。
辛普森公式是一种更精确的数值积分方法，它的原理是利用二次多项式对函数曲线进行近似。
具体步骤如下：1. 将积分区间[a, b]等分为2n个小区间，每个小区间的宽度为h = (b - a) / (2n)。
2. 计算函数在每个小区间的端点处的函数值，得到f(a), f(a+h), f(a+2h), ..., f(b)。
3. 使用辛普森公式计算每个小区间的面积：Area = (h/3) * (f(i) + 4f(i+1) + f(i+2))，其中i为当前小区间的序号。
4. 将每个小区间的面积相加得到最终的积分值：Integral = Area1 + Area2 + ... + Area2n。
辛普森公式相比梯形公式在积分计算中具有更高的精度。
它通过利用二次多项式对函数曲线进行逼近，可以更准确地估计出积分值。
但是，辛普森公式的缺点是需要计算更多的函数值，导致计算量增加。
利用复化梯形公式、复化simpson 公式计算积分
教师签名：
年月日
上图是利用复化simpson公式所画出的误差。其中：红线是计算误差，‘＊’号是实际误差。
注：纵轴是0.0001。
0.5208、0.1029、0.0326、0.0133、0.0064、0.0035、0.0020、0.0013、0.0008是n值分别为2到10的实际误差，-0.0083是计算误差。
成
绩
评
注：由于被积函数是x.^4，它的四阶倒数是24，所以它的积分余项是：(-(b-a)/180*((b-a)/2).^4*24)
实
验
结
果
分
析
及
心
得
体
会
上图是利用复化梯形公式所画出的误差。其中：红线是计算误差，‘＊’号是实际误差。-0.0017是计算误差。0.0417、0.0185、0.0104、0.0067 0.0046、0.0034、0.0026、0.0021、0.0017是n值分别为2到10的实际误差。
x=linspace(0,1,9);
plot(x,A,'r',x,R,'*')％将计算误差与实际误差用图像画出来
注：由于被积函数是x.^2，它的二阶倒数为2，所以积分余项为：(-(b-a)/12*h.^ 2*2)
实
验
原
理
（
算
法
流
程
图
或
者
含
注
释
的
源
代
码
）
利用复化simpson公式的程序代码如下：
同样首先建立被积函数的函数文件：
f=x.^2;％首先建立被积函数，以便于计算真实值。
a=0;％积分下线
b=1;％积分上线
复合梯形积分和复合Simpson积分计算数值积分
实验五一、实验名称复合梯形积分和复合Simpson 积分计算数值积分二、实验目的与要求：实验目的: 掌握复合梯形积分和复合Simpson 积分算法。
实验要求：1.给出复合梯形积分和复合Simpson 积分算法思路，2.用C 语言实现算法，运行环境为Microsoft VisualC++。
三、算法思路：我们把整个积分区间[a,b]分成n 个子区间[xi,xi+1]，i=0，1,2,…,n，其中x0=a ，xn+1=b 。
这样求定积分问题就分解为求和问题：⎰∑⎰=-==b a n i x x i i dx x f dx x f S 11)()(当这n+1个结点为等距结点时，即n a b h ih a x i /)(-=+=，其中，i=0,1,2,…,n ，复化梯形公式的形式是∑=-+=ni i i n x f x f h S 11)]()([2 算法：input n0.0←Sfor i=1 to n do ))()((21i i x f x f h S S ++←- end dooutput S如果n 还是一个偶数，则复合Simpson 积分的形式是∑=--++=2/121222)]()(4)([3n i i i i n x f x f x f h S 算法：input n0.0←Sfor i=1 to n/2 do ))()(4)((321222i i i x f x f x f h S S +++←-- end dooutput S四、实验题目：五、问题的解：编写程序（程序见后面附录），输出结果如下：为了便于看清数值积分结果与原函数积分实际结果的差异。
我在运行程序时故意计算了一下原函数积分的实际结果。
分析并比较得到的数据可以看出，当k 越来越大时，数值积分的结果越来越靠近原函数积分实际结果，并且复合Simpson 积分的结果更迅速地靠近原函数积分实际结果，这是有原因的，从两种方法的误差项即可看出。
复合梯形公式与复合辛普森公式对比
SHANGHAI JIAO TONG UNIVERSITY题目名称：复合梯形公式与复合辛普森公式对比学生姓名:学生学号:班级:学院(系):目录1.概述 (3)2.问题提出 (4)3.算法推导 (5)4.算法框图 (6)4.1复合梯形公式算法流程图 (6)4.2 复合辛普森公式算法流程图 (7)5.MATLAB源程序 (8)6.结论与展望 (9)图表目录图4-1 复合梯形公式算法流程图 (6)图4-2 复合辛普森公式算法流程图 (7)图6-1 MATLAB计算结果 (9)表2-1函数计算结果表 (4)1.概述梯形求积公式和辛普森求积公式分别是牛顿-科斯特公式中n=1和n=2时的情形。
其中梯形求积公式可表示为由于牛顿-科斯特公式在n≥8时不具有稳定性，故不可能通过提高阶的方法来提高求积精度。
为了提高精度通常可把积分区间分成若干子区间（通常是等分），再在每个子区间上用低阶求积公式。
这种方法称为复合求积法。
本文主要讨论复合梯形公式和复合辛普森公式在同一数学问题中的应用。
首先给出了复合梯形公式和复合辛普森公式的推导过程以及其余项的表达形式，然后用流程图的形式介绍算法思路，再运用MATLAB编写代码计算结果，最后对结果进行对比讨论。
希望通过两个算法在同一个算例中的应用对比，更好的理解和掌握复合梯形公式和复合辛普森公式的适用范围和适用条件。
并且能够熟悉MATLAB编程求解问题的流程，掌握编程化的思想方法。
同时对两种方法的计算结果对比分析，讨论两种求积方法的计算精度。
2.问题提出对于函数f(x)=sinxx给出的函数表如下，试用复合梯形公式和复合辛普森公式计算积分I=∫sinxx dx1。
表 2-1函数计算结果表3. 算法推导3.1复合梯形公式根据梯形公式，将区间[a,b]划分为n 等份，分点x k =a +kℎ,h =b−a n,k =0,1,…,n ，在每个子区间[x k ,x k+1](k =0,1,…,n −1)上采用梯形公式，则得：记则T n 为复合梯形公式。
梯形多步法和辛普森积分
2020/4/23
华南师范大学数学科学学院 谢骊玲
12
例7.2 对于梯形公式，h=1
1 f (x)dx 1 ( f (0) f (1)) 1 (1.00000 0.72159) 0.86079
0
2
2
对于辛普森公式，h=1/2
1 f (x)dx 1/ 2 ( f (0) 4 f (1) f (1)) 1 (1.00000 4(1.55152) 0.72159) 1.32128
2
几个简单的数值积分公式
左/中/右矩形公式
b
I a f (x)dx f (a)(b a)
I b f (x)dx f ( a b)(b a)
a
2
b
I a f (x)dx f (b)(b a)
梯形公式
I b f (x)dx b a [ f (a) f (b)]
a
2
2020/4/23
左矩形公式梯形公式中矩形公式右矩形公式谢丽玲华南师范大学数学科学学院2020年8月5日积分导论数值积分的目的是通过计算有限个采样点上fx的值来逼近区间a和b中fx的定积分71定义为ax0
第7章 数值积分
2020/4/23
华南师范大学数学科学学院 谢骊玲
1
数值积分问题
数值积分是工程师和科学家经常使用的基本工具， 用来计算无法解析求解的定积分的近似解
组合辛普森公式
1 f (x)dx 1/ 4 ( f (0) 4 f (1) 2 f (1) 4 f ( 3) f (1))
0
3
4
2
4
1 (1.00000 4(1.65534) 2(1.55152) 4(1.06666) 0.72159) 1.30938 12
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2.0 0
0
1.5
0
f ( x)dx
2(0.5) (7(1.00000) 32(1.55152) 12(0.72159) 32(0.93765) 7(1.13390)) 2.29444 45
华南师范大学数学科学学院 谢骊玲
2013-7-29
2
2
1.5
1.5
1
1
0.5 3.5
华南师范大学数学科学学院 谢骊玲
例7.2 对于梯形公式，h=1
1
0
1 1 f ( x)dx ( f (0) f (1)) (1.00000 0.72159) 0.86079 2 2
对于辛普森公式，h=1/2
1
0
f ( x)dx
1/ 2 1 1 ( f (0) 4 f ( ) f (1)) (1.00000 4(1.55152) 0.72159) 1.32128 3 2 6
0.5
0
1.0
0.5 f ( x)dx (1.00000 1.55152) 0.63788 2
0.5 f ( x)dx (1.00000 4(1.55152) 0.72159) 1.32128 3
f ( x)dx 3(0.5) (1.00000 3(1.55152) 3(0.72159) 0.93765) 1.64193 8
组合梯形公式
1
0
f ( x)dx
1/ 4 1 1 3 ( f (0) 2 f ( ) 2 f ( ) 2 f ( ) f (1)) 2 4 2 4 1 (1.00000 2(1.65534) 2(1.55152) 2(1.06666) 0.72159) 1.28358 8
积分简介
数值积分的目的是，通过在有限个采样点上计算 f (x)的值来逼近 f (x)在区间[a,b]上的定积分
定义7.1 设a=x0<x1<…<xM=b. 称形如
Q[ f ] wk f ( xk ) w0 f ( x0 ) w1 f ( x1 ) wM f ( xM )
I f ( x)dx f (b)(b a )
a
梯形公式
I
b a
ba f ( x)dx [ f (a) f (b)] 2
华南师范大学数学科学学院 谢骊玲
2013-7-29
2
2
1.5
1.5
1
1
0.5 3.5
0.5 3.5
3
0.5
1
1.5
-0.5
2
3
2.5
0.5
3
1
3.5
故截断误差的一般形式为E[ f ]=K f (n+1)(c)，其中K是一个 合理选择的常数，n为精度
2013-7-29
注意：积分公式的数值精度定义没有指定积分区间 华南师范大学数学科学学院 谢骊玲
基于多项式插值的面积公式
M 通过M+1个等距点{( xk , f ( xk ))}k 0 存在唯一的 次数小于等于M的多项式PM(x)。当用该多 项式来近似[a,b]上的f (x)时，PM(x)的积分 就近似等于f (x)的积分，这类公式称为牛 顿－科特斯公式。当使用采样点x0=a和 xM=b时，称为闭型牛顿－科特斯公式
x3
x0
3h 3h5 (4) f ( x)dx ( f 0 3 f1 3 f 2 f3 ) f ( ) 8 80
布尔公式的精度为n=5，如果f∈C6[a,b]，则
2013-7-29
x4
x0
2h 8h7 (6) f ( x)dx (7 f 0华南师范大学数学科学学院f3 谢骊玲4 ) 32 f1 12 f 2 32 7 f f ( ) 45 945
组合辛普森公式
1
0
f ( x)dx
1/ 4 1 1 3 ( f (0) 4 f ( ) 2 f ( ) 4 f ( ) f (1)) 3 4 2 4 1 (1.00000 4(1.65534) 2(1.55152) 4(1.06666) 0.72159) 1.30938 12
x4
x0
利用N－C公式求数值积分
例7.1 函数f(x)=1+e-xsin(4x)，等距面积节点为x0=0.0, x1=0.5,x2=1.0,x3=1.5,x4=2.0,对应的函数值为f0=1.00000, f1=1.55152, f2=0.72159, f3=0.93765, f4=1.13390,h=0.5
2013-7-29
华南师范大学数学科学学院 谢骊玲
组合公式
例7.3 在区间[0,1]上取相同的步长h=1/4，进行5次函数求值
组合梯形公式
x4
x0
f ( x)dx f ( x)dx f ( x)dx f ( x) dx f ( x) dx
x0 x1 x2 x3
x1
0.5 3.5
0.5 3
1
1.5
-0.5
2 3
2.5
0.5
3
1
3.5
1.5
2
-0.5 2.5
[x0,x1]上y=P1(x)的梯形积分公式
-0.5 2.5
[x0,x2]上y=P2(x)的辛普森积分公式
-1 2
-1 2
-1.5 1.5
-1.5 1.5
-2 1
-2 1
-2.5 0.5
-2.5 0.5
0.5
布尔公式的结果 该定积分的真解
2013-7-29
(5)
5 0
t3 dt 4.8998922 t e 1
本章的目的是推导数值积分的基本原理
华南师范大学数学科学学院 谢骊玲
2013-7-29
几个简单的数值积分公式
左/中/右矩形公式
I f ( x)dx f ( a)(b a )
a b
I
b
a b
ab f ( x)dx f ( )(b a ) 2
x2
x3
x4
h h h h ( f 0 f1 ) ( f1 f 2 ) ( f 2 f3 ) ( f3 f 4 ) 2 2 2 2 h ( f 0 2 f1 2 f 2 2 f3 f 4 ) 2
组合辛普森公式 x
x4
0
f ( x)dx f ( x)dx f ( x )dx
华南师范大学数学科学学院 谢骊玲
2013-7-29
5
5
4
例7.2
该定积分的真解为
4
3
1
0
21e 4cos(4) sin(4) f ( x)dx 1.3082506046426 17e
3 2
2
1 5
1 5
0.5 4 -1 3 -2 2 -3 1 -4
1
1.5 4
2
0.5
2.5
1
[0,1]上y=P1(x)的梯形积分公式
对于辛普森3/8公式，h=1/3
1
0
f ( x)dx
3(1/ 3) 1 2 ( f (0) 3 f ( ) 3 f ( ) f (1)) 8 3 3 1 (1.00000 3(1.69642) 3(1.23447) 0.72159) 1.31440 8
对于布尔公式，h=1/4
2013-7-29
华南师范大学数学科学学院 谢骊玲
闭型牛顿－科特斯面积公式
定理7.1 设xk=x0+kh为等距节点，且fk=f(xk)。前4个闭型N－C 面积公式为 x1 h （梯形公式） x0 f ( x)dx 2 ( f0 f1 )
x2
x0
h f ( x)dx ( f 0 4 f1 f 2 ) 3
1
1.5
-0.5
2
2.5
0.5
3
1
3.5
1.5
2
-0.5
[x0,x3]上y=P3(x)的辛普森3/8积分公式
2013-7-29
-0.5
[x0,x4]上y=P4(x)的布尔积分公式
-1
-1 华南师范大学数学科学学院 谢骊玲
N－C公式的精度
推论7.1 设f(x)充分可微，则N－C面积公式的E[f]包含一个高 阶的导数项。 梯形公式的精度为n=1，如果f∈C2[a,b]，则
（辛普森公式）
2013-7-29
x3
x0
3h f ( x)dx ( f0 3 f1 3 f 2 f3 ) （辛普森3/8公式） 8
2h f ( x)dx (7 f 0 32 f1 12 f 2 32 f3 7 f 4 ) （布尔公式） 45
华南师范大学数学科学学院 谢骊玲
步长的选择
因为各个公式所需节点个数不同，如果固 定求积区间[a,b]的端点，则对不同公式要 采用不同的步长。梯形公式、辛普森公式、 辛普森3/8公式和布尔公式的步长分别为 h=b-a,h=(b-a)/2,h=(b-a)/3和h=(b-a)/4
例7.2 分别将区间[0,1]作1、2、3、4等分
2013-7-29
第7章 数值积分
2013-7-29
华南师范大学数学科学学院 谢骊玲
数值积分问题
数值积分是工程师和科学家经常使用的基本工具， 用来计算无法解析求解的定积分的近似解
如： ( x)