计算方法第4章数值积分.docVIP

第四章 数值积分 ( Numerical Integration and Differentiation) 记号： ——积分， ——求积公式， 误差： ， ( Quadrature Formula ) 节点：， 求积系数：； 4.1 内插求积, Newton-Cotes公式 利用插值多项式： 积分 ; 例如,插值多项式取Lagrange形式 , 便有 及 此类通过在节点 处满足插值条件的插值多项式导出的求积公式称为内插求积公式。其中： 求积系数 ， 误差 显然，当函数为不超过次的多项式时，必有.由此可见,若函数为不超过某次的多项式，便有求积公式的误差，即. 定义：若求积公式对任何不超过m次的多项式成立：，而对次多项式等式不再成立，则称该求积公式的代数精度为m (事实上,可对一般的函数与其相应的近似公式定义代数精度). 代数精度——何时没有误差。 本节讨论的公式的节点为等距分布： ——称为Newton-Cotes 公式 4.1.1 Newton-Cotes 公式 1）n=1，梯形(Trapezoidal)公式 ， 误差：利用积分中值定理（参见附录I，定理I.6：积分中值定理） 2）n=2 Simpson公式， 误差 注：该误差公式不能简单地对仅用此三点的插值多项式的余项公式积分获得。 3）n=4 Cotes 公式 4.1.2 复化求积公式(Composite Numerical Integration) 梯形公式、Simpson公式的 优点：简单，容易估计误差； 缺点：大区间 误差大 改进：利用积分的可加性 ， 1）复化梯形公式 误差：若连续，由连续函数的介值定理(见例I.2)，及， 结论：代数精度 1 ； 计算精度 —— 一般，若有误差，误差与步长的关联程度； （ 若，则 ） 2）复化Simpson公式 误差：若连续，与复化梯形公式误差的推导一样，有 结论：代数精度 3 计算精度 ； 4.1.3 变步长积分法——步长的选取 通常，对数值积分有一定的误差要求：给定误差界，使； 1) 复化梯形积分法，由误差公式 设在变化不大，即将以上两式相除，即得 由此，得梯形积分的误差估计： ； 从而，由可计算的值 估计； 由于这是由约等式导出的，因此通常取 估计数值积分是否达到误差要求。 2) Simpson积分公式，由误差公式 设在变化不大，即将以上两式相除，即得 ； 从而，由 ； 3) 变步长 比较步长之比为2:1的复化梯形积分公式： ， ， 可得（理解）: 由此可见，容易通过计算步长对分以后两个不同的复化梯形积分公式计算值的差，估计复化梯形积分公式的误差，由此确定是否达到所要求的精度。 4.1.4 Romberg积分 由梯形积分的误差估计： ，设想将误差部分“归还”给梯形积分，应可以得到积分的一个更精确的近似值。通过计算，有： = = ； 注意，梯形积分公式 的精度是，而通过适当的组合，可得到Simpson公式，它的精度是，精度得到很大的提高；同时代数精度也得到提高。这种方法并没有增加函数计算量，仅增加了极少的代数运算，就能大大地提高计算精度，这类方法通常被称为“外推法”。 用类似的方法，将Simpson公式适当地组合，可得到精度为，代数精度为5的 Cotes公式；再将Cotes公式适当地组合，可得到精度为，代数精度为7的 Romberg公式： ； 在计算时，可列出如下的表进行计算： 一般, Romberg 方法用到 为止，尽管理论上还可以继续进行： 以取得精度更高的公式。但我们可以从两方面衡量该公式的实际意义。 1、将此公式改为：，由于，因此可将此式认作是对的修正，而修正值是；当较大时，与已是积分的比较准确的近似值，因此这部分的修正对的进一步改进影响一般比较小； 2、注意到公式的误差：，公式的误差：，公 式的误差：，公式的误差需要用到f的高阶导数，而f的高阶导数的性态是很难确定的（如同插值多项式余项估计）， 因此，一般不再继续进行此类外推（Romberg）方法，仅用到 为止。 例 积分： 例： 积分： （效果并不明显） 由此，也可从此体会到，Romberg 方法的效果（准确性），与被积函数的高阶导数的性态有关。 4.1.5 待定系数法 例：求确定以下公式的系数 使之具有尽可能高的代数精度，求其代数精度，并求其误差： 解：取插值多项式，先求差商表： 由此，通过对该插值多项式积分，得；因