數值积分的计算方法论文.docVIP

HYPERLINK \l 目录摘要 摘 要 本文应用插值积分法和逼近论的思想，简单重述了推导Newton-Cotes公式和Gauss-Legendre求积公式的过程，以及这两个公式的系数、精度等问题。并以这两种数值积分的求解方法为基础，应用quad、guass函数编写具体Matlab程序，通过计算机软件计算出所给题目的近似数值积分。对二者所得的结果进行比较，从而研究了用Newton-Cotes和Gauss-Legendre公式求积分的方法和二者的精确度问题。得知，这两种求积公式所得的结果在精度上的确存在差异，结合理论部分更加充分地说明了，n相同时Gauss-Legendre公式比Newton-Cotes公式具有更高的代数精度，但当代数精度相同时，二者计算的结果仍存在细微的差异。 关键字：插值积分、Newton-Cotes公式、Gauss-Legendre公式 数值积分 HYPERLINK \l 目录第1章 理论依据 逼近论——构造一个简单函数p(x)近似表示f(x),然后对 p(x)求积分得到 f(x)的积分的近似值。基于插值原理，推导出数值积分的基本公式。 §1插值求积公式 为了用数值方法求，对被积函数f(x)在给定的n+1个节点上作Lagrange插值，用插值函数Pn(x)代替f(x)，就可用I（Pn(x)）构造求积公式，近似地计算定积分I(f(x))。 §2Newton—Cotes公式 §2.1Newton—Cotes公式的推导 当§1.1插值求积公式的插值节点为等距节点时，就得到Newton—Cotes公式。 将区间[a,b]n等分，，n+1个节点为 xk=a+kh (k=0,1,…,n) 在节点上对f(x)的Lagrange插值多项式是： 用Pn(x)代替f(x)构造求积公式： 记Akn=ab 作代换x=a+th带入上式，变为： 其中：Ck(n)= 这个积分是有理多项式积分，它与被积函数f(x)和区间[a,b]无关。只要确定n就能计算出系数Ck 于是得到称为Newton—Cotes公式的求积公式： （1-2） 其中Ckn称为Newton— 表1 Newton—Cotes系数 n C 1 1/2 1/2 2 1/6 4/6 1/6 3 1/8 3/8 3/8 1/8 4 7/90 32/90 12/90 32/90 7/90 5 19/288 25/96 25/144 25/144 25/90 19/288 6 41/840 9/35 9/280 34/105 9/280 9/35 41/840 §2.2Newton—Cotes公式误差和稳定性 在积分公式中用插值多项式Pn(x)代替f(x)的插值误差是 因此，Newton—Cotes公式的截断误差是 （1-3） 讨论舍入误差对计算结果产生的影响，设（1-2）式近似计算 其中计算函数值f(xn)有误差值εn（k=0,1,2, …,n）。在（1-2）式中令fx≡1,pnx=1k=0n 如果Ckn皆为正，并设ε=max0≤k≤n|εk|，则|e 但当n8时，Ckn将出现负数，这时，数值计算的稳定性不能保证，所以节点超过8时Newton—Cotes公式 当n为偶数时，Newton—Cotes积分公式具有n+1次代数精度。 §2.3经典Newton—Cotes公式 当n=4,5点公式称为经典Newton—Cotes公式 其中xk=a+kb-a4 §3 Gauss-Legendre求积公式 在积分区间[a,b]内对积分节点不作限制，不取等距，积分节点和求积系数都作为待定未知量。通过适当选择节点和求积系数，能构造更有效的高精度求积公式。 §3.1计算n阶求积公式 若有m次代数精度，对(k=0,1,…)应有 i=0 而。 §3.2 Gauss求积公式的基本原理 更一般形式： （2-1） 为权函数，设0，且在[a,b]上可积，构造n阶求积公式： （2-2） 积分点使得（2-2）式达到2n+1次代数精度，则积分点称为Gauss点，（2-2）式称为Gauss求积公式。 §3.3 Gauss-Legendre求积公式 求积分，权数=1， 其中（i=0，1，…，n）是n+1阶Legendre多项式的零点，求积系数为： （i=0，1，…，n） 具体Gauss-Legendre公式的插值节点和系数见表2（其中n为插值节点个数，xi为积分点，Ai 表二Gauss-Legendre公式的插值节点和系数 n x A n x A 0 0 2 0 0.56888888889 1 ± 1 5 ± 0.1713244924 2 ± 0.555555556 ± 0.3607615370 0 0. ± 0 3 ± 0.3478548451 6 ± 0.1294849662