800 likes | 1.13k Views
MÍNIMOS CUADRADOS. A mxn. A x = b Sistema Inconsistente. A x = b consistente. b está en C A. A x = b inconsistente. b no está en C A. b. A x* es un vector del espacio columna C A. A x* = proy C A b A x* = b* Sistema Consistente. b – A x* mínima.
E N D
A mxn A x = b Sistema Inconsistente
A x = b consistente b está en CA A x = b inconsistente b no está en CA
Ax* es un vector del espacio columna CA A x* = proy CAb A x* = b* Sistema Consistente
x* es una solución de A x* = b* x* es una solución de aproximación de A x = b
A fin de encontrar x* a partir de A x* = b* podríamos partir de A x* = proyCAb
Existe una mejor manera de conseguirlo ( b – A x* ) es ortogonal a cada vector de CA
Por lo tanto, ( b – A x* ) es ortogonal a cada vector columnade A
c1 . ( b – A x* ) = 0 c2 . ( b – A x* ) = 0 c1T . ( b – A x* ) = 0 c2T . ( b – A x* ) = 0
AT . ( b – A x* ) = 0 AT b – ATA x* = 0 ATA x* = AT b
A mxn y b en Rm A x = b siempre tiene al menos una solución por mínimos cuadrados ( o por aproximación ) x*
x* es una solución por mínimos cuadrados de A x = b si y sólo si x* es una solución de las ecuaciones normales ATA x* = AT b
A tiene columnas LI si y sólo siATA es InvertibleEn este casola solución de aproximación de A x = b es única y está dada por
x* = ( AT A )-1ATb seudoinversa de A
x - y = 0 x + y = 0 SEL y = 1 Inconsistente A x = b =
Columnas de A LIATA = Invertiblex* única solución por aproximación
x* = ( AT A )-1ATb x* =
x* solución por mínimos cuadrados de A x = b b – A x* error de mínimos cuadrados
= b – A x* vector de error de mínimos cuadrados
vector de error de mínimos cuadrados 1 - 2 3 = b - A x* =
error de mínimos cuadrados =b – A x* 0,2721655
Observar las ecuacionesdel sistema x - y = 0 0 – ( 0 -1/3 ) = 1/31 x + y = 0 0 – ( 0 +1/3 ) = -1/3 2 y = 1 1 – ( 0 +1/3 ) = 2/3 3
Columnas de A LDATA No es Invertible las ecuaciones normales ATA x* = AT b tienen un número infinito de soluciones
Buscaremos entonces la solución x* de menor longitud (la más cercana al origen)
APLICACIONES Ajuste de Curvas por Mínimos Cuadrados
Encontrar la recta que da el mejor ajuste para los puntos (1,4) , (-2,5), (3,-1) y (4,1) y = b + mx
4 = b + m 5 = b - 2m -1 = b + 3m 1 = b + 4m Sistema Inconsistente
= Columnas de A LIATA Invertiblex* única solución por aproximación
x* = ( AT A )-1ATy x* = y = 3,57 – 0,88 x
vector de error de mínimos cuadrados = _ _ A x* y =
error de mínimos cuadrados =b – A x* 2,579224
Observar la primera ecuacióndel sistema 4 = b + m 4 = 3,57 + (- 0,88) 4 – 2,69 = 1,31= 1 (primer componente del vector )
Encontrar el mejor ajuste cuadrático para los puntos(1,4) , (-2,5), (3,-1) y (4,1) y = a + bx + cx2
4 = a + b + c 5 = a - 2b + 4c -1 = a + 3b + 9c 1 = a + 4b + 16c Sistema Inconsistente
= Columnas de A LIATA Invertiblex* única solución por aproximación
x* = ( AT A )-1ATy x* = y = 3,75 – 0,81 x – 0,04 x2
vector de error de mínimos cuadrados _ = _ y x* A =
error de mínimos cuadrados =b – A x* 2,5244009
Observar la primera ecuacióndel sistema 4 = a + b + c 4 = 3,75 +(- 0,81)+(- 0,04) 4 – 2,9 = 1,1= 1 (primer componente del vector )