首页 > 新闻中心

亚博app|数值积分方式的总结(从简单矩形积分到龙贝格积分、自适应积分、高斯积分等)

2021-02-21 20:18:13  来源:网络整理 作者:佚名 点击量: 89

各种数值积分方式总结(龙贝格积分、自适应积分、高斯积分等)

本文整理了各类类型的数值积分,从简单的矩形积分、辛普森积分到高精度的自适应积分、龙贝格积分和高斯积分;对于多重积分(高维积分),本文先简单介绍一下蒙特卡洛方式(适用于积分维数小于4的积分),后续会发贴专门针对多重积分方式做介绍。

(如有疏忽,欢迎见谅,谢谢~)

1 (一重积分)常用的数值积分方式

常用积分方式包括矩形积分、辛普森积分、自适应积分、龙贝格积分和高斯积分等;其中自适应积分、龙贝格积分和高斯积分属于高精度积分。

1.1 牛顿-科茨(Newton-Cotes)积分公式

牛顿-科茨(Newton-Cotes)积分公式:

在这里插入图片描述

包含了梯形法则、辛普森法则、布尔法则等。

1.1.1 梯形法则(2点积分)

梯形法则积分估算公式:

其中,

在这里插入图片描述

在这里插入图片描述

于是得:

在这里插入图片描述

梯形公式要估算两个点的函数值f(a)和f(b)。所以说是两点积分

在这里插入图片描述

梯形法则的偏差为:

在这里插入图片描述

1.1.2-3 辛普森法则(3点积分和4点积分)

辛普森法则分为辛普森1/3法则和辛普森3/8法则,就是用二次/三次配准曲线来取代直线进行积分。

等积法公式_积化和差与和差化积公式_三角函数积化和差和差化积公式

在这里插入图片描述

1.1.2 辛普森1/3法则(3点积分)

辛普森1/3法则法则积分估算公式:

在这里插入图片描述

1.1.3 辛普森3/8法则(4点积分)

辛普森3/8法则法则积分估算公式:

在这里插入图片描述

1.1.4 布尔法则(5点积分)

布尔法则与辛普森法则类似,只是更高阶的方程积分等积法公式,其公式为:

在这里插入图片描述

1.1.5牛顿-科茨积分公式总结

在这里插入图片描述

注意:1.对于3点积分与4点积分具有相同的偏差阶;5点积分与6点积分具有相同的偏差阶。对于点数更多的积分也一样。因此优先使用质数个点的积分。

2.对于点数小于等于9的牛顿-科茨积分,求积公式的稳定性得不到保证(具体缘由可以找数值积分的数瞧瞧,本文源自《现代数值积分》第5章 数值积分与数值微分),积分不能保证收敛,因此实际估算中通常不采用高阶牛顿-科茨积分公式。(处于效率考虑,一般也极少用超过5个点的牛顿-科茨积分公式)。

(穿插)积分改进思路1

牛顿-科茨积分公式思想:根据积分点构造近似的配准方程,进而借助方程积分表示原积分。

对于高次的方程配准,会出现龙格现象(对于高次的方程配准,插值方程会出现不收敛的现象,成为龙格现象,想深入了解可以去查询相关资料。并且次数越高,越容易出现不收敛现象)。

针对龙格现象,更好的选择切比雪夫节点来进行插补(想深入了解可以去查询相关资料,本文部份源自《现代数值积分》第3章 多项式配准与样条配准)。

在这里插入图片描述

由于切比雪夫节点的特殊性,对于高阶配准的可以收到较好的疗效,因此,对于“牛顿-科茨积分公式总结”中提及的点数小于等于9的牛顿-科茨积分,可以采用切比雪夫估算节点,然后构造配准方程,再估算积分(此方式实现和推广比较麻烦,此处不做扩充,感兴趣的可以自己下去推论实现)。

1.2 复合积分(复合矩形积分、复合辛普森积分等)

三角函数积化和差和差化积公式_积化和差与和差化积公式_等积法公式

上述主要介绍的是单个区间的积分。为了估算更准确,我们可以将区间界定为若干个小的区间,然后对每位新村间进行积分,然后对各个新村间的积分求和。由于对于单个新村间的积分可分为梯形法则,辛普森法则,和布尔法则,因此复合积分的每位新村间的积分可以是这种方式中的任意一种。

下面主要讲一下 复合矩形积分和复合辛普森积分。

1.2.1 复合矩形积分

小区间的界定又分为等长区间和不等长区间两种界定方式。

不等长区间,复合矩形积分估算公式:

在这里插入图片描述

区间不等长情况应用场景太有限。

等长区间,复合矩形积分估算公式:

在这里插入图片描述

等长区间,复合矩形积分估算偏差:

在这里插入图片描述

1.2.2 复合辛普森积分

同样的,复合辛普森积分区间的界定也可以分为等长区间和不等长区间两种界定方式。

由于不等长区间应用极少,因此直接介绍等长区间,复合辛普森1/3积分估算公式:

在这里插入图片描述

等长区间,复合辛普森1/3积分估算偏差:

在这里插入图片描述

由于3点积分(辛普森1/3积分)与4点积分(辛普森3/8积分)具有相同的偏差阶,因此通常优先使用复合辛普森1/3积分,复合辛普森3/8积分应用很少,因此本文不做介绍。感兴趣的朋友可以下去自己推一下复合辛普森3/8积分估算公式,应该也很简单。

(穿插)积分改进思路2

可以看出,将区间界定为更小的区间可以有效的增强积分估算的精度。那么我们是不是可以将积分区间不断地细分,从而至得到更精确的积分结果?理论上是可以的!!

方法1:(区间折半法)

等积法公式_积化和差与和差化积公式_三角函数积化和差和差化积公式

a.计算2个区间的积分,用复合辛普森1/3法则(当然,也可以用梯形法则或则布尔法则等都可以,但是推荐复合辛普森1/3法则);

b.将两个区间等分为4个区间,用复合辛普森1/3法则估算四个子区间的积分和;

c.比较a和b两步估算的结果相对误差是否接近一个太小的数:(Ia - Ib)/Ib

方法2:改进方法1
	a.对每个子区间一分为二后,比较两个区间的积分值与一个区间的积分值是否相近;
	b.如果相近则返回该值;否则对每个子区间执行a。(递归实现)

上面说,理论上是可以的!!理论上是可以的!!但是实际不行!!!

因为计算机浮点数运算是存在偏差的。当区间数目很大的时侯,浮点数估算的舍入误差(不理解的可以自行去学习了解一下)变得很大,会限制积分的提升。下面就是一个挺好的反例。

在这里插入图片描述

1.3常用的高精度积分

上述区间二分的思想是很有用的。在龙贝格积分和自适应积分中就会找到它的影子。

1.3.1 龙贝格积分

龙贝格积分的思想,及推论过程:

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

类推:联合两个偏差为O(h^4)的积分,改进后可以得到一个偏差为O(h ^ 6)的积分:

在这里插入图片描述

由此形成龙贝格积分:

在这里插入图片描述

在这里插入图片描述

积化和差与和差化积公式_等积法公式_三角函数积化和差和差化积公式

(穿插)积分改进思路3

龙贝格积分是用矩形公式推论得到的,那么我们是不是可以用辛普森公式做类似的推论,也能抵达一套新的公式呢?答案是:理论上是可以的!!

由辛普森1/3公式推论下来的加入偏差修正的积分计算式:

在这里插入图片描述

那么问题就来了。

我们的计算机抒发的double值是有范围的,而对于辛普森外,采用递归,递归8次就是16^8 = 2 ^ 32。计算机递归估算的结果很容易就形成溢出了。除非采用多精度浮点数估算才行,而采用多精度浮点数可能存在一些估算效率的问题。。。更深入的问题我就不继续深入剖析了(有兴趣的可以试试)。。。好了,总之来说,用辛普森推论下来的方式不是很实用。还是矩形龙贝格积分精典管用。

1.3.2 自适应积分

龙贝格积分中有两点:

a。每次就会将所有的步长折半;(那么我们能不能只对估算偏差较大的那种区间进行折半成两个子区间再积分呢?)

b。辛普森公式推论的方式,由于计算机double表示范围有限没用上。

下面还会在自适应积分中派上用场。自适应积分思想:(递归方式)

step1.辛普森1/3法则估算区间[a,b]的积分;

step2.区间步长折半,复合辛普森1/3法则估算区间[a,b]的积分;

如果两次偏差大于给定的误差限,停止执行,返回估算结果

在这里插入图片描述

如果两次偏差小于给定的误差限,对两个子区间分别执行step1和step2。

伪代码:

在这里插入图片描述

1.3.3 高斯积分

高斯积分应用比较广泛的是高斯-勒让德公式,此外还有高斯-切比雪夫公式;以及区间【0,正无穷】的广义积分 高斯-拉盖尔公式;以及区间【负无穷,正无穷】的广义积分 高斯-埃尔米特公式。

三角函数积化和差和差化积公式_积化和差与和差化积公式_等积法公式

高斯-勒让德公式

介绍高斯积分的文档就好多了,本文不细讲推论的过程,给出估算方式:

高斯-勒让德

多点高斯-勒让德公式:

在这里插入图片描述

在这里插入图片描述

高斯勒让德公式的区间为【-1,1】.对于任意区间为【a,b】;需要做一下简单的变量代换,

在这里插入图片描述

高斯积分偏差:

高斯勒让德公式的区间为【-1,1】,误差:

在这里插入图片描述

对于任意区间为【a等积法公式,b】,误差:

在这里插入图片描述

在这里插入图片描述

由上式继续往下推论,可以得出,当区间宽度【b-a】> 2(n+1)时,误差会发散。因此对于高斯点为n的高斯积分,积分区间宽度不要小于2(n+1),否则积分结果太可能会不确切。

此外,高斯积分的偏差还与被积函数的光顺性有关,被积函数光顺性越差,误差越大。详细可参考《现代数值积分》第5章 数值积分与数值微分,P139页)。

高斯积分的优点:高斯积分时给定节点数下代数精度最高的求积公式。

高斯积分的不足:高斯积分每次改变积分点的个数,所以的积分点的函数值都须要重新估算。

高斯-切比雪夫公式

高斯点为切比雪夫节点。

(未完待续。。。)

广义积分:高斯-拉盖尔公式 广义积分:高斯-埃尔米特公式 2 多重积分方式 二重积分的辛普森公式 二重积分的蒙特卡洛方式 多重(高维)积分的蒙特卡洛方式

亚博app|老王

新闻中心
新闻分类1
新闻分类2
支持与帮
隐私策略
反馈邮箱
admin@admin.com
产品中心
产品分类1
产品分类2
联系我们
公司地址:江苏XXX市XX区XX路88号
苏ICP备xxxxxxxx号
Copyright © 2012-2019 老王机械公司/老王有限公司 版权所有
亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app亚博app