基于复分析的生成函数法求解递推关系的疑难咨询
Hey there! 你提到的直接用高阶导数求生成函数系数的思路虽然理论上可行,但实际操作起来会非常繁琐——对于有理函数形式的生成函数,部分分式分解+几何级数展开才是生成函数法求解递推关系的标准高效路径。咱们一步步拆解这个问题:
直接计算$f(z)$的n阶导数在$z=0$处的值会涉及复杂的莱布尼茨法则展开,每一步的计算量都会随着n增大而急剧增加,完全没必要——生成函数法的核心优势就是通过代数变形把复杂的生成函数转化为我们熟悉的、能直接展开成幂级数的形式,从而快速提取系数。
1. 对生成函数的分母做因式分解
首先解二次方程$z^2 + z - 1 = 0$,得到两个根:
$$r_1 = \frac{-1 + \sqrt{5}}{2}, \quad r_2 = \frac{-1 - \sqrt{5}}{2}$$
因此分母可以写成$(z - r_1)(z - r_2)$,生成函数变为:
$$f(z) = \frac{-z}{(z - r_1)(z - r_2)}$$
2. 执行部分分式分解
我们把$f(z)$拆成两个简单分式的和:
$$f(z) = \frac{A}{z - r_1} + \frac{B}{z - r_2}$$
通过代入根的值求解系数A和B:
- 代入$z = r_1$:$-r_1 = A(r_1 - r_2) \implies A = \frac{-r_1}{r_1 - r_2}$
- 代入$z = r_2$:$-r_2 = B(r_2 - r_1) \implies B = \frac{r_2}{r_1 - r_2}$
计算$r_1 - r_2 = \sqrt{5}$,代入后得到:
$$A = \frac{1 - \sqrt{5}}{2\sqrt{5}}, \quad B = \frac{-1 - \sqrt{5}}{2\sqrt{5}}$$
3. 转化为几何级数的标准形式
为了方便展开成幂级数,我们把分式改写为$\frac{1}{1 - x}$的形式(因为$\sum_{n=0}^\infty x^n = \frac{1}{1 - x}$,当$|x| < 1$时成立):
$$\frac{A}{z - r_1} = \frac{-A}{r_1} \cdot \frac{1}{1 - \frac{z}{r_1}}, \quad \frac{B}{z - r_2} = \frac{-B}{r_2} \cdot \frac{1}{1 - \frac{z}{r_2}}$$
代入A、B、$r_1$、$r_2$的值并化简后,可以得到更简洁的形式:
$$f(z) = \frac{1}{\sqrt{5}} \left( \frac{1}{1 - \phi z} - \frac{1}{1 - \psi z} \right)$$
其中$\phi = \frac{1 + \sqrt{5}}{2}$(黄金分割比),$\psi = \frac{1 - \sqrt{5}}{2}$。
4. 展开幂级数并提取系数
利用几何级数的展开式,把两个分式分别展开:
$$\frac{1}{1 - \phi z} = \sum_{n=0}^\infty (\phi z)^n = \sum_{n=0}^\infty \phi^n z^n$$
$$\frac{1}{1 - \psi z} = \sum_{n=0}^\infty (\psi z)^n = \sum_{n=0}^\infty \psi^n z^n$$
代入回$f(z)$的表达式:
$$f(z) = \frac{1}{\sqrt{5}} \sum_{n=0}^\infty (\phi^n - \psi^n) z^n$$
对比生成函数的定义$f(z) = \sum_{n=0}^\infty c_n z^n$,直接对应系数就得到斐波那契数列的显式公式:
$$c_n = \frac{\phi^n - \psi^n}{\sqrt{5}}$$
如果坚持要用$c_n = \frac{1}{n!} f^{(n)}(0)$,可以利用留数定理或者有理函数的n阶导数公式,但过程会繁琐很多。比如对于有理函数$f(z) = \frac{P(z)}{Q(z)}$,你可以通过将$f(z)$展开为洛朗级数,或者反复使用莱布尼茨法则对分子分母求导,但这完全是绕远路,远不如部分分式分解高效。
内容的提问来源于stack exchange,提问作者Shreyas B.

