Algoritmo para Sistemas de Ecuaciones
Algoritmo para Sistemas de Ecuaciones
Algoritmo para Sistemas de Ecuaciones
Resumen
1. Introducción
117
118 Daniel Gómez-Garcı́a y Humberto Madrid de la Vega
5 x1 + 4 x 2 = 2
.
2 x1 + 3 x 2 = 1
Solución
• Primera etapa, de eliminación:
La intención es generar una matriz triangular superior de elementos
enteros, a partir de la matriz aumentada del sistema. Aplicando la
eliminación Gaussiana [8, p. 3], se forman las matrices equivalentes
5 4 2 5 4 2 5 4 2
∼ 2 2 = 7 1 .
2 3 1 0 3− (4) 1 − (2) 0
5 5 5 5
Una manera de convertir estos números fraccionarios en números en-
teros consiste en multiplicarlos por el denominador de los multipli-
cadores Gaussianos:
⎡ ⎤
5 4
5 4 2 5 4 2 2
∼ ∼ ⎣ 5 4 5 2 ⎦
2 3 1 0 5(3) − 2(4) 5(1) − 2(2) 0
2 1
2 3
5 4 2
= .
0 7 1
Sistemas de ecuaciones lineales con coeficientes enteros 119
5 4 2
5 4 5 2
2 3 1 con los cálculos |A| = = 7, = 1.
2 3 2 1
7 1
y para operar en esta etapa con aritmética entera, basta con encontrar
múltiplos de x, y = |A| x, —los numeradores de la regla de Cramer
que como se sabe son enteros cuando la matriz aumentada original es
entera— obtenidos por sustitución regresiva en el sistema modificado
5 4 7(2) 5 4 14
=
0 7 7(1) 0 7 7
3 x1 + 2 x 2 − x3 = 1
x1 − x2 + 2 x 3 = 0 .
4 x1 − 2 x 2 + 3 x 3 = 0
120 Daniel Gómez-Garcı́a y Humberto Madrid de la Vega
Solución
• Primera etapa, de eliminación:
Se forman las matrices equivalentes de elementos enteros:
⎡ ⎤ ⎡ ⎤ ⎡ ⎤
3 2 −1 1 3 2 −1 1 3 2 −1 1
⎣ 1 −1 ⎦ ⎣
2 0 ∼ 0 −5 7 −1 ⎦ ∼ ⎣ 0 −5 7 −1 ⎦ ,
4 −2 3 0 0 −14 13 −4 0 0 33 6
⎡ ⎤
1
1 ⎣ ⎦
sistema es x = 5 . Como era de esperarse, los numeradores de
11
2
Cramer son números enteros y la solución es racional.
En el proceso de sustitución, también se ha usado aritmética entera.
Una forma de organizar los cálculos, cuando se hacen con papel y
lápiz, consiste en preparar un arreglo tabular compacto que contiene
las cuentas de la primera etapa y la columna modificada. La sustitu-
ción regresiva se hace directamente para encontrar los numeradores de
Cramer y luego la solución:
A b |A| b
11 − 2(5) + 1(2) 1
3 2 −1 1 11 y1 = = 1, x1 =
3 11
1 −1 2 0
4 −2 3 0
−11 − 7(2) 5
−5 7 −1 −11 y2 = = 5, x2 =
−5 11
−14 13 −4
22 2
|A| = 11 2 22 y3 = = 2, x3 =
11 11
En los numeradores
de esta solución se
A b |A| b han anotado directa-
1
3 2 −1 1 11 x1 = mente los resultados
11
1 −1 2 0 de la sustitución re-
4 −2 3 0 gresiva, y.
5
−5 7 −1 −11 x2 =
11
−14 13 −4 Los denominadores
2
|A| = 11 2 22 x3 = de la solución com-
11 parten el valor de
|A| = 11.
122 Daniel Gómez-Garcı́a y Humberto Madrid de la Vega
2 x1 + x2 − 2 x3 + 2 x4 = 1
3 x1 + x2 − x3 − x4 = 2
.
5 x1 + 2 x3 + x4 = 3
x1 − x 2 + 3 x3 + 4 x4 = −1
Solución
A partir de ésta, se obtienen ma-
⎡ ⎤ trices equivalentes, cuyos elemen-
2 1 −2 2 1 tos se calculan con determinantes
⎢ 3 1 −1 −1 2 ⎥ de segundo orden, divididos en-
⎢ ⎥
⎣ 5 0 2 1 3 ⎦ tre su pivote anterior (excepto en
1 −1 3 4 −1 el primer paso. Optativamente, el
pivote anterior inicial se define
como 1.)
⎡ ⎤
2 1 −2 2 1
⎢ 2 −2 2 ⎥
⎢ 0 2 1 2 2 1 ⎥
⎢ 3 −1 3 −1 ⎥
⎢ 3 1 3 2 ⎥
⎢ ⎥
∼⎢
⎢ 0
2 1
2 −2
2 2
2 1 ⎥=
⎥
⎢ 5 0 5 2 5 1 5 3 ⎥
⎢ ⎥
⎢ ⎥
⎣ 2 1 2 −2 2 2 2 1 ⎦
0
1 −1 1 3 1 4 1 −1
⎡ ⎤
2 1 −2 2 1
⎢ 0 −1 4 −8 1 ⎥
=⎢
⎣ 0
⎥
−5 14 −8 1 ⎦
0 −3 8 6 −3
⎡ ⎤
2 1 −2 2 1
⎢ ⎥
⎢ 0 −1 4 −8 1 ⎥
⎢ ⎥
⎢ ⎥
⎢ −1 4 −1 −8 −1 1 ⎥
⎢ ⎥
⎢ −5 14 −5 −8 −5 1 ⎥
∼⎢ ⎥=
⎢ 0 0 ⎥
⎢ 2 2 2 ⎥
⎢ ⎥
⎢ −1 4 −1 −8 −1 1 ⎥
⎢ ⎥
⎣ −3 8 −3 6 −3 −3 ⎦
0 0
2 2 2
Sistemas de ecuaciones lineales con coeficientes enteros 123
⎡ ⎤
2 1 −2 2 1
⎢ 0 −1 4 −8 1 ⎥
=⎢⎣ 0
⎥
0 3 −16 2 ⎦
0 0 2 −15 3
⎡ ⎤
2 1 −2 2 1
⎢ ⎥ ⎡ ⎤
⎢ 0 −1 4 −8 1 ⎥ 1 −2
⎢ ⎥ 2 2 1
⎢ ⎥ ⎢
⎢ −16 ⎥ ⎢ 0 −1 4 −8 1 ⎥
⎥.
∼⎢ 0 0 3 2 ⎥=⎣
⎢ ⎥ 0 0 3 −16 2 ⎦
⎢ 3 −16 3 2 ⎥
⎢ ⎥ 0 0 0 13 −5
⎣ 2 −15 2 3 ⎦
0 0 0
−1 −1
Los cálculos de las matrices equivalentes anteriores, añadiendo el
vector del lado derecho modificado para obtener los numeradores de la
regla de Cramer, toman la forma alternativa más compacta:
A b |A| b
13 − (−45) + 2(−18) − 2(−5)
2 1 −2 2 1 13 y1 = = 16
2
3 1 −1 −1 2
5 0 2 1 3
1 −1 3 4 −1
13 − 4(−18) + 8(−5)
−1 4 −8 1 13 y2 = = −45
−1
−5 14 −8 1
−3 8 6 −3
26 + 16(−5)
3 −16 2 26 y3 = = −18
3
2 −15 3
−65
|A| = 13 −5 −65 y4 = = −5
13
⎡ ⎤
16
⎢ ⎥
⎢ ⎥
⎢ −45 ⎥
1 ⎢⎢ ⎥
x= ⎥.
⎢
13 ⎢ ⎥
−18 ⎥
⎢ ⎥
⎣ ⎦
−5
1 1 1 1
x1 + 2 x2 + 3 x3 + 4 x4 = 4
1 1 1 1
2 x1 + 3 x2 + 4 x3 + 5 x4 = 0
1 1 1 1
.
3 x1 + 4 x2 + 5 x3 + 6 x4 = 0
1 1 1 1
4 x1 + 5 x2 + 6 x3 + 7 x4 = 0
124 Daniel Gómez-Garcı́a y Humberto Madrid de la Vega
Solución
Para operar con el algoritmo, primero se multiplican los renglones por
12, 60, 60 y 420, respectivamente, resultando el arreglo de elementos
enteros A∗ x = b∗ :
A∗ b∗ |A∗ | b∗
12
12 6 4 3 3 9 x1 = =4
3
30 20 15 12 0
20 15 12 10 0
105 84 70 60 0
−90
60 60 54 −90 −270 x2 = = −30
3
60 64 60 −60
378 420 405 −315
180
20 30 150 450 x3 = = 60
3
210 324 1260
−105
|A∗ | = 3 −105 −315 x4 = = −35
3
3 Cuando las incógnitas x tienen
|A| = |H4 | = valores enteros, pueden obte-
12 · 60 · 60 · 420
nerse del lado derecho b, sin
1 necesidad de usar su valor
|A| = .
6 048 000 modificado |A|b.
y ası́ sucesivamente.
Para mayor claridad se han enmarcado los elementos pivotales.
Cuadro 1: Eliminaciones
Gaussiana para B = A.
y del algoritmo
(1) (0) (1) (4) (3) (4)
Observe que: bij × a00 = aij , . . . , bij × a33 = aij .
B A
7 3 1 4 7 3 1 4
3 4 2 1 (0) 3 4 2 1
× a00 =
1 2 2 1 1 2 2 1
1 1 2 2 1 1 2 2
2,7143 1,5714 −0,7143 19 11 −5
(1)
1,5714 1,8571 0,4286 × a11 = 11 13 3
0,5714 1,8571 1,4286 4 13 10
0,9474 0,8421 (2) 18 16
× a22 =
1,5263 1,5789 29 30
(3)
0,2222 × a33 = 4
(0) (1) (2) (3) (4)
donde: a00 = 1, a11 = 7, a22 = 19, a33 = 18, a44 = 4.
7 3 1 4
7 3 1
(3) (4) 3 4 2 1
δ43 = 3 4 2 = 29,
δ44 = = 4 = |A|.
1 1 2 1 2 2 1
1 1 2 2
(p) (p)
δij = aij (1b)
p = 1, 2, . . . , n − 1, i, j > p :
(p) (p)
(p+1) (p) bip bpj (p)
bij = bij − (p)
, bip = 0
bpp
p = 1, 2, . . . , n − 1, i, j > p :
a11 · · · a1p a1j
.. .. ..
. . .
ap1 · · · app apj
(p+1)
ai1 · · · aip aij δij
(p+1)
bij = = .
a11 · · · a1p (1) (2)
b11 b22 · · · bpp
(p)
.. ..
. .
ap1 · · · app
Sistemas de ecuaciones lineales con coeficientes enteros 129
p = 1, 2, . . . , n − 1, i, j > p :
a11 · · · a1p a1j
.. .. ..
(p+1) . . . = b(1) b(2) · · · b(p) · b(p+1) .
δij = 11 22 pp ij (2a)
ap1 · · · app apj
ai1 · · · aip aij
Referencias
[1] D. Carlson, Ch. Johnson, D. Lay, A. Porter, Editores. “Linear al-
gebra gems assets for undergraduate mathematics”. The Mathe-
matical Association of America, MAA Notes No. 59, Washington.,
2002.