几种数值积分方法的对比分析.pdfVIP

2014年第 6期 信 息通 信 2014 (总第 138期) INFORM_AT10N C0～ 仉 NICAT10NS (Sum．No138) 几种数值积分方法的对比分析 樊 一园，黄伟伟，柳 华 ，景 杨，吝维军 (扬州大学数学科学学院。江苏 扬州 225002) 摘要：首先简单介绍了复化牛顿一科特斯积分、高斯积分、自适应积分和龙贝格积分的基本思想，然后从节点选取、误差 及运算速度对这四种积分方法进行了对比分析，最终得 出了四种积分方法的适用条件及优缺点。 关键词：数值积分；Matlab；-i~差 中图分类号：P312．4 文献标识码：A 文章编号：1673．1131(2014)06．0020．02 O引言 (J1)=9[f(a)+2∑，(口+砌)+，(易)】 (4) 在数学分析中，给定函数的定积分的计算不总是可行的， 许多定积分不能用已知的积分公式得到精确值。因此，数值积 ()= (JjI) (5) 分在生产生活中显得尤为重要，借助于Matlab等软件，数值积 (：)一()2m() 分可以快速而有效地计算复杂的积分，能够以简单的方法求解 +l()=—— — ——～ ( =1,2，…)(6) 具体数值问题，但不同的数值积分方法在耗时、精度等方面不 1一(二 同，有时更出现数值不稳定现象，因此有必要对不同积分方法 、2 进行对比分析。 1．4 自适应积分 本文针对 目前较为常用的复化牛顿一科特斯积分 (New． 自适应的求积方法有不同种形式，本文仅引用Simpson求 ton．Cotes)、高斯积分(Gauss)、自适应辛普森积分(AdaptSim— 积公式导出的一种 自适应的求积形式。基本思想如下： pson)和龙贝格积分 (Romberg)这四种积分方法进行对比研 在区间[a，b]上利用辛普森求积公式算出近似结果S，若f 究。由文 [1】我们知道，牛顿．科特斯求积 公式当n≥8时，稳 (x)在 [a，b]有4阶连续导数，则存在芎。使得 定性得不到保证，故常利用复化的方法来提高精度。由文 [2】 知取n=4时，有较好的积分效果，因此本文取n=4(即复化布 r )ax=S一 () (7) h：—b-a (8) 尔求积公式)进行对 比研究。而在高斯积分公式中，以权函 — 数p(x)一l尤为常见，故本文选取p(x) 1(即高斯勒让德积 将区间[a，b]分为两个相等的子区间，在每个子区间上利用 分)进行对 比研究。 辛普森求积公式得近似积分结果S。、S：，则 l算法的简介 r，(x)ax=St+S2-1h5．，4()() (9) 1．1复化布尔积分 将积分区间[a，b】等分为N个子区间[ 一4，4]，子区 设，‘’() f ‘’()贝0误差估计公式为 间长度h= ，在每个子区间上用阶数n=4时的牛顿科特斯 R(f)= l·+ 2一 l (1o) 公式 (即布尔公式)求积分，则有 若 R(Oe，则停止计算，S+S 为输出结果；否则，对每个 r )dx嘉 ㈩+32毫 3)+l2妻 …)+… 子区间再重复上述步骤。 3