牛顿插值公式看起来非常复杂,除非您引入正确的符号。如果符号正确,它看起来就像泰勒级数。这种表示法不仅更简单、更容易记住,而且还暗示了扩展。
我们需要的符号分为两部分。首先,我们需要前向差分算子 Δ 定义为
及其扩展 Δ k通过将 Δ 应用于函数k次来定义。
我们需要的另一个符号是下降幂,用指数下方的小条表示:
Δ 算子类似于微分算子D ,下降幂类似于普通幂。使用这种表示法,牛顿插值公式基于
看起来类似于基于a的幂级数的泰勒公式
牛顿公式适用于多项式,无限和实际上是有限和,因为对于超出某个点的所有k, Δ k f ( a ) 均为 0。
牛顿公式是泰勒公式的离散模拟,因为它仅使用离散点处的f值,即整数处(移位a )并且仅涉及有限运算:有限差分不涉及极限。
收敛
正如我今天早上在@AlgebraFact上提到的,
将有限级数写为具有有限数量的非零系数的无限级数通常很有用。
这消除了显式跟踪术语数量的需要,并且可以建议下一步做什么。
将牛顿公式写成无限级数使我们不必写下线性插值的一个版本、二次插值的另一个版本、三次插值的另一个版本等等。(当您刚开始接触时,写出这些特殊情况是一个很好的练习。主题,但请记住接下来的无限级数版本。)
至于建议下一步做什么,很自然地会探索如果无限级数确实是无限的,即如果f不是多项式,会发生什么。级数在什么情况下收敛?如果它确实收敛到某个值,那么它是否一定在每个x处收敛到f ( x ) ?
示例f ( x ) = sin(π x ) 表明牛顿定理并不总是成立,因为对于该函数,当a = 0 时,牛顿定理右侧的级数完全相同为零,因为所有 Δ k项为零。但卡尔森定理 [1] 本质上是说,对于沿虚轴增长慢于 sin(π x ) 的整个函数,牛顿定理中的级数收敛于f 。
说一个函数是“完整的”意味着它在整个复平面上是解析的。这意味着上面的泰勒级数必须处处收敛才能使牛顿级数收敛[2]。
相关帖子
[1] 卡尔森没有e 。不要与卡尔森关于傅立叶级数收敛性的定理混淆。
[2] 卡尔森原始定理要求f是完整的。后来的改进表明, f在开右半平面上解析且在闭右半平面上连续就足够了。
后离散泰勒级数首次出现在John D. Cook上。
原文: https://www.johndcook.com/blog/2025/01/14/discrete-taylor-series/