FACTORIZACIÓN QR

Martha C. Moreno

26 de octubre de 2011

Departamento de Matemáticas
Universidad Nacional de Colombia

MCMP FACTORIZACIÓN QR
Factorización QR
En muchas aplicaciones de las matrices en diferentes ciencias, se
presentan con frecuencia cálculos matriciales complejos, la idea es
encontrar procedimientos que los simplifiquen aún cuando estos
sean realizados por métodos informáticos. Uno de tales
procedimientos está basado en la aplicación del proceso de
Gram-Schmidt y es conocido como la Factorización QR.

MCMP FACTORIZACIÓN QR
Factorización QR
En muchas aplicaciones de las matrices en diferentes ciencias, se
presentan con frecuencia cálculos matriciales complejos, la idea es
encontrar procedimientos que los simplifiquen aún cuando estos
sean realizados por métodos informáticos. Uno de tales
procedimientos está basado en la aplicación del proceso de
Gram-Schmidt y es conocido como la Factorización QR.

Teorema
Si A es una matriz m × n con columnas l.i, entonces A puede
factorizarse en la forma:

A = QR

En la que Q es una matriz con columnas ortonormales y R es una

matriz triangular superior no singular.

MCMP FACTORIZACIÓN QR
Sea A una matriz m × n con columnas l.i

A = (C1 C2 ··· Cn )

MCMP FACTORIZACIÓN QR
Sea A una matriz m × n con columnas l.i

A = (C1 C2 ··· Cn )

{C1 , C2 , · · · , Cn } es una base del espacio columna de A (espacio

generado por las columnas de A).

MCMP FACTORIZACIÓN QR
Sea A una matriz m × n con columnas l.i

A = (C1 C2 ··· Cn )

{C1 , C2 , · · · , Cn } es una base del espacio columna de A (espacio

generado por las columnas de A).

Aplicando el proceso de Gram-Schmidt podemos construir a partir

de esta una base ortonormal {Q1 , Q2 , · · · , Qn } del espacio
columna.

MCMP FACTORIZACIÓN QR
Sea A una matriz m × n con columnas l.i

A = (C1 C2 ··· Cn )

{C1 , C2 , · · · , Cn } es una base del espacio columna de A (espacio

generado por las columnas de A).

Aplicando el proceso de Gram-Schmidt podemos construir a partir

de esta una base ortonormal {Q1 , Q2 , · · · , Qn } del espacio
columna.
Donde:
vi
Qi = kvi k

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1
v2 = C2 − proyv1 C2

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1
v2 = C2 − proyv1 C2
v3 = C3 − proyv1 C3 − proyv2 C3

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1
v2 = C2 − proyv1 C2
v3 = C3 − proyv1 C3 − proyv2 C3
..
.

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1
v2 = C2 − proyv1 C2
v3 = C3 − proyv1 C3 − proyv2 C3
..
.
vn = Cn − proyv1 Cn − proyv2 Cn − · · · − proyvn−1 Cn

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1
v2 = C2 − proyv1 C2
v3 = C3 − proyv1 C3 − proyv2 C3
..
.
vn = Cn − proyv1 Cn − proyv2 Cn − · · · − proyvn−1 Cn

{v1 , v2 , · · · , vn } es una base ortogonal.

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1
v2 = C2 − proyv1 C2
v3 = C3 − proyv1 C3 − proyv2 C3
..
.
vn = Cn − proyv1 Cn − proyv2 Cn − · · · − proyvn−1 Cn

{v1 , v2 , · · · , vn } es una base ortogonal.

La normalizamos: { kvv11 k , kvv22 k , · · · , kvvnn k } y

MCMP FACTORIZACIÓN QR
La base de la que partimos es: {C1 , C2 , · · · , Cn }

v1 = C 1
v2 = C2 − proyv1 C2
v3 = C3 − proyv1 C3 − proyv2 C3
..
.
vn = Cn − proyv1 Cn − proyv2 Cn − · · · − proyvn−1 Cn

{v1 , v2 , · · · , vn } es una base ortogonal.

La normalizamos: { kvv11 k , kvv22 k , · · · , kvvnn k } y

Qi = kvvii k

MCMP FACTORIZACIÓN QR
Como {Q1 , Q2 , · · · , Qn } es una base del espacio columna,
entonces cada columna de A se puede expresar como una c.l de sus
elementos, es decir:

MCMP FACTORIZACIÓN QR
Como {Q1 , Q2 , · · · , Qn } es una base del espacio columna,
entonces cada columna de A se puede expresar como una c.l de sus
elementos, es decir:

C1 = r11 Q1 + r21 Q2 + · · · + rn1 Qn

MCMP FACTORIZACIÓN QR
Como {Q1 , Q2 , · · · , Qn } es una base del espacio columna,
entonces cada columna de A se puede expresar como una c.l de sus
elementos, es decir:

C1 = r11 Q1 + r21 Q2 + · · · + rn1 Qn

C2 = r12 Q1 + r22 Q2 + · · · + rn2 Qn
..
.

MCMP FACTORIZACIÓN QR
Como {Q1 , Q2 , · · · , Qn } es una base del espacio columna,
entonces cada columna de A se puede expresar como una c.l de sus
elementos, es decir:

C1 = r11 Q1 + r21 Q2 + · · · + rn1 Qn

C2 = r12 Q1 + r22 Q2 + · · · + rn2 Qn
..
.
Cn = r1n Q1 + r2n Q2 + · · · + rnn Qn

MCMP FACTORIZACIÓN QR
Como {Q1 , Q2 , · · · , Qn } es una base del espacio columna,
entonces cada columna de A se puede expresar como una c.l de sus
elementos, es decir:

C1 = r11 Q1 + r21 Q2 + · · · + rn1 Qn

C2 = r12 Q1 + r22 Q2 + · · · + rn2 Qn
..
.
Cn = r1n Q1 + r2n Q2 + · · · + rnn Qn

De donde:

MCMP FACTORIZACIÓN QR
MCMP FACTORIZACIÓN QR
A = (C1 C2 ··· Cn )

MCMP FACTORIZACIÓN QR
A = (C1 C2 ··· Cn )  
r11 r12 · · · r1n
 r21 r22 · · · r2n 
A = (Q1 Q2 ··· Qn )  .. .. .. ..
 

 . . . . 
rn1 rn2 · · · rnn

MCMP FACTORIZACIÓN QR
A = (C1 C2 ··· Cn )  
r11 r12 · · · r1n
 r21 r22 · · · r2n 
A = (Q1 Q2 ··· Qn )  .. .. .. ..
 

 . . . . 
rn1 rn2 · · · rnn
A=Q R

MCMP FACTORIZACIÓN QR
A = (C1 C2 ··· Cn )  
r11 r12 · · · r1n
 r21 r22 · · · r2n 
A = (Q1 Q2 ··· Qn )  .. .. .. ..
 

 . . . . 
rn1 rn2 · · · rnn
A=Q R
Como la base es ortonormal entonces:
 
C1 • Q 1 C2 • Q 1 · · · C n • Q1
 C1 • Q 2 C2 • Q 2 · · · C n • Q2 
R= .. .. .. ..
 

 . . . . 
C1 • Q n C2 • Q n · · · C n • Qn

MCMP FACTORIZACIÓN QR
Como Ci • Qj = 0 si i < j, pues cada Qj es ortogonal a
C1 , C2 , · · · , Cj−1 , entonces:

MCMP FACTORIZACIÓN QR
Como Ci • Qj = 0 si i < j, pues cada Qj es ortogonal a
C1 , C2 , · · · , Cj−1 , entonces:
 
C1 • Q 1 C2 • Q 1 · · · C n • Q1
 0 C2 • Q2 · · · C n • Q2 
R = .. .. .. ..
 

 . . . . 
0 0 ··· C n • Qn

MCMP FACTORIZACIÓN QR
Como Ci • Qj = 0 si i < j, pues cada Qj es ortogonal a
C1 , C2 , · · · , Cj−1 , entonces:
 
C1 • Q 1 C2 • Q 1 · · · C n • Q1
 0 C2 • Q2 · · · C n • Q2 
R = .. .. .. ..
 

 . . . . 
0 0 ··· C n • Qn

Matriz Triangular superior y no singular ¿Porqué?

MCMP FACTORIZACIÓN QR
Como Ci • Qj = 0 si i < j, pues cada Qj es ortogonal a
C1 , C2 , · · · , Cj−1 , entonces:
 
C1 • Q 1 C2 • Q 1 · · · C n • Q1
 0 C2 • Q2 · · · C n • Q2 
R = .. .. .. ..
 

 . . . . 
0 0 ··· C n • Qn

Matriz Triangular superior y no singular ¿Porqué?

Note que: rij = Qi • Cj

MCMP FACTORIZACIÓN QR
Nota
Recuerde que:
Para poder aplicar el proceso de factorización QR, la matriz
debe tener columnas linealmente independientes.

MCMP FACTORIZACIÓN QR
Nota
Recuerde que:
Para poder aplicar el proceso de factorización QR, la matriz
debe tener columnas linealmente independientes.
Es posible factorizar matrices no cuadradas.

MCMP FACTORIZACIÓN QR
Nota
Recuerde que:
Para poder aplicar el proceso de factorización QR, la matriz
debe tener columnas linealmente independientes.
Es posible factorizar matrices no cuadradas.

Teorema
Si A es cuadrada no singular, entonces es factorizable QR, y la
matriz Q es ortogonal.

MCMP FACTORIZACIÓN QR
Ejemplo
 
1 2
Sea A =
−1 3

MCMP FACTORIZACIÓN QR
Ejemplo
 
1 2
Sea A =
−1 3
   
1 2
Como las columnas y son l.i, entonces la matriz
−1 3
es factorizable QR.

MCMP FACTORIZACIÓN QR
Ejemplo
 
1 2
Sea A =
−1 3
   
1 2
Como las columnas y son l.i, entonces la matriz
−1 3
es factorizable QR.

Construimosuna base
 ortonormal
 a partir de la base:
1 2
B = {C1 = , C2 = }
−1 3

MCMP FACTORIZACIÓN QR
  
1
v1 = C 1 =
−1

MCMP FACTORIZACIÓN QR
  
1
v1 = C 1 =
−1
5
     
2 −1 1 2
v2 = C2 − proyv1 C2 = − 2 = 5
3 −1 2

MCMP FACTORIZACIÓN QR
  
1
v1 = C 1 =
−1
5
     
2 −1 1 2
v2 = C2 − proyv1 C2 = − 2 = 5
3 −1 2

Normalizando:
 
1 !
√1
 
−1 2
Q1 = √
2
= −1
√
2

MCMP FACTORIZACIÓN QR
  
1
v1 = C 1 =
−1
5
     
2 −1 1 2
v2 = C2 − proyv1 C2 = − 2 = 5
3 −1 2

Normalizando:
 
1 !
√1
 
−1 2
Q1 = √
2
= −1
√
  2
5
 2  !
5 √1
2 2
Q2 = 5
√
2
= √1
2
2

MCMP FACTORIZACIÓN QR
Entonces:
!
√1 √1
Q= 2 2
−1
√ √1
2 2

MCMP FACTORIZACIÓN QR
Entonces:
!
√1 √1
Q= 2 2
−1
√ √1
2 2
  √ √
− 2
!
Q1 • C 1 Q1 • C 2 2 2
R= = √
Q2 • C 1 Q2 • C 2 0 5 2
2

MCMP FACTORIZACIÓN QR
Entonces:
!
√1 √1
Q= 2 2
−1
√ √1
2 2
  √ √
− 2
!
Q1 • C 1 Q1 • C 2 2 2
R= = √
Q2 • C 1 Q2 • C 2 0 5 2
2
! √ √ !
√1 √1 2 − 2
A= 2 2 2
√
−1
√ √1 5 2
2 2 0 2

MCMP FACTORIZACIÓN QR