自适应辛普森(Simpson)积分及二重积分
概述
辛普森积分法利用拉格朗日插值法来近似拟合原函数,然后通过对拟合函数的积分来简化近似。
近似公式:
∫ a b f ( x ) d x ≈ b − a 6 ( f ( a ) + 4 f ( a + b 2 ) + f ( b ) \int_{a}^{b}{f(x)}\operatorname d{x} \approx \frac{b-a}6(f(a)+4f(\frac{a+b}2)+f(b) ∫abf(x)dx≈6b−a(f(a)+4f(2a+b)+f(b)
拉格朗日插值
拉格朗日插值利用二次曲线(抛物线)来拟合。
在求通过点集 ( x 1 , y 1 ) , ( x 2 , y 2 ) , ( x 3 , y 3 ) , . . . , ( x n , y n ) {(x_1,y_1),(x_2,y_2),(x_3,y_3),...,(x_n,y_n)} (x1,y1),(x2,y2),(x3,y3),...,(xn,yn)的曲线上,拉格朗日利用多条二次曲线来拟合。首先二次曲线 g = f ( x ) g=f(x) g=f(x)满足: f i ( x j ) { 1 i = j 0 i = ̸ j {f_i(x_j)\left\{\begin{array}{lc}1&i=j\\0&i =\not j\end{array}\right.} fi(xj){ 10i=ji≠j。因此 y n f n ( x ) y_nf_n(x) ynfn(x)在 x n x_n xn处,取值为 y n y_n yn,其余点取值为0。
拉格朗日通过对 f i ( x ) f_i(x) fi(x)的构造如下:
f i ( x ) = ∏ j = ̸ i 1 ⩽ j ⩽ n ( x − x j ) ( x i − x j ) f_i(x)=\prod_{j=\not i}^{1\leqslant j \leqslant n}\frac{(x-x_j)}{(x_i-x_j)} fi(x)=j≠i∏1⩽j⩽n(xi−xj)(x−xj)
最终我们得到:
f ( x ) = ∑ i = 1 n y i f i ( x ) f(x)=\sum_{i=1}^{n}y_if_i(x) f(x)=i=