Scribd
Importer
29 vues1 page

TP Interpolation et Intégration en Mécanique

Ce document présente un TP sur l'interpolation et l'intégration dans le cadre d'un cours de mécanique à Sorbonne-Université. Les étudiants doivent réaliser une approximation polynômiale à partir de points de mesure, calculer des valeurs numériques, et tracer des graphiques à l'aide de gnuplot. Des méthodes d'intégration comme les trapèzes et Simpson sont également explorées, avec une analyse de l'erreur en fonction du pas d'espace.

Transféré par

Perez Jerome
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
29 vues1 page

TP Interpolation et Intégration en Mécanique

Ce document présente un TP sur l'interpolation et l'intégration dans le cadre d'un cours de mécanique à Sorbonne-Université. Les étudiants doivent réaliser une approximation polynômiale à partir de points de mesure, calculer des valeurs numériques, et tracer des graphiques à l'aide de gnuplot. Des méthodes d'intégration comme les trapèzes et Simpson sont également explorées, avec une analyse de l'erreur en fonction du pas d'espace.

Transféré par

Perez Jerome
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Sorbonne-Université, Faculté des Sciences et Ingénierie

Licence de Mécanique Année universitaire 2023–2024

LU2ME003 : Méthodes mathématiques et numériques pour la mécanique 1

TP1 - Interpolation / intégration

On se donne n = 3. Les n + 1 points de mesure fj d’une grandeur f (x) sont donnés en n + 1


abscisses xj (0 Æ j Æ n = 3) et notés dans le tableau suivant (cf. TD 2) :
j 0 1 2 3
xj -2 0 4 6
fj 3 5 8 5
L’objectif de ce TP est de réaliser une approximation polynômiale P (z) de f (z) dans le
domaine [≠2, 6] et ensuite de calculer numériquement l’intégrale du polynôme interpolé.

0 Créer un fichier [Link] contenant les mesures sur deux colonnes : xj , fj et


tracer les points sur une figure.
1. Soit ztest œ [≠2, 6] quelconque. Utiliser les éléments du tableau (xj , fj ) donné ci-
dessus pour implémenter le calcul de la valeur du polynôme Q interpolé de Lagrange
R
ÿ n Ÿ (ztest ≠ xj )
n
Ptest pour cette valeur de ztest , donnée par Ptest = fi a b.
i=0 j=0,j”=i (xi ≠ xj )
Pour cela utiliser 2 boucles imbriquées. Affichez à l’écran la valeur numérique obtenue.
2. On souhaite construire le polynôme P point par point. Pour cela, on discrétise le
domaine [≠2, 6] en m sous-intervalles, avec m = 10 : on pose h = (xn ≠ x0 )/m le
pas d’espace, puis on introduit les abscisses zk = x0 + k ú h, pour k = 0,.., m. On
calcule un vecteur P à m + 1 composantes dont chaque composante Pk est la valeur
du polynôme de Lagrange pour l’abscisse zk . Calculer et afficher à l’écran les m + 1
couples de valeurs zk , Pk (utiliser une boucle extérieure pour varier k).
3. Créer un fichier [Link] contenant deux colonnes : zk , Pk . Utiliser gnuplot pour
tracer sur une même figure les points (zk , Pk ) et les quatre points de mesure (ou de
collocation) initiaux.
s6
4. Calculer et afficher à l’écran une valeur approchée Ih de l’intégrale I(x) = ≠2 P (z)dz
par la méthode composite des trapèzes, en utilisant les m + 1 valeurs de P obtenues
à la question 2. Calculer et afficher à l’écran l’erreur E entre la valeur approchée
précédente et la valeur exacte de l’intégrale.
5. Faire de même qu’à la question précédente par la méthode composite de Simpson.
6. ,. Faire varier la valeur de m pour que h = 0.1, 0.01, 0.001, 0.0001 et 0.00001.
Il est possible pour cela de conserver le même programme en modifiant juste le pa-
ramétrage. Créer un fichier [Link] à deux colonnes : h et E(h). Tracer l’erreur
E(h) en fonction du pas d’espace h en échelle log-log, et vérifier l’ordre de l’erreur
pour la méthode composite des trapèzes et celle de Simpson.

Vous aimerez peut-être aussi