1. 程式人生 > 其它 >第六章 數值積分與數值微分

第六章 數值積分與數值微分

6.1 積分與數值積分

6.1.1 定積分簡介

給定有界函式\(f(x)\)以及區間\([a,b]\),任取一組分點

\(\qquad a=x_0<x_1<x_2<\cdots<x_{n-1}<x_n=b\)

把區間\([a,b]\)分成n個小區間\([x_i,x_{i+1}],i=0:n-1\),再任取\(\xi_i\in[x_i,x_{i+1}]\),令

\(\qquad R_n=\sum\limits_{i=0}^{n-1}f(\xi_i)\Delta x_i,\quad\Delta x_i=x_{i+1}-x_i,\)

\(\lambda=\max\limits_{0\le i\le n-1}\{\Delta x_i\}\)

,如果不論\([a,b]\)怎麼分,不論\(\xi_i\)如何選取,只要\(\lambda\rightarrow0\),和式的極限都存在,則把它稱為函式\(f(x)\)\([a,b]\)上的定積分,記為

\(\qquad I=\int_a^bf(x)\mbox{d}x=\lim\limits_{\lambda\rightarrow 0}\sum\limits_{i=0}^{n-1}f(\xi_i)\Delta x_i\)

和式\(R_n=\sum\limits_{i=0}^{n-1}f(\xi_i)\Delta x_i\)稱為Riemann和。

通常來說,函式\(f(x)\)只要滿足連續、分段連續、單調這三者之一,定積分總存在。

1)Riemann和能否近似計算定積分

答案是否定的,為得到較好的近似效果,\(\xi_i\)\(x_i\)都要針對性選取,且n通常會比較大,計算成本較高與初心相悖。

2)牛頓-萊布尼茨公式等解析方法的侷限性

一些函式找不到初等函式表示的原函式,\(f(x)\)還可能是一個表函式,根本不知道具體的表示式。

6.1.2 數值積分

2)數值求積分的思路1

有效利用插值多項式進行數值求積

注意到\(p_n(x)=y_0l_0(x)+\cdots+y_nl_n(x),\)

\(\int_{a}^bp_n(x)\mbox{d}x=\sum\limits_{k=0}^{n}y_k\int_{a}^bl_k(x)\mbox{d}x=\sum\limits_{k=0}^nA_ky_k,\)

其中\(A_k=\int_{a}^bl_k(x)\mbox{d}x\),這裡\(A_k\)只依賴節點,不依賴於被積函式,具有一定可重複性。

3)數值求積分的思路2