斐波那契数列与黄金分割比的关系

前言

众所周知,斐波那契数列的递推公式为:

fi={0i=01i=1fi1+fi2i2f_i = \begin{cases} 0 & i=0 \\ 1 & i = 1 \\ f_{i-1}+f_{i-2} & i \ge 2 \end{cases}

为什么会把斐波那契数列与黄金分割比联系再一起呢?其实和斐波那契的通项公式有关。

斐波那契数列的通项公式

我们考虑使用生成函数推导斐波那契数列的通项公式(别问我为什么,问就是博客正好写了,搬运一下)。

关于生成函数
生成函数是指以某个序列为系数的多项式,其形式为:

F(x)=i=0aixiF(x) = \sum_{i=0} a_i x^i

其中 xx 的存在并没有什么意义,所以又被称为形式幂级函数

我们定义斐波那契数列的生成函数 F(x)F(x) 为:

F(x)=i0fixiF(x) = \sum_{i \ge 0} f_i x^i

我们对这个式子进行变换

F(x)=i0fixi=f0+xf1+i2fixi\begin{aligned} F(x) &= \sum_{i \ge 0} f_i x^i \\ &= f_0 + x f_1 + \sum_{i \ge 2} f_i x^i \\ \end{aligned}

i2i \ge 2 时,fi=fi1+fi2f_i = f_{i-1} + f_{i - 2},所以:

F(x)=f0+xf1+i2(fi1+fi2)xi=f0+xf1+i2fi1xi+i2fi2xi=f0+xf1+i2fi1xi+i2fi2xi=f0+xf1+xi2fi1xi1+x2i2fi2xi2=f0+xf1+xi1fixi+x2i1fixi\begin{aligned} F(x) &= f_0 + x f_1 + \sum_{i \ge 2} (f_{i-1} + f_{i-2}) x^i \\ &= f_0 + x f_1 + \sum_{i \ge 2} f_{i-1} x^i + \sum_{i \ge 2} f_{i-2} x^i \\ &= f_0 + x f_1 + \sum_{i \ge 2} f_{i-1} x^i + \sum_{i \ge 2} f_{i-2} x^i \\ &= f_0 + x f_1 + x \sum_{i \ge 2} f_{i-1} x^{i-1} + x^2 \sum_{i \ge 2} f_{i-2} x^{i-2} \\ &= f_0 + x f_1 + x \sum_{i \ge 1} f_{i} x^{i} + x^2 \sum_{i \ge 1} f_{i} x^i \end{aligned}

然后由 F(x)=i0fixi\displaystyle F(x) = \sum_{i \ge 0} f_i x^i 且当 i=0i=0 时,f0x0=0f_0 x^0 = 0 可得:

F(x)=f0+xf1+xi1fixi+x2i1fxiF(x)=f0+xf1+xF(x)+x2F(x) f0=0,f1=1 F(x)=x+xF(x)+x2F(x)\begin{aligned} & F(x) = f_0 + x f_1 + x \sum_{i \ge 1} f_{i} x^i + x^2 \sum_{i \ge 1} f_{} x^i \\ & F(x) = f_0 + x f_1 + x F(x) + x^2 F(x) \\ \because \ & f_0 = 0, f_1 = 1 \\ \therefore \ & F(x) = x + x F(x) + x^2 F(x) \\ \end{aligned}

合并同类项,系数化为一,得:

(1xx2)F(x)=x F(x)=x1xx2\begin{aligned} & (1 - x - x^2) F(x) = x \\ \therefore \ & F(x) = \frac{x}{1 - x - x^2} \end{aligned}

这样,我们就得到了 F(x)F(x)封闭形式

对于这个式子,我们考虑使用待定系数法将其部分分式分解:

A1ax+B1bx=x1xx2AAbx+BaBx1(a+b)x+abx2=x1xx2(A+B)(Ab+aB)x1(a+b)x+abx2=x1xx2\begin{aligned} \frac{A}{1-ax} + \frac{B}{1-bx} &= \frac{x}{1 - x - x^2} \\ \\ \frac{A - Abx + B - aBx}{1 - (a + b)x + abx^2} &= \frac{x}{1 - x - x^2} \\ \\ \frac{(A+B) - (Ab + aB)x}{1 - (a + b)x + abx^2} &= \frac{x}{1 - x - x^2} \\ \end{aligned}

观察两边的系数,得到一个方程组:

{A+B=0Ab+aB=1a+b=1ab=1{A=15B=15a=1+52b=152\begin{cases} A + B = 0 \\ Ab + aB = -1 \\ a + b = 1 \\ ab = -1 \\ \end{cases} \Rightarrow \begin{cases} A = \frac{1}{\sqrt{5}} \\ B = - \frac{1}{\sqrt{5}} \\ a = \frac{1 + \sqrt{5}}{2} \\ b = \frac{1 - \sqrt{5}}{2} \\ \end{cases}

接下来我们要介绍一个前置知识——几何级数(等比数列求和)。

关于几何级数
几何级数是指形式为

n=0rn\sum_{n=0}^{\infty} r^n

的级数。其中 rr 称为公比(其实就是等比数列)。r<1|r| < 1 时,级数收敛,其和为:

n=1rn=11r\sum_{n=1}^{\infty} r^n = \frac{1}{1 - r}

这个就是我们接下来展开要用到的公式。

我们运用几何级数将上面最初的表达式展开:

A1ax=Ai=0(ax)i=Ai=0aixiB1bx=Bi=0(bx)i=Bi=0bixi\begin{aligned} & \frac{A}{1-ax} = A\sum_{i=0}^{\infty} (ax)^i = A\sum_{i=0}^{\infty} a^i x^i \\ & \frac{B}{1-bx} = B\sum_{i=0}^{\infty} (bx)^i = B\sum_{i=0}^{\infty} b^i x^i \\ \end{aligned}

收敛性证明

用几何级数展开公式时有一个前提:r<1|r| < 1,而上面两个式子的 rr 分别为 axaxbxbx
其中 a=1+52,b=152a = \frac{1 + \sqrt{5}}{2}, b = \frac{1 - \sqrt{5}}{2},所以:

ax=1+52xbx=152x\begin{aligned} & ax = \frac{1 + \sqrt{5}}{2} x \\ & bx = \frac{1 - \sqrt{5}}{2} x \end{aligned}

对于 a,b,xa,b,x 存在限制:

{1+52x<1152x<1{1+52x<1152x<1\begin{cases} |\frac{1 + \sqrt{5}}{2} x| < 1 \\ \\ |\frac{1 - \sqrt{5}}{2} x| < 1 \end{cases} \Rightarrow \begin{cases} |\frac{1 + \sqrt{5}}{2}| |x| < 1 \\ \\ |\frac{1 - \sqrt{5}}{2}| |x| < 1 \end{cases}

所以

{x<512x<5+12\begin{cases} |x| < |\frac{\sqrt{5} - 1}{2}| \\ \\ |x| < |-\frac{\sqrt{5} + 1}{2}| \end{cases}

x<min(512,5+12)=512|x| < \min\Big(|\frac{\sqrt{5} - 1}{2}|, |-\frac{\sqrt{5} + 1}{2}|\Big) = \frac{\sqrt{5} - 1}{2}

因为 F(x)F(x) 定义为形式幂级数,所以我们仅关注 xix^i 的系数,而不关注 xx 的具体的值,只要由满足条件的 xx 存在,就可以用上面的方式对原式展开。

用上面两个式子等量代换,得:

F(x)=A(i=0aixi)+B(i=0bixi)=15(i=0(1+52)ixi)15(i=0(152)xi)=15(i=0(1+52)ixii=0(152)ixi)=15i=0((1+52)ixi(152)ixi)=15i=0xi((1+52)i(152)i)=i=0(1+52)i(152)i5xi\begin{aligned} F(x) &= A \Bigg(\sum_{i=0}^{\infty} a^i x^i \Bigg) + B \Bigg( \sum_{i=0}^{\infty} b^i x^i \Bigg) \\ &= \textcolor{red}{\frac{1}{\sqrt{5}}} \Bigg( \sum_{i=0}^{\infty} \bigg( \frac{1 + \sqrt{5}}{2} \bigg)^i x^i \Bigg) - \textcolor{red}{\frac{1}{\sqrt{5}}} \Bigg( \sum_{i=0}^{\infty} \bigg( \frac{1 - \sqrt{5}}{2} \bigg) x^i \Bigg) \\ &= \frac{1}{\sqrt{5}} \Bigg( \sum_{\textcolor{blue}{i=0}}^{\infty} \bigg( \frac{1 + \sqrt{5}}{2} \bigg)^i x^i - \sum_{\textcolor{blue}{i=0}}^{\infty} \bigg( \frac{1 - \sqrt{5}}{2} \bigg)^i x^i \Bigg) \\ &= \frac{1}{\sqrt{5}} \sum_{i=0}^{\infty} \bigg( \Big(\frac{1 + \sqrt{5}}{2} \Big)^i x^i - \Big( \frac{1-\sqrt{5}}{2} \Big)^i x^i \bigg)\\ &= \frac{1}{\sqrt{5}} \sum_{i=0}^{\infty} \textcolor{green}{x^i} \bigg( \Big(\frac{1 + \sqrt{5}}{2} \Big)^i - \Big( \frac{1-\sqrt{5}}{2} \Big)^i \bigg) \\ &= \sum_{i=0}^{\infty} \frac{\big(\frac{1 + \sqrt{5}}{2} \big)^i - \big( \frac{1-\sqrt{5}}{2} \big)^i}{\sqrt{5}} x^i \\ \end{aligned}

对比 F(x)F(x) 最初的的定义:

F(x)=i=0fixiF(x) = \sum_{i=0}^{\infty} f_i x^i

得斐波那契数列的通项公式为:

fi=(1+52)i(152)i5f_i = \frac{\big(\frac{1 + \sqrt{5}}{2} \big)^i - \big(\frac{1 - \sqrt{5}}{2} \big)^i}{\sqrt{5}}

其中就有黄金分割比 ϕ=512\phi = \frac{\sqrt{5} - 1}{2} 的相反数

更直观的体现

我们已经发现了斐波那契数列的通项公式中的黄金分割比,其实黄金分割比在斐波那契数列中还有更直观的体现。

我们取斐波那契数列相邻两项的比值:

fifi+1=(1+52)i(152)i(1+52)i+1(152)i+1=21+5(1+52)i+1215(152)i+1(1+52)i+1(152)i+1=512(1+52)i+1(5125)(152)i+1(1+52)i+1(152)i+1=512(1+52)i+1512(152)i+15(152)i+1(1+52)i+1(152)i+1=5125(152)i+1(1+52)i+1(152)i+1=5125(3+52)i+11\begin{aligned} \frac{f_i}{f_{i+1}} &= \frac{\big ( \frac{1+\sqrt{5}}{2} \big)^i - \big ( \frac{1 - \sqrt{5}}{2} \big)^i}{\big ( \frac{1+\sqrt{5}}{2} \big)^{i+1} - \big ( \frac{1 - \sqrt{5}}{2} \big)^{i+1}} \\ &= \frac{ \frac{2}{1+\sqrt{5}} \big( \frac{1 + \sqrt{5}}{2} \big)^{\textcolor{red}{i+1}} - \frac{2}{1-\sqrt{5}} \big( \frac{1 - \sqrt{5}}{2} \big)^{\textcolor{red}{i+1}} }{ \big(\frac{1+\sqrt{5}}{2} \big)^{i+1} - \big ( \frac{1 - \sqrt{5}}{2} \big)^{i+1}} \\ &= \frac{ \frac{\sqrt{5} - 1}{2} \big( \frac{1 + \sqrt{5}}{2} \big)^{i+1} - \big(\frac{\sqrt{5} - 1}{2} - \sqrt{5}\big) \big( \frac{1 - \sqrt{5}}{2} \big)^{i+1} }{ \big(\frac{1+\sqrt{5}}{2} \big)^{i+1} - \big ( \frac{1 - \sqrt{5}}{2} \big)^{i+1}} \\ &= \frac{ \frac{\sqrt{5} - 1}{2} \textcolor{red}{\big( \frac{1 + \sqrt{5}}{2} \big)^{i+1}} - \frac{\sqrt{5} - 1}{2} \textcolor{red}{\big( \frac{1 - \sqrt{5}}{2} \big)^{i+1}} - \sqrt{5}\big( \frac{1 - \sqrt{5}}{2} \big)^{i+1}}{ \textcolor{red}{\big(\frac{1+\sqrt{5}}{2} \big)^{i+1}} - \textcolor{red}{\big ( \frac{1 - \sqrt{5}}{2} \big)^{i+1}}} \\ &= \frac{\sqrt{5} - 1}{2} - \frac{\sqrt{5}\big( \frac{1 - \sqrt{5}}{2} \big)^{i+1}}{\big(\frac{1+\sqrt{5}}{2} \big)^{i+1} - \big ( \frac{1 - \sqrt{5}}{2} \big)^{i+1} } \\ &= \frac{\sqrt{5} - 1}{2} - \textcolor{red}{\frac{\sqrt{5}}{ \big( - \frac{3 + \sqrt{5}}{2} \big)^{i+1} - 1 }} \\ \end{aligned}

ii \to \infty 时,后面的分式的值趋向于 00,即:

limififi+1=512\lim_{i \to \infty} \frac{f_i}{f_{i+1}} = \frac{\sqrt{5} - 1}{2}

所以,斐波那契数列相邻两项的比值趋向于黄金分割比 ϕ\phi