Suite de Fibonacci#

Leonardo Fibonacci (circa 1175 – 1250), mathématicien italien, proposa comme un exercice mathématique amusant le modèle suivant pour décrire la croissance naturelle d’une population de lapins. On suppose qu’on dispose d’un couple de lapins tout juste nés, maintenus isolés de leurs congénères. On suppose de plus que chaque mois, chaque couple âgé d’au moins deux mois engendre un nouveau couple de lapins. On ne tient pas compte de la mortalité naturelle des lapins. Le nombre \(C_n\) de couples de lapins suit par conséquent une relation de récurrence,

\[ C_0 = C_1 = 1 \ \ \ , \ \ \ C_n = C_{n-1}+C_{n-2} \]

La suite d’entiers \(C_n\) est connue comme la suite de Fibonacci. On commence par inclure les paquets numpy et matplotlib

Hide code cell source 
import numpy as np
import matplotlib.pyplot as plt

  1. Créez un tableau ns de 50 entiers de 0 à 49 inclus.

Hide code cell source 
N = 50
ns = np.arange(50)

  1. En utilisant une boucle for, la valeur des deux premiers termes \(C_0\) et \(C_1\), et la relation de récurrence, calculez les 50 premières valeurs de cette suite, et stockez-les dans un tableau numpy C.

Hide code cell source 
C = np.empty(N, dtype = 'int64')
C[0] = 1
C[1] = 1
for i in range(2, N):
    C[i] = C[i - 1] + C[i - 2]

  1. Représentez l’évolution du nombre de couples de lapins en fonction du temps (exprimé en mois), pendant 50 mois. Pensez à ajouter un titre, à indiquer les axes et les unités.

Hide code cell source 
plt.plot(ns, C, 'b+-')
plt.title("Évolution du nombre de couples de lapins")
plt.xlabel("Temps [mois]")
plt.ylabel("Couples de lapins")
plt.show()
../../../_images/f7ff7fd0f4d22715c1a0106a3a948aa2106c3dcc0961c90cb06362259487119a.png

  1. Ajoutez dans votre programme le calcul de la suite des rapports \(r_n\) entre deux termes successifs de la suite de Fibonacci (pour \(n\) allant de 1 à 49 inclus) : $\( r_n = \frac{C_{n + 1}}{C_n} \)$

Hide code cell source 
# Avec boucle :
r = np.empty((N - 1,))
for i in range(N - 1):
    r[i] = C[i + 1] / C[i]

# Ou sans boucle :
r = np.array(C[1:], dtype = 'float') / C[:-1] # Si on veut être sûr des conversions de type

# Ou plus simplement :
r = C[1:] / C[:-1]

  1. Tracez le graphe de ce rapport \(r_n\) en fonction du mois \(n\), pour \(n\) allant de 1 à 49 inclus.

Hide code cell source 
plt.plot(ns[1:], r, 'b+-')
plt.title("Suite de Fibonacci : rapport entre termes successifs")
plt.xlabel("Temps $n$ [mois]")
plt.ylabel("$r_n = C_n/C_{n-1}$")
plt.show()
../../../_images/c21891d021b6d5f446fcb7ae45f6f0ead072c11e1c2035a63abd6ada600ac867.png

On peut montrer que le rapport \(r_n\) a pour limite le nombre d’or \(\phi = \frac{1}{2} \left(1+\sqrt{5}\right)\). Pour des \(n\) suffisamment grands, la suite \(C_n\) tend donc vers une suite géométrique de la forme \(D_n = A \phi^n\).

  1. Tracez les 50 premiers termes de la suite de Fibonacci sur un graphe semi-logarithmique (échelle log selon les ordonnées). La suite de Fibonacci tend-elle bien vers une suite géométrique ?

Hide code cell source 
plt.semilogy(ns, C, 'b+-')
plt.title("Suite de Fibonacci")
plt.xlabel("Temps [mois]")
plt.ylabel("Couples de lapins")
plt.show()
../../../_images/8bbd4deb0e9e7846d60b316ba6611541e4199272d573bd41d7f866c8e2af9f19.png

Comme \(\displaystyle C_n\underset{n\to\infty}{\sim} D_n\) et \(\displaystyle \lim_{n\to\infty}C_n = +\infty\), on a : $\( \log(C_n) \underset{n\to\infty}{\sim} \log(D_n) = \log A + n \log \phi \)$

  1. Calculez la suite \(D_n\) et stockez-la dans un tableau numpy.

Hide code cell source 
D = np.log(C)