基于复化辛卜生求积公式的变步长求积算法
1.辛卜生求积公式
对于 I ∗ = ∫ a b f ( x ) d x I^*=\int_{a}^{b}{f(x)dx} I∗=∫abf(x)dx ,差值型求积公式就是构造多项式 P n ( x ) P_n(x) Pn(x),使 I ∗ ≈ I = I ∗ = ∫ a b f ( x ) d x I^*\approx I= I^*=\int_{a}^{b}{f(x)dx} I∗≈I=I∗=∫abf(x)dx。
如果以 a , c = a + b 2 , b a,c=\frac{a+b}{2},b a,c=2a+b,b为三个差值节点,构造二次插值多项式
P 2 ( X ) = ( x − c ) ( x − b ) ( a − c ) ( a − b ) f ( a ) + ( x − a ) ( x − b ) ( c − a ) ( c − b ) f ( c ) + ( x − a ) ( x − c ) ( b − a ) ( b − c ) ) f ( b ) P_2(X)=\frac{(x-c)(x-b)}{(a-c)(a-b)}f(a)+\frac{(x-a)(x-b)}{(c-a)(c-b)}f(c)+\frac{(x-a)(x-c)}{(b-a)(b-c))}f(b) P2(X)=(a−c)(a−b)(x−c)(x−b)f(a)+(c−a)(c−b)(x−a)(x