TD 2

L’objectif de ce TP est d’approfondir la récursivité et les calculs de complexité simple pour certaines fonctions.

Rappel Rendre le TP

Le TD doit être rendu sur Moodle avant le début du TD suivant. En cas de TD non rendu, vous aurez automatiquement 0, sauf cas exceptionnel.

Le rendu doit être un dossier compressé au format zip contenant:

• Un fichier contenant les réponses
• Un fichier par exercice du TD

Chaque fonction présente dans le TD doit contenir une documentation et un nom clair.

Les TP peuvent être rendus en binômes et le nom du binôme doit apparaître clairement dans le fichier réponses. Si une réponse est faite à plus que deux, vous devez l’indiquer dans le fichier réponse. Si ce n’est pas indiqué, cela sera considéré comme du plagiat.

Si vous utilisez du code sur internet, vous devez indiquer la source dans le code source et dans le fichier réponse. Si ce n’est pas indiqué, cela sera considéré comme du plagiat.

Exercice 1

Un étudiant de L2 MIASHS un rendu un TD en oubliant de respecter les consignes! Il a très mal nommé ses fonctions et ne leur à pas mis de documentations.

Afin qu’il n’ait pas 0, renommez les fonctions et écrivez pour lui la documentation.

Exercice 2

La suite de Fibonacci \((f_n)\) est définie par récurrence

\[ \begin{cases} f_0 = f_1 = 1 \\ f_{n+2} = f_{n+1} + f_n \end{cases} \]

1. Écrire une fonction fibonacci récursive qui prend en argument un entier \(n\) et retourne \(f_n\), le \(n\)ème terme de la suite de Fibonacci.

2. Analyser la complexité de la fonction.

3. Écrire une fonction récursive de complexité linéaire qui prend en argument un entier \(n\) et retourne le couple \((f_n, f_{n-1})\). En déduire une fonction fibonacci2 qui prend en argument un entier \(n\) et retourne \(f_n\).

4. (Difficile) Mémoisation

Écrire un décorateur (une fonction de fonction, voir ici) qui prend en argument:

• Une fonction \(f\) qui prend en argument un entier
• Un dictionnaire \(D\)

et qui retourne une nouvelle fonction \(nf\) qui prend en argument un entier et:

• si cet entier est une clef de \(D\) retourne la valeur associé
• sinon exécute la fonction \(f\) et stocke le résultat dans \(D\) et le retourne.

Décorez la fonction fibonacci de la première question. Quelle la complexité de la fonction décorée? Quelle est la complexité en espace?

Exercice 3:

Voici la spécification du problème decoupe

1. Écrire une version itérative de la fonction decoupe

2. Écrire une version récursive de la fonction decoupe

3. Faire l’analyse de complexité de chaque implémentation

Exercice 4: Les ensembles de Mandelbrot

Pour \(c\in\mathbb{C}\), on définit la suite complexe

\[ \begin{cases} z_0 = 0\\ z_{n+1} = z_n^2 + c \end{cases} \]

Cette suite peut converger ou diverger en fonction de la valeur de \(c\). On sait qu’elle diverge quand son module est supérieur à 2.

1. Écrire une fonction récursive mandelbrot qui prend en argument le nombre complexe \(c\) et \(n\) et retourne \(z_n\) ou \(2\) si le module de \(z_n\) est supérieur à \(2\). Écrire une fonction mandelbrot_module qui prends en argument un nombre complexe \(c\) et un entier \(n\) et retourne le module de \(z_n\) si celui ci est plus petit que \(2\) et \(2\) sinon.

2. Écrire une fonction fenetre qui prend en argument un nombre complexe \(O\) (l’origine), un float \(D\) (le diametre), et un entier \(n\) (nombre de points) et retourne la liste de liste \(L\) tel que pour tout \(a,b< n\) on a \(L[a][b] = O + (a*\frac{D}{n} - \frac{D}{2}) + (b*\frac{D}{n} - \frac{D}{2})j\)

3. À l’aide d’une ou plusieurs application de map, construisez une fonction qui prends une liste de liste de complexe et applique la fonction Mandelbrot sur chaque point.

4. Pour dessiner une matrice en Python, on peut utiliser le module matplotlib et numpy comme suit:

import matplotlib.pyplot as plt
import numpy as np

L = [
[0, 0, 0, 0, 0, 0, 0],
[0, 1, 1, 1, 1, 1, 0],
[0, 1, 0, 1, 0, 1, 0],
[0, 1, 1, 1, 1, 1, 0],
[0, 0, 1, 0, 1, 0, 0],
[0, 0, 1, 1, 1, 0, 0],
[0, 0, 0, 1, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 0]
]

M = np.array(L) # Tranforme la liste de liste en tableau numpy
plt.imshow(M)
plt.axis("off") # enlève les axes, inutile ici
plt.show() # dessine le resultat

Les couleurs peuvent être changées (voir ici)

Dessiner les ensembles de Mandelbrot:

• centrés en \(0\), de diamètre 3 avec un nombre d’itération de 20

• centrés en -0.06783611264225832 + 0.6617460391250546j de diamètre 0.001, en faisant varier l’itération entre 33 et 160, on obtient les images suivantes (assemblées en un gif):

$ convert -delay 10 -loop 0 *.png mandelbrot.gif

/bin/sh: 1: $: not found

5. (Bonus) Ajoutez des jolies images de fractales (voire des images animés) et le code qui les génèrent. Les plus jolies images seront ajoutés pour la postérité au sujet du TD (avec un remerciement).

Compiled the: dim. 07 janv. 2024 23:19:18 CET