Trabajo2 Pilarsaldivia

Informe Taller Grupal N°2
Álgebra Lineal para Ingeniería
Resolución de un sistema de ecuaciones lineales, en el sentido de mínimos

cuadrados, mediante factorización QR
Integrantes: Antonio Bottai Saez

Ignacio Caamaño San Martín.
Matías Iván Godoy Poblete.
Nicolás Antonio Gómez Rivas.
Pilar Francisca Saldivia Saldivia.
Profesor: Juan Pablo Concha Rosales.
Valdivia, 03 de Noviembre del 2017.

1.3 Ejercicios:
1. De la expresión (1) deduzca la expresión (2).
- De la expresión (1) se desarrolla la expresión a partir de las propiedades de producto interno
𝑚𝑖𝑛𝑥 ‖𝐴𝑥 − 𝑏‖2 = 〈𝐴𝑥 − 𝑏, 𝐴𝑥 − 𝑏〉 / Producto interno
〈𝐴𝑥, 𝐴𝑥 − 𝑏〉 + 〈−𝑏, 𝐴𝑥 − 𝑏〉 / Aplicación producto interno
〈𝐴𝑥, 𝐴𝑥 − 𝑏〉 − 〈𝑏. 𝐴𝑥 − 𝑏〉 / Distribución
〈𝐴𝑥, 𝐴𝑥〉 − 2〈𝐴𝑥, 𝑏〉 + 〈𝑏, 𝑏〉 / Desarrollo cuadrado de binomio producto interno
〈𝐴𝑥, 𝐴𝑥〉 − 2𝑏 𝑡 𝐴𝑥 + 〈𝑏, 𝑏〉 / Desarrollo por propiedad
𝑡 𝑡
(𝐴𝑥) 𝐴𝑥 − 2𝑏 𝐴𝑥 + 𝑏 𝑏 / 𝑡
𝐴𝑡 𝑋 𝑡 𝐴𝑥 − 2𝑏 𝑡 𝐴𝑥 + 𝑏 𝑡 𝑏/ (2)
⇒ 𝑚𝑖𝑛𝑥 ‖𝐴𝑥 − 𝑏‖2 = 𝐴𝑡 𝑋 𝑡 𝐴𝑥 − 2𝑏 𝑡 𝐴𝑥 + 𝑏 𝑡 𝑏
Luego de métodos algebraico y por propiedades de producto interno comprobamos que de la

expresión 𝑚𝑖𝑛𝑥 ‖𝐴𝑥 − 𝑏‖2 (1) efectivamente llegamos 𝐴𝑡 𝑋 𝑡 𝐴𝑥 − 2𝑏 𝑡 𝐴𝑥 + 𝑏 𝑡 𝑏 (2)
2. De la expresión (3) deduzca la expresión (4).

- De la expresión (3) se desarrolla a partir de las propiedades que se mencionan
𝐴𝑡 𝐴𝑥 = 𝐴𝑡 𝑏 (3)
1
𝐴𝑡 𝑄𝑅𝑥 = 𝐴𝑡 / Multiplicar por 𝐴𝑡 𝑄 por derecha
1
𝑅𝑥 = 𝐴𝑡 𝑄 𝐴𝑡 𝑏
𝑅𝑥 = (𝐴𝑡 𝑄)−1 𝐴𝑡 𝑏
𝑅𝑥 = (𝐴𝑡 )−1 𝑄 −1 𝐴𝑡 𝑏
𝑅𝑥 = 𝐴𝑄 −1 𝐴−1 𝑏 / Propiedad una matriz A es ortogonal si 𝐴𝑡 = 𝐴−1 y (𝐴𝑡 )−1 = 𝐴
1
𝑄 −1 𝐴 = 𝑄 𝑡 𝐴/Multiplicar por 𝐴 por la derecha
1
𝑄 −1 𝐴 𝐴 = 𝑄 𝑡
𝑄 −1 𝐴𝐴−1 = 𝑄 𝑡
𝑅𝑥 = 𝑄 𝑡 𝑏
Por lo tanto, se demuestra que a partir de la expresión 𝐴𝑡 𝐴𝑥 = 𝐴𝑡 𝑏 (3) se puede deducir 𝑅𝑥 = 𝑄 𝑡 𝑏 (4)
3. Expliquen para qué es el proceso de ortogonalización de Gram-Schmidt y las condiciones que se
deben cumplir para poder utilizarlo. Investigue los comandos para realizar el proceso de
ortogonalización de Gram-Schmidt y para hacer la factorización QR, en el software wxMaxima.
De un pequeño ejemplo de su uso.
R: Ortogonalización de Gram-Schmidt es parte fundamental del proceso de factorización QR, ya que
una vez realizado este, de manera correcta se puede proceder a realizar el resto del procedimiento que
se menciona en el marco teórico de este taller, en el inciso 1.1, se menciona que para realizar la
factorización QR se necesitan de un conjunto de vectores ortonormales los cuales se obtienen de este
tipo de ortogonalización. Como condición las columnas de A, deben ser linealmente independientes.
Para realizar el proceso de ortogonalización de Gram-Schmidt en Wxmáxima, se carga la primera
función “eigen” el orden de los comandos es el siguiente:
1.- load (“eigen”)
2.- A: matrix [(V1), (V2), (V3)]
3.- y: gramschmidt (A)
4.- map (innerproduct, [y[1], y[2], y[3]], [y[2], y[3], y[1]])
Para realizar la factorización QR, se carga la primera función “lapack” el orden de los comandos el
siguiente:
1.- load (“lapack”)
2.- [Q,R]:dgeqrf(A)
Se adjunta el ejemplo de uso de los comandos recién presentados:

- Ortogonalización de Gram-Schmidt:
Comando 1 y 2 Load(“eigen”) para cargar la función, y se ingresan las coordenadas de las matrices
previamente corroborado que estas son linealmente independientes.
Se ingresa el segundo comando, en el cual se realiza el proceso de ortogonalización de la matriz A

recién ingresada:
A partir de ello se obtiene y1, y2 e y3, como muestra la siguiente imagen
La normalización se comprueba realizando producto interno entre estos vectores para ello se utiliza el
ultimo comando, como se muestra en la teoría este conjunto es ortogonal si su producto interno da como
resultado [0,0,0]
Así se prueba el proceso de ortogonalización de un vector, utilizando WxMaxima.
- Factorización QR:
Comando 1 y 2 se carga la función (“lapack”) y se ingresa el comando 2 el cual nos entrega lo que sería
QyR
A partir de esto obtenemos que:
Q:
R:
Todos los procedimientos se realizaron a partir de los conceptos entregados en el archivo teórico más
lo investigado sobre comandos para utilizar WxMáxima.
4.La matriz de Hilbert de orden n se define como . Realice una
factorización QR de la matriz de Hilbert de orden 5, además indique cual es la base ortonormal
encontrada por el proceso de ortogonalización de Gram-Schmidt.
Se ingresa la definición de la matriz de Hilbert al software Wxmáxima, además de la condición de
que esta debe ser una matriz de orden 5. Se adjunta este procedimiento en lo siguiente:
Se realiza la ortogonalización de Gram-Schmidt y posteriormente su normalización, la

ortogonalización se realiza utilizando el mismo comando que en el ejercicio anterior, lo que entrega
los vectores, que en el siguiente paso se deben normalizar
Se procede a normalizar cada uno de estos vectores, para ello se ingresa por separado en el siguiente
orden:
Siendo los vectores A0, A2, A3, A4, A5 los ingresados anteriormente, se procede a normalizarlos
cada uno por separado:
La siguiente matriz se forma a partir de los vectores X1,X2, X3,X4,X5 ESTAS SON LAS BASES
ORTONORMALES y como resultado generan Qt ya que se ordena por filas.
A esta última matriz (Qt ), se ingresa el comando “transpose” para transponer esta matriz es decir que
las filas se cambian por columnas, como muestra las siguiente imagen lo que se obtiene es Q.
Para comprobar estos resultados, como dice la teoría QtQ = Id, lo cual se cumple como se muestra en
lo siguiente
Como ultimo procedimiento se debe obtener la matriz R, realizamos desde la expresión QtH = R
Esta última es la matriz R, asi se concluye la factorización QR.

5. Considere las siguientes matrices:
(a) Determine el sistema normal asociado a Ax = b.
Se tiene que 𝐴𝑥 = 𝑏 queda expresada de la siguiente manera
Al realizar la multiplicación correspondiente obtiene

Por consiguiente se tiene la minimizacion de ésta.
Para determinar el sistema normal 𝐴𝑥 = 𝑏, se debe cumplir lo siguiente

𝐴𝑡 𝐴𝑥 = 𝐴𝑡 𝑏
Donde A= At= b= ; x=
Entonces se tiene que:
Se desarrolla el sistema normal asociado a B.

(b) Obtenga la descomposición QR de la matriz A.
Se ingresa a WxMaxima, las coordenadas de la matriz A, luego de ello se ingresa el comando
que se utilizó en la pregunta 3 sobre factorización QR, obteniéndose lo siguiente
Se carga la función “lapack” + comando [Q.R]:dgeqrf(A)
Según lo obtenido no se puede realizar la factorización QR ya que la matriz R, no es diagonal

superior, la cual es su característica, por ello la matriz A no se puede factorizar
(c) Encuentre la solución del sistema en el sentido de los mínimos cuadrados (si es que existe)
usando la descomposición hallada en el inciso anterior.
No se puede realizar el desarrollo para encontrar la solución al sistema en el sentido de
mínimos cuadrados a partir del inciso anterior, ya que la matriz A no es posible factorizarla,
por factorización QR, como se menciona en la pregunta anterior
Los cálculos se realizaron por escrito y en WxMaxima.

Trabajo2 Pilarsaldivia

Cargado por

Copyright:

Formatos disponibles

Trabajo2 Pilarsaldivia

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Trabajo2 Pilarsaldivia

Cargado por

Copyright:

Formatos disponibles

Informe Taller Grupal N°2

Álgebra Lineal para Ingeniería

Resolución de un sistema de ecuaciones lineales, en el sentido de mínimos

Integrantes: Antonio Bottai Saez

Profesor: Juan Pablo Concha Rosales.

Valdivia, 03 de Noviembre del 2017.

Luego de métodos algebraico y por propiedades de producto interno comprobamos que de la

2. De la expresión (3) deduzca la expresión (4).

Se adjunta el ejemplo de uso de los comandos recién presentados:

Se ingresa el segundo comando, en el cual se realiza el proceso de ortogonalización de la matriz A

Así se prueba el proceso de ortogonalización de un vector, utilizando WxMaxima.

Se realiza la ortogonalización de Gram-Schmidt y posteriormente su normalización, la

Esta última es la matriz R, asi se concluye la factorización QR.

(a) Determine el sistema normal asociado a Ax = b.

Se tiene que 𝐴𝑥 = 𝑏 queda expresada de la siguiente manera

Al realizar la multiplicación correspondiente obtiene

Para determinar el sistema normal 𝐴𝑥 = 𝑏, se debe cumplir lo siguiente

Entonces se tiene que:

Se desarrolla el sistema normal asociado a B.

Se carga la función “lapack” + comando [Q.R]:dgeqrf(A)

Según lo obtenido no se puede realizar la factorización QR ya que la matriz R, no es diagonal

Los cálculos se realizaron por escrito y en WxMaxima.

También podría gustarte