(微积分的数值计算方法.docVIP

第七章 微积分的数值计算方法 7.1 微积分计算存在的问题/数值积分的基本概念 微分计算问题 求函数的导数(微分)，原则上没有问题。当然，这是指所求函数为连续形式且导数存在的情形。但如果函数一表格形式给出，要求函数在某点的导数值；或者是希望某点的导数值只用其附近离散点上的函数值近似地表示，这就是新问题了，它称为微分的数值计算，或称为数值微分。 2．定积分计算问题 计算函数在上的定积分 当被积函数的原函数能用有限形式给出时，可用积分基本公式来计算： 然而，问题在于：① 的原函数或者很难找到，或者根本不存在；②可能给出一个函数表；③仅仅知道是某个无穷级数的和或某个微分方程的解等等。这就迫使人们不得不寻求定积分的近似计算，也称数值积分。 3．数值积分的基本形式 数值积分的基本做法是构造形式如下的近似公式 (7.1.1) 或记成 (7.1.2) 和 分别成为上的的数值求积公式及其余项(截断误差)，和分别称为求积节点和求积系数(求积系数与被积函数无关)。 这种求积公式的特点是把求积过(极限过程)程转化为乘法与加法的代数运算。构造这种求积公式需要做的工作是：确定节点及系数，估计余项以及讨论的算法设计及其数值稳定性。 4．插值型求积公式 如何构造求积公式呢？基本的技术是用被积函数的Lagrange插值多项式近似代替，也即对上指定的个节点及相应的函数值，作 代入(7.1.2)式等号左边有 或写成形如(7.1.2)式的一般形式： (7.1.4) 其中 (7.1.5) (7.1.6) 称(7.1.4)为插值型求积公式。下面将要介绍的几种实用求积公式无不都是插值型公式的某种具体形式。 由上述定义，用余项公式可以衡量数值求积公式的精确程度。不过，由余项公式(7.1.6)可见，如果为次数的多项式，则中有,故，从而公式(7.1.4)成为等式 ，这就是说，当被积函数为次数的多项式时，其相应的插值型求积公式不是近似公式，而是准确公式。据此，人们引入了另一个衡量数值求积公式近似程度好坏的“代数精度” 概念。 5. 代数精度 定义7.1.1 如果数值求积公式，当是次多项式(可表示为)时，均有，也即有，便称求积公式至少具有次代数精度；如果当是次多项式(可表示为)时，只能是，也即有 ，便称求积公式具有 次代数精度。 定理 7.1.1 有个节点的插值型求积公式(7.1.4)，至少具有次代数精度。 例7.1.1 确定节点和系数，使得下列形式的求积公式 具有3次代数精度。 解 利用代数精度定义 取 令 计算得 即得至少具有3次代数精度的求积公式为 验证上式两端对，有 左边=2/5 ，右边=1/3 两者不相等，故可知所得公式正好为3次代数精度。 7.2 ?Newton-Cotes型求积公式 构造具体形式的插值型求积公式(6.1.4)的第一个想法是考虑求积节点为等距节点的情况(这时公式可能比较简单，而且节点也就定下来了)，这就是Newton-Cotes(牛顿-柯特斯公式)，而Newton-Cotes型公式当节点只有两点和三点时，就是熟知的梯形公式和Simpson(辛普生)公式。 Newton-Cotes 求积公式 考虑积分区间上的节点为等距节点 (7.2.1 ) 由(7.1.5)式可得插值型求积系数 (通过变换) 记 (7.2.2) 并称它为Cotes系数，从而得等距节点的插值型的阶Newton-Cotes求积公式，。(7.2.3) 显然，对不同的以及,按公式(7.2.2)算出，按公式可构造出不同阶的Newton-Cotes公式。 下面是Cotes系数表 （7.2.1） 1 1/2 1/2 2 1/6 4/6 1/6 3 1/8 3/8 3/8 1/8 4 7/90 16/45 2/15 16/45 7/90 …… …… …… …… …… …… 8 989/28350 5888/28350 -928/28350 10496/28350 -4540/28350 …… 989/28350 梯形公式/Simpson公式 当，即两个节点，取积分区间两个端点为节点时，由上表(7.