Incertitude due linterpolation sur

une table de donnes numriques

FLORIAN PLATEL
www.metgen.org

1. Formulation mathmatique du problme

Soient x1 et x2 deux rels. On note f, la fonction
traduisant le phnomne physique que lon tudie.
Celle-ci est suppose continue sur [x1, x2]. Cette
fonction nest connue que de manire discrte aux
points x1 et x2 par le tableau 1.

y2 + I(y2)

P2
y2 - I(y2)

y1 + I(y1)
x

f(x)

x1

y1

x2

y2

P1
y1 - I(y1)

Tableau 1
Table dtalonnage donnant f

x1 - I(x1) x1 + I(x1)

x2 - I(x2) x2 + I(x2)

Fig. 1. - Exemple de fonction passant par P1 et P2

Dans le cas le plus gnral, lon a construit cette table

partir de plusieurs expriences : on a gnr plusieurs
valeurs de xi , i {1, 2} et lon a mesur les valeurs de
yi correspondantes. Ainsi, xi (resp. yi) est affecte de
lincertitude I(xi) (resp. I(yi)).

2. Interpolation linaire
Soient la droite passant par les points (x1, y1) et
(x2, y2) et x [x1, x2]. On suppose dautre part que x
est affecte de lincertitude I(x).

Posons :
Posons :
Pi = [xi - I(xi), xi + I(xi)] x [yi - I(yi), yi + I(yi)] ,
avec i {1, 2}.
En notant
lensemble des fonctions continues
passant par les pavs P1 et P2 , il vient alors,
f (fig. 1).

x =

x x1
.
x 2 x1

(1)

Dans ces conditions , la valeur interpole est :

y = x (y2 - y1) + y1 .

INCERTITUDE DUE A LINTERPOLATION SUR UNE TABLE DE DONNEES NUMERIQUES

(2)

2.1. Calcul de lincertitude type sur y

Puisque x [0, 1] , on a

En supposant les variables de lquation (2) non

corrles, lincertitude type sur y est calcule par la loi
de propagation des variances :

donc : u 2 ( x )

y
u ( y ) =
x
2

u 2 ( x ) +

yi
i =1

u 2 ( x )

Px2

ce qui conduit :
u 2 ( y ) = ( y 2 y1 ) 2 u 2 ( x )
+ (1 x ) u ( y1 ) +
2

2x

u ( y2 )
2

(3)

u 2 ( x) +
x i

i =1
2

u ( y)

u 2 ( x i ) ,

1
( x 2 x1 )

u m2 ( y i )

u ( y)

(6)

Py
+
P
x

u m2 ( x i ) + u 2 ( x) .

(7)

u m2 ( y i )

Py
+
Px

u m2 ( x i ) + u 2 ( x) ,

avec :
um(xi) = Max {u(x
{ u(x
u(x;2)}u(x
; 2) } ;
1) 1;)
um(yi) = Max {u(y
{ u(y
)
;
)
u(y
;

)}
u(y
; 2) } ;
1 1
2
Px = x2 - x1 ;
Py = y2 - y1 .

u 2 ( x)
+ (1 x ) 2 u 2 ( x1 )

Rsultat :n 1 :
Rsum
Lincertitude type lie linterpolation linaire
est majore de la faon suivante :
2

ce qui conduit :
u 2 ( x ) =

um2 ( xi ) + u 2 ( x) .

2.4. Expression finale du rsultat

En supposant les variables de lquation (1) non

corrles, lincertitude type sur y est calcule par la loi
de propagation des variances :

u ( x ) = x
x

}.

En combinant (5) et (6), il vient :

2.2. Calcul de lincertitude type sur x

um2 ( xi ) + u 2 ( x)

En posant Px = x2 - x1 , le pas selon x , lexpression

prcdente scrit :

u 2 ( y i ) ,

1
( x2 x1 )2

(1 x ) 2 + 2x 1 , et

(4)

+ 2x u 2 ( x 2 )

2.3. Simplification des expressions

3. Incertitude due lassimilation une droite

2.3.1. Simplification de u(y)

Lassimilation de la fonction f la droite introduit

une erreur. La fonction ntant connue que de manire
discrte, il nest en principe pas possible de dterminer
directement lerreur. Nanmoins, on peut quand mme
la majorer.

En posant : um(yi) = Max {u(y1) ; u(y2)} ,

il vient daprs lexpression (3) :

u 2 ( y ) (1 x ) 2 + 2x u m2 ( y i ) + ( y 2 y1 ) 2 u 2 ( x )

3.1. Existence dun majorant

Puisque x [0, 1] , il vient (1 x ) +

2

donc : u ( y )
2

u m2 ( y i )

2x

1 , et

+ ( y 2 y1 ) u ( x ) .
2

En posant Py = y2 - y1 , le pas selon y , lexpression

prcdente scrit :
u 2 ( y ) u m2 ( y i ) + Py2 u 2 ( x ) .
(5)
2.3.2. Simplification de u(x)
En posant : um(xi) = Max {u(x1) ; u(x2)} ,
il vient daprs lexpression (4) :
1
(1 x ) 2 + 2x u m2 ( x i )
u 2 ( x )
( x 2 x1 ) 2

{(

+ u 2 ( x)

METGEN DOSSIER METROLOGIE N 1

Soit lquation de la droite . Pour x [x1, x2] ,

on dfinit (x) = (x) - f(x) qui est la fonction donnant
lerreur commise en assimilant f la droite .
Par hypothse, la fonction tant continue sur
lintervalle [ x1, x2] qui est un compact de
, on
dduit que ([x1, x2]) est aussi un compact de : ce
qui signifie quil existe c [x1, x2] tel que
x [x1 , x 2 ] ( x) (c) .

3.2. Calcul dun majorant

En conservant les notations prcdentes, lon
obtient ( c) = 0 . Posons (x) = K(x - x1) + f(x1) ,

DM1 Rv. 5

f ( x 2 ) f ( x1 )
. Il vient alors f( c) = K .
x 2 x1
Soit : x [x1, x2]. Par hypothse, f est de classe C2 sur
[x1, x2] et drivable lordre 3 sur ] x1, x2[ . La formule
de Lagrange applique entre x et x1 conduit dire quil
existe ]0, 1[ tel que :
avec

K=

x [x1 , x 2 ]

.
( x1 x) 2
f " ( x1 + ( x1 x))
2

M=

On a alors successivement :
( x) = K ( x x1 ) + f ( x1 ) f ( x1 ) ( x x1 ) f ' ( x)
+

Px2

f"
8 [Sup
x,x ]

( x)

Rsultat n 2 :
La valeur absolue de lerreur commise en
assimilant f la droite sur lintervalle [ x1, x2]
est majore par :

f ( x1 ) = f ( x) + ( x1 x) f ' ( x)
+

Donc :

Px2

f"
8 [Sup
x,x ]
1

( x x1 )
f " ( x1 + ( x1 x))
2
2

Remarque : la principale difficult est maintenant de majorer | f | .

Dans le cas le plus gnral, on ne peut quvaluer f partir de trois
couples dobservations ( xi, yi).

puis :

( x) = ( x x1 ) K f ' ( x)
+

( x x1 )
f " ( x1 + ( x1 x))
2

.
3.3. Calcul dun majorant de |f|

Cette formule applique x = c conduit dire quil

existe 1 ]0, 1[ tel que :

Remarque : pour un intervalle suffisamment petit , cette

hypothse semble lgitime.

(c x1 ) 2
(c ) =
f " ( x1 + 1 ( x1 c)) ,
2

Soient p 2[X] et E = {(xi, yi) , i {1 ; 2 ; 3}} le

graphe de la fonction f. Pour x , on note :
p(x) = a0 +a1x + a2x2 . Le triplet (a0, a1, a2) vrifie le
systme dquations :

et donc :
(c )

(c x1 ) 2
Sup f "
2
[x , x ]
1

Comme (c) est un majorant de sur [x1, x2] :

x [x1 , x 2 ]

( x )

(c x1 ) 2
Sup f "
2
[x , x ]
1

(8)

( x)

(c x 2 )
Sup f "
2
[x , x ]
2

y1 = a 0 + a1 x1 + a 2 x12

2
y 2 = a 0 + a1 x 2 + a 2 x 2 .

2
y 3 = a 0 + a1 x 3 + a 2 x 3

De mme en remplaant x1 par x2 dans la dmonstration

prcdente, on trouve :
x [x1 , x 2 ]

A ce stade, on fait lhypothse que f peut tre

modlise par un polynme de [X] de degr 2.

(9)

En multipliant membre membre les inquations (8) et

(9), lon obtient :
x [x1 , x 2 ]
2
[( x)] [(c x1 ) (c x 2 )] Sup
4

[ x1 , x2 ]

f "

Sur lintervalle [ x1, x2] , la fonction [(c - x1) (c -x2)]2 est

x + x2
(x x1 )4 .
maximale pour c = 1
et vaut 2
2
16

Posons :
1 x1
2 = 1 x2
1 x3

y1
y2
y3

1 x1

x12

et = 1 x 2
1 x3

x 22

x32

Avec ces notations a2 = 2 / . Dj est un

dterminant de Vandermonde dordre 3 et vaut :
= (x3 -x2) (x3 - x1) (x2 - x1). On calcule ensuite 2
par la rgle de Sarrus :
2 = (x2 - x1) y3 + (x1 - x3) y2 +(x3 - x2) y1 .
On en dduit :
a2 =

( x 2 x1 ) y 3 + ( x1 x 3 ) y 2 + ( x 3 x 2 ) y1
.
( x 3 x 2 ) ( x 3 x1 ) ( x 2 x1 )

INCERTITUDE DUE A LINTERPOLATION SUR UNE TABLE DE DONNEES NUMERIQUES

Comme p( x) = 2 a2 , on retiendra :
f " = 2 a2

5. Application numrique

(11)

Rsultat n 3 :
La valeur absolue de lerreur commise en
assimilant f la droite D sur lintervalle [ x1, x2]
est majore par :
M=

Px2
a2 ,
4

avec :
a2 =

( x 2 x1 ) y 3 + ( x1 x 3 ) y 2 + ( x 3 x 2 ) y1
( x 3 x 2 ) ( x 3 x1 ) ( x 2 x1 )

et : (x1, y1), (x2, y2), (x3, y3) les points

dtalonnage donnant la fonction f.

Soit mesurer la temprature en un point dun bain

thermostatique. Pour ce faire, on utilise le thermomtre
numrique de rsolution 0,001 C dont on a reproduit
sur le tableau 2 le relev de mesures extrait de son
certificat dtalonnage.

Valeur lue
(L)
C

Valeur
talon (T)
C

Correction
(C)
C

Incertitude
(sur C)
C

Immersion
mm

-80,173

-80,129

0,044

0,009

380

-0,055

0,000

0,055

0,005

380

100,010

100,102

0,092

0,011

380

199,854

199,988

0,133

0,011

380

298,117

298,307

0,190

0,016

380

399,017

399,264

0,246

0,017

380

Tableau 2
Table dtalonnage du thermomtre

3.4. Calcul de lincertitude type

A priori, toute valeur de lerreur est quiprobable, ce
qui conduit considrer une distribution rectangulaire
de variance :

(2 M )2
v=
12

Dautre part, daprs le rapport de caractrisation,

lcart type exprimental du thermomtre 0 C vaut
0,002 C et 0,004 C 100 C.

do :
u droite =

Lindication brute du thermomtre a pour valeur

moyenne 50,203 C et pour cart type exprimental
0,003 C.

Px2 a 2

4 3

(12)

Remarque : on voit bien ce stade lintrt deffectuer ltalonna ge

de ses instruments de mesure en interne puisque lon dispose alors
de la rptabilit aux points dtalonnage. Cette information nest en
gnral pas contenue dans le certificat dtalonnage, et pourtant, elle
est ncessaire pour pouvoir lexploiter en dehors des points
dtalonnage

Quelle est la correction c a appliquer au thermomtre ?

Quelle est lincertitude type u(c) associe ?

4. Rsultat final
Daprs (7) et (12) :
Py
u 2 = u m2 ( y i ) +
Px

P4 a2
u m2 ( xi ) + u 2 ( x) + x 2 .

48

En suivant le mode opratoire figurant en annexe 1, on

trouve :

Cest dire :
u=

u m2 ( y i )

Py
+
Px

P4 a2
u m2 ( x i ) + u 2 ( x) + x 2 .

48

METGEN DOSSIER METROLOGIE N 1

c = 0,074 C et u(c) = 0,005 8 C .

DM1 Rv. 5

Annexe 1
Mode opratoire Calcul dune correction dtalonnage et de son incertitude
On a la table dtalonnage suivante (donne issue du certificat dtalonnage) :
Valeur appareil
x1
x2
x3

Valeur talon
y1
y2
y3

Correction
c1
c2
c3

Incertitude
I1
I2
I3

Lappareil mesure x (compris entre x1 et x2). Lcart type exprimental sur x vaut S(x).
La rsolution de lappareil est note r.
Remarque : certains appareils peuvent changer automatiquement de rsolution en fonction des valeurs mesures. Si un changement se produit entre
les points (x1, x2, x3) la mthode ne pourra pas tre applique car a priori la fonction dfinie par ci = f(xi) nest pas continue.

1.

Calcul des incertitudes lmentaires :

1.1. Lincertitude type sur c1 note u(c1) est obtenue partir de lincertitude donne dans le certificat dtalonnage.
Cette dernire ayant t calcule en prenant un coefficient dlargissement k = 2, on obtient u(c1) en divisant I1 par
deux : u(c1) = I1 / 2
1.2. Lincertitude type sur c2 note u(c2) est obtenue partir de lincertitude donne dans le certificat dtalonnage.
Cette dernire ayant t calcule en prenant un coefficient dlargissement k = 2, on obtient u(c2) en divisant I2 par
deux : u(c2) = I2 / 2
1.3. Lincertitude type sur x1 note u(x1) est obtenue partir de la rsolution et dune rptabilit effectue la valeur
x1 qui vaut S(x1) :
u ( x1 ) =

r2
+ S 2 ( x1 )
12

1.4. Lincertitude type sur x2 note u(x2) est obtenue partir de la rsolution et dune rptabilit effectue la valeur
x2 qui vaut S(x2) :
u( x 2 ) =

r2
+ S 2 ( x2 )
12

1.5. Lincertitude type sur x note u(x) est obtenue partir de la rsolution et de lcart type exprimentale S(x) et
vaut :
u ( x) =

r2
+ S 2 ( x)
12

1.6. um(ci) = Max { u(c1), u(c2) } et um(xi) = Max { u(x1), u(x2) }

1.7. Px = x2 - x1 et Pc = c2 - c1
1.8. a2 =
2.

( x2 x1 ) c3 + ( x1 x3 ) c2 + ( x3 x2 ) c1
( x3 x2 ) ( x3 x1 ) ( x2 x1)

Rsultats

2.1. La correction ajouter la valeur x est : c =

Pc
( x x1 ) + c1
Px

2.2. Lincertitude type sur la correction est : u (c) =

u m2 (c i )

P
+ c
Px

P4 a2
u m2 ( x i ) + u 2 ( x) + x 2
48

INCERTITUDE DUE A LINTERPOLATION SUR UNE TABLE DE DONNEES NUMERIQUES

Annexe 2
Mode opratoire Fiche de calcul

Fiche de calcul N
Correction dtalonnage
1.

Informations gnrales
N talon :

Certificat dtalonnage N :

Unit :

Date :
Rapport de caractrisation N :
Date :

2.

Donnes dtalonnage

2.2. Points utiliss pour linterpolation (issus du certificat dtalonnage)

Point
n
1

Valeur lue

Valeur talon

Correction

Incertitude

2
2.3. Point supplmentaire (issu du certificat dtalonnage)
Point
n
3
3.

Valeur lue

Valeur talon

Correction

Informations complmentaires sur ltalon

3.1. Rsolution de ltalon

r=

3.2. Ecarts types de ltalon (issus du rapport de caractrisation)

4.

5.

Au point n 1

S=

Au point n 2

S=

Valeur mesure par ltalon

Mesure

x=

Ecart type

S=

Rsultat

5.1. Calculs intermdiaires

u(c1) =

u(c2) =

u(x1) =

u(x2) =

um(ci) =

um(xi) =

Px =

Pc =

u(x) =

a2 =

5.2. Rsultat final

Correction
Incertitude type

METGEN DOSSIER METROLOGIE N 1

c=
u(c) =

DM1 Rv. 5