13.07.2015 Views

Simulation numérique et expérimentale du comportement ...

Simulation numérique et expérimentale du comportement ...

Simulation numérique et expérimentale du comportement ...

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

IV. Choix <strong>et</strong> identification des lois de <strong>comportement</strong>Algorithmes de Gauss-Newton <strong>et</strong> de Levenberg-MarquardtCes algorithmes sont particulièrement adaptés aux fonctionnelles coûts moindres carrés <strong>du</strong>type :2Équation 31 J ( x) = j k( x)Nk=1exp( u )EFoù par exemple j ( x) u ( x)k= .kkLa particularité de ces fonctionnels coûts réside dans le fait que l’on connaît la forme de leursdérivées premières (Équation 32) <strong>et</strong> secondes (Équation 33):JxÉquation 32 (,xJ ) = = 2i ijNkjkk= 1 xi2 JxxÉquation 33 ( H )ij= = 2ijNk=1$"jk# xijxkj+jk2 jkxxij!Alors, en supposant que le deuxième terme de l'Équation 33 est négligeable devant le premier(ce qui est le cas quand on s'approche de l’optimum), on peut écrire :Équation 342 Jxxij- 2Nk=1$"jk# xijxkj!Il est intéressant à ce niveau d’intro<strong>du</strong>ire la matrice de sensibilité ou matrice jacobiennedéfinie par :Équation 353 j11 x111 j2A = 1 x11 ...1jN112x1j1x2j2x2...jNx2............j10x.n.j2.x.1... .j.N.xN./On peut ainsi exprimer le gradient <strong>et</strong> le hessien par :Équation 36,cJ = 2 AjÉquation 37H= 2 ATAavec j = [ j ,........ ] T1j N87

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!