主题
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 方法。


另请参见

后向差分, 贝塞尔有限差分公式, 导数, 差分方程, 差分表, 埃弗雷特公式, 有限元方法, 前向差分, 高斯后向公式, 高斯前向公式, 插值, 杰克逊差分扇, 牛顿后向差分公式, 牛顿-科茨公式, 牛顿差商插值公式, 牛顿前向差分公式, 商差表, 递推方程, 斯蒂芬森公式, 斯特林有限差分公式, 影子微积分

使用 Wolfram|Alpha 探索

参考文献

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.

在 Wolfram|Alpha 中引用

有限差分

引用为

Weisstein, Eric W. "有限差分." 来自 MathWorld--Wolfram 网络资源。 https://mathworld.net.cn/FiniteDifference.html

主题分类