数值分析课件第章数值积分与数值微分.docVIP

第4章 数值积分与数值微分 1 数值积分的基本概念 实际问题当中常常需要计算定积分。在微积分中，我们熟知，牛顿—莱布尼兹公式是计算定积分的一种有效工具，在理论和实际计算上有很大作用。对定积分，若在区间上连续，且的原函数为，则可计算定积分 似乎问题已经解决，其实不然。如 1）是由测量或数值计算以数据表形式给出时，Newton-Leibnitz公式无法应用。 2）许多形式上很简单的函数，例如 等等，它们的原函数不能用初等函数的有限形式表示。 3）即使有些被积函数的原函数能通过初等函数的有限形式表示，但应用牛顿—莱布尼兹公式计算，仍涉及大量的数值计算，还不如应用数值积分的方法来得方便，既节省工作量，又满足精度的要求。例如下列积分 对于上述这些情况，都要求建立定积分的近似计算方法—数值积分法。 1.1 数值求积分的基本思想 根据以上所述，数值求积公式应该避免用原函数表示，而由被积函数的值决定。由积分中值定理：对，存在，有 表明，定积分所表示的曲边梯形的面积等于底为而高为的矩形面积（图4-1）。问题在于点的具体位置一般是不知道的，因而难以准确算出。我们将称为区间上的平均高度。这样，只要对平均高度提供一种算法，相应地便获得一种数值求积分方法。 如果我们用两端的算术平均作为平均高度的近似值，这样导出的求积公式 (1.1) 便是我们所熟悉的梯形公式（图4-2）。而如果改用区间中点的“高度”近似地取代平均高度，则可导出所谓中矩形公式（简称矩形公式） (1.2) 更一般地，我们可以在区间上适当选取某些节点，然后用加权平均得到平均高度的近似值，这样构造出的求积公式具有下列形式： (1.3) 式中称为求积节点；成为求积系数，亦称伴随节点的权。权仅仅与节点的选取有关，而不依赖于被积函数的具体形式。 这类由积分区间上的某些点上处的函数值的线性组合作为定积分的近似值的求积公式通常称为机械求积公式，它避免了Newton-Leibnitz公式寻求原函数的困难。对于求积公式(1.3)，关键在于确定节点和相应的系数。 1.2 代数精度的概念 由Weierstrass定理可知，对闭区间上任意的连续函数，都可用多项式一致逼近。一般说来，多项式的次数越高，逼近程度越好。这样，如果求积公式对阶多项式精确成立，那么求积公式的误差仅来源于阶多项式对连续函数的逼近误差。因此自然有如下的定义 定义1 如果某个求积公式对于次数不超过的多项式均准确地成立，但对于次多项式就不准确成立，则称该求积公式具有次代数精度。 例1 判断求积公式 的代数精度。 解 记 因为 所以求积公式具有5次代数精度。 例2给定形如的求积公式，试确定系数，使公式具有尽可能高的代数精度。 解 求积公式中有三个参数，因此至少对 应精确成立，即 当时，得 当时，得 当时，得 解得，于是有 当时，，而上式右端为，故公式对不精确成立，其代数精度为2。 1.3插值型的求积公式 最直接自然的一种想法是用在上的插值多项式代替。由于代数多项式的原函数是容易求出的，我们以在上的积分值作为所求积分的近似值，即 这样得到的求积分公式称为插值型求积公式。通常采用Lagrange插值。 设上有个互异节点，的次Lagrange插值多项式为 其中，插值型求积公式为 (1.4) 其中。可看出，仅由积分区间与插值节点确定，与被积函数的形式无关。求积公式(1.4)的截断误差为 (1.5) 定义2 求积公式 如其系数，则称此求积公式为插值型求积公式。 定理1 形如(1.3)的求积公式至少有次代数精度的充分必要条件是插值型的。 证明 如果求积公式(1.3)是插值型的，由公式(1.5)可知，对于次数不超过的多项式，其余项等于零，因而这时求积公式至少具有次代数精度。 反之，如果求积公式(1.3)至少具有次代数精度，那么对于插值基函数应准确成立，并注意到，即有 所以求积公式(1.3)是插值型的。 1.4 广义皮亚诺定理及求积公式的余项 余项公式(1.5)在实际应用中很不方便。为了推导出统一的只包含函数导数的余项公式，需要引进皮亚诺定理。 先引进线性泛函的概念。 给定一个函数集合（函数空间），如果对于这个集合中的任何一个函数，都有一个确定的实数与之对应，则说给定了这个函数集合上的一个泛函。 设是定义在线性空间上的泛函，对函数集合中的函数，如果对于任何常数有 则称此泛函是线性泛函。 函数的函数值、积分值、导数值都是线性泛函。 广义皮亚诺（Peano）定理：设求积公式(1.5)的余项是空间上的线性泛函，且的代数精度为，那么对任意，有 (1.6) 其中 (1.7) 这里，其中 是区间上的任意点。 证明 设是以为节点的插值多项式，则 故 将皮亚诺定理应用于求积余项公式(1