数值积分的若干问题研究.pdfVIP

第一章绪论 1．1 引言 设，是区间【o，6】上某一给定的可积函数，现在要求计算定积分．ef(z)dt． 微积分教课书中有许多计算这类定积分的方法，或是借助于原函数，或是不经 过原函数，而是借助于各种各样的，大部分是技巧性的方法来计算．当然，对 于莱些可以求出精确表达式，但是我们又不熟悉的积分，我们可以参考积分 表．然而，随着电子计算机的高速发展，人们往往不满足于现有的积分表，更 何况现有的这些积分表或多或少存在一些数学上的或者印刷上的错误．Klerer 和Grossmma曾经专门研究了八种十分流行的数学用表，发觉表中的差错率几 乎超过百分之五． 不管是借助于原函数还是查找积分表，用这些方法只能解决很狭隘的一类 积分，在它的范围外通常采用各种近似计算的方法．在工程计算中，由于许多函 数的不定积分无法用简单函数表达出来，甚至函数本身都无法详尽地描述，而 代之以表格的形式给出一些离散点上的函数值．在上述这些情况下，理论上的 积分式就可能完全用不上了，我们必须采用数值积分．数值积分是研究如何求 出一个积分的数值，确切地说是积分的近似值．为此，我们定义任何一个能够近 似定积分Cf(z)dt的数值表达式，通常称为数值积分公式或求积公式． 当然，可积函数的种类是极其多种多样的．若考虑函数的近似积分的某个 完全确定的方法，则不可能对所有一般可积函数预先指出逼近的信值，这种估 值显然等于无穷大．例如，考虑梯形公式．显然，可以构造在区间【n，6】端点处 等于零的函数而使它在这区间上的定积分大于任何预先给定的数．按梯形公 式逼近这样函数的定积分所得的误差就等于这个积分本身，从而可以任意的 大f511．所以，数值积分的基本问题就是针对某些函数类，选择合适的求积节点 和求积系数，使得求积公式尽可能好地逼近．Cf(x)dt．正如前面所说的，在许 多情况下，我们所面l临的问题是要从实验数据求出积分，这时理论上的积分式 就可能完全用不上了．这样，寻求基于给定信息的求积公式就变得十分必要和 有意义了． 另外，我们遇到的函数可能不光滑，甚至并不连续，有时函数还会发生振荡 第一章 绪论 2 或者在某点无界．在这些情况下，我们必须考虑一些特殊的处理手段，以适应不 同的函数类．若积分式中的被积函数在求积区间上无界或积分域无界，则这种 积分就是大家熟知的反常积分(也称之为奇异积分或非正常积分)．我们知道， 微分方程边值问题的解都可以用积分的形式表示，而且也常会遇到一些被积函 数在积分区域内具有奇异性的积分．另外，奇异积分方程在力学、空气动力学及 数学物理等学科中的应用非常广泛．而上述问题中，对奇异积分的处理是数值 计算实现的关键．由此可见，如何计算奇异积分的近似值是一类很重要的问题． 1．2最佳求积 在介绍最佳求积公式之前，我们先来看看一些经典的求积公式． 1．2．1经典的求积公式 数值积分，即用被积函数值的一种线性组合(偶尔也采用非线性组合)来 逼近所求积分： ，6 n 札捌(，)：=／f(t)dt≈∑Ai，(瓯)， (1．1) 。o 1商1 在(1．1)式中，z1，茁2，…，％称为求积公式的节点，A1，A2，…，k称为求积系 数，或者称为相应于这些求积节点的n个“权”．有时，(1．1)式最后项还可能出 现被积函数的导数．我们将与函数值有关的信息称为Lagrange信息，与导数值 有关的信息，则称之为Hermite信息． 所以，只要确定节点≈和权九，求积公式也就确定了．一种很经典也最简 单的方法，是取区间【a，6】上的等距节点作为求积节点，即Xi=o+(i--1)^，h= 项式： 胁)：=善Ⅱ蓑mt)121… 因此，我们可以用丑。，q(厶)作为％棚(，)的一个近似，这样得到的求积公式 ％捌(厶)称为Newton—Cotes求积公式．这类求积公式最大的优点就是表达式简 单，而且当n=2时，得到梯形公式；当礼=3时，就是Simpson公式或抛物型 第一章 绪论