在数值分析中,数值积分是计算定积分数值的方法和理论。在数学分析中,给定函数的定积分的计算不总是可行的。许多定积分不能用已知的积分公式得到精确值。数值积分是利用黎曼积分和积分中值等数学定义和定理,用数值逼近的方法近似计算给定的定积分值。借助于电子计算设备,数值积分可以快速而有效地计算复杂的积分,能够以简单的方法求解具体数值问题,但数值积分的难点在于计算时间有时会过长,有时会出现数值不稳定现象,需要较强的理论支撑。

黎曼积分(Riemann integral)

在实数分析中,由黎曼创立的黎曼积分(Riemann integral)首次对函数在给定区间上的积分给出了一个精确定义。对于一在区间上之给定非负函数,我们想要确定所代表的曲线与坐标轴所夹图形的面积,作为曲线与坐标轴所夹面积的黎曼积分。黎曼积分的核心思想就是试图通过无限逼近来确定这个积分值。如函数取负值,则相应的面积值亦取负值。

积分中值定理(Mean value theorem of integrals)

积分中值定理揭示了一种将积分化为函数值,或者是将复杂函数的积分化为简单函数的积分的方法,若函数f(x) 在 闭区间 [a, b]上连续,则在积分区间 [a, b]上至少存在一个点 ξ,使下式成立

Integral(f(x)) on [a, b] = f(ξ)(b – a)

其中,a、b、ξ满足:a≤ξ≤b

数值积分的必要性

数值积分的必要性源自计算函数的原函数的困难性。利用原函数计算定积分的方法建立在牛顿-莱布尼兹公式之上。然而,原函数可以用初等函数表示的函数为数不多,大部分的可积函数的积分无法用初等函数表示,甚至没有解析表达式(“积不出来”的函数)。例如常见的正态分布函数的原函数就无法用初等函数表示。

不仅如此,在很多实际应用中,可能只能知道积分函数在某些特定点的取值,或者积分函数可能是某个微分方程的解,这些都是无法用求原函数的方法计算函数的积分。另外,当积分区域是曲面、三维形体以至于高维流形时,牛顿-莱布尼兹公式也不再适用,因此只能使用数值积分计算函数的近似值。

矩形法

矩形法是一种计算定积分近似值的方法,其思想是求若干个矩形的面积之和,这些矩形的高由函数值来决定。将积分区间[a, b] 划分为n个长度相等的子区间,每个子区间的长度为(a-b)/n 。这些矩形左上角、右上角或顶边中点在被积函数上。这样,这些矩形的面积之和就约等于定积分的近似值。

由函数上的点为矩形的左上角、右上角或顶边中点来决定,又分别被称为下(左)矩形公式、上(右)矩形公式和中矩形公式。当 n 逐渐扩大时,此近似值更加准确。矩形法的计算本质上是与黎曼积分的定义相吻合的。上述的点无论取哪个值,最终和式的值都将趋近于定积分的值。

梯形法

为了计算出更加准确的定积分,采用梯形代替矩形计算定积分近似值,其思想是求若干个梯形的面积之和,这些梯形的长短边高由函数值来决定。这些梯形左上角和右上角在被积函数上。这样,这些梯形的面积之和就约等于定积分的近似值。

辛普森法(Simpson’s rule)

矩形法和梯形法都是用直线线段拟合函数曲线的方法,另一种形式是采用曲线段拟合函数,实现近似逼近的数值积分方法。辛普森法(Simpson’s rule)是以二次曲线逼近的方式取代矩形或梯形积分公式,以求得定积分的数值近似解。

一般插值方法

另一种数值积分的思路是用一个容易计算积分而又与原来的函数“相近”的函数来代替原来的函数。这里的“相近”是指两者在积分区间上定积分的值比较接近。最自然的想法是采用多项式函数。比如说,给定一个函数后,在积分区间中对原来的函数进行拉格朗日插值。得到拉格朗日插值多项式以后,计算这个多项式的积分。

拉格朗日插值(Lagrange Interpolation)

拉格朗日插值是一种多项式插值方法,可以找到一个多项式,其恰好在积分区间中取的各个点取到给定函数的值。这样的多项式称为拉格朗日(插值)多项式。

数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。对于给定的n+1个点,对应于它们的次数不超过n的拉格朗日多项式有且只有一个。

牛顿-科茨公式(Newton-Cotes rule / Newton-Cotes formula)

牛顿-柯特斯公式(Newton-Cotes rule / Newton-Cotes formula)是以拉格朗日多项式插值的一般方法。梯形法则和辛普森法则便是牛顿-柯特斯公式的特例情况。

由于该拉格朗日多项式的系数都是常数,所以积函数的系数都是常数。这种方法缺点是对于次数较高的多项式而有很大误差(龙格现象),不如高斯积分法。

龙格现象(Runge Phenomenon)

在数值分析领域中, 龙格现象是用高阶多项式进行多项式插值时所出现的问题。

在某些高阶多项式等距点xi 进行插值,那么插值结果就会出现震荡。可以证明,在多项式的阶数增高时插值误差甚至会趋向无限大。

解决龙格现象的办法是使用切比雪夫节点代替等距点可以减小震荡,在这种情况下,随着多项式阶次的增加最大误差逐渐减小。这个现象表明高阶多项式通常不适合用于插值。使用分段多项式样条可以避免这个问题。如果要减小插值误差,那么可以增加构成样条的多项式的数目,而不必是增加多项式的阶次。第一类切比雪夫多项式的根(即切比雪夫节点)可以用于多项式插值。相应的插值多项式能最大限度地降低龙格现象,并且提供多项式在连续函数的最佳一致逼近。

代数精度评估

数值积分方法的代数精度用于衡量原函数和数值积分结果两者的逼近程度。若E(f)=0对f(x)=x^k(k=0,1,…,d)精确成立,而当f(x)=x^(d+1)时不再是精确等式,则说求积公式的代数精度是d。根据K.外尔斯特拉斯的多项式逼近定理,就一般的连续函数而言,d越大E(f)越小,因此可以用代数精度的高低说明数值积分公式的优劣

梯形公式和中矩形公式具有1次代数精度,辛普森公式有3次代数精度

最后修改日期: 2018年10月20日

作者

留言

撰写回覆或留言

发布留言必须填写的电子邮件地址不会公开。