Résoudre un sytème au moindres carrés ...

(Cliquez-ici pour accéder à la version originale de cette discussion avec couleurs et images)







Posted by: totor

Bonjour,

J'ai une question sur la résolution d'un système au sens des moindres carrés. D'abord, ce que je comprend c'est que quand on résoud un système au sens des moindres carrés c'est que la solution n'existe pas, on essaye d'en trouver une qui sois "au mieux", c'est à dire \min_w || Aw - g||^2. En général on résoud un système qui a plus d'équation que d'inconnues, par exemple ( en deux dimensions) :
<br />
\begin{array}{cc}<br />
au(x,y) + bv(x,y) &amp;=&amp; g_1 \\<br />
cu(x,y) + dv(x,y) &amp;=&amp; g_2<br />
\end{array}<br />
ou a,b,c,d sont des paramètres connus, g = [g_1,g_2]^T est le vecteur connu et w = [u(x,y), v(x,y)]^T sont le vecteur solution inconnu que l'on cherche à déterminer. En fait, ca donne sous forme matricielle :
<br />
Aw = g<br />
avec
<br />
A=<br />
\begin{pmatrix}<br />
a &amp; b\\<br />
c &amp; d <br />
\end{pmatrix}<br />
et la solution se trouve avec les equations normales :
<br />
w = inv(A^TA)A^Tg<br />

Bref, ma question c'est comment on fait quand on cherche à minimiser un truc avec des dérivées ...?..???? par exemple si on a le système avec des dérivées de la solution :
<br />
\begin{array}{cc}<br />
au(x,y) + bv(x,y) + \frac{\partial u(x,y)}{\partial x} + \frac{\partial u(x,y)}{\partial y} &amp;=&amp; g_1 \\<br />
cu(x,y) + dv(x,y) + \frac{\partial u(x,y)}{\partial x} + \frac{\partial u(x,y)}{\partial y} &amp;=&amp; g_2<br />
\end{array}<br />

Quelle technique utiliser ??? Sauriez vous comment faire à ce moment là..?


merci beaucoup par avance!!

(ps : je suis nul en maths! mais ne le dites à personne :p)



Posted by: totor

up ...











-