数值积分-计算方法讲课.docVIP

数值积分 理论依据 逼近论——构造一个简单函数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)构造求积公式： 记，(k=0,1,…,n) 作代换x=a+th带入上式，变为： 其中： (k=0,1,…,n) （1-1） 这个积分是有理多项式积分，它与被积函数f(x)和区间[a,b]无关。只要确定n就能计算出系数。 于是得到称为Newton—Cotes公式的求积公式： （1-2） 其中称为Newton—Cotes系数。如表1所示。 表1 Newton—Cotes系数 n 1 1/2 1/2 2 1/6 4/6 1/6 3 1/8 3