Calcul d'incertitude

Un article de Wikipédia, l'encyclopédie libre.

Le calcul d'incertitude permet d'évaluer correctement les erreurs qui se produisent lors de mesures liées à la vérification d'une relation entre différentes grandeurs physiques. Les instruments de mesure n'étant pas de précision infinie, les mesures faites pendant une expérience ne sont pas exactes. Il faut donc évaluer ces incertitudes pour répondre à la question : « la relation n'est pas vérifiée exactement parce qu'elle est fausse ou parce que les mesures sont incertaines ? » On en déduit des marges d'erreurs, en dehors desquelles la relation sera invalidée. Cela fait partie intégrante de la méthode scientifique.

Sommaire

1 Utilisation des différentielles totales exactes
2 Utilisation de calculatrices

[modifier] Utilisation des différentielles totales exactes

Une loi physique s'exprime par une relation algébrique entre un certain nombre de grandeurs mesurables.

[modifier] Exemples simples : surface et volume

Le calcul de la surface d'un rectangle de côtés L et l :

$S=L \cdot l$

devient lorsque les côtés deviennent L+dL et l+dl:

$S(L+dL,l+dl) = (L+dL) \cdot(l+dl) = L\cdot l + L \cdot dl +l\cdot dL + dl\cdot dL$

Donc la variation de la surface dS peut s'écrire :

$dS = (L+dL)\cdot (l+dl) - L \cdot l = L\cdot dl +l\cdot dL + dL\cdot dl$

que l'on approche par :

$dS = L\cdot dl +l\cdot dL$ car dL.dl est négligeable.

Noter que

$dS = \frac{\partial (L\cdot l) }{\partial L}dL+\frac{\partial (L\cdot l) }{\partial l}dl$

car

$\frac{\partial (L\cdot l) }{\partial L}= l ; \frac{\partial (L\cdot l) }{\partial l}=L$

d'où

$dS = \frac{\partial S(L,l) }{\partial L}dL+\frac{\partial S(L,l) }{\partial l}dl$

De même la variation de volume d'une boîte de côtés x, y, z de volume V=xyz :

$V(x+dx,y+dy,z+dz) = (x+dx)\cdot(y+dy)\cdot(z+dz)$

$= x\cdot y\cdot z +dx\cdot y \cdot z+x\cdot dy\cdot z + x\cdot y\cdot dz + x\cdot dy\cdot dz + y\cdot dx\cdot dz + z\cdot dx\cdot dy + dx\cdot dy\cdot dz$

peut s'écrire

$dV =V(x+dx,y+dy,z+dz) - x\cdot y\cdot z =$

$= dx\cdot y \cdot z+x\cdot dy\cdot z + x\cdot y\cdot dz + x\cdot dy\cdot dz + y\cdot dx\cdot dz + z\cdot dx\cdot dy + dx\cdot dy\cdot dz$

que l'on approche par :

$dV = y\cdot z\cdot dx +z\cdot x\cdot dy + x\cdot y\cdot dz$

Noter que

$dV = yz dx +zx dy + xy dz = \frac{\partial (xyz) }{\partial x}dx+\frac{\partial (xyz) }{\partial y}dy+\frac{\partial (xyz) }{\partial z}dz$

car

$\frac{\partial (xyz) }{\partial x}= yz ; \frac{\partial (xyz) }{\partial y}=xz ; \frac{\partial (xyz) }{\partial z}=xy$

et donc

$dV = \frac{\partial (xyz) }{\partial x}dx+\frac{\partial (xyz) }{\partial y}dy+\frac{\partial (xyz) }{\partial z}dz= \frac{\partial V(x,y,z) }{\partial x}dx+\frac{\partial V(x,y,z) }{\partial y}dy+\frac{\partial V(x,y,z) }{\partial z}dz$

[modifier] La variation d'une fonction f(x,y,z)

Et plus généralement, pour le calcul de la variation d'une fonction f(x,y,z).

$\frac{\partial f(x,y,z) }{\partial x}$ = dérivée partielle par rapport à x

$d f(x,y,z) = \frac{\partial f(x,y,z) }{\partial x}dx+\frac{\partial f(x,y,z) }{\partial y}dy+\frac{\partial f(x,y,z)}{\partial z}dz$

[modifier] Loi des gaz parfaits

Prenons par exemple la loi des gaz parfaits reliant :

P : la pression du gaz
V : le volume occupé par le gaz
n : la quantité de gaz en moles (1 mole = 6,022 10²³ molécules)
R : la constante des gaz parfaits = 8,314 J.K^-1.mol^-1
T : la température absolue du gaz, en kelvin.

$P =\frac{n \times R \times T}{V}$ exprime la pression en fonction de n, R, T et V.

Écrivons sa différentielle :

$dP (T,R,n,V)=P (T+dT,R+dR,n+dn,V+dV)-P (T,R,n,V) = \frac{n \times R}{V} dT + \frac{n \times T}{V} dR + \frac{ R \times T}{V}dn - \frac{n \times R \times T}{V^2}dV$ .

la variation la plus grande s'obtiendra lorsque les 4 termes ci-dessus s'ajouteront :

$\delta P =\frac{n \times R}{V}\delta T + \frac{n \times T}{V}\delta R +\frac{ R \times T}{V}\delta n + \frac{n \times R \times T}{V^2}\delta V$

donne l'erreur absolue sur P déduite du calcul de P à partir de la connaissance des erreurs sur T, R, n et V.

Dans ce cas particulier, on a :

$\frac{dP (T,R,n,V)}{P} = \frac{P (T+dT,R+dR,n+dn,V+dV)-P (T,R,n,V)}{P}$ .

$\frac{dP (T,R,n,V)}{P} = \frac{ \frac{n \times R}{V} dT + \frac{n \times T}{V} dR +\frac{ R \times T}{V} dn - \frac{n \times R \times T}{V^2} dV }{P}= \frac{dT}{T} + \frac{dR}{R} +\frac{dn}{n} - \frac{dV}{V}$ .

et donc dans l'absolu :

$\frac{\delta P}{P} =\frac{\delta T}{T} + \frac{\delta R}{R} +\frac{\delta n}{n} + \frac{\delta V}{V}$ .

On peut aussi utiliser la différentielle logarithmique :

$P =\frac{n \times R \times T}{V}$ .

Donc

$\ln(P) =\ln(n) +\ln(R) +\ln(T) -\ln(V)\,$ .

En dérivant, on obtient :

$\frac{dP}{P} =\frac{dT}{T} + \frac{dR}{R} +\frac{dn}{n} - \frac{dV}{V}$ .

Cette méthode plus rapide s'applique lorsqu'on cherche à faire la différentielle d'une fonction, quotient ou produit de plusieurs variables.

Les incertitudes relatives s'ajoutent lorsque l'on a un produit de variables et ce résultat est remarquable car il est facile à retenir : les incertitudes relatives s'ajoutent lorsque la formule ne comporte que des produit (au sens large : une division est un produit par l'inverse).

[modifier] Utilisation de calculatrices

Ce qui vient d'être fait peut être fait par calcul direct avec une calculatrice ou un tableur (sur ordinateur):

Utilisation de graphes et de barres d'erreurs

Reprenons l'exemple de l'étude des gaz parfaits. Si l'on trace P en fonction de 1/V, on obtiendra théoriquement une droite passant par l'origine, avec comme pente RnT : y=(RnT).x.

n et T étant maintenus constants (l'enceinte ou cellule de mesure contenant le gaz étant sans fuite et thermostatée avec T connu à 0,2%), P étant mesuré, en utilisant un manomètre, avec 5% d'erreur relative, et V étant mesuré avec 2% d'erreur relative, pour chaque point de mesure expérimentale (P,1/V), on trace des barres d'erreurs représentant l'erreur absolue.

Un programme de « fit » ou d'ajustage de courbe, basé sur l'idée de minorer la distance de la droite (ou courbe) à tous les points expérimentaux, permet de tracer la droite théorique et de calculer sa pente nRT avec un coefficient de confiance r² proche de l'unité, si le fit est bon. On utilise la "méthode des moindres carrés" : le programme utilisé somme les distances entre la droite et chaque point, le minimum de cette somme correspondant à la meilleure droite de régression.

Dans le cas de figure ci dessus, on obtient ainsi nRT= 2.54 (1 ± 0.07) Joule

Ceci permet de dire que à n et T constants, l'expérience confirme que PV est constant à 7% près pour le gaz étudié et que pour améliorer ce résultat, il faut mesurer V à mieux que 5%.