Dérivation et intégration : Méthode de base
Contents
La page ci-présente existe en version notebook téléchargeable grâce au bouton (choisir le format
.ipynb
). On rappelle qu’l faut ensuite l’enregistrer dans un répertoire adéquat sur votre ordinateur (capa_num
par exemple dans votre répertoire personnel) puis lancer Jupyter Notebook depuis Anaconda pour accéder au notebook, le modifier et exécutez les cellules de code adéquates.
2.3. Dérivation et intégration : Méthode de base#
On se propose de mettre en place les calculs d’intégration et de dérivation sur des cas simples basés sur l’électrocinétique. On considère un circuit RC soumis à un échelon de tension
avec
Pour la suite, on prendra
L’intensité qui circule à travers le condensateur (et donc aussi la résistance) est donné par
La puissance dissipée à un instant
par effet Joule est donné par . L’énergie dissipée par effet Joule entre et est donc:
Dans ce cas, simple, on peut calculer analytiquement
2.3.1. Calcul de l’intensité.#
Exercice A:
Importer les bibliothèques scientifiques
numpy
etmatplotlib.pyplot
.Ecrire deux instructions pour créer deux vecteurs
tk
etuk
contenant respectivement les instantséquirépartis entre 0 et 5s et les valeurs à ces instants. On prendra 1000 valeurs pour chaque vecteur. Ecrire une boucle avec son initialisation qui calcule les
par un calcul dérivée à gauche avec un pas puis les ajoute à une liste ik
.Transformer la liste
ik
en un vecteurnumpy
.Tracer sur un même graphique
obtenu par dérivation numérique et obtenu par calcul analytique puis jouer sur le pas pour voir si la dérivation est encore acceptable si on augmente .
On rappelle que analytiquement, on trouve:
Exercice B:
Reprendre les instructions précédentes et réaliser le calcul de la dérivation numérique par intégration à droite puis centrée. Comparer graphiquement à la solution analytique.
2.3.2. Calcul de l’énergie dissipée#
Exercice C:
Ecrire une instruction pour créer un vecteur
Pk
contenant les puissancesdissipées par effet Joule dans R à chaque instant. Ecrire une boucle avec son initialisation qui calcule les
, énergie dissipée entre et , par la méthode des rectangles à gauche avec un pas puis les ajoute à une liste Ek
. Pensez à utiliser une récurrence.Transformer la liste
Ek
en un vecteurnumpy
.Tracer sur un même graphique
obtenu par dérivation numérique et obtenu par calcul analytique puis jouer sur le pas pour voir si la dérivation est encore acceptable si on augmente .
On rappelle que analytiquement, on trouve:
Exercice D:
Reprendre les instructions précédentes et réaliser le calcul de l’intégration numérique par intégration à droite, centrée puis par les trapèzes. Comparer graphiquement à la solution analytique.
2.3.3. A partir de données expérimentales.#
Les données expérimentales sont celle déjà utilisée précédemment et accessible par ce lien. On rappelle que la méthode d’importation d’un fichier est disponible ici.
Note
Si votre navigateur ouvre le fichier au lieu de le télécharger, faites un clic droit sur le lien et choisir Enregistrer la cible du lien
.
Exercice E:
Après avoir télécharger le fichier et l’avoir placé dans le bon répertoire, ouvrez le pour lire l’en-tête et la structure du fichier. Ecrire alors deux instructions : la première importe le fichier et crée deux vecteurs
tk
etukb
contenant les temps de mesure et les mesures de tension en bits et la seconde instruction qui transforme la listeukb
en une listeuk
contenant les valeurs de tension en volt. S’aider de la ligne de commande proposée ci-dessous et de la description du fichier.Ecrire une boucle avec son initialisation qui calcule les
par dérivation numérique (cf. remarque après) et qui le stocke dans une liste ik
puis transformer la liste en un vecteur.Ecrire une boucle avec son initialisation qui calcule les
, énergie dissipée entre et par intégration numérique et qui le stocke dans une liste Ek
puis transformer la liste en un vecteur. On utilisera la méthode des rectangles (cf. remarque après).Tracer trois graphiques:
sur trois graphiques différents (on utilisera soit trois fenêtres graphiques, soit un graphique multiple si on se sent à l’aise avec ce dernier).
Pour l’import, on rappelle la commande permettant de séparer les colonne en plusieurs variables:
temps, ucb = np.loadtxt('chemin_vers_fichier', skiprows=nlignes_a_sauter, delimiter="delimiteur de colonne", unpack=True)
Pour les méthodes numériques, on ne peut choisir ici le pas d’intégration car on ne connaît que des valeurs discrètes.
Dérivation numérique: le pas est imposé par les écarts entre les
. De plus, comme précisé précédemment, vous devez utiliser une dérivée à droite en et une dérivée à gauche en . Pour les instants intermédiaires, il est demandé d’utiliser une dérivée centrée.Intégration numérique: le pas est à nouveau imposé par les
. On utilisera la méthode des trapèzes pour l’intégration.