主題
Search

有限差分


有限差分是導數的離散 аналог(模擬)。函式 f_p 的有限前向差分定義為

 Deltaf_p=f_(p+1)-f_p,
(1)

有限後向差分定義為

 del f_p=f_p-f_(p-1).
(2)

前向有限差分在 Wolfram 語言中實現為DifferenceDelta[f, i]。

如果值以間距 h 製表,則使用符號

 f_p=f(x_0+ph)=f(x)
(3)

k前向差分然後可以寫成 Delta^kf_p,類似地,k後向差分可以寫成 del ^kf_p

然而,當 f_p 被視為連續函式 f(x) 的離散化時,有限差分有時寫為

Deltaf(x)=f(x+1/2)-f(x-1/2)
(4)
=2AdjustmentBox[I, BoxMargins -> {{0.101266, -0.101266}, {0.375, -0.375}}, BoxBaselineShift -> -0.375]AdjustmentBox[I, BoxMargins -> {{0, 0}, {-0.25, 0.25}}, BoxBaselineShift -> 0.25](x)*f(x),
(5)

其中 * 表示卷積AdjustmentBox[I, BoxMargins -> {{0.101266, -0.101266}, {0.375, -0.375}}, BoxBaselineShift -> -0.375]AdjustmentBox[I, BoxMargins -> {{0, 0}, {-0.25, 0.25}}, BoxBaselineShift -> 0.25](x) 是奇脈衝對。因此,有限差分運算元可以寫成

 Delta^~=2AdjustmentBox[I, BoxMargins -> {{0.101266, -0.101266}, {0.375, -0.375}}, BoxBaselineShift -> -0.375]AdjustmentBox[I, BoxMargins -> {{0, 0}, {-0.25, 0.25}}, BoxBaselineShift -> 0.25]*.
(6)

n具有常數 n 階有限差分。例如,取 n=3 並製作一個差分表

 x; 1; 2; 3; 4; 5x^3; 1; 8; 27; 64; 125Delta; 7; 19; 37; 61Delta^2; 12; 18; 24Delta^3; 6; 6Delta^4; 0.
(7)

Delta^3 列是常數 6。

有限差分公式對於外推有限量的資料以嘗試找到通項非常有用。具體來說,如果函式 f(n) 僅在少數離散值 n=0, 1, 2, ... 處已知,並且希望確定 f 的解析形式,則如果 f 被假定為多項式函式,則可以使用以下過程。將感興趣的序列中的第 n 個值表示為 a_n。然後將 b_n 定義為前向差分 Delta_n=a_(n+1)-a_n,將 c_n 定義為二階前向差分 Delta_n^2=b_(n+1)-b_n,等等,構建如下表

 a_0=f(0)  a_1=f(1)  a_2=f(2)  ...  a_p=f(p)
(8)
   b_0=a_1-a_0  b_1=a_2-a_1  ...  b_(p-1)=a_p-a_(p-1)
(9)
 c_0=b_1-b_0  ...  ...
(10)
 ...
(11)

繼續計算 d_0e_0 等,直到獲得 0 值。然後,給出值 a_n多項式函式由下式給出

f(n)=sum_(k=0)^(p)alpha_k(n; k)
(12)
=a_0+b_0n+(c_0n(n-1))/2+(d_0n(n-1)(n-2))/(2·3)+....
(13)

當使用符號 Delta_0=a_0Delta_0^2=b_0 等時,這個優美的方程稱為牛頓前向差分公式。為了看一個具體的例子,考慮一個序列,其前幾個值為 1, 19, 143, 607, 1789, 4211 和 8539。然後差分表由下式給出

 1  19  143  607  1789  4211  8539
18  124  464  1182  2422  4328
106  340  718  1240  1906
234  378  522  666
144  144  144
0  0
(14)

讀取每行中的第一個數字得到 a_0=1b_0=18c_0=106d_0=234e_0=144。將這些代入得到方程

f(n)=1+18n+1/2106n(n-1)+1/6234n(n-1)(n-2)+1/(24)144n(n-1)(n-2)(n-3)
(15)
=6n^4+3n^3+2n^2+7n+1,
(16)

它確實完全符合原始資料。

導數的公式由下式給出

f^'(x_0)=1/h[f(x_0+h)-f(x_0)]-1/2hf^('')(xi)
(17)
=1/(2h)[-3f(x_0)+4f(x_0+h)-f(x_0+2h)]+1/3h^2f^((3))(xi)
(18)
=1/(2h)[f(x_0+h)-f(x_0-h)]-1/6h^2f^((3))(xi)
(19)
=1/(12h)(f_(-2)-8f_(-1)+8f_1-f_2)+1/(30)h^4f^((5))(xi)
(20)
=1/(12h)(-25f_0+48f_1-36f_2+16f_3-3f_4)+1/5h^4f^((5))(xi)
(21)
f^('')(x_0)=1/(h^2)(f_(-1)-2f_0+f_1)-1/(12)h^2f^((4))(xi)
(22)
=1/(h^2)(f_0-2f_1+f_2)+1/6h^2f^((4))(xi_1)-hf^((3))(xi_2)
(23)
f^((3))(x_0)=1/(h^3)(f_3-3f_2+3f_1-f_0)+O(h)
(24)
=1/(2h^3)(f_2-2f_1+2f_(-1)-f_(-2))+O(h^2)
(25)
f^((4))(x_0)=1/(h^4)(f_4-4f_3+6f_2-4f_1+f_0)+O(h)
(26)
=1/(h^4)(f_2-4f_1+6f_0-4f_(-1)+f_(-2))+O(h^2)
(27)

(Beyer 1987, pp. 449-451; Zwillinger 1995, p. 705)。

有限差分積分的公式

 int_(x_0)^(x_n)f(x)dx=hint_0^nf_pdp
(28)

由 Beyer (1987, pp. 455-456) 給出。

有限差分導致差分方程,即微分方程的有限 аналог(模擬)。事實上,影子微積分顯示了連續函式的許多優雅 аналог(模擬)的著名恆等式。偏微分方程的常用有限差分格式包括所謂的 Crank-Nicolson、Du Fort-Frankel 和 Laasonen 方法。


另請參見

後向差分, 貝塞爾有限差分公式, 導數, 差分方程, 差分表, 埃弗雷特公式, 有限元方法, 前向差分, 高斯後向公式, 高斯前向公式, 插值, 傑克遜差分扇, 牛頓後向差分公式, 牛頓-科茨公式, 牛頓差商插值公式, 牛頓前向差分公式, 商差表, 遞推方程, 斯蒂芬森公式, 斯特林有限差分公式, 影子微積分

使用 探索

參考文獻

Abramowitz, M. 和 Stegun, I. A. (編). "差異." §25.1 in 數學函式手冊,包含公式、圖表和數學表格,第 9 次印刷。 紐約: Dover, pp. 877-878, 1972.Beyer, W. H. CRC 標準數學表格,第 28 版。 Boca Raton, FL: CRC Press, pp. 429-515, 1987.Boole, G. 和 Moulton, J. F. 有限差分微積分專著,第 2 版修訂版。 紐約: Dover, 1960.Conway, J. H. 和 Guy, R. K. "牛頓有用的微小公式." In 數字之書。 紐約: Springer-Verlag, pp. 81-83, 1996.Fornberg, B. "有限差分公式中權重的計算." SIAM Rev. 40, 685-691, 1998.Iyanaga, S. 和 Kawada, Y. (編). "插值." 附錄 A,表 21 in 數學百科全書。 Cambridge, MA: MIT Press, pp. 1482-1483, 1980.Jordan, C. 有限差分微積分,第 3 版。 紐約: Chelsea, 1965.Levy, H. 和 Lessman, F. 有限差分方程。 紐約: Dover, 1992.Milne-Thomson, L. M. 有限差分微積分。 倫敦: Macmillan, 1951.Richardson, C. H. 有限差分微積分導論。 紐約: Van Nostrand, 1954.Spiegel, M. 有限差分和微分方程微積分。 紐約: McGraw-Hill, 1971.Stirling, J. Methodus differentialis, sive tractatus de summation et interpolation serierum infinitarium。 倫敦, 1730. Holliday, J. 的英文翻譯 微分方法:無限級數求和與插值專著。 1749.Tweddle, C. 詹姆斯·斯特林:生平和著作草圖及其科學通訊。 牛津,英格蘭: 牛津大學出版社, pp. 30-45, 1922.Weisstein, E. W. "關於有限差分方程的書籍." http://www.ericweisstein.com/encyclopedias/books/FiniteDifferenceEquations.html.Zwillinger, D. (編). "差分方程" 和 "數值微分." §3.9 和 8.3.2 in CRC 標準數學表格和公式。 Boca Raton, FL: CRC Press, pp. 228-235 和 705-705, 1995.

在 中引用

有限差分

引用為

Weisstein, Eric W. "有限差分." 來自 網路資源。 https://mathworld.tw/FiniteDifference.html

主題分類