numpy计算斐波那契数列

import numpy as np

# 第一种方法
F = np.matrix([[1, 1], [1, 0]])
print(F)
‘‘‘
[[1 1]
 [1 0]]
‘‘‘
print((F ** 7)[0, 0])
# 21
print((F ** 2))
‘‘‘
[[2 1]
 [1 1]]
‘‘‘
print((F ** 3))
‘‘‘
[[3 2]
 [2 1]]
‘‘‘
print((F ** 4))
‘‘‘
[[5 3]
 [3 2]]
‘‘‘

# 第二种方法
# 利用黄金分割公式(比奈公式), 加上取整函数, 计算斐波那契数列
n = np.arange(1, 9)
print(n)
# [1 2 3 4 5 6 7 8]
sqrt5 = np.sqrt(5)  # 取开方
print(sqrt5)
# 2.2360679775
phi = (1 + sqrt5) / 2   # 计算黄金分割率
print(phi)
# 1.61803398875
fibonacci = np.rint((phi**n - (-1/phi)**n)/sqrt5)   # rint函数, 对浮点数取整, 但不改变浮点数类型

print(fibonacci)
# [  1.   1.   2.   3.   5.   8.  13.  21.]