ELIMINACIÓN DE
GAUSS
César Guerra
Eliminación Gaussiana
Operaciones elementales de fila en una matriz
Multiplicar un renglón por un número distinto de cero.
Sumar el múltiplo de un reglón a otro renglón.
Intercambiar dos renglones.
Ejemplo:
Sea el sistema de ecuaciones:
2 x1 + 4 x2 + 6 x3 =
18
4 x1 + 5 x2 + 6 x3 =
24
3 x1 + x2 − 2 x3 =
4
Encontrar la solución del sistema.
Eliminación de Gauss
Llevar la matriz de coeficientes aumentada s una matriz por filas escalonada.
2 4 6 18 1 2 3 9 R →R −4 R 1 2 3 9
R1 → 2 R1
1
R3 → R3 −3 R1
2 2 1
4 5 6 24 → 4 5 6 24 → 0 −3 − 6 −12
3 1 −2 4 3 1 −2 4 0 −5 −11 −23
1 2 3 9 1 2 3 9
R2 →− 13 R2 R3 → R3 +5 R2
→ 0 1
2 4 → 0 1 2 4
0 −5 −11 −23 0 0 −1 −3
1 2 3 9
Y luego realizar sustituciones R3 →− R3
→ 0 1 2 4
hacia atrás para hallar las 0 0 1 3
soluciones.
Eliminación de Gauss Jordan
Llevar la matriz de coeficientes aumentada a una matriz por filas escalonada
reducida:
2 4 6 18 1 2 3 9 R →R −4 R 1 2 3 9
R1 → 12 R1 R3 → R3 −3 R1
2 2 1
4 5 6 24 → 4 5 6 24 → 0 −3 − 6 −12
3 1 −2 4 3 1 −2 4 0 −5 −11 −23
1 2 3 9 R → R − 2 R 1 0 −1 1
R2 →− 13 R2 R13 → R13 +5 R22
→ 0 1
2 4 → 0 1 2 4
0 −5 −11 −23 0 0 −1 −3
1 0 −1 1 R → R + R 1 0 0 4
R3 →− R3 R2 → R2 − 2 R3
1 1 3
→ 0 1 2 4 → 0 1 0 −2
0 0 1 3 0 0 1 3
La solución viene dada por x1=4, x2=-2, x3=3
Eliminación de Gauss
Mediante la eliminación de Gauss o Gauss-Jordan se puede resolver un
sistema lineal para varios vectores b y para la inversa de A.
Considérese el conjunto de sistemas
A ⋅=
x1 b1 A ⋅=
x2 b2 A ⋅=
x3 b3 A ⋅=
Y 1
Las soluciones se obtienen a partir de la matriz aumentada:
[ A | b1 | b 2 | b3 | I ]
Después de la eliminación de Gauss o Gauss-Jordan tenemos:
∆ sup | b 1 | b 2 | b 3 | Y
[I | x1 | x2 | x3 | Y ]
y luego sustitución hacia atrás
Fase de eliminación
En general, consideramos el caso de n ecuaciones con n incógnitas y partimos
de la matriz aumentada de coeficientes:
a11 a12 a1n b1
a a22 a2 n b2
A = 21
an1 an 2 ann bn
Hay (n-1) pasos de eliminación hacia adelante. Eliminamos el primer
coeficiente de la fila i≥2 (y realizamos la misma operación para j≥2):
A′ ( i,:) =
0 ai 2 − ai1 a12 ain − ai1 a1n bi − ai1 b1
a a a
11 11 11
A′ ( i,:) = [ 0 ai′2 ain′ bi′]
Fase de eliminación
Después de esta transformación tenemos:
a11 a12 a1n b1
0 ′ a2′ n
a22 b2′
A' =
0 an′ 2 ann
′ bn′
Y repetimos el procedimiento, eliminando el primer coeficiente de la fila i≥3
(y realizando la misma operación para j≥3):
A′′ ( i,:) =0 0 ai′3 − aa22i′2 a23
′
′ ain′ − aa22i′′2 a2′ n bi′ − aa22i′′2 b2′
A′′ ( i,:) = [ 0 0 ai′′3 ain′′ bi′′]
Fase de eliminación
Siguiendo estas transformaciones (n-1) pasos:
a11 a12 a1n b1
0 ′
a22 a2′ n b2′
A= 0 0 ′′
a33 a3′′n b3′′
0 0 0 ( n −1)
ann bn( n −1)
Donde los coeficientes con i, j ≥ k+1 cumplen la relación:
( −)k 1
a
aij( k ) aij( k − ) − ik( k −1) akj( k − ) ¿te convence?
1 1
=
akk
Fase de substitución
La fase de substitución consiste en resolver las ecuaciones resultantes
de la fase de eliminación, comenzando con la última ecuación:
bn( n 1)
−
xn = ( n −1)
ann
y sustituyendo hacia atrás:
n
( i −1)
bi − ∑ aij(i −1) x j
x=
i
j = i +1
( i −1)
i
= {n − 1, ,1}
a ii
Determinante de una Matriz Cuadrada
Si un múltiplo de una fila de A se le suma otra fila de A, dando una
nueva matriz B, entonces: det(A) = det(B)
El determinante de una matriz triangular A está dada por:
n
det(A) = a11 × a12 × ann = ∏ aii
i =1
Entonces, usando solo eliminación de Gauss, sin reducir los pivotes a la
unidad, tendríamos:
An×n ⇒ ∆ sup,n×n
det(=
A) det(∆ sup,n×n )