5 Método de mínimos cuadrados

Sea AX=bAX=b un sistema de ecuaciones lineales con AMm×n(K)A\in M_{m\times n}(K). Dado X~Kn\tilde{X}\in K^{n}, llamamos error cometido al tomar X~\tilde{X} como solución del sistema a ϵX~=AX~b0\epsilon_{\tilde{X}}=\left\lVert A\tilde{X}-b\right\rVert\geq 0.

Remark 5.1.

Si X~\tilde{X} es solucion del sistema, ϵX~=0\epsilon_{\tilde{X}}=0.

Nuestro objetivo, dado un sistema AX=bAX=b, es hallar X~\tilde{X} tal que eX~e_{\tilde{X}} sea lo más pequeño posible.

Sea SS el subespacio generado por las columnas de AA. Como vv, tomamos el término independiente. Buscamos X~\tilde{X} tal que AX~A\tilde{X} sea pS(b)p_{S}(b), ya que es el vector de SS que “mas se parece a bb”, es decir, si X~\tilde{X} es tal que AX~=pS(b)A\tilde{X}=p_{S}(b) entonces ϵX~=AX~b=pS(b)b\epsilon_{\tilde{X}}=\left\lVert A\tilde{X}-b\right\rVert=\left\lVert p_{S}(b)% -b\right\rVert es minimo.

Como hallamos X~\tilde{X} tal que AX~=pS(b)A\tilde{X}=p_{S}(b)?

Necesitamos hallar X~\tilde{X} tal que bAX~Sb-A\tilde{X}\in S^{\perp}. Por lo tanto, hay que imponer que ,,0=A(n),bAX~,\ldots,0=\langle A^{(n)},b-A\tilde{X}\rangle.

0\displaystyle 0
=A(1),bAX~=(A(1))t¯(bAX~)\displaystyle{}=\langle A^{(1)},b-A\tilde{X}\rangle=\overline{(A^{(1)})^{t}}% \cdot(b-A\tilde{X})
0\displaystyle 0
=A(2),bAX~=(A(2))t¯(bAX~)\displaystyle{}=\langle A^{(2)},b-A\tilde{X}\rangle=\overline{(A^{(2)})^{t}}% \cdot(b-A\tilde{X})
\displaystyle{}\vdots
0\displaystyle 0
=A(n),bAX~=(A(n))t¯(bAX~)\displaystyle{}=\langle A^{(n)},b-A\tilde{X}\rangle=\overline{(A^{(n)})^{t}}% \cdot(b-A\tilde{X})

es decir, At¯(bAX~)=(00)\overline{A^{t}}\cdot(b-A\tilde{X})=\begin{pmatrix}0\\ \vdots\\ 0\\ \end{pmatrix}, es decir At¯AX~=At¯b\overline{A^{t}}\cdot A\tilde{X}=\overline{A^{t}}b sistema auxiliar S.C. (puede ser determinado o indeterminado, pero ϵX~\epsilon_{\tilde{X}} es unico).

La solucion de este sistema auxiliar es la solucion de minimos cuadrados del sistema original.