**问题:**如何求斐波那契数列的通项?
notes:斐波那契数列:后面一项是前面两项的求和。
eg:1,1,2,3,5,8,13,21。。。
**解答:**我们将斐波那契数列的问题写成矩阵形式,即:
[
x
k
+
2
x
k
+
1
]
=
[
1
1
1
0
]
∗
[
x
k
+
1
x
k
]
\begin{bmatrix} x_{k+2} \\ x_{k+1} \end{bmatrix}= \begin{bmatrix} 1&1 \\ 1 & 0 \end{bmatrix} *\begin{bmatrix} x_{k+1} \\ x_{k} \end{bmatrix}
[xk+2xk+1]=[1110]∗[xk+1xk]
前一天的例子说明了,我们要求一个通向的长期稳态解,需要考察最大特征值对应的特征向量。这里,我们对上述矩阵求特征值和特征向量,得到:
l
1
=
1.61803399
l
2
=
−
0.61803399
l_1=1.61803399 ~~~~~~~~l_2=-0.61803399
l1=1.61803399 l2=−0.61803399
l
1
和
l
2
l_1 和l_2
l1和l2表示矩阵的特征值,而通解的形式可以写作:
[
x
k
+
2
x
k
+
1
]
=
c
1
∗
l
1
k
∗
v
1
+
c
2
∗
l
2
k
∗
v
2
\begin{bmatrix} x_{k+2} \\ x_{k+1} \end{bmatrix}=c_1*l_1^k*v_1+c_2*l_2^k*v_2
[xk+2xk+1]=c1∗l1k∗v1+c2∗l2k∗v2
其中
v
1
和
v
2
v_1和v_2
v1和v2是他们各自的特征向量,
c
1
和
c
2
c_1和c_2
c1和c2是他们的系数。而最终由最大的特征值即1.618这一项决定,通过计算可得
c
1
=
1
/
5
c_1=1/\sqrt{5}
c1=1/5。
答案:
于是,我们直接计算 1.61 8 k / 5 1.618^k/\sqrt{5} 1.618k/5,来得到斐波那契数列的通解。
以下是,我用电脑算出的前二十项:
0.7236067977499789
1.1708203932499368
1.8944271909999157
3.065247584249853
4.959674775249769
8.024922359499623
12.984597134749393
21.009519494249016
33.99411662899841
55.00363612324743
88.99775275224584
144.0013888754933
232.9991416277391
377.0005305032324
609.9996721309716
987.0002026342041
1596.9998747651757
2584.00007739938
4180.999952164556