數值積分

在數值分析中，數值積分（英語：Numerical integration）是計算定積分數值的方法和理論。在數學分析中，給定函數的定積分的計算不總是可行的。許多定積分不能用已知的積分公式得到精確值。數值積分是利用黎曼積分等數學定義，用數值逼近的方法近似計算給定的定積分值。藉助於電子計算設備，數值積分可以快速而有效地計算複雜的積分。

必要性

數值積分的必要性源自計算函數的原函數的困難性。利用原函數計算定積分的方法建立在牛頓-萊布尼茲公式之上。然而，原函數可以用初等函數表示的函數為數不多，大部分的可積函數的積分無法用初等函數表示，甚至無法有解析表達式。例如常見的正態分佈函數：

f(x)=e^{-{\frac {x^{2}}{2}}}

的原函數就無法用初等函數表示。

不僅如此，在很多實際應用中，只能知道積分函數在某些特定點的取值，比如天氣測量中的氣溫、濕度、氣壓等，醫學測量中的血壓、濃度等等。另外，積分函數有可能是某個微分方程的解。由於很多微分方程只能數值求解，因此只能知道函數在某些點上的取值。這時是無法用求原函數的方法計算函數的積分的。

另外，當積分區域是曲面、三維形體以至於高維流形時，牛頓-萊布尼茲公式不再適用，只能使用更廣泛的格林公式或斯托克斯公式，以轉化為較低維數上的積分，但只能用於少數情況。因此，只能使用數值積分計算函數的近似值。

矩形法

用一系列矩形的和來逼近積分的精確值。

矩形法是一種常見的數值積分方法，用來計算一維定積分的近似值。矩形法的主要思想是將積分區間 $I=[a,b]$ 分割成許多足夠小的分區間的總和： $I=[x_{0},x_{1})\cup [x_{1},x_{2})\cup \cdots \cup [x_{n-1},x_{n}]$ ， $a=x_{0}<x_{1}<x_{2}<\cdots <x_{n-1}<x_{n}=b$ ，使得能夠假設積分函數 $f$ 在各個小區間 $[x_{i},x_{i+1})$ 上的取值變化不大。這時，可以在每個分區間上取一個代表性的點 $c_{i}\in [x_{i-1},x_{i})$ （稱為節點），並將分區間的長度乘以積分函數在這一點上的值，以近似得到函數在這一段小區間上的積分。直觀上來看，就是取一個矩形，用它的面積來代替積分函數的曲線在這一小段區間上圍出來的曲邊梯形的面積。總體上，將所有這樣的矩形面積加起來（這個和稱為黎曼和），就近似地等於函數在這個區間上的定積分。

\int _{I}f(t)dt\approx \mathbf {I} _{\mbox{appr}}(c_{1},x_{1},c_{2},x_{2},\cdots ,c_{n},x_{n})=\sum _{i=0}^{n-1}(x_{i+1}-x_{i})f(c_{i})

根據黎曼積分的定義，只要區間被分得足夠精細，那麼這樣的分割所得到的黎曼和會無限趨近於函數的積分。

公式

根據每個小區間中節點的選取方式，可以得到不同的數值積分公式。

上矩形公式：取每個小區間中的「最高點」（ $f$ 的最大值或上確界）作為節點。

\mathbf {I} _{\mbox{appr}}(c_{1},x_{1},c_{2},x_{2},\cdots ,c_{n},x_{n})=\sum _{i=0}^{n-1}(x_{i+1}-x_{i})\sup(f_{[x_{i-1},x_{i})})

下矩形公式：取每個小區間中的「最低點」（ $f$ 的最小值或下確界）作為節點。

\mathbf {I} _{\mbox{appr}}(c_{1},x_{1},c_{2},x_{2},\cdots ,c_{n},x_{n})=\sum _{i=0}^{n-1}(x_{i+1}-x_{i})\inf(f_{[x_{i-1},x_{i})})

中矩形公式：取每個小區間中央的一點作為節點。

\mathbf {I} _{\mbox{appr}}(c_{1},x_{1},c_{2},x_{2},\cdots ,c_{n},x_{n})=\sum _{i=0}^{n-1}(x_{i+1}-x_{i})f({\frac {x_{i-1}+x_{i}}{2}})

插值法

另一種數值積分的思路是用一個容易計算積分而又與原來的函數「相近」的函數來代替原來的函數。這裏的「相近」是指兩者在積分區間上定積分的值比較接近。最自然的想法是採用多項式函數。比如說，給定一個函數 $f$ 後，在積分區間 $I=[a,b]$ 中取 $a=x_{0}<x_{1}<x_{2}<\cdots <x_{n-1}<x_{n}=b$ ，就可以對原來的函數進行拉格朗日插值。得到拉格朗日插值多項式以後，計算這個多項式的積分。

\int _{I}f(t)dt\approx \mathbf {I} _{\mbox{appr}}(x_{1},x_{2},\cdots ,x_{n})=\sum _{i=0}^{n}f(x_{i})L_{i}

其中 $L_{i}$ 是拉格朗日插值的基本多項式。

牛頓-柯特斯公式

牛頓-柯特斯公式是一種插值型公式。假設 $I=[a,b]$ 中取 $x_{k}=a+k{\frac {b-a}{n}}$ ,可以寫成：

\int _{I}f(t)dt\approx \mathbf {I} _{\mbox{appr}}(x_{1},x_{2},\cdots ,x_{n})=(b-a)\sum _{k=0}^{n}C_{k}^{(n)}f(x_{k})

其中的 $\ C_{k}^{(n)}={\frac {1}{n}}\int _{0}^{n}\prod _{k\neq j}{\frac {t-j}{k-j}}dt={\frac {(-1)^{n-k}}{n\cdot k!(n-k)!}}\int _{0}^{n}\prod _{k\neq j}(t-j)dt$

如果 $n=1$ ，那麼牛頓-柯特斯公式就變成梯形公式（取每個小區間兩端點，做成梯形，梯形的值也和矩形一樣，趨於原來的函數的積分）。

\mathbf {I} _{\mbox{appr}}(x_{1},x_{2},\cdots ,x_{n})=\sum _{i=0}^{n-1}(x_{i+1}-x_{i})({\frac {f(x_{i+1})+f(x_{i})}{2}})

參見

參考來源

李慶揚. 《数值分析（第四版）》. 華中科技大學出版社. 2006. ISBN 7-560-93742-X.