【数值分析】复化积分公式

对于积分：

只要找到被积公式的原函数F(x)，利用牛顿莱普利兹公式有：

但是，实际使用这种求积分的方法往往是有困难的，因为大量的被积函数的原函数是不能用初等函数表示的；另外，当f(x)是由测量或数值计算给出的一张数据表时，牛顿莱普利兹公式也无法直接运用，因此有必要研究积分的数值计算问题。

对于一些理论的推导，大家可以看看维基百科，下面我主要给出牛顿-科特斯公式在n=1(梯形求积公式)、n=2（辛普森公式）的情况，并通过代码实现。

梯形公式：

辛普森公式：

应用高阶牛顿-科特斯公式计算积分时，会出现数值不稳定的情况，而低阶公式往往因为积分步长过大使得离散误差变大，因此，为了提高求积公式的精度，可以把积分区间分成若干个子区间，在每个子区间上使用低阶求积公式，然后将结果加起来，这种方法称为复化求积法。

复化梯形公式

将区间[a,b]划分为n等分，步长为h=(b-a)/h，节点为

，在每个子区间

使用梯形公式得：

复化辛普森公式

根据复化梯形公式的推导，同理可得复化辛普森公式为：

下面我们通过实例来实现复化梯形公式和复化辛普森公式：

对于函数f(x)=sin(x)/x，试用复化梯形公式和复化辛普森公式计算函数f(x)在[0，1]上的积分。

具体的程序实现如下：

#include<stdio.h>
#include<math.h>
double Function(double x)//所要计算积分的函数f(x)
{if(x==0)//sin(x)/x在0处的取值为1return 1;elsereturn sin(x)/x;
}
//复化梯形公式
double Trapz(double a,double b,int n)
{double h=(b-a)/n;double T=0;for(int i=1;i<n;i++){T=T+Function(a+i*h);}T*=2;T=(Function(a)+Function(b)+T)*h/2;return T;
}
//复化辛普森公式
double MulripleSimpson(double a,double b,int n)
{double h=(b-a)/n;double T=0;for(int i=0;i<n;i++){T=T+Function(a+i*h)+4*Function(a+(i+0.5)*h)+Function(a+(i+1)*h);}T=T*h/6;return T;
}
void main()
{printf("使用复化梯形公式可得：%f\n",Trapz(0,1,8));printf("使用复化辛普森公式可得：%f\n",MulripleSimpson(0,1,4));
}

运行结果如下图：

结果分析：

比较复化梯形公式和复化辛普森公式两种方法的运行结果，我们发现复化辛普森公式与准确值0.9460831更加接近，复化梯形公式只有2位有效数字，而复化辛普森公式有6为有效数字。

原文：http://blog.csdn.net/tengweitw/article/details/43311685

作者：nineheadedbird

【数值分析】复化积分公式相关推荐

c语言simpson积分计算方法,数值分析复化Simpson积分公式和复化梯形积分公式计算积分的通用程序...
数值分析复化Simpson积分公式和复化梯形积分公式计算积分的通用程序数值分析第五次程序作业 PB09001057 孙琪 [问题] 分别编写用复化Simpson积分公式和复化梯形积分公式计算积分的通 ...
数值分析复化梯形公式matlab,数值分析复化梯形公式,复化Simpson公式MATLAB程序
<数值分析复化梯形公式,复化Simpson公式MATLAB程序>由会员分享,可在线阅读,更多相关<数值分析复化梯形公式,复化Simpson公式MATLAB程序(1页珍藏版)>请 ...
数值分析复化求积matlab,MATLAB数值分析实验二(复合梯形、辛普森和龙贝格求积,以及二重积分计算等)...
1.理解如何在计算机上使用数值方法计算定积近似值; 2.学会复合梯形.复合Simpson和龙贝格求积分公式的编程与应用. 3.探索二重积分在矩形区域的数值积分方法. 佛山科学技术学院实验报告 ...
数值分析复化求积matlab,数值分析：复化公式求解定积分
复化公式求解定积分用熟悉的计算机语言编程上机完成 (1)用复化梯形公式.复化Simpson公式和复化Cotes公式计算积分,自己设置不同精度要求,对结果进行比较分析. (2)用Romberg积分法计 ...
[计算机数值分析]复化梯形的递推公式的变步长算法求积分
如下图变步长的算法框图,其中T1和T2分别代表二分前后的积分值. 问题描述: 运行示例: 源码: #include<iostream> #include<cmath> doub ...
数值分析复化求积matlab,数值分析实验指导-7积分
内江师范学院数值分析实验报告册编制张莉审定牟廉明专业: 班级:级学号: 姓名: 数学与信息科学学院 2013年9月说明一.学生在做实验之前必须要准备实验,主要包括预 ...
【matlab】数值积分公式的程序实现
(一)专题实验(Newton-Cotes积分公式) 1.编写[a,b]上梯形积分公式.Simpson积分公式. 2.利用自己编写的程序计算定积分,计算一下数值解和精确解之间差的绝对值. 梯形积分: f ...
数值积分-龙贝格(Romberg)积分
数值积分在工程上是个比较有用的数学工具.在工程上有很多数学问题,看似简单,计算所用的数学公式不算复杂,但是求解起来却很困难,很难获得解析解的公式,这个时候就需要用到数值求解的办法来获取满足工程需要的近 ...
MATlAB运用——数值积分
实验一分别利用变步长复化梯形公式.变步长复化Simpson公式和复化Guass-Legendre I型公式计算下列式子,要求绝对误差限为0.5∗10−70.5*10^{-7}0.5∗10−7,并比较 ...

【数值分析】复化积分公式

【数值分析】复化积分公式相关推荐

最新文章

热门文章