Buenas,
Me surgió una duda más bien teórica al estar haciendo el Ej7 del Prático 5.
Estoy computando la matriz ortonormal y triangular superior para una matriz 
con el método de Householder en octave, y si bien consigo computar matrices con estas características (llamemosles 
y 
) al hacer 
esperaba conseguir la matriz original a menos de errores muy pequeños (orden de 
, con 
) pero veo que realmente no consigo eso. Tengo diferencias de hasta 10 unidades en una matriz simple que hice (tamaño 4*3, con |entradas|<10).
Entonces vi el código de Householder que esta colgado en EVA para ver que tan esperable eran estos resultados y vi que las matrices que se consiguen allí (
, 
) al multiplicarlas tampoco devuelven exacto la matriz original .
Mis preguntas específicas son:
¿El método de Householder nos debería dar la factorizacion 
tal que 
?
¿Qué tanto error puede tener este método?
Saludos!
