Valentine Nadeau
Kym N’diaye
MAE S6 E1
TP n°2 Métrologie Qualité
Simulation de Monte Carlo
Table des matières
I. Introduction .................................................................................................................................... 1
II. Mesures .......................................................................................................................................... 1
III. Monte-Carlo : codage en python.................................................................................................... 4
IV. Conclusion .................................................................................................................................. 6
I. Introduction
Le but de ce TP est d’effectuer une simulation de la méthode Monte-Carlo. Cela sera fait à travers le
calcul de la constante gravitationnelle g. Nous commencerons par déterminer cette constante de
manière expérimentale, grâce à un pendule oscillant. Puis, nous mettrons en place un programme
python afin de traiter nos mesures et de mettre en place la méthode Monte-Carlo
La méthode Monte Carlo est une technique permettant de simuler des distributions de probabilités.
En métrologie, cette méthode est utilisée pour estimer des incertitudes.
𝐿
Rappelons que la période d’une oscillation correspond à 𝑇 = 2𝜋√ où :
𝑔
- T correspond à la période (s) ;
- L à la longueur (m) ;
- g est la constante gravitationnelle terrestre, g = 9,81 N/kg (ou m/s²)
II. Mesures
Nous cherchons à déterminer la constante gravitationnelle de manière expérimentale. Pour ce faire,
nous sommes munis de :
- Un écrou
- Un fil de pêche
- Une règle
- Un rapporteur d’angle
- Une potence avec une pince
- Un chronomètre
1
Valentine Nadeau
Kym N’diaye
MAE S6 E1
Figure 1 : schéma de notre montage
Protocole de mesure :
1. Le boulon est déjà fixé au fil de pêche, qui lui est attaché à la potence avec un pince ;
2. Mesurer la longueur du fil ;
3. Avec le rapporteur d’angle, déplacer l’écrou de 20° maximum ;
4. Sans appliquer de force, lâcher l’écrou ;
5. Démarrer le chronomètre et compter 10 périodes d’oscillation de l’écrou, noter le temps
mesuré ;
6. Répéter cette mesure 3 fois ;
7. Répéter les étapes 2 à 6 pour 10 longueurs en tout.
On obtient donc 30 mesures (présentées dans le tableau 1 ci-dessous) de périodes pour différentes
longueurs. La formule définit en introduction nous permet de calculer la valeur de g :
𝐿
𝑔=
𝑇 2
(2𝜋)
Cette formule nous permet de calculer gexpérimental. La moyenne de nos 30 mesure est de 9,96 N/kg ce
qui présente un écart relatif de 1,5% par rapport à la valeur théorique.
2
Valentine Nadeau
Kym N’diaye
MAE S6 E1
nombre de
L(m) t(s) T(s) T^2 (s) T/2PI g (N/kg)
périodes
0,1 10 6,12 0,612 0,37 0,10 10,54
1 0,1 10 6,17 0,617 0,38 0,10 10,37
0,1 10 6,18 0,618 0,38 0,10 10,34
0,25 10 9,95 0,995 0,99 0,16 9,97
2 0,25 10 9,84 0,984 0,97 0,16 10,19
0,25 10 9,98 0,998 1,00 0,16 9,91
0,355 10 11,74 1,174 1,38 0,19 10,17
3 0,355 10 11,86 1,186 1,41 0,19 9,96
0,355 10 11,39 1,139 1,30 0,18 10,80
0,415 10 12,77 1,277 1,63 0,20 10,05
4 0,415 10 13,06 1,306 1,71 0,21 9,61
0,415 10 12,99 1,299 1,69 0,21 9,71
0,55 10 14,57 1,457 2,12 0,23 10,23
5 0,55 10 14,9 1,49 2,22 0,24 9,78
0,55 10 14,96 1,496 2,24 0,24 9,70
0,625 10 15,83 1,583 2,51 0,25 9,85
6 0,625 10 15,95 1,595 2,54 0,25 9,70
0,625 10 15,73 1,573 2,47 0,25 9,97
0,81 10 18,2 1,82 3,31 0,29 9,65
7 0,81 10 18,23 1,823 3,32 0,29 9,62
0,81 10 17,98 1,798 3,23 0,29 9,89
0,85 10 18,45 1,845 3,40 0,29 9,86
8 0,85 10 18,5 1,85 3,42 0,29 9,80
0,85 10 18,25 1,825 3,33 0,29 10,08
1 10 19,94 1,994 3,98 0,32 9,93
9 1 10 20,39 2,039 4,16 0,32 9,50
1 10 19,72 1,972 3,89 0,31 10,15
1,315 10 23,06 2,306 5,32 0,37 9,76
10 1,315 10 23,09 2,309 5,33 0,37 9,74
1,315 10 22,96 2,296 5,27 0,37 9,85
Tableau 1 : mesures d'oscillation à longueurs variables
𝐿 4𝜋2 𝑇2 4𝜋2
Puisque, 𝑇 2 = 4𝜋 2 𝑔 → 𝑔
𝐿 → 𝐿
= 𝑔
on peut tracer l’évolution du carré de la période en
fonction de la longueur, la pente de cette courbe nous permettra de calculer g.
3
Valentine Nadeau
Kym N’diaye
MAE S6 E1
𝑇2 4𝜋2 4𝜋2
On a 𝐿
= 4,0574 − 0,035 = 𝑔
, donc, 𝑔 = 4,0574−0,035
= 9,8146 N/kg
Ce résultats présente un écart relatif de 0,08% par rapport à la valeur théorique de 9,80665 N/Kg.
III. Monte-Carlo : codage en python
Nous avons pour objectif de mettre en place un programme python permettant de tracer l’évolution
de la période d’oscillation en fonction de la longueur du fils de pêche. Nous souhaitons ensuite générer
un nuage de points autour de chaque mesure selon la méthode Monte-Carlo, cela nous permettra de
déterminer une valeur de g plus précise.
Notre code est présenté ci-dessous, nous nous sommes aidés du site StackOverflow ainsi que des
indications de l’enseignante pour le compléter.
Nous commençons par définir nos listes de mesures en x (longueur) et en y (période d’oscillation). Ces
valeurs correspondent à nos mesures acquises dans la partie précédente du TP.
4
Valentine Nadeau
Kym N’diaye
MAE S6 E1
Puis, nous générons 50 points répartis de manière aléatoire autour de chaque point de mesure (à ±0,5).
Puis, nous déterminons l’équation de la droite f(x)=y. Le coefficient directeur et l’ordonnée à l’origine
sont utilisés pour calculer g (comme dans la partie précédente).
Chaque valeur de g calculée est ensuite moyennée. Finalement, on trace le graphique et la valeur
moyenne de g est affichée en sortie.
Figure 2 : Graphique tracé par le code python
5
Valentine Nadeau
Kym N’diaye
MAE S6 E1
Le résultat obtenu de 9,78341 N/kg présente un écart relatif de 0,2% par rapport à la valeur théorique
de la constante gravitationnelle. Ce résultat est moins précis que ce que nous avons obtenu
expérimentalement, cela ne devrait pas être le cas.
IV. Conclusion
Ce TP nous a permis de déterminer la constante de gravitation de manière expérimentale ainsi
qu’à travers la méthode de Monte-Carlo. Bien que nos résultats aient été moins concluants à l’aide de
cette méthode statistique, nous avons pu comprendre en quoi elle est plus efficace pour estimer des
incertitudes suivant une distribution gaussienne.
