目錄
- 1. 引言
- 幾個常用積分公式及其複合公式
1. 引言
高數中計算積分思路基本是牛頓萊布尼茲法:
\[I[f]=\int_{a}^{b}f(x)\mathrm{d}x=F(b)-F(a),
\]
\[F^{\prime}(x)=f(x).
\]
實際計算中,原函式一般無法求出.給不出解析解,只能求出數值解.
設在區間 [a,b]( 不妨先設 a,b 為有限數 ) 上 ,\(f(x) ≈ P_n (x), P_n (x)\) 為某個較“簡單”的函式 , 則顯然有
\[\int_a^bf(x)\mathrm{d}x\approx\int_a^bP_n(x)\mathrm{d}x.
\]
如果\(\operatorname*{max}_{a\leqslant x\leqslant b}|f(x)-P_{n}(x)|\leqslant\varepsilon\),則誤差估計:
\[\left|\int_a^bf(x)\mathrm{d}x-\int_a^bP_n(x)\mathrm{d}x\right|\leqslant(b-a)\varepsilon.
\]
幾個常用積分公式及其複合公式
- 中點公式
對f(x),使用\(f(\frac{a+b}{2})\)近似代替.有:
\[\int_a^bf(x)\mathrm{d}x\approx\int_a^bf\left(\frac{a+b}{2}\right)\mathrm{d}x=(b-a)f\left(\frac{a+b}{2}\right).
\]
誤差估計:
\[\int_a^bf(x)\mathrm{d}x-(b-a)f\left(\frac{a+b}{2}\right)=\frac{1}{12}(b-a)^3f''(\xi).
\]
- 梯形公式
拉格朗日插值多項式\(L_n(x)\):
\[L_n(x)=\sum_{j=0}^ny_jl_j(x)=\sum_{j=0}^ny_j\prod_{\substack{i=0\\i\neq j}}^n\frac{x-x_i}{x_j-x_i}.
\]
\(n=1\)時,\(L_1(x)=l_0(x)y_0+l_1(x)y_1\),用\(L_1(x)\) 近似代替 f(x) 稱為線性插值 , 公式(3.9)稱為線性插值多項式或一次插值多項式.即:
\[L_1(x)=\frac{x-x_1}{x_0-x_1}y_0+\frac{x-x_0}{x_1-x_0}y_1.
\]
\(n=2\)時,\(L_2(x)=l_0(x)y_0+l_1(x)y_1+l_2(x)y_2\),用\(L_2(x)\)近似代替 f(x) 稱為二次插值或拋物線插值 , 稱式 (3.10) 為二次插值多項式
\[L_2(x)=\frac{(x-x_1)(x-x_2)}{(x_0-x_1)(x_0-x_2)}y_0+\frac{(x-x_0)(x-x_2)}{(x_1-x_0)(x_1-x_2)}y_1+\frac{(x-x_0)(x-x_1)}{(x_2-x_0)(x_2-x_1)}y_2
\]
基於\(x=a,x=b\)兩節點構造線性插值函式\(L_1(x)\),近似代替原函式\(f(x)\),得到梯形公式.
\[L_1(x)=\frac{x-b}{a-b}f(a)+\frac{x-a}{b-a}f(b)
\]
\[\begin{aligned}
\int_{a}^{b}f(x)\mathrm{d}x& \approx\int_{a}^{b}L_{1}(x)\mathrm{d}x=\int_{a}^{b}\left[{\frac{x-b}{a-b}}f(a)+{\frac{x-a}{b-a}}f(b)\right]\mathrm{d}x \\
&=\frac{1}{2}(b-a)\left[f(a)+f(b)\right].
\end{aligned}\]
誤差估計:
\[\int_a^bf(x)\mathrm{d}x-\frac{1}{2}(b-a)\left[f(a)+f(b)\right]=-\frac{1}{12}(b-a)^3f''(\xi), \xi\in(a,b)
\]
- 辛普森 (Simpson) 公式(拋物型公式)
\(\text{若 }f(x)\text{ 用透過節點 }x_0=a, x_1=\frac{a+b}{2}, x_2=b\text{ 的二次插值多項式 }L_2(x)\text{ 代替}\)
\[f(x)\approx L_2(x)=\frac{(x-x_1)(x-x_2)}{(x_0-x_1)(x_0-x_2)}f(x_0)+\frac{(x-x_0)(x-x_2)}{(x_1-x_0)(x_1-x_2)}f(x_1)+\frac{(x-x_0)(x-x_1)}{(x_2-x_0)(x_2-x_1)}f(x_2)
\]
可以得到積分公式:
\[\int_a^bf(x)\mathrm{d}x\approx\int_a^bL_2(x)\mathrm{d}x=\frac16(b-a)\left[f(a)+4f\left(\frac{a+b}2\right)+f(b)\right].
\]
誤差估計:
\[\int_a^bf(x)\mathrm dx-\frac{1}{6}(b-a)\left[f(a)+4f\left(\frac{a+b}{2}\right)+f(b)\right]=-\frac{(b-a)^5}{2880}f^{(4)}(\xi),\quad\xi\in(a,b).
\]