复化科特斯公式matlab_【原创】牛顿-柯特斯数值积分公式及其MATLAB的实现

一、数值积分基本公式

数值求积基本通用公式如下

Eqn1.gif (1.63 KB, 下载次数: 40)

2009-11-20 23:23 上传

xk：求积节点

Ak：求积系数，与f(x)无关

数值积分要做的就是确定上式中的节点xk和系数Ak。可以证明当求积系数Ak全为正时，上述数值积分计算过程是稳定。

二、插值型数值积分公式

对f(x)给定的n+1个节点进行Lagrange多项式插值，故

Eqn2.gif (2.95 KB, 下载次数: 32)

2009-11-20 23:23 上传

即求积系数为

Eqn3.gif (3.29 KB, 下载次数: 32)

2009-11-20 23:23 上传

三、牛顿-柯特斯数值积分公式

当求积节点在[a,b]等间距分布时，插值型积分公式(先使用Lagrange对节点进行多项式插值，再计算求积系数，最后求积分值)称为Newton-Cotes积分公式。

由于Newton-Cotes积分是通过Lagrange多项式插值变化而来的，我们都知道高次多项式插值会出现Runge振荡现象，因此会导致高阶Newton-Cotes公式不稳定。

Newton-Cotes积分公式的求积系数为

Eqn4.gif (3.38 KB, 下载次数: 31)