{ "cells": [ { "cell_type": "markdown", "id": "3aefd380", "metadata": {}, "source": [ "
\n", "\n", "Ce fichier reflète les énoncés des exercices d'un chapitre de la deuxième édition du livre *[Programmation Python par la pratique](https://www.dunod.com/sciences-techniques/programmation-python-par-pratique-problemes-et-exercices-corriges-0)*. Tous les énoncés peuvent être téléchargés au format [Jupyter Notebook](https://python.guillod.org/livre/) ou exécutés directement en ligne sur [mybinder](https://mybinder.org/v2/gh/guillod/livre-python/HEAD). Les corrigés sont disponibles dans le livre en version papier (ISBN [9782100862122](https://www.dunod.com/sciences-techniques/programmation-python-par-pratique-problemes-et-exercices-corriges-0)) ou numérique (ISBN [9782100875153](https://www.dunod.com/sciences-techniques/programmation-python-par-pratique-problemes-et-exercices-corriges-0)) publiés chez Dunod. Ce fichier reflète les exercices tels que publiés dans cet ouvrage et diffère quelque peu des exercices présentés à la page [python.guillod.org](https://python.guillod.org/)." ] }, { "cell_type": "markdown", "id": "524ec2c5", "metadata": {}, "source": [ "# 4 Représentations graphiques" ] }, { "cell_type": "markdown", "id": "650d8931", "metadata": {}, "source": [ "
\n", "\n", "Représenter graphiquement des résultats issus de calculs numériques ou symboliques est indispensable pour les analyser ou les interpréter. Le module Matplotlib permet de faire des représentations graphiques très variées. La documentation de Matplotlib est disponible [ici](https://matplotlib.org/users/).\n", "\n", "Pour l'utiliser, il est d'usage de l'importer ainsi:" ] }, { "cell_type": "code", "execution_count": 1, "id": "5d5b9cb4", "metadata": { "collapsed": false }, "outputs": [], "source": [ "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "id": "8c154bef", "metadata": {}, "source": [ "et souvant NumPy est également très utile:" ] }, { "cell_type": "code", "execution_count": 2, "id": "52857805", "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np" ] }, { "cell_type": "markdown", "id": "9f135bf3", "metadata": {}, "source": [ "**Concepts abordés:**\n", "\n", "* représentation graphique unidimensionnelle\n", "\n", "* représentations graphiques bidimensionnelle et tridimensionnelle\n", "\n", "* diagramme araignée\n", "\n", "* diagramme de bifurcation\n", "\n", "* optimisation par parallélisation" ] }, { "cell_type": "markdown", "id": "bb938740", "metadata": {}, "source": [ "# Exercice 4.1: Représentations graphiques" ] }, { "cell_type": "markdown", "id": "976216fa", "metadata": {}, "source": [ "Par exemple la fonction `plot` peut être utilisée pour représenter la fonction $x^2$:" ] }, { "cell_type": "code", "execution_count": 3, "id": "48a5b0ad", "metadata": { "collapsed": false }, "outputs": [], "source": [ "x = np.linspace(0,1,50)\n", "y = x**2\n", "plt.plot(x,y)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "1cdd3bb0", "metadata": {}, "source": [ "Afin de définir une jolie figure pouvant être exportée, la syntaxe est la suivante:" ] }, { "cell_type": "code", "execution_count": 4, "id": "7b0a1eec", "metadata": { "collapsed": false }, "outputs": [], "source": [ "plt.figure(figsize=(8,5)) # taille de la figure (en inches)\n", "plt.title(r'Graphique de $x^2$') # titre de la figure (du code LaTeX peut être inclus)\n", "plt.xlabel(r'$x$') # titre de l'axe horizontal\n", "plt.ylabel(r'$y$') # titre de l'axe vertical\n", "plt.plot(x, y, marker='o', label=r\"$x^2$\") # légende\n", "plt.legend() # affiche la légende\n", "plt.savefig(\"test.pdf\") # exporte la figure en PDF\n", "plt.savefig(\"test.png\", dpi=100) # exporte la figure en PNG\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "7a3278d9", "metadata": {}, "source": [ "**a)**\n", "Représenter graphiquement sur la même figure les fonctions $\\sin(kx)$ et $\\cos(kx)$ pour $k=1,2,3$ sur $x\\in[0,2\\pi]$. Faire en sorte que les graduations sur l'axe horizontal soient tous les $\\frac{\\pi}{2}$:\n", "\n", "
\n", "\n", "**Indication.**\n", "Utiliser la fonction `xticks` décrite [ici](https://matplotlib.org/api/_as_gen/matplotlib.pyplot.xticks.html)." ] }, { "cell_type": "markdown", "id": "c7f72dfe", "metadata": {}, "source": [ "**b)**\n", "Regarder l'aide de la fonction `imshow` et l'utiliser pour représenter graphiquement une matrice de nombres aléatoires dans $[0,1]$ de taille $10\\times10$:\n", "\n", "
" ] }, { "cell_type": "markdown", "id": "dc0b1745", "metadata": {}, "source": [ "**c)**\n", "Représenter graphiquement la fonction $f(x,y) = \\frac{-y}{5} + e^{-x^2-y^2}$ pour $x\\in[-3,3]$ et $y\\in[-3,3]$ en densité et avec des courbes de niveau:\n", "\n", "
" ] }, { "cell_type": "markdown", "id": "3609c3dd", "metadata": {}, "source": [ "# Exercice 4.2: Chaos déterministe" ] }, { "cell_type": "markdown", "id": "c518c936", "metadata": {}, "source": [ "Le but est d'étudier un modèle extrêmement élémentaire d'évolution\n", "d'une population. Malgré son caractère simpliste, ce modèle présente\n", "de nombreuses propriétés très intéressantes, en particulier chaotiques.\n", "\n", "La proportion d'une population à un temps discret $i\\in\\mathbb{N}$\n", "est notée $x_{i}\\in[0,1]$. L'évolution de la population est donnée\n", "par $x_{i+1}=f_{\\!\\mu}(x_{i})$ où $f_{\\!\\mu}:[0,1]\\to[0,1]$ est la fonction\n", "définie par:" ] }, { "cell_type": "markdown", "id": "9e02130f", "metadata": {}, "source": [ "$$\n", "f_{\\!\\mu}(x)=\\mu x(1-x) \\,.\n", "$$" ] }, { "cell_type": "markdown", "id": "e64f960c", "metadata": {}, "source": [ "Le paramètre $\\mu\\in[0,4]$ décrit la croissance de la population.\n", "L'application $f_{\\!\\mu}$ est appelée l'application logistique\n", "de paramètre $\\mu$.\n", "\n", "**a)**\n", "Définir la fonction $f_{\\!\\mu}$ en Python comme $\\mathtt{f}(\\mathtt{x},\\mathtt{mu})=f_{\\!\\mu}(x)$ et la représenter graphiquement pour plusieurs valeurs\n", "de $\\mu$." ] }, { "cell_type": "markdown", "id": "db134552", "metadata": {}, "source": [ "**b)**\n", "Pourquoi la valeur de $\\mu$ ne peut-elle pas être supérieure à 4 dans le modèle ?" ] }, { "cell_type": "markdown", "id": "8607c1ed", "metadata": {}, "source": [ "**c)**\n", "Définir une fonction qui prend comme argument une valeur de $\\mu$,\n", "une donnée initiale $x_{0}\\in[0,1]$ et un nombre $n\\in\\mathbb{N}$\n", "et retourne la liste $S_{\\mu}^{n}=(x_{0},x_{1},x_{2},\\dots,x_{n})$.\n", "Modifier cette fonction pour qu'elle puisse prendre un paramètre\n", "optionnel $m\\in\\mathbb{N}$ et retourne la liste $S_{\\mu}^{m,n}=(x_{m},x_{m+1},\\dots,x_{n})$,\n", "c'est-à-dire retire les $m-1$ premiers éléments à $S_{\\mu}^{n}$." ] }, { "cell_type": "markdown", "id": "54e401e0", "metadata": {}, "source": [ "**d)**\n", "Tester cette fonction en faisant la représentation graphique de la\n", "liste $S_{\\mu}^{n}$ pour différentes valeurs des paramètres $x_{0}$\n", "et $\\mu$. Observer les différents comportements de la suite." ] }, { "cell_type": "markdown", "id": "164f3439", "metadata": {}, "source": [ "**Diagramme araignée.**\n", "Une façon d'étudier plus spécifiquement ce qui se passe lorsque $\\mu$\n", "varie est de faire un diagramme dit araignée qui consiste à relier\n", "par des droites les points:" ] }, { "cell_type": "markdown", "id": "b9a93058", "metadata": {}, "source": [ "$$\n", "\\bigl\\{(x_{0},0),(x_{0},x_{1}),(x_{1},x_{1}),(x_{1},x_{2}),(x_{2},x_{2}),\\dots,(x_{n},x_{n}),(x_{n},x_{n+1})\\bigr\\} .\n", "$$" ] }, { "cell_type": "markdown", "id": "391d0695", "metadata": {}, "source": [ "**e)**\n", "Définir une fonction qui retourne la liste des points nécessaire pour\n", "construire le diagramme araignée.\n", "\n", "**Indication.**\n", "Dans le but d'utiliser ensuite Matplotlib, il est parfois\n", "plus simple de représenter une liste de points $\\bigl\\{(x_{0},y_{0}),(x_{1},y_{1}),\\dots,(x_{n},y_{n})\\bigr\\}$\n", "comme le vecteur des abscisses et le vecteur des ordonnées, *i.e.*\n", "$(x_{0},x_{1},\\dots,x_{n})$ et $(y_{0},y_{1},\\dots,y_{n})$." ] }, { "cell_type": "markdown", "id": "a1849925", "metadata": {}, "source": [ "**f)**\n", "Définir une fonction qui trace le graphe de la fonction $f_{\\!\\mu}$, le\n", "graphe de la fonction identité, ainsi que les segments reliant les\n", "points de la liste précédente." ] }, { "cell_type": "markdown", "id": "39457be4", "metadata": {}, "source": [ "**g)**\n", "En expérimentant, étudier qualitativement les effets des paramètres\n", "$\\mu$ et $x_{0}$. Décrire le comportement observé lorsque $\\mu$ augmente." ] }, { "cell_type": "markdown", "id": "36cc5c59", "metadata": {}, "source": [ "**Diagramme de bifurcation.**\n", "Les expérimentations précédentes laissent penser que le comportement\n", "de la suite $x_{i}$ en temps grand (*i.e.* lorsque $i$ est\n", "grand) est indépendant du choix de la condition initiale $x_{0}$\n", "mais dépend beaucoup de la valeur du paramètre $\\mu$. Le but de cette\n", "section est de représenter graphiquement pour chaque valeur de $\\mu$\n", "l'ensemble des points $S_{\\mu}^{m,n}=(x_{m},x_{m+1},\\dots,x_{n})$,\n", "\\emph{i.e.} en mettant $\\mu$ sur l'axe horizontal et toutes les valeurs\n", "de $S_{\\mu}^{m,n}$ sur l'axe vertical.\n", "Un bon choix de paramètres pour nettoyer le diagramme et garder seulement\n", "le comportement en temps long du système est $n=200$ et $m=100$\n", "par exemple.\n", "\n", "**h)**\n", "Définir une fonction qui pour une liste $L$ de valeurs de $\\mu$ donnée,\n", "une donnée initiale $x_{0}\\in[0,1]$ et des entiers $m,n\\in\\mathbb{N}$\n", "retourne la liste des points $\\bigl\\{(\\mu,x):x\\in S_{\\mu}^{m,n}\\;\\text{pour}\\;\\mu\\in L\\bigr\\}$." ] }, { "cell_type": "markdown", "id": "5a12d0d1", "metadata": {}, "source": [ "**i)**\n", "Définir une fonction qui représente graphiquement cette liste de points.\n", "Il est suggéré de prendre pour $L$ une liste de 1 000 valeurs dans\n", "l'intervalle $[0,4]$.\n", "\n", "**Indication.**\n", "Utiliser la fonction `scatter` de\n", "Matplotlib avec les paramètres `s=1` et `edgecolor='none'`." ] }, { "cell_type": "markdown", "id": "eaabc064", "metadata": {}, "source": [ "**j)**\n", "Interpréter le diagramme obtenu, en particulier ce qu'il dit sur le\n", "comportement en temps long du système. Déterminer approximativement\n", "pour quelles valeurs de $\\mu$ le système :\n", "\n", "* possède zéro comme unique point fixe;\n", "\n", "* possède un unique point fixe non nul;\n", "\n", "* oscille entre deux valeurs distinctes (cycle de longueur deux);\n", "\n", "* oscille entre quatre valeurs distinctes (cycle de longueur quatre);\n", "\n", "* oscille entre trois valeurs distinctes (cycle de longueur trois).\n", "\n", "**Représentation de l'attracteur.**\n", "Pour des valeurs de $\\mu$ proches de 4, les valeurs de la population\n", "$x_{i}$ ont l'air d'être plus ou moins aléatoires. Pourtant, le système\n", "est purement déterministe dans le sens où pour une valeur initiale\n", "$x_{0}$ donnée, la population $x_{i}$ est définie sans aléa. Ce comportement\n", "*a priori* aléatoire est appelé chaos déterministe.\n", "Pour bien remarquer que les points $x_{i}$ ne sont pas déterminés de\n", "manière aléatoire, le but est de représenter graphiquement les\n", "points $(x_{n},x_{n+1})$ et de voir que $x_{n+1}$ n'est pas du tout\n", "aléatoire par rapport à $x_{n}$.\n", "\n", "**k)**\n", "Pour chaque valeur de $\\mu$ fixée, définir une fonction qui retourne\n", "la liste de points:" ] }, { "cell_type": "markdown", "id": "b7956391", "metadata": {}, "source": [ "$$\n", "\\bigl\\{(x_{m},x_{m+1}),(x_{m+1},x_{m+2}),\\dots,(x_{n},x_{n+1})\\bigr\\} \\,.\n", "$$" ] }, { "cell_type": "markdown", "id": "3270ae68", "metadata": {}, "source": [ "**l)**\n", "Représenter graphiquement ces points pour différentes valeurs de $\\mu$.\n", "Par exemple $n=5 000$ et $m=100$ est un bon choix de paramètres." ] }, { "cell_type": "markdown", "id": "d1b84094", "metadata": {}, "source": [ "**m)**\n", "Comment serait le graphique précédent si chaque $x_{i}$ était tiré\n", "aléatoirement dans l'intervalle $[0,1]$ indépendamment de $x_{i-1}$\n", "?" ] }, { "cell_type": "markdown", "id": "ed3ccb3c", "metadata": {}, "source": [ "# Exercice 4.3: Ensemble de Mandelbrot" ] }, { "cell_type": "markdown", "id": "51595371", "metadata": {}, "source": [ "L'ensemble de Mandelbrot est défini comme l'ensemble des points $c\\in\\mathbb{C}$ pour lesquels la suite de nombres complexes définis récursivement par $z_0=0$ et" ] }, { "cell_type": "markdown", "id": "d13486e7", "metadata": {}, "source": [ "$$\n", "z_{n+1} = z_n^2+c \\,,\n", "$$" ] }, { "cell_type": "markdown", "id": "500812f7", "metadata": {}, "source": [ "est bornée.\n", "Il est possible de montrer que $c\\in\\mathbb{C}$ est dans l'ensemble de Mandelbrot si et seulement si $|z_n|\\leq2$ pour tout entier $n$.\n", "\n", "**a)**\n", "Écrire une fonction `mandelbrot(c)` qui vérifie si le point $c\\in\\mathbb{C}$ est dans l'ensemble de Mandelbrot de manière approximative en testant les cent premières itérations." ] }, { "cell_type": "markdown", "id": "c14b4a61", "metadata": {}, "source": [ "**b)**\n", "Tester la fonction précédente avec $c=0$ et $c=1+i$. Qu'est-il attendu théoriquement ?\n", "\n", "**Indication.**\n", "Par exemple, le nombre complexe $2+3i$ se définit en Python par `2+3j`." ] }, { "cell_type": "markdown", "id": "79af17ba", "metadata": {}, "source": [ "**c)**\n", "Écrire une fonction `mandelbrot_set(N)` qui génère un tableau de taille $N\\times N$ représentant l'ensemble $c \\in \\big\\{x+iy: x\\in[-2,2] \\; \\text{et} \\; y\\in[-2,2]\\big\\}$ et qui retourne un tableau de booléens de taille $N\\times N$ déterminant si le point associé est dans l'ensemble de Mandelbrot ou non." ] }, { "cell_type": "markdown", "id": "6b24728e", "metadata": {}, "source": [ "**d)**\n", "En utilisant la fonction précédente avec $N=100$, représenter graphiquement avec `imshow` une approximation de l'ensemble des points appartenant à l'ensemble de Mandelbrot." ] }, { "cell_type": "markdown", "id": "500decbe", "metadata": {}, "source": [ "**e)**\n", "En adaptant les fonctions précédentes, représenter graphiquement le logarithme du nombre d'itérations nécessaires avant de ne plus satisfaire $|z_n|\\leq 2$ au lieu d'un booléen:\n", "\n", "
" ] }, { "cell_type": "markdown", "id": "c4e327a4", "metadata": {}, "source": [ "**f)**\n", "!! La méthode précédente a le désavantage de procéder séquentiellement au calcul pour chaque valeur de $c$, ce qui rend l'évaluation assez lente. Proposer une nouvelle implémentation permettant de calculer parallèlement toutes les valeurs en utilisant les indexages NumPy." ] }, { "cell_type": "markdown", "id": "1797f8e0", "metadata": {}, "source": [ "# Exercice 4.4: ! Représentations graphiques avancées" ] }, { "cell_type": "markdown", "id": "49c00b76", "metadata": {}, "source": [ "Le but de cet exercice est de découvrir un panel des possibilités offertes par Matplotlib.\n", "\n", "**a)**\n", "Tracer les lignes de courant du champ de vecteur de l'oscillateur de Van der Pol:" ] }, { "cell_type": "markdown", "id": "793ced19", "metadata": {}, "source": [ "$$\n", "\\begin{pmatrix}y\\\\ \n", "-x+\\mu(1-x^{2})y\n", "\\end{pmatrix}\n", "$$" ] }, { "cell_type": "markdown", "id": "858b093d", "metadata": {}, "source": [ "pour différentes valeurs de $\\mu\\in\\mathbb{R}$." ] }, { "cell_type": "markdown", "id": "7a147cf5", "metadata": {}, "source": [ "**b)**\n", "Représenter graphiquement la courbe paramétrique:" ] }, { "cell_type": "markdown", "id": "e5f9ef98", "metadata": {}, "source": [ "$$\n", "\\begin{pmatrix}\\left(1+t^{2}\\right)\\sin(2\\pi t)\\\\ \n", "\\left(1+t^{2}\\right)\\cos(2\\pi t)\\\\ \n", "t\n", "\\end{pmatrix}\n", "$$" ] }, { "cell_type": "markdown", "id": "20744c60", "metadata": {}, "source": [ "pour $t\\in[-2,2]$." ] }, { "cell_type": "markdown", "id": "2e8607b7", "metadata": {}, "source": [ "**c)**\n", "Représenter la fonction de deux variables:" ] }, { "cell_type": "markdown", "id": "29743d33", "metadata": {}, "source": [ "$$\n", "f(x,y)=\\sin\\left(\\sqrt{x^{2}+y^{2}}\\right)\n", "$$" ] }, { "cell_type": "markdown", "id": "33a07398", "metadata": {}, "source": [ "en trois dimensions pour $x\\in[-5,5]$ et $y\\in[-5,5]$." ] }, { "cell_type": "markdown", "id": "3489e30d", "metadata": {}, "source": [ "**d)**\n", "!! Représenter le ruban de Möbius donné comme surface paramétrique :" ] }, { "cell_type": "markdown", "id": "172e49d9", "metadata": {}, "source": [ "$$\n", "\\begin{pmatrix}\\left(3+v\\cos(u/2)\\right)\\cos u\\\\ \n", "\\left(3+v\\cos(u/2)\\right)\\sin u\\\\ \n", "v\\sin(u/2)\n", "\\end{pmatrix}\n", "$$" ] }, { "cell_type": "markdown", "id": "85fa798c", "metadata": {}, "source": [ "pour $u\\in[0,2\\pi]$ et $v\\in[-1,1]$." ] }, { "cell_type": "markdown", "id": "8e0663fa", "metadata": {}, "source": [ "**e)**\n", "!! Regarder les exemples disponibles [ici](https://matplotlib.org/tutorials/introductory/sample_plots.html) et en choisir deux à comprendre et à modifier." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 5 }