算法学习：插值型求积公式

算法学习：插值型求积公式

牛顿-柯斯特(Newton-Cotes)求积公式

定义

牛顿-柯斯特(Newton-Cotes)求积公式是插值型求积公式的特殊形式
在插值求积公式

$\int_a^bf(x)dx \approx \int_a^bP(x)dx=\sum\limits_{k=0}^n A_kf(x_k)$

中所取节点是等距时称为牛顿-柯斯特公式
其中插值多项式

$P(x)=\sum\limits_{k=0}^n\ell_k(x)f(x_k)$

求积系数

$A_k=\int_a^b\ell_k(x)dx$

这里$\ell_k(x)$指的是插值基函数，即有

$A_k=\int_a^b\ell_k(x)dx=\int_a^b\prod\limits_{j=0,j\neq k}^n\frac{x-x_j}{x_k-x_j}dx$

推导

在$[a,b]$区间内设置等距的插值基点$a=x_0<x_1\cdots <x_n=b$，设节点步长为$h=\frac{b-a}{n},x_k=a+kh,k\in \{0,1,\cdots ,n\}$
积分作变量替换$x=a+th$
$x_k-x_j=(k-j)h,x-x_j=(t-j)h,dx=hdt$
(xk−x0)⋯(xk−xk−1)(xk−xk+1)⋯(xk−x