Fibonacci数列的一般形式——二阶常系数齐次差分方程的解法

原创已于 2022-12-28 23:44:53 修改 · 2.1k 阅读

6 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#算法 #差分方程 #数列通项

于 2022-09-28 01:18:45 首次发布

理论学习同时被 2 个专栏收录

46 篇文章

订阅专栏

数学基础

6 篇文章

订阅专栏

本文介绍了一种解决二阶线性差分方程的方法，通过构造等比数列，推导出方程的通项公式，并给出了Fibonacci数列的特例解析。

1. 初等解法

对于如下二阶线性差分方程：
$b\neq0\tag{1-1}$
我们希望构造出含 $f (n)$ 和 $f (n - 1)$ 的等比数列。
用待定系数法，两边同时减 $k f (n - 1)$ :
$f(n)-kf(n-1)=(a-k)f(n-1)+bf(n-2)\tag{1-2}$
令 $g (n) = f (n) - k f (n - 1)$ ,我们希望构造出
$\begin{align*} g(n)&=(a-k)g(n-1)\\ &=(a-k)(f(n-1)-kf(n-2))\\ \tag{1-3} \end{align*}$ ,则有 $- k (a - k) = b$ ,进而有
$k^2-ak-b=0\tag{1-4}$
该方程称为特征方程。
根据求根公式，得到两个解 $k1=a+a2+4b2,k2=a−a2+4b2k_1=\frac{a+\sqrt{a^2+4b}}{2},k_2=\frac{a-\sqrt{a^2+4b}}{2}$

根据韦达定理：
$\begin{align*} k_1+k_2&=a\\ k_1k_2&=-b\\ \tag{1-5} \end{align*}$
令 $k=k_1$ ，则有
$\begin{align*} g(n)&=f(n)-k_1f(n-1)\\ &=(a-k_1)g(n-1)\\ &=k_2^{n-2}g(2) \tag{1-6} \end{align*}$
也即有
$f(n)=k_1f(n-1)+k_2^{n-2}g(2)\tag{1-7}$
对上式进行展开：
$\begin{align*} f(n)&=k_1f(n-1)+k_2^{n-2}g(2)\\ &=k_1(k_1f(n-2)+k_2^{n-3}g(2))+k_2^{n-2}g(2)\\ &=k_1^2f(n-2)+k_1k_2^{n-3}g(2)+k_2^{n-2}g(2)\\ &=k_1^{n-1}f(1)+k_1^{n-2}k_2^{0}g(2)+k_1^{n-3}k_2^{1}g(2)+\cdots+k_1^{0}k_2^{n-2}g(2)\\ \tag{1-8} \end{align*}$
若 $k_1=k_2=k$ ，则
$\begin{align*} f(n)&=k^{n-1}f(1)+g(2)(n-1)k^{n-2}\\ &=(\frac{f(1)}{k}-\frac{g(2)}{k^2}+\frac{g(2)}{k^2}n)k^n \tag{1-9} \end{align*}$
令：
$\begin{align*} c_1&=\frac{f(1)}{k}-\frac{g(2)}{k^2}\\ c_2&=\frac{g(2)}{k^2}\\ \tag{1-10} \end{align*}$
则：
$f(n)=c_1k_1^n+c_2k_2^n\tag{1-11}$

若 $k1≠k2k_1\neq k_2$ ，则
$\begin{align*} f(n)&=f(1)k_1^{n-1}+g(2)\frac{k_1^{n-1}-k_2^{n-1}}{k_1-k_2}\\ &=\frac{1}{k_1}(f(1)+\frac{g(2)}{k_1-k_2})k_1^n+\frac{1}{k_2}(-\frac{g(2)}{k_1-k_2})k_2^n\\ \tag{1-12} \end{align*}$
令
$\begin{align*} c_1&=\frac{1}{k_1}(f(1)+\frac{g(2)}{k_1-k_2})\\ c_2&=\frac{1}{k_2}(-\frac{g(2)}{k_1-k_2})\\ \tag{1-13} \end{align*}$
则
$f(n)=c_1k_1^n+c_2k_2^n\tag{1-14}$

1.1. 特例之Fibonacci数列

当 $,b=1f(1)=1,\ f(2)=1,\ a=1,\ ,b=1$ 时，就是大家熟知的Fibonacci数列的递推式。此时特征方程为:
$k^2-k-1=0\tag{1-15}$
两个解分别为： $k2=1−52k_1=\frac{1+\sqrt{5}}{2},\ k_2=\frac{1-\sqrt{5}}{2}$ , 辅助变量 $g (2) = 1 - k$
代入(1-13)和(1-14)得通项公式为：
$f(n)=\frac{1}{\sqrt{5}}\left(\frac{1+\sqrt{5}}{2}\right)^n-\frac{1}{\sqrt{5}}\left(\frac{1-\sqrt{5}}{2}\right)^n\tag{1-16}$

2. 高等解法

（待补充）