《具体数学》学习笔记: 3.和式(下)

《具体数学》学习笔记: 3.和式(下)

序言:
本文将记录一个全新的数学知识—有限微积分，并将之与我们熟知的无限微积分作比对。

1. 无限微积分
高等数学或是微积分课程中，对于一个函数微分，借助微分算子(operator) $\mathrm{d}$ 表达为:
$$\mathrm{d}f(x)=\underset{h\to 0}{\lim }\frac{f(x+h)-f(x)}{h}$$
算子(operator) $\mathrm{d}$ 的作用是在一个函数上给出新的函数。如果 $f$ 是从实数到实数的一个适当光滑的曲线，那么 $\mathrm{d}f$亦然。
反过来，我们还能够对函数作积分(逆微分)，借助积分算子(operator) $\int$，并且可以与微分进行联系：
$$g(x)=\mathrm{d}f(x)\mathrm{i}\mathrm{f}\mathrm{a}\mathrm{n}\mathrm{d}\mathrm{o}\mathrm{n}\mathrm{l}\mathrm{y}\mathrm{i}\mathrm{f}\int g(x)\mathrm{d}x=f(x)+C$$
这里，$\int g(x)\mathrm{d}x$ 是 $g(x)$ 的不定积分，它是导数等于 $g(x)$ 的一个函数类。另外，我们还有一个重要的微积分基本定理(牛顿-莱布尼兹公式)：
$${\int }_a^{b}g(x)\mathrm{d}x=f(x){|}_a^b=f(b)-f(a)$$
以解决定积分的计算。

2. 有限微积分
对于整数，我们无法直接套用传统的无限微积分，但是我们也希望能够参考这种思路。由此，数学家们发展出了与传统无限微积分类似的“有限微积分”。定义差分算子 $\Delta$：
$$\Delta f(x)=\frac{f(x+h)-f(x)}{h}{|}_{h=1}=f(x+1)-f(x)$$
对于整数而言， $h$ 至小为$1$，因此 $h=1$ 已经是所能达到的最接近的“极限”。在无限微积分中，我们有一个很简洁的关系式：
$$\mathrm{d}(x^m)=m{x}^{m-1}$$
而我们定义的 $\Delta$ 算子，似乎却不行
$$\Delta (x^2)=(x+1{)}^2-x^2=2x+1̸=2(x+1)$$
不过，对于另一种“$m$次幂”，却有着非凡的效果。定义
$$x^{\underline{m}}=\stackrel{m}{\overbrace{x(x-1)\cdots (x-m+1)}},m\ge 0$$
注意 $x$ 的上标 $m$ 中有一个下划线，它表示等式右边这 $m$ 个因子是阶梯般地”向下再向下“。相应地，还有“向上再向上”的定义法
$$x^{\overline{m}}=\stackrel{m}{\overbrace{x(x+1)\cdots (x+m-1)}},m\ge 0$$
当 $m=0$时，我们有 $x^{\underline{0}}=x^{\overline{0}}=1$，因为不含任何因子的乘积定义为1。我们一般称 $x^{\underline{m}}$ 为下降阶乘幂(falling factorial power)，$x^{\frac{}{m}}$为上升阶乘幂(rising factorial power)。特别的，$n!=n^{\underline{n}}=1^{\overline{n}}$。
现在，来尝试对下降阶乘幂求差分：
$$\begin{array}{rl}\Delta x^{\underline{m}}& =(x+1{)}^{\underline{m}}-x^{\underline{m}}\\ & =(x+1)x\cdots (x-m+2)-x\cdots (x-m+2)(x-m+1)\\ & =mx(x-1)\cdots (x-m+2)\\ & =m{x}^{\underline{m-1}}\end{array}$$
恰好地，这正能与无限微积分公式 $\mathrm{d}(x^m)=m{x}^{m-1}$ 相媲美。

继续，我们来考虑 $\Delta$ 的逆运算， 便是我们的求和算子 $\sum$，类似于不定积分，我们可以导出不定和式：
$$g(x)=\Delta f(x)\mathrm{i}\mathrm{f}\mathrm{a}\mathrm{n}\mathrm{d}\mathrm{o}\mathrm{n}\mathrm{l}\mathrm{y}\mathrm{i}\mathrm{f}\sum g(x)\delta x=f(x)+C$$
不定积分中，$C$ 是任意常数，而不定和式中的 $C$ 则满足 $p(x+1)=p(x)$ 的任意一个 $p(x)$。此外，我们还能模仿出类似于定积分的确定和式:
$$\sum _a^{b}g(x)\delta x=f(x){|}_a^b=f(b)-f(a)$$

使用数学归纳法使上式更直观化，并理解其意义：

1. 命 $b=a$，得到
   $$\sum _a^{a}g(x)\delta x=f(a)-f(a)=0$$
2. 命 $b=a+1$，得到
   $$\sum _a^{a+1}g(x)\delta x=f(a+1)-f(a)=g(a)$$
3. 命 $b=b+1$，得到
   $$\begin{array}{rl}\sum _a^{b+1}g(x)\delta x-\sum _a^{b}g(x)\delta x& =(f(b+1)-f(a))+(f(b)-f(a))\\ & =g(b)\end{array}$$

因此，${\sum }_a^{b}g(x)\delta x={\sum }_a^{b-1}g(x)$。与无限积分不同的是，这里我求和时将抛去上限。

至此，我们可以对求和有一个新的认识了。首先，刚才已经证明 $\Delta x^{\underline{m}}=m{x}^{\underline{m-1}}$，那么反过来便有
$$\sum _{0\le k\le n}{k}^{\underline{m}}=\frac{k^{\underline{m+1}}}{m+1}{|}_0^{n+1}=\frac{(n+1{)}^{\underline{m+1}}}{m+1}$$

1. 命 $m=1$，得到
   $$\begin{array}{rl}\sum _{0\le k\le n}{k}^{\underline{1}}& =\sum _{0\le k\le n}k\\ & =\frac{(n+1{)}^{\underline{2}}}{2}\\ & =\frac{(n+1)n}{2}\end{array}$$
   这就是连续自然和公式，进一步地
2. 命 $m=2$，根据 $k^2=k^{\underline{2}}+k^{\underline{1}}$
   $$\begin{array}{rl}\sum _{0\le k\le n}{k}^2& =\sum _{0\le k\le n}{k}^{\underline{2}}+\sum _{0\le k\le n}{k}^{\underline{1}}\\ & =\frac{(n+1{)}^{\underline{3}}}{3}+\frac{(n+1{)}^{\underline{2}}}{2}\\ & =\frac{1}{3}\left[(n+1)n(n-1)+\frac{3}{2}(n+1)n\right]\\ & =\frac{1}{3}(n+1)(n+\frac{1}{2})n\end{array}$$

我们用新的方法再一次求得了连续自然数的平方和公式。此外，连续自然数的 $n$ 阶和也能以此递推。

这里只记录了正指数的下降幂，关于负指数下降幂的相关讨论参考《具体数学(第三版)》的 2.6 节后半部分。