Factorizacin QR
Factorizacin QR
En lgebra lineal, la descomposicin o factorizacin QR de una matriz es una descomposicin de la misma como producto de una matriz ortogonal por una triangular superior. La descomposicin QR es la base del algoritmo QR utilizado para el clculo de los vectores y valores propios de una matriz.
Definicin
La descomposicin QR de una matriz cuadrada real A es donde Q es una matriz ortogonal (QTQ = I ) y R es una matriz triangular superior.
Clculo de la descomposicin QR
Mediante el mtodo de ortogonalizacin de Gram-Schmidt
Recurriendo al mtodo de ortogonalizacin de Gram-Schmidt, con las columnas de A como los vectores a procesar. . Entonces
Naturalmente, utilizamos los ais de A para obtener:
Ahora estas ecuaciones pueden ser escritas en forma matricial de esta manera:
:::::::::
El producto de cada fila con cada columa de las matrices de arriba, nos da la respectiva columna de A con la que comenzamos y, por tanto, dada la matriz A, la hemos factorizado en una matriz ortogonal Q (la matriz de eks), aplicando el proceso de Gram-Schmidt, y la matriz resultante triangular superior es R. Alternativamente, la matriz puede calcularse de la siguiente manera:
Factorizacin QR Recordemos que: Entonces, tenemos
Note que Ejemplo Si se considera la descomposicin de
, entonces
Se busca la matriz ortogonal
tal que
Por lo que calculamos
mediante Gram-Schmidt como sigue:
Por lo tanto, tenemos
Mediante el uso de reflexiones de Householder
Una transformacin de Householder o reflexin de Householder es una transformacin que refleja el espacio con respecto a un plano determinado. Esta propiedad se puede utilizar para realizar la transformacin QR de una matriz si tenemos en cuenta que es posible elegir la matriz de Householder de manera que un vector elegido quede con una nica componente no nula tras ser transformado (es decir, premultiplicando por la matriz de Householder). Grficamente, esto significa que es posible reflejar el vector elegido respecto de un plano de forma que el reflejo quede sobre uno de los ejes de la base cartesiana. La manera de elegir el plano de reflexin y formar la matriz de Householder asociada es el siguiente: Sea un vector columna arbitrario m-dimensional tal que || || = ||, donde es un escalar; (si el algoritmo se implementa utilizando aritmtica de coma flotante, entonces debe adoptar el signo contrario que 1 para evitar prdida de precisin). Entonces, siendo el vector (1,0,...,0)T, y |||| la norma eucldea, se define:
Factorizacin QR es un vector unitario perpendicular al plano de reflexin elegido. dicho plano. Esta propiedad se puede usar para transformar gradualmente los vectores columna de una matriz A de dimensiones m por n en una matriz triangular superior. En primer lugar se multiplica A con la matriz de Householder Q1 que obtenemos al elegir como vector la primera columna de la matriz. Esto proporciona una matriz QA con ceros en la primera columna (excepto el elemento de la primera fila). es una matriz de Householder asociada a
el procedimiento se puede repetir para A (que se obtiene de A eliminando la primera fila y columna), obteniendo as una matriz de Householder Q2. Hay que tener en cuenta que Q2 es menor que Q1. Para conseguir que esta matriz opere con Q1A en lugar de A es necesario expandirla hacia arriba a la izquierda, completando con un uno en la diagonal, o en general:
Tras repetir el proceso
veces, donde
es una matriz triangular superior. De forma que tomando
es una descomposicin QR de la matriz
Este mtodo tiene una estabilidad numrica mayor que la del mtodo de Gram-Schmidt descrito arriba. Una pequea variacin de este mtodo se utiliza para obtener matrices semejantes con la forma de Hessenberg, muy tiles en el clculo de autovalores por acelerar la convergencia del algoritmo QR reduciendo as enormemente su coste computacional. Ejemplo Vamos a calcular la descomposicin de la matriz
En primer lugar necesitamos encontrar una reflexin que transforme la primera columna de la matriz A, vector , en usando la expresin,
en nuestro caso : y Por lo tanto
Factorizacin QR
, entonces
Ahora observamos:
con lo que ya casi tenemos una matriz triangular. Slo necesitamos hacer cero en el elemento (3,2). Tomando la submatriz bajo el (1, 1) y aplicando de nuevo el proceso a
Mediante el mismo mtodo que antes obtenemos la matriz de Householder
Finalmente obtenemos
La matriz Q es ortogonal y R es triangular superior, de forma que A = QR es la descomposicin QR buscada.
Mediante rotaciones de Givens
Las descomposiciones QR tambin puden calcularse utilizando una serie de rotaciones de Givens. Cada rotacin anula (hace cero) un elemento en la subdiagonal de la matriz, formando de este modo la matriz R. La concatenacin de todas las rotaciones de Givens realizadas, forma la matriz ortogonal Q. En la prctica, las rotaciones de Givens no se utilizan en la actualidad para construir una matriz completa y realizar un producto de matrices. En su lugar, se utiliza un procedimiento de rotacin de Givens, que es equivalente a la multiplicacin reducida de matrices de Givens, sin el trabajo extra de manejar los elementos reducidos. El procedimiento de rotacin de Givens es til en situaciones donde slo pocos elementos fuera de la diagonal necesitan ser anulados y es ms fcil de paralelizar que las transformaciones de Householder.
Factorizacin QR Ejemplo Calculemos la descomposicin de
Primero, necesitamos formar una matriz de rotacin tal que hagamos cero el elemento ms inferior a la izquierda, . Construimos esta matriz empleando el mtodo de la rotacin de Givens y llamamos la matriz resultante . Rotamos primero el vector , representndolo a lo largo del eje X. Este vector forma un ngulo :
. Creamos la matriz ortogonal de rotacin de Givens,
Y el resultado de
tiene ahora un cero en el
elemento.
Procedemos anlogamente con las matrices de Givens y , formando una matriz triangular todas las matrices de Givens descomposicin QR es .
, que hacen cero los elementos subdiagonales es formada a partir del producto en cadena de , y la
. La matriz ortogonal
. Luego tenemos
Relacin con el determinante
Es posible utilizar la descomposicin QR para encontrar el valor absoluto del determinante de una matriz. Suponiendo que una matriz se descompone segn . Entonces se tiene
Puesto que Q es unitaria,
. Por tanto,
donde
son los valores de la diagonal de R.
Fuentes y contribuyentes del artculo
Fuentes y contribuyentes del artculo
Factorizacin QR Fuente: http://es.wikipedia.org/w/index.php?oldid=71922683 Contribuyentes: GermanX, Gonchibolso12, Juan Mayordomo, Juankvillegas, Kender00, Nasil, Newtonfn, Plux, Rdaneel, Republicanito, Xosema, Yeza, 33 ediciones annimas
Licencia
Creative Commons Attribution-Share Alike 3.0 //creativecommons.org/licenses/by-sa/3.0/
