Propagation des erreurs - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs est disponible ici.

Une mesure est toujours entachée d'erreur. Lorsqu'une valeur mesurée est utilisée dans une formule, il faut savoir estimer l'erreur induite sur le résultat de la formule. On parle de propagation d'erreur.

Approches pragmatiques

Report des extrêmes dans le calcul

La première solution consiste à effectuer les calculs avec les extrêmes de l'intervalle d'erreur. Si la mesure a pour valeur

a ± Δa

alors la " valeur réelle " est supposée être dans l'intervalle [aa;aa]. On calcule donc ici

y1 = ƒ(aa)
y2 = ƒ(aa)

et, selon l'ordre de y1 et de y2, on prend [y1;y2] ou [y2;y1] comme intervalle d'erreur.

Cette méthode n'est valable que si la loi est monotone (c'est-à-dire croissante ou décroissante) sur l'intervalle [aa;aa].

Estimation à partir de la dérivée

Une manière simple, utilisée fréquemment en physique, consiste à utiliser un développement limité du premier ordre, c'est-à-dire à remplacer la loi ƒ par sa " tangente " pour estimer l'erreur. On estime ainsi l'erreur avec une loi uniforme (linéaire) et simple.

On a  :

ƒ(x) = ƒ(a) + ƒ '(a)·(x-a) + o(x)

o(x) est une fonction qui " tend vite " vers 0. Si l'on remplace x par a + Δa, on a alors

ƒ(a + Δa) = ƒ(a) + ƒ '(a)·Δa + o(a + Δa)

On peut donc estimer

Δy ≈ ƒ '(a) · Δa

Cette erreur est sous-estimée si l'on a une loi convexe.

Approche mathématique

Notations

par exemple, si x est un vecteur (x1, x2,…, xn), alors

\partial_i y(\underline{x}) = \frac{\partial y}{\partial x_i}(\underline{x})

Formules

Une fonction de variables aléatoires

y = y (\underline{x})

est elle-même une variable aléatoire. Si les erreurs sont petites, la variance du développement limité de y au premier ordre autour des valeurs moyennes μ des x est une bonne estimation de la variance de y :

y(x) = y(\mu) + \sum \left(x_i - \mu_i \right) \partial_i y

On néglige les termes d'ordre supérieur dans l'expansion, il vient :

V(y(x)) = \sum_i\sum_j \partial_i y \partial_j y V_{ij}(x)

Si les x sont indépendantes

V(y(x)) = \sum_i \left( \partial_i y \right)^2 V_{ii}(x)

Applications

Mesure d'une résistance

Une application pratique est la mesure expérimentale d'une résistance R à partir de la chute de tension U entre ses bornes et du courant I. La résistance est décrite par la loi d'Ohm.

R = \frac{U}{I}

Nous avons

\langle R \rangle = \frac{\langle U \rangle}{\langle I \rangle}

\partial_U R = \frac{1}{\langle I \rangle} = \frac{\langle R \rangle}{\langle U \rangle} et \partial_I R = - \frac{\langle U \rangle}{\langle I \rangle^2} = - \frac{\langle R \rangle}{\langle I \rangle}

Il vient

\frac{V_R}{\langle R \rangle^2} = \frac{V_U}{\langle U \rangle^2} + \frac{V_I}{\langle I \rangle^2}

Dans ce cas simple, l'erreur relative sur R correspond à la moyenne géométrique des erreurs relatives sur U et I.

Utilisation des différentielles totales exactes

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

La pression en fonction de n, R, T et V s'exprime par

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

écrivons sa différentielle :

dP =\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

Si l'on " remplace " des variations élémentaires de variables dx par les erreurs sur les variables δx, on obtient :

\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

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

Autres exemples simples :

S = Ll et S + dS = (L + dL)(l + dl) = Ll + Ldl + ldL + dldL
peut s'écrire
dS = ((L + dL)(l + dl) − Ll) = Ldl + ldL + dLdl
que l'on approxime par
dS = Ldl + ldL
  • le calcul d'un volume V = x·y·z
V(x + dx,y + dy,z + dz) = (x + dx)(y + dy)(z + dz) = xyz + dxyz + xdyz + xydz + xdydz + ydxdz + zdxdy + dxdydz
peut s'écrire
dV = yzdx + zxdy + xydz + dxdydz
que l'on approxime par dV = yzdx + zxdy + xydz
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
rappel: \frac{\partial (xyz) }{\partial x}= yz ; \frac{\partial (xyz) }{\partial y}=xz ; \frac{\partial (xyz) }{\partial z}=xy
  • et plus généralement pour le calcul de la variation d'une fonction ƒ(x,y,z).
si \frac{\partial f(x,y,z) }{\partial x} est la 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
Page générée en 0.112 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales
Version anglaise | Version allemande | Version espagnole | Version portugaise