北京交通大学(数字分析研究生课程)5数值积分与数值微分.doc

第5章 数值积分与数值微分方法1 基本概念梯形公式中矩形公式 则上式为一个数值求积公式.称为求积系数，称为求积节点;而称 为求积余项或求积公式的截断误差。从定义可以看到，数值求积公式依赖于求积节点个数n、求积节点和求积系数，这三个量有一个发生变化，则产生不同的求积公式.定义1 若求积公式对于次数不超过的多项式准确成立，而对于次多项式不准确成立，则称该求积公式具有次代数精度为. 一般，一个求积公式的代数精度越大，则该求积公式越好.确定代数精度的方法依次取代入公式并验证是否成立.若第一个使不成立的值为，则对应的代数精度为.例1确定求积公式的代数精度.解 取代入求积公式有易验证，但，故本题求积公式代数精度为3.例 2确定下面求积公式的参数A，B，C，使它具有尽可能高的代数精度，并指出相应的代数精度.解 本题要先求出具体的求积公式，然后再判断所求公式的代数精度。公式有3个待定参数，故利用3个条件得到的3个等式关系就可以解决求出具体求积公式的问题.依次取代入求积公式并取等号，有解之得故所求的求积公式为为确定其代数精度，再取代入求出的公式继续计算，有，故所求的求积公式具有二次代数精度.插值型求积公式 考虑关于个节点的Lagrange插值多项式与的余项，有这里两边取积分,有记 则有 若舍去，得求积公式(求积系数)该公式是插值型求积公式。插值型求积公式的求积余项当为次数不超过次的多项式时，有，对应的. 因此个节点的插值型求积公式的代数精度至少为若求积公式的代数精度至少是，则该公式是插值型求积公式.2. Newton-Cotes求积公式点的Newton-Cotes公式将求积节点取为a,b上的等距节点做积分变量变换 有 记称为Cotes系数.求积公式 称为Newton-Cotes求积公式.易验证2 点的Newton-Cotes公式 这正是我们熟悉的梯形公式.3点的Newton-Cotes公式为 称它为Simpson公式.例1 试分别用梯形公式和Simpson公式计算解 用梯形公式计算，有用Simpson公式计算，有梯形公式与Simpson公式的余项 梯形公式余项为利用积分中值定理可有梯形公式余项 Simpson公式的余项 部分Cotes系数n12345678 当较大时Cotes系数会出现负数，此时Newton-Cotes不具有数值稳定性，因而一般不用较大的Newton-Cotes公式来做计算.3 复化求积公式1）复化梯形公式取等距节点将积分区间a,b n等分,在每个小区间上用梯形公式做近似计算，就有得求积公式-复化梯形公式 复化梯形公式的余项 记故复化梯形公式的求积余项 2) 复化Simpson公式取等距节点将积分区间a,b n等分,在每个小区间上用Simpson公式做近似计算，再累加起来就有式中，得复化Simpson公式 复化Simpson公式的余项记有复化Simpson公式的求积余项例1 分别用复化梯形公式和复化Simpson公式计算，要求误差不超过.解 数值计算结果列表，其中代表求积余项.N复化梯形公式复化Simpson公式2-17.389 2595.32-11.592 840-0.47822-13.336 0231.27-11.984 94423-12.382 1620.312-12.064 20924-12.148 004-12.069 95125-12.089 74226-12.075 19427-12.071 55828-12.070 649本题积分的准确值为，可见复化梯形公式和复化Simpson公式能求出精度较高的解。例2 考虑用复化Simpson公式计算要使误差小于，那么求积区间0,1应分成多少个子区间？以此计算积分近似值。解 复化Si