Álgebra y Geometría

Tema 1. Matrices

Rodrigo García Manzanas

Ruth Carballo Fidalgo
DPTO. DE MATEMÁTICA APLICADA Y CIENCIAS DE LA
COMPUTACIÓN
Este tema se publica bajo Licencia:
Creative Commons BY-NC-SA 4.0
C ONTENIDOS

1 Introducción

2 Operaciones con matrices

3 Determinantes

4 Matriz inversa

5 Matrices elementales

6 Formas escalonada y reducida

7 Factorización de matrices

G1954: Álgebra y Geometría Matrices 1/37

 Introducción

M ATRICES

D EFINICIÓN
Una matriz no es más que un conjunto de elementos (normalmente números) dispuestos en
filas y columnas

···
 
a11 a12 a1n
 a21 a22 ··· a2n 
A = {aij } =  . .. .. .. 
 .
. . . .

am1 am2 ··· amn

El índice i hace referencia a las filas y el j a las columnas. Una matriz A con m filas y n
columnas se denota como Am×n y se dice que es de tamaño m × n. Su dimensión es igual
al número de elementos que contiene: m × n

Ejemplos:
 
  1 2
1 2 3
A2×3 = A3×2 = 3 4
4 5 6
5 6
tamaño 2 × 3, dimensión 6
tamaño 3 × 2, dimensión 6
G1954: Álgebra y Geometría Matrices 2/37
 Introducción

T IPOS DE MATRICES

Atendiendo a su tamaño:
Matriz rectangular: m 6= n
 
Matriz fila: m = 1, 1 2 3
 
1
 
Matriz columna: n = 1, 
2

3
!
1 2
Matriz cuadrada: m = n (orden), A2×2 = A2 =
3 4
Atendiendo a sus elementos:
Matriz real: aij ∈ R, ∀ i, j
Matriz compleja: ∃ algún aij ∈ C
Matriz nula: aij = 0, ∀ i, j
G1954: Álgebra y Geometría Matrices 3/37
 Introducción

T IPOS DE MATRICES

Matrices cuadradas
2 0 0

Matriz diagonal 0 7 → en rojo: diagonal principal
0
0 0 −3
1 0 0
Matriz identidad I3 = 0 1 0 → unos en la diagonal principal
001
2 1 0
Matriz triangular superior 0 −3 3
0 0 0
2 0 0
Matriz triangular inferior 4 −3 0
3 1 0
0 1 0
Matriz estrictamente triangular (superior) 0 0 3
000
5 1 3
Matriz simétrica 1 8 2 ⇔ aij = aji ∀ i 6= j
329
 
0 1 −3
Matriz antisimétrica −1 0 2 ⇔ aij = −aji ∀ i 6= j; aii = 0 ∀ i
3 −2 0

En las matrices   se llama traza a la suma de los elementos de la diagonal principal.

cuadradas,
1 2
Por ejemplo, tr =1+4=5
3 4
G1954: Álgebra y Geometría Matrices 4/37
 Operaciones con matrices

O PERACIONES BÁSICAS CON MATRICES

Sean A y B matrices del mismo tamaño

A y B son iguales si y sólo si aij = bij , ∀ i, j
   
5 1 5 1
4 −2 = 4 −2
−1 3 −1 3

A ± B = C ⇔ cij = aij ± bij , ∀ i, j

     
5 1 2 1 7 2
4 −2 + −2 −2 = 2 −4
−1 3 4 −6 3 −3

α escalar ∈ R ⇒ αA = {αaij }, ∀ i, j
   
5 1 −20 −4
−4 4 −2 = −16 8
−1 3 4 −12

G1954: Álgebra y Geometría Matrices 5/37

 Operaciones con matrices

S UMA DE MATRICES

Sean A, B, C, O ∈ Mm×n , siendo O la matriz nula

Propiedades:
Asociativa
A + (B + C) = (A + B) + C
Conmutativa
A+B =B+A
Elemento neutro
A+O =O+A=A
Elemento opuesto
A + (−A) = O

G1954: Álgebra y Geometría Matrices 6/37

 Operaciones con matrices

P RODUCTO DE UN ESCALAR POR UNA MATRIZ

Sean α, β ∈ R y A, B ∈ Mm×n

Propiedades:
Asociativa respecto del producto por escalares
(αβ)A = α(βA)
Conmutativa
αA = Aα
Distributiva respecto de la suma de matrices
α(A + B) = αA + αB
Distributiva respecto de la suma de escalares
(α + β)A = αA + βA

G1954: Álgebra y Geometría Matrices 7/37

 Operaciones con matrices

P RODUCTO DE MATRICES

Para que dos matrices A y B se puedan multiplicar, el número de columnas

de A tiene que coincidir con el de filas de B. En estas circunstancias, se
define el producto AB como otra matriz C que tiene tantas filas como A y
columnas como B
(m × n) · (n × p) = (m × p)
| {z } | {z } | {z }
tamaño de A tamaño de B tamaño de C

Los elementos de C se calculan del siguiente modo:

X
AB = C = {cij } | cij = aik bkj , ∀i, j
k
Ejemplo:  
  1 0
2 3 1 −1 2 =
0 −1 −2
−2 3
   
2 · 1 + 3 · (−1) + 1 · (−2) 2·0+3·2+1·3 −3 9
= =
0 · 1 + (−1) · (−1) + (−2) · (−2) 0 · 0 + (−1) · 2 + (−2) · 3 5 −8

G1954: Álgebra y Geometría Matrices 8/37

 Operaciones con matrices

P RODUCTO DE MATRICES

Sean A, B, C ∈ Mm×n

Propiedades (I):

Asociativa
A(BC) = (AB)C

Distributiva respecto de la suma de matrices

A(B + C) = AB + AC
(B + C)A = BA + CA

El producto de matrices no siempre es conmutativo. De hecho, en

general, AB 6= BA
Ejemplo: Calcula AB y BA para las siguientes matrices:
34

62


A= 12 B= 32

Cuando dos matrices verifican que AB = BA, se dicen conmutativas. Si verifican que
AB = −BA, se dicen anticonmutativas

G1954: Álgebra y Geometría Matrices 9/37

 Operaciones con matrices

P RODUCTO DE MATRICES

Sean A, B, C ∈ Mm×n y O la matriz nula

Propiedades (II):

El producto de matrices tiene divisores de O

En general, AB = O ; A = O y/o B = O
Ejemplo: Calcula AB para las siguientes matrices:
11

−3 −2


A= 11 B= 3 2

El producto de matrices no verifica la propiedad de simplificación

En general, AB = AC ; B = C
Ejemplo: Calcula AB y AC para las siguientes matrices:
2 −1
11

−3 −2



A= 11 B= 3 2
C= −2 1

G1954: Álgebra y Geometría Matrices 10/37

 Operaciones con matrices

P RODUCTO DE MATRICES

Sean A, B, C ∈ Mnxn (matrices cuadradas)

Propiedades:
El producto de dos matrices triangulares, ambas superiores (inferiores), es otra matriz
triangular superior (inferior)
         
121 211 2 2 10 100 200 2 00
030 004 = 0 0 12 120 300 = 8 00
001 001 00 1 233 121 16 6 3

El producto de dos matrices diagonales es otra matriz diagonal que se obtiene

multiplicando, elemento a elemento, las dos diagonales
    
100 300 300
020 020 = 040 ←− el producto de matrices diagonales es conmutativo
003 001 003

Para la matriz identidad, se verifica Am×n In = Am×n , Im Am×n = Am×n

102

1 0 0 102





023 010 = 023
10 212 = 212
001 01 143 143

Potencia de una matriz cuadrada

Ak = AA · · · AA
| {z }
k veces
Por convenio: A0 = I

G1954: Álgebra y Geometría Matrices 11/37

 Operaciones con matrices

M ATRIZ TRANSPUESTA

Dada una matriz cualquiera A, se llama transpuesta (At ) a la matriz

que resulta de cambiar ordenadamente las filas por las columnas
1 4
A = ( 14 25 36 ) → At = 25
1 36
A= 2 → At = ( 1 2 3 )
3

Propiedades:
(At )t = A
(λA)t = λAt , ∀ λ ∈ R
(A ± B)t = At ± B t
(AB)t = B t At
A simétrica ⇔ A = At
A antisimétrica ⇔ −A = At
G1954: Álgebra y Geometría Matrices 12/37
 Operaciones con matrices

M ATRIZ TRANSPUESTA

T EOREMAS
1 Dada una matriz cuadrada A, A + At es una matriz simétrica
2 Dada una matriz cuadrada A, A − At es una matriz antisimétrica
3 Toda matriz cuadrada A se puede expresar de forma única como
suma de una matriz simétrica y otra antisimétrica
4 Dada una matriz cualquiera A, AAt y At A son matrices simétricas

Ejercicio:
1 Demuestra los cuatro teoremas anteriores
2 1


2 Expresa la matriz −3 5 como suma de una matriz simétrica y otra
antisimétrica
2 3 
3 Halla una matriz simétrica a partir de 1 0
4 −1

G1954: Álgebra y Geometría Matrices 13/37

 Determinantes

D ETERMINANTES

A toda matriz cuadrada A se le puede asociar un número llamado

determinante, que se representa como |A|. Por ejemplo,
23 = −10
41

Cálculo del determinante:

Dependiendo del tamaño y estructura de la matriz, se emplean
distintos métodos:
Método de Sarrus
Método de los adjuntos
Método de los pivotes

G1954: Álgebra y Geometría Matrices 14/37

 Determinantes

D ETERMINANTES

Método de Sarrus
Suele emplearse para el cálculo de determinantes de orden dos o tres.

Ejercicio: Calcula los siguientes determinantes

2 −1 1
2 −3
0 1 −2
1 1
−2 0 3
G1954: Álgebra y Geometría Matrices 15/37
 Determinantes

D ETERMINANTES
Método de los adjuntos (o cofactores)
Suele emplearse para el cálculo de determinantes de orden superior a tres.

D EFINICIONES
Menor: Dada una matriz cuadrada A de orden n, se llama menor (o menor complementario)
del elemento aij al determinante de orden n − 1 que resulta de suprimir la fila i y la columna
j en A. Lo llamaremos mij . Adjunto (del elemento aij ) → Aij = (−1)i+j mij

A partir de estas definiciones se puede calcular el determinante de orden n como la suma de

los productos de los elementos de una fila o columna por sus respectivos adjuntos.
desarrollo por fila desarrollo por columna

1 0 2 0
1 2 0 1
Ejercicio: Calcula el determinante
−1 1 4 −1
3 −1 −3 −2
G1954: Álgebra y Geometría Matrices 16/37
 Determinantes

D ETERMINANTES

Método de los pivotes

T EOREMA
Si a los elementos de una fila o columna se suman los correspondientes de otras
paralelas multiplicados por un número, el valor del determinante no varía

Basándonos en esta propiedad, podemos ir operando para obtener un determinante

cuyo valor sea el mismo pero en el que todos los elementos salvo uno de una fila o
columna sean nulos. A partir de ahí, se puede aplicar fácilmente el método de los
adjuntos.

Ejercicio: Calcula el siguiente determinante:

1 2 1 2 1
0 0 1 1 1
1 1 0 0 0
0 0 1 1 2
1 2 2 1 1
Sugerencia: Primera operación: F3 − F1 , segunda operación: F5 − F1
G1954: Álgebra y Geometría Matrices 17/37
 Determinantes

D ETERMINANTES
Propiedades:
Si todos los elementos de una fila o columna son nulos, el determinante es
nulo
Si una fila o columna es combinación lineal de las demás, el determinante es
nulo. Lógicamente, si hay dos filas o columnas iguales, el determinante es nulo
Si intercambiamos dos filas o columnas, el valor del determinante cambia de
signo
Para multiplicar un número por un determinante, se multiplica dicho número
por los elementos de una (y sólo una) fila o columna cualquiera
El determinante de una matriz triangular es igual al producto de los elementos
de la diagonal
|At | = |A|
∀ λ ∈ R, |λA| = λn |A|, donde n es el orden de la matriz A
|AB| = |A||B|
1
|A−1 | = |A|
, donde A−1 es la matriz inversa (la veremos a continuación)

G1954: Álgebra y Geometría Matrices 18/37

 Matriz inversa

M ATRIZ INVERSA

Dada una matriz cuadrada A decimos que tiene inversa, A−1 , si

AA−1 = A−1 A = I

Matriz regular (o invertible): Aquella que tiene inversa

Matriz singular: No tiene inversa

Propiedades:
Si existe, A−1 es única
(A−1 )−1 = A
(AB)−1 = B −1 A−1
1 −1
(λA)−1 = λ
A , ∀ λ 6= 0 ∈ R
(At )−1 = (A−1 )t

G1954: Álgebra y Geometría Matrices 19/37

 Matriz inversa

C ÁLCULO DE LA MATRIZ INVERSA

D EFINICIÓN
Matriz adjunta: Dada una matriz cualquiera A, se llama adjunta, Adj(A),
a la matriz que resulta de sustituir cada elemento de A por su
correspondientes adjunto

A partir de la matriz adjunta se calcula la inversa como

(Adj(A))t
A−1 = , A tiene inversa si y sólo si |A| 6= 0
|A|
Ejercicio: Calcula la inversa de la siguiente matriz:
 
1 −1 −3
2 −2 1 
−1 3 −1

G1954: Álgebra y Geometría Matrices 20/37

 Matriz inversa

M ATRICES ORTOGONALES

Una matriz cuadrada A es ortogonal si

AAt = At A = I

Propiedades:
La inversa de una matriz ortogonal es su transpuesta
At = A−1 , y es ortogonal también
El determinante de una matriz ortogonal es 1 o −1
El producto de matrices ortogonales es ortogonal

Ejercicio: Demuestra las tres propiedades anteriores.

G1954: Álgebra y Geometría Matrices 21/37

 Matrices elementales

M ATRICES ELEMENTALES

Se llama matriz elemental a una matriz cuadrada que resulta de

efectuar una operación elemental sobre una fila o columna en la
matriz identidad.
O PERACIONES ELEMENTALES

Tipo 1: Cambiar entre sí dos filas o columnas

Tipo 2: Multiplicar una fila o columna por un número real λ 6= 0
Tipo 3: Sumar a una fila o columna, otra fila o columna
multiplicada por un número real λ 6= 0

Ejemplos de matrices elementales:

  F (−3)  1   F  

F1,2 0 1
100 2 0 0 100 3,2 (2) 100
10 −→ 010 −→ 0 −3 0 010 −→ 010
01 10 001 0 0 1 001 021

En este curso, trabajaremos siempre con operaciones elementales

por filas. Representaremos las matrices elementales con la letra E
G1954: Álgebra y Geometría Matrices 22/37
 Matrices elementales

M ATRICES ELEMENTALES

T EOREMA
Si en una matriz A efectuamos una operación elemental por filas la
matriz que obtenemos es EA, donde E es la matriz elemental
resultante de efectuar la misma operación elemental (sobre la matriz
identidad)
Ejemplo:
1 0 −1
2 
Partimos de la matriz A = 1 −2 1 −2
3 −1 0 1

1 0 −1 2 1 0 −1
2 F  
2,1 (−2)
1 −2 1 −2 −→ −3 −4 1 0
3 −1 0 1 3 −1 0 1
 F  
100 2,1 (−2) 1 00
I3 = 010 −→ −2 1 0 =E
001 0 01
0 −1 2 1 0 −1
 1 00
2 1
  
EA = −2 1 0 1 −2 1 −2 = −3 −4 1 0
0 01 3 −1 0 1 3 −1 0 1

G1954: Álgebra y Geometría Matrices 23/37

 Matrices elementales

M ATRICES ELEMENTALES

O PERACIONES ELEMENTALES INVERSAS

Son aquellas que anulan la acción de una operación elemental

Operación elemental (por filas) Operación elemental inversa (por filas)

Cambiar la fila i por la j Cambiar la fila j por la i
1
Multiplicar una fila por λ 6= 0 Multiplicar una fila por λ
Sumar a la fila i, la j mutiplicada por λ 6= 0 Sumar a la fila i, la j mutiplicada por −λ

Ejemplos:
12 0  F   12 0
 12 0
 12 F 12 
1,3 0 0 −1 F2 (2) 0 2,3 (−3) 0
23 1 −→ 23 1 23 1 −→ 46 2 23 1 −→ 23 4
0 0 −1 12 0 0 0 −1 0 0 −1 0 0 −1 0 0 −1
1
 F 12  12 F 12  12 F 12 
0 0 −1 3,1 0 0 2( 2 ) 0 0 2,3 (3) 0
23 1 −→ 23 1 46 2 −→ 23 1 23 4 −→ 23 1
12 0 0 0 −1 0 0 −1 0 0 −1 0 0 −1 0 0 −1

G1954: Álgebra y Geometría Matrices 24/37

 Matrices elementales

M ATRICES ELEMENTALES

T EOREMA
Cuando aplicamos una operación elemental sobre I, obtenemos una
matriz elemental E. Cuando aplicamos una operación elemental
inversa sobre I, obtenemos la inversa de E, es decir, E −1 . Por tanto,
toda matriz elemental E tiene inversa, que es también una matriz
elemental
Ejercicios:
1 Dadas las matrices elementales 3 × 3 que se obtienen de realizar las
operaciones elementales F1,3 , F2 (2) y F2,3 (−3), halla sus matrices inversas
1 −1


2 Dada la matriz A = 0 2 :
halla las matrices elementales tales que E2 E1 A = I
halla las matrices E1−1 y E2−1
escribe A como producto de matrices elementales
escribe A−1 como producto de matrices elementales
G1954: Álgebra y Geometría Matrices 25/37
 Matrices elementales

M ATRICES EQUIVALENTES

T EOREMA
Si partiendo de una matriz A podemos llegar a otra B mediante
operaciones elementales (por filas), también podremos volver a A
desde B realizando las operaciones elementales inversas en orden
inverso. Se dice entonces que A y B son equivalentes (por filas)

Ek Ek−1 · · · E2 E1 A = B ⇔ A = E1−1 E2−1 · · · Ek−1

−1
Ek−1 B

Ejercicio:
Demuestra que las matrices A y B son equivalentes por filas
0 −1
2  1 2 3 
A= 1 2 3 B= 2 0 −1
5 −2 1 4 −4 −2

Sugerencia: Aplícale a A las operaciones elementales F1,2 y F3,1 (−1)

G1954: Álgebra y Geometría Matrices 26/37

 Matrices elementales

C ÁLCULO DE LA INVERSA MEDIANTE OPERACIONES ELEMENTALES

T EOREMA
Las operaciones elementales que nos sirven para convertir una
matriz cuadrada A en I, efectuadas sobre I, nos dan A−1

Ejercicio:
Calcula, mediante operaciones elementales, la inversa de la siguiente
matriz (si existe)
 −1 1 1 
A= 0 10
1 01

Procedimiento:
 Comienza formando
 la matriz (A|I), en este caso
−1 1 1 1 0 0
 0 1 0 0 1 0 . A partir de ahí, aplica las operaciones elementales
1 0 1 0 0 1
1
necesarias: por ejemplo F1 (−1), F3,1 (−1), F3,2 (−1), F3 ( 2 ), F1,3 (1), F1,2 (1). Acabarás
−1
llegando a la matriz (I|A )

G1954: Álgebra y Geometría Matrices 27/37

 Formas escalonada y reducida

F ORMAS ESCALONADA Y REDUCIDA DE UNA MATRIZ

F ORMA ESCALONADA

Se llama forma escalonada por filas de una matriz A a otra matriz que se
obtiene a partir de A mediante operaciones elementales y que verifica:

Si tiene filas cuyos elementos son todos nulos, están en las filas
inferiores

El primer elemento distinto de cero de una fila (empezando por la

izquierda) se llama pivote, y a su columna, columna pivotal

Dadas dos filas sucesivas, el pivote de la segunda fila está más a la

derecha que el de la primera

Ejercicio: Di si son formas escalonadas o no las siguientes

145 1067
1045 256 023






12 000 0312 2350 000 001 0248
03 021 0000 0001 007 002 0013
000 0052

G1954: Álgebra y Geometría Matrices 28/37

 Formas escalonada y reducida

F ORMAS ESCALONADA Y REDUCIDA DE UNA MATRIZ

En la práctica, para obtener una forma escalonada de una matriz dada se utiliza la
eliminación gaussiana:
1 Buscamos en la primera columna un elemento distinto de cero, que llevaremos a la
primera fila. Este elemento será el primer pivote. A continuación, mediante operaciones
elementales, haremos ceros por debajo de él
2 Buscamos en la segunda columna un elemento distinto de cero en la segunda o demás
filas inferiores. Operamos hasta tener un segundo pivote en esa fila, y hacemos ceros
por debajo de él
3 Seguimos recorriendo el resto de columnas hasta no encontrar más pivotes

Ejemplo: Halla una forma escalonada de la siguiente matriz

 
1 2 1 2 1 0
1 2 4 −1 4 3 
−2
2 4 5 1 5

1 3 1 2 6 0

Sugerencia: Puedes aplicar, en este orden, las operaciones elementales F2,1 (−1), F3,1 (−2), F4,1 (−1), F4,2 ,
F4,3 (−1)

G1954: Álgebra y Geometría Matrices 29/37

 Formas escalonada y reducida

F ORMAS ESCALONADA Y REDUCIDA DE UNA MATRIZ

F ORMA REDUCIDA

Se llama forma escalonada reducida (o simplemente reducida) por

filas de la matriz A a toda matriz escalonada obtenida mediante
operaciones elementales por filas sobre A en la que los pivotes son
1 y los demás elementos de la columna pivotal 0.

Ejercicio: Di si son formas escalonadas reducidas o no las siguientes

1 0 3 1 3 0 0 2 5 6
012
000
0010
0001
014
003
( 10 20 01 34 )

Nota: Dependiento del orden de actuación y las operaciones realizadas sobre la matriz
original, se pueden obtener diferentes formas escalonadas (de hecho, infinitas). Sin embargo,
la forma reducida es única.

G1954: Álgebra y Geometría Matrices 30/37

 Formas escalonada y reducida

R ANGO

D EFINICIÓN
El rango de una matriz es el número de filas con algún elemento
distinto de cero que hay en cualquier forma escalonada por filas, o, lo
que es lo mismo, el número de columnas pivotales. Como veremos
más adelante, esta definición equivale a decir que el rango de una
matriz es el número de filas (o columnas) que son linealmente
independientes. Por último, también se puede ver el rango como el
orden de la mayor submatriz cuyo determinante sea distinto de cero.

Ejercicio: ¿Cuál es el rango de las siguientes matrices?

 
  1 2 1 −1 1 0    
1 −1 0 5 1 2 3 2 4 9
0 −1 4 3
0 3 0 0 5 10 15 3 6 7
0 0 0 0
0 0

G1954: Álgebra y Geometría Matrices 31/37

 Factorización de matrices

FACTORIZACIÓN LU
T EOREMA
Dada una matriz A, se puede encontrar mediante la aplicación de operaciones elementales
del tipo 3 (quedarían excluídas el intercambio de filas y la multiplicación de una fila por un
escalar ) una factorización de la forma A = LU , donde L es una matriz triangular inferior
con unos en la diagonal y U es una matriz triangular superior, en la que la primera fila coincide
con la primera fila de A  
  a1,1 a1,2 ··· a1,n
1 0 ··· 0
0 ∗ ··· ∗
∗ 1 ··· 0  
  U = . . . .
L=  .. . . . . . . . 
. . .
 . . . .
. . . .
0 0 ··· ∗
∗ ∗ ··· 1
Si A es invertible, la factorización LU es única (en caso contrario, podría no serlo). Esta
factorización es útil para el cálculo de determinantes e inversas (en el caso de matrices A
cuadradas), así como para la resolución de sistemas de ecuaciones lineales (lo veremos más
adelante)

Razonamiento: Hemos visto que, a partir de A, se puede llegar a una matriz escalonada U aplicando operaciones
elementales, por lo que:
−1
Ek Ek−1 · · · E2 E1 A = U ⇒ A = (Ek Ek−1 · · · E2 E1 ) U
| {z }
−1 −1 −1 −1
L=E E ···E E
1 2 k−1 k

Por tanto, en teoría el problema se reduce al cálculo de una serie de matrices elementales inversas
G1954: Álgebra y Geometría Matrices 32/37
 Factorización de matrices

FACTORIZACIÓN LU
Sin embargo, en la práctica no hace falta calcular estas matrices, sino que se emplea el
algoritmo de los multiplicadores cambiados de signo, que se ilustra a continuación. Nótese que
en este algoritmo L es inicialmente la matriz identidad.
Ejemplo: Halla la factorización LU de la siguiente matriz
 
−2 4 −1
A= 4 −5 4 
−6 −3 −14
Proceso de construcción de la matriz U :
−2 4 −1
 −2 4 −1

F2,1 (2) F3,1 (−3)
A= 4 −5 4 −
−−−−−
→ 0 3 2 −−−−−−−→
−6 −3 −14 −6 −3 −14
−2 4 −1
 −2 4 −1

F3,2 (5)
0 3 2 −
−−−−−
→ 0 3 2 =U
0 −15 −11 0 0 −1
Proceso de construcción de la matriz L:
1 0 0
  1 0 0
  1 0 0

F2,1 (2) F3,1 (−3) F3,2 (5)
I = 0 1 0 −
−−−−−
→ −2 1 0 −−−−−−−→ −2 1 0 −
−−−−−
→
0 0 1 0 0 1 3 0 1
 1 0 0

−2 1 0 =L
3 −5 1

G1954: Álgebra y Geometría Matrices 33/37

 Factorización de matrices

FACTORIZACIÓN LU
Hay matrices que no se pueden factorizar mediante el método que acabamos de ver. Bastaría
por ejemplo con que el elemento a11 de A fuese un 0. Para solventar este tipo de situaciones
se cuenta con una generalización de la factorización A = LU , la P A = LU (que no es
única), en la que P es una matriz de permutaciones. A continuación se ilustra el
procedimiento para obtener una factorización P A = LU (en el que sólo se permiten
operaciones elementales de los tipos 1 y 3) para un caso en el que se realiza un único
intercambio de filas. Nótese que, inicialmente, tanto P como L son la matriz identidad.
Ejemplo: Halla una factorización P A = LU de la siguiente matriz
0 0 2

A= 3 5 −6
1 4 −2
Proceso de construcción de la matriz U :
0 0 2
 1 4 −2
 1 4 −2

F1,3 F2,1 (−3)
A= 3 5 −6 −
−−−
→ 3 5 −6 −−−−−−−→ 0 −7 0 =U
1 4 −2 0 0 2 0 0 2
Proceso de construcción de la matriz P :
1 0 0
 0 0 1
 0 0 1

F1,3 F2,1 (−3)
I = 0 1 0 −
−−−
→ 0 1 0 −−−−−−−→ 0 1 0 =P
0 0 1 1 0 0 1 0 0
Proceso de construcción de la matriz L:
1 0 0
 1 0 0
 1 0 0

F1,3 F2,1 (−3)
I = 0 1 0 −
−−−
→ 0 1 0 −−−−−−−→ 3 1 0 =L
0 0 1 0 0 1 0 0 1

G1954: Álgebra y Geometría Matrices 34/37

 Factorización de matrices

FACTORIZACIÓN LU

Veamos a continuación un caso en el que se realizan dos intercambios de filas para

llegar a una factorización P A = LU .
Ejemplo: Halla una factorización P A = LU de la siguiente matriz
 
0 0 2
A = −1 5 −2
3 6 7
Proceso de construcción de la matriz U :
 0 0 2
 −1 5 −2
 −1 5 −2

F1,2 F3,1 (3) F2,3
A= −1 5 −2 −
−−−
→ 0 0 2 −
−−−−−
→ 0 0 2 −
−−−
→
3 6 7 3 6 7 0 21 1
−1 5 −2

0 21 1 =U
0 0 2
Proceso de construcción de la matriz P :
1 0 0
 0 1 0
 0 1 0
 0 1 0

F1,2 F3,1 (3) F2,3
I = 0 1 0 −
−−−
→ 1 0 0 −
−−−−−
→ 1 0 0 −
−−−
→ 0 0 1 =P
0 0 1 0 0 1 0 0 1 1 0 0
Proceso de construcción de la matriz L:
1 0 0
 1 0 0
  1 0 0
  1 0 0

F1,2 F3,1 (3) F2,3
I = 0 1 0 −
−−−
→ 0 1 0 −
−−−−−
→ 0 1 0 −
−−−
→ −3 1 0 =L
0 0 1 0 0 1 −3 0 1 0 0 1

G1954: Álgebra y Geometría Matrices 35/37

 Factorización de matrices

FACTORIZACIÓN DE C HOLESKY

T EOREMA
Dada una matriz real A, simétrica y definida positiva*, se puede
encontrar una matriz triangular inferior L tal que A = LLt
*Criterio de Sylvester: Una matriz simétrica se dice definida positiva cuando todos sus menores principales (incluyendo su
propio determinante) son positivos.
 
1 2
Ejemplo: Halla una factorización de Cholesky de la matriz A =
2 8
 
l11 0
Se trata de buscar una matriz L = tal que
l21 l22
      
1 2 l11 0 l11 l21 2
l11 l11 l21
= = 2 + l2
. La factorización de
2 8 l21 l22 0 l22 l21 l11 l21 22
Cholesky no es única. En este caso, una de las soluciones posibles que obtendríamos
igualando
 términos
 sería  {l11
 = 1, l21 = 2, l22 = ±2}. Se puede comprobar, por ejemplo,
1 2 1 0 1 2
que = . Otra opción igualmente válida sería
2 8 2 2 0 2
| {z } | {z }
L Lt
{l11 = −1, l21 = −2, l22 = ±2}
G1954: Álgebra y Geometría Matrices 36/37
 Factorización de matrices

FACTORIZACIÓN DE C HOLESKY

En general, los elementos de la matriz triangular inferior L, lij , se pueden

obtener de acuerdo a las siguientes fórmulas, donde aij son los términos
de la matriz A:
 v
 u j−1
X
u
2

 ljj = ± ajj − ljk

 t

L: k=1
j−1
 !
 X
1
 lij = ljj aij − lik ljk , i > j



k=1

Ejercicio: Halla una factorización de Cholesky de la matriz

 
4 2 1
A = 2 2 0
1 0 3

Al igual que la factorización LU , la factorización de Cholesky resulta útil para el cálculo de

determinantes, inversas y para la resolución de sistemas de ecuaciones lineales (lo veremos
más adelante)
G1954: Álgebra y Geometría Matrices 37/37