斐波那契数列

最新推荐文章于 2022-10-23 17:20:08 发布

zixiaqian

最新推荐文章于 2022-10-23 17:20:08 发布

阅读量2.3k

点赞数

文章标签：生物算法编程

转http://zh.wikipedia.org/wiki/%E6%96%90%E6%B3%A2%E9%82%A3%E5%A5%91%E6%95%B0%E5%88%97

表达式

为求得斐波那契数列的一般表达式，可以借助线性代数的方法。高中的初等数学知识也能求出。

高中的初等数学知识解法

已知

$a 1 = 1$
$a 2 = 1$
$a n = a n - 1 + a n - 2$

1首先构建等比数列

设 $a n + α a n - 1 = β(a n - 1 + α a n - 2)$
化简得
$a n = (β - α) a n - 1 + αβ a n - 2$
比较系数可得:
$/begin{cases} /beta-/alpha=1 // /alpha/beta=1 /end{cases}$
不妨设 $β > 0α > 0$
解得：

$/begin{cases} /alpha=/dfrac{/sqrt{5}-1}{2} // /beta=/dfrac{/sqrt{5}+1}{2} /end{cases}$
所以有 $a n + α a n - 1 = β(a n - 1 + α a n - 2)$ 即{ $a n + α a n - 1$ }为等比数列。

2求出数列{ $a n + α a n - 1$ }

有以上可得：
$/begin{align}a_{n+1} +/alpha a_{n} &=(a_2+/alpha a_1)/beta^{n-1}// & = (1+/alpha)/beta^{n-1}// & =/beta^n // /end{align}$

变形得： $/frac{a_{n+1}}{/beta^{n+1}}+/frac{/alpha}{/beta}/frac{a_{n}}{/beta^{n}}=/frac{1}{/beta}$ 令 $b_n=/frac{a_n}{/beta^n}$

3求数列{ $b n$ }进而得到{ $a n$ }

$b_{n+1}+/frac{/alpha}{/beta}b_{n}=/frac{1}{/beta}$
设 $b_{n+1}+/lambda=-/frac{/alpha}{/beta} (b_{n}+/lambda)$ 解得 $/lambda=-/frac{1}{/alpha+/beta}$ 故数列 $b n + λ$ 为等比数列
即 $b_n+/lambda=/left(-/frac{/alpha}{/beta}/right)^{n-1}/left(b_1+/lambda/right)$ 而 $b_1=/frac{a_1}{/beta}=/frac{1}{/beta}$ 故有 $b_n+/lambda=/left(-/frac{/alpha}{/beta}/right)^{n-1}/left(/frac{1}{/beta}+/lambda/right)$
又有 $/begin{cases} /alpha=/dfrac{/sqrt{5}-1}{2} // /beta=/dfrac{/sqrt{5}+1}{2} /end{cases}$ 和 $b_n=/frac{a_n}{/beta^n}$
可得 $a_{n}=/frac{/sqrt{5}}{5} /cdot /left[/left(/frac{1 + /sqrt{5}}{2}/right)^{n} - /left(/frac{1 - /sqrt{5}}{2}/right)^{n}/right]$

得出 $a n$ 表达式

$a_{n}=/frac{/sqrt{5}}{5} /cdot /left[/left(/frac{1 + /sqrt{5}}{2}/right)^{n} - /left(/frac{1 - /sqrt{5}}{2}/right)^{n}/right]$

线性代数解法

1 首先构建一个矩阵方程

设J_n为第n个月新出生的兔子数量，A_n为这一月份的兔子数量。

${J_{n+1}/choose A_{n+1}} = /begin{pmatrix}0&1//1&1/end{pmatrix} /cdot {J_n/choose A_{n}},$

上式表达了两个月之间，兔子数目之间的关系。而要求的是，A_n+1的表达式。

2 求矩阵的特征值: $λ$

行列式：- $λ$ *(1- $λ$ )-1*1= $λ$ ²- $λ$ -1

当行列式的值为0，解得 $λ 1$ = $/frac{1}{2} (1 + /sqrt{5})$ 或 $λ 2$ = $/frac{1}{2} (1 - /sqrt{5})$

3特征向量

将两个特征值代入

$/begin{pmatrix}0&1//1&1/end{pmatrix}-/lambda /cdot E) /cdot/vec x = 0$

求特征向量 $/vec x$ 得

$/vec x_1$ = $/begin{pmatrix} 1///frac{1}{2} (1 + /sqrt{5})/end{pmatrix}$

$/vec x_2$ = $/begin{pmatrix} 1///frac{1}{2} (1 - /sqrt{5})/end{pmatrix}$

4 分解首向量

第一个月的情况是兔子一对，新生0对。

${J_{1}/choose A_{1}} = /begin{pmatrix}0//1/end{pmatrix}$

将它分解为用特征向量表示。

$/begin{pmatrix}0//1/end{pmatrix}=/frac{1}{/sqrt{5}} /cdot /begin{pmatrix}1///frac{1}{2} (1 + /sqrt{5})/end{pmatrix}-/frac{1}{/sqrt{5}} /cdot /begin{pmatrix}1///frac{1}{2} (1 - /sqrt{5})/end{pmatrix}$ （4）

5用数学归纳法证明

从

${J_{n+1}/choose A_{n+1}} = /begin{pmatrix}0&1//1&1/end{pmatrix} /cdot {J_n/choose A_{n}}$ = $/lambda /cdot {J_n/choose A_{n}}$

可得

${J_{n+1}/choose A_{n+1}} = /begin{pmatrix}0&1//1&1/end{pmatrix}^n /cdot {J_{1}/choose A_{1}} =/lambda^n /cdot {J_{1}/choose A_{1}}$ （5）

6 化简矩阵方程

将（4）代入（5）

${J_{n+1}/choose A_{n+1}} = /lambda^n /cdot [/frac{1}{/sqrt{5}} /cdot /begin{pmatrix}1///frac{1}{2} (1 + /sqrt{5})/end{pmatrix}-/frac{1}{/sqrt{5}} /cdot /begin{pmatrix}1///frac{1}{2} (1 - /sqrt{5})/end{pmatrix}]$

根据 3

${J_{n+1}/choose A_{n+1}} = /frac{1}{/sqrt{5}} /cdot /lambda_1^n /cdot /begin{pmatrix}1///frac{1}{2} (1 + /sqrt{5})/end{pmatrix}- /frac{1}{/sqrt{5}} /cdot /lambda_2^n/cdot /begin{pmatrix}1///frac{1}{2} (1 - /sqrt{5})/end{pmatrix}$

7 求A的表达式

现在在6的基础上，可以很快求出A_n+1 的表达式，将两个特征值代入 6 中

$A_{n+1}=/frac{1}{/sqrt{5}} /cdot /lambda_1^{n+1} - /frac{1}{/sqrt{5}} /cdot /lambda_2^{n+1}$

$A_{n+1}=/frac{1}{/sqrt{5}} /cdot (/lambda_1^{n+1} - /lambda_2^{n+1})$

$A_{n+1}=/frac{1}{/sqrt{5}} /cdot ((/frac{1}{2} (1 + /sqrt{5}))^{n+1} - (/frac{1}{2} (1 - /sqrt{5}))^{n+1})$ (7)

(7)即为A_n+1 的表达式

近似值

$F_n /approx /frac{1}{/sqrt{5}} a^n = /frac{1}{/sqrt{5}} /cdot ( /frac{1}{2} (1 + /sqrt{5}) )^n /approx 0.223606798 /cdot 3.236067977^n$

用计算机求解

可通过编程观察斐波那契数列。分为两类问题，一种已知数列中的某一项，求序数。第二种是已知序数，求该项的值。

可通过递归的算法解决此两个问题。

和黄金分割的关系

开普勒发现两个斐波那契数的比会趋近黄金分割：

$/frac {f_{n+1}}{f_n} /approx a = /frac{1}{2} (1 + /sqrt{5}) = /Phi /approx 1{,}618{...}$

斐波那契数亦可以用连分数来表示：

$/frac{1}{1} = 1 /qquad /frac{2}{1} = 1+/frac{1}{1} /qquad /frac{3}{2} = 1+/frac{1}{1+ /frac{1}{1}} /qquad /frac{5}{3} = 1+/frac{1}{1+ /frac{1}{1+ /frac{1}{1}}} /qquad /frac{8}{5} = 1+/frac{1}{1+ /frac{1}{1+ /frac{1}{1+ /frac{1}{1}}}}$

$F_n = /frac{1}{/sqrt{5}} /left/{ /left( /frac{1+/sqrt{5}}{2} /right)^n - /left( /frac{1-/sqrt{5}}{2} /right)^n /right/} = {/phi^n /over /sqrt{5}} - {(1-/phi)^n /over /sqrt{5}}$

而黄金分割数亦可以用无限连分数表示：

$/Phi = 1+/frac{1}{1+ /frac{1}{1+ /frac{1}{1+ /frac{1}{1+ ...}}}}$

和自然的关系

许多的生物构成都和斐波那契数列有正相关。例如人体从肚脐至头顶之距离和从肚脐至脚底之距趋近于 $/boldsymbol{/lim_{n /to /infty}}/frac{/digamma_n}{/digamma_(n-1)}$ 向日葵的种子螺旋排列99%是 $/digamma_n$ 。

恒等式

证明以下的恒等式有很多方法。以下会用组合论述来证明。 $F n$ 可以表示成用多个1和多个2相加令其和等于<mat 不失一般性，我们假设n ≥ 1。 $F n + 1$ 是计算了将1和2加到n的方法的数目。若第一个被加数是1，有 $F n$ 种方法来完成对n-1的计算；若第一个被加数是2，有F(n-1)来完成对n-2的计算。因此，共有 $F n + F n - 1$ 种方法来计算n的值。