Vincent Trinquet IVP 2 Avril 2007 - EIVP
Vincent Trinquet IVP 2 Avril 2007 - EIVP
Vincent Trinquet IVP 2 Avril 2007 - EIVP
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
cad<br />
y1= a0 + a1 (-1)+ a2 (-1) + a12 (-1)(-1)<br />
y2= a0 + a1(1) + a2 (-1) + a12 (1) (-1)<br />
y3= a0 + a1 (-1)+ a2 (1) + a12 (-1) (1)<br />
y1= a0 + a1 (1) + a2 (1) + a12 (1) (1)<br />
cad<br />
y1= a0 - a1 - a2 + a12<br />
y2= a0 + a1 - a2 - a12<br />
y3= a0 - a1 + a2 - a12<br />
y4= a0 + a1 + a2 + a12<br />
y1, y2, y3 , y4 sont mesurés expérimentalement donc connus.<br />
On en déduit a0, a1, a2, a12 grâce à ce système de 4 équations à 4 inconnues<br />
Son expression matricielle est :<br />
⎛ y1<br />
⎞ ⎛1<br />
⎜ ⎟ ⎜<br />
⎜ y2<br />
⎟ ⎜1<br />
⎜ ⎟ =<br />
y ⎜<br />
3 1<br />
⎜ ⎟ ⎜<br />
⎜ ⎟ ⎜<br />
⎝ y4<br />
⎠ ⎝1<br />
−1<br />
+ 1<br />
−1<br />
+ 1<br />
−1<br />
−1<br />
+ 1<br />
+ 1<br />
+ 1⎞⎛a<br />
⎟⎜<br />
−1⎟⎜<br />
a<br />
−1⎟⎜<br />
a<br />
⎟⎜<br />
+ 1⎟⎜<br />
⎠⎝a<br />
⎛ y1<br />
⎞ ⎛a<br />
⎜ ⎟ ⎜<br />
⎜ y2<br />
⎟ ⎜a<br />
y=X a avec y= ⎜ ⎟ , a =<br />
y<br />
⎜<br />
3 a<br />
⎜ ⎟ ⎜<br />
⎜ ⎟<br />
⎝ y ⎜<br />
4 ⎠ ⎝a<br />
⎛1<br />
⎜<br />
⎜1<br />
et X= ⎜1<br />
⎜<br />
⎝1<br />
−1<br />
+ 1<br />
−1<br />
+ 1<br />
_1<br />
−1<br />
+ 1<br />
+ 1<br />
+ 1⎞<br />
⎟<br />
−1⎟<br />
−1⎟<br />
⎟<br />
+ 1<br />
⎟<br />
⎠<br />
0<br />
1<br />
2<br />
12<br />
0<br />
1<br />
2<br />
⎞<br />
⎟<br />
⎟<br />
⎟<br />
⎟<br />
⎠<br />
12<br />
⎞<br />
⎟<br />
⎟<br />
⎟<br />
⎟<br />
⎠<br />
On remarque que, par construction,<br />
- la première colonne de la matrice X n’est composée que de 1<br />
- la 2 ème colonne de la matrice X est composée de signes de x1<br />
- la 3 ème colonne de la matrice X est composée de signes de x2<br />
- la 4 ème colonne de la matrice X est composée de signes de x1x2<br />
Ceci est dû au fait que ymod =a0 + a1 x1+ a2 x2+ a12 x1 x2<br />
En renversant le système, on trouve a :<br />
a=X -1 y avec X -1 =matrice inverse de X<br />
Comme X est une matrice d’Hadammard, X –1 peut être obtenue très facilement (cf annexe).<br />
4