UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ

FACULTAD DE INGENIERÍA MECÁNICA

SOLUCIÓN DE SISTEMA DE ECUACIONES LINEALES
1. INTRODUCCIÓN
Por sistema de ecuaciones lineales se entiende un conjunto de ecuaciones que deben
resolverse simultáneamente y que presentan la siguiente estructura:

𝒂𝟏𝟏 𝒙𝟏 + 𝒂𝟏𝟐 𝒙𝟐 + 𝒂𝟏𝟑 𝒙𝟑 + ⋯ + 𝒂𝟏𝑵 𝒙𝑵 = 𝒃𝟏

𝒂𝟐𝟏 𝒙𝟏 + 𝒂𝟐𝟐 𝒙𝟐 + 𝒂𝟐𝟑 𝒙𝟑 + ⋯ + 𝒂𝟐𝑵 𝒙𝑵 = 𝒃𝟐
𝒂𝟑𝟏 𝒙𝟏 + 𝒂𝟑𝟐 𝒙𝟐 + 𝒂𝟑𝟑 𝒙𝟑 + ⋯ + 𝒂𝟑𝑵 𝒙𝑵 = 𝒃𝟑
. . . . . (1)
. . . . .
. . . . .
𝒂𝑴𝟏 𝒙𝟏 + 𝒂𝑴𝟐 𝒙𝟐 + 𝒂𝑴𝟑 𝒙𝟑 + ⋯ + 𝒂𝑴𝑵 𝒙𝑵 = 𝒃𝑴

Este sistema de M ecuaciones algebraicas lineales con N incógnitas puede escribirse

en forma matricial como:

𝑨𝑴×𝑵 𝑿 = 𝑩

𝒂𝟏𝟏 𝒂𝟏𝟐 𝒂𝟏𝟑 ⋯ 𝒂𝟏𝑵 𝒙𝟏𝟏 𝒃𝟏𝟏

𝒂𝟐𝟏 𝒂𝟐𝟐 𝒂𝟐𝟑 ⋯ 𝒂𝟐𝑵 𝒙𝟏𝟐 𝒃𝟏𝟐
⋮ ⋮ ⋮ ⋮
𝑨= 𝑿= 𝑩= (𝟐)
⋮ ⋮ ⋮ ⋮ ⋮ ⋮
𝒂𝑴𝟏 𝒂𝑴𝟐 𝒂𝑴𝟑 ⋯ 𝒂𝑴𝑵 ⋮ ⋮
[ ] [𝒙𝟏𝑵 ] [𝒃𝟏𝑴 ]
La matriz de coeficientes A se llama matriz del sistema. La matriz formada por A, a la
que se le ha agregado el vector de términos independientes B como última columna,
se le llama la matriz ampliada o matriz aumentada del sistema de ecuaciones, que se
representa por [A | B] y X es el vector de incógnitas.

Antes de proceder a resolver un sistema de ecuaciones es necesario determinar si

dicho sistema tiene o no solución y, en caso de tenerla, cuántas posibles soluciones
tiene. A continuación se presentan las diversas alternativas:

MG. ARTURO H. GAMARRA MORENO

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA MECÁNICA

Si el vector de términos independientes B del sistema dado en (2) es diferente de cero

se dice que el sistema de ecuaciones es no homogéneo y en caso contrario el
sistema es homogéneo.

Sistema compatible o consistente. Es aquél que tiene solución y en este caso se

cumple que (Teorema de Rouché-Frobenius):

𝒓𝒂𝒏𝒈𝒐[𝑨] = 𝒓𝒂𝒏𝒈𝒐[𝑨 | 𝑩]
El rango de una matriz es el número de columnas linealmente independientes.
También es el orden de determinante no nulo de mayor orden que puede obtenerse
de esa matriz.

Sistema incompatible o inconsistente. Es aquél que no tiene solución y cumple la

relación:

𝒓𝒂𝒏𝒈𝒐[𝑨] < 𝑟𝑎𝑛𝑔𝑜[𝑨 | 𝑩]

Sistema determinado. Es un sistema compatible que presenta solución única y en
este caso se verifica que:

𝒓𝒂𝒏𝒈𝒐[𝑨] = 𝑵ú𝒎𝒆𝒓𝒐 𝒅𝒆 𝒊𝒏𝒄ó𝒈𝒏𝒊𝒕𝒂𝒔

Un sistema homogéneo que es determinado tiene únicamente la solución trivial X=0.

Un sistema compatible que presenta infinidad de soluciones se conoce como sistema

indeterminado y se caracteriza por:

𝒓𝒂𝒏𝒈𝒐[𝑨] < 𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑖𝑛𝑐ó𝑔𝑛𝑖𝑡𝑎𝑠

METODOS DIRECTOS E ITERATIVOS PARA RESOLVER SISTEMAS DE

ECUACIONES LINEALES

MÉTODO DE GAUSS-JORDAN

El método de Gauss-Jordan utiliza operaciones con matrices para resolver sistema de

ecuaciones de “n” número de variables. Para aplicar este método solo hay que recordar
que cada operación que se realice se aplicará a toda la fila o a toda la columna en su
caso.

Ejemplo 1): Utilizar el método de Gauss-Jordan hasta que |∈𝒂 | < 𝟏% para aproximar la
solución del siguiente sistema de ecuaciones:

MG. ARTURO H. GAMARRA MORENO

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA MECÁNICA
𝒙 + 𝟐𝒚 + 𝒛 = 𝟑
𝟐𝒙 + 𝟓𝒚 − 𝒛 = 𝟒
𝟑𝒙 − 𝟐𝒚 − 𝒛 = 𝟐

Paso 1: Se forma la matriz aumentada:

𝟏 𝟐 𝟏⋮ 𝟑
[𝟐 𝟓 −𝟏⋮ 𝟒]
𝟑 −𝟐 −𝟏⋮ 𝟐

Nota: El objetivo del método es lograr formar una matriz identidad de esta forma:

𝟏 𝟎 𝟎⋮ 𝒂
[𝟎 𝟏 𝟎⋮ 𝒃]
𝟎 𝟎 𝟎⋮ 𝒄

Donde el sistema tiene la siguiente solución:

𝒙 = 𝒂; 𝒚 = 𝒃; 𝒛=𝒄

Paso 2: Como se busca obtener una diagonal de “1” en la primera fila, ya se tiene un
número 1. Nuestro objetivo ahora será hacer obtener ceros debajo de este número “1”:

Identificamos fila, columna y elemento pivote:

Elemento
pivote 𝟏 𝟐 𝟏⋮ 𝟑 Fila pivote
[𝟐 𝟓 −𝟏⋮ 𝟒]
𝟑 −𝟐 −𝟏⋮ 𝟐

Columna pivote

Ejemplo 1): Utilizar el método de Gauss-Jordan para hallar la solución del siguiente
sistema de ecuaciones:

𝟐𝒙 + 𝟒𝒚 − 𝒛 = −𝟏𝟒
−𝟏𝟎𝒙 + 𝟔𝒚 + 𝟑𝒛 = −𝟓𝟔
𝟓𝒙 − 𝟏𝟐𝒚 + 𝟕𝒛 = 𝟓𝟔

MG. ARTURO H. GAMARRA MORENO

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA MECÁNICA
MÉTODO DE GAUSS-SEIDEL

Es un método iterativo para resolver sistemas de ecuaciones lineales

𝑨𝒙=𝒃

𝒂𝟏𝟏 𝒂𝟏𝟐 𝒂𝟏𝟑 𝒙𝟏 𝒃𝟏

[𝒂𝟐𝟏 𝒂𝟐𝟐 𝒂𝟐𝟑 ] [𝒙𝟐 ] = [𝒃𝟐 ]
𝒂𝟑𝟏 𝒂𝟑𝟐 𝒂𝟑𝟑 𝒙𝟑 𝒃𝟑

Un sistema de ecuaciones puede ser o no ser convergente.

 Comprobar convergencia:
Un sistema de ecuaciones converge si la matriz A es diagonalmente dominante

𝒂𝟏𝟏 𝒂𝟏𝟐 𝒂𝟏𝟑

[𝒂𝟐𝟏 𝒂𝟐𝟐 𝒂𝟐𝟑 ]
𝒂𝟑𝟏 𝒂𝟑𝟐 𝒂𝟑𝟑

Esto se comprueba con las siguientes condiciones:

|𝒂𝟏𝟏 | ≥ |𝒂𝟏𝟐 | + |𝒂𝟏𝟑 |

|𝒂𝟐𝟐 | ≥ |𝒂𝟐𝟏 | + |𝒂𝟐𝟑 |

|𝒂𝟑𝟑 | ≥ |𝒂𝟑𝟏 | + |𝒂𝟑𝟐 |

Si se cumple lo anterior se tiene una matriz diagonalmente dominante y se puede

aplicar este método, si no es así intercambiando las filas se puede hacer que la matriz
sea diagonalmente dominante.

 Se tiene como ecuaciones del sistema:

𝑬𝒄. 𝟏 → 𝒂𝟏𝟏 𝒙𝟏 + 𝒂𝟏𝟐 𝒙𝟐 + 𝒂𝟏𝟑 𝒙𝟑 = 𝒃𝟏

𝑬𝒄. 𝟐 → 𝒂𝟐𝟏 𝒙𝟏 + 𝒂𝟐𝟐 𝒙𝟐 + 𝒂𝟐𝟑 𝒙𝟑 = 𝒃𝟐
𝑬𝒄. 𝟑 → 𝒂𝟑𝟏 𝒙𝟏 + 𝒂𝟑𝟐 𝒙𝟐 + 𝒂𝟑𝟑 𝒙𝟑 = 𝒃𝟑

Despejamos 𝒙𝟏 de la Ec. 1; 𝒙𝟐 de la Ec. 2 y 𝒙𝟑 de la Ec. 3 respectivamente:

𝒃𝟏 − 𝒂𝟏𝟐 𝒙𝟐 − 𝒂𝟏𝟑 𝒙𝟑
𝒙𝟏 =
𝒂𝟏𝟏

𝒃𝟐 − 𝒂𝟐𝟏 𝒙𝟏 − 𝒂𝟐𝟑 𝒙𝟑
𝒙𝟐 =
𝒂𝟐𝟐

𝒃𝟑 − 𝒂𝟑𝟏 𝒙𝟏 − 𝒂𝟑𝟐 𝒙𝟐
𝒙𝟑 =
𝒂𝟑𝟑

MG. ARTURO H. GAMARRA MORENO

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA MECÁNICA

 Comenzamos con las iteraciones:

Como todo método iterativo tenemos nuestros valores iniciales:

𝒙𝟏 = 𝟎; 𝒙𝟐 = 𝟎; 𝒙𝟑 = 𝟎

Obteniendo el error porcentual en cada una de las iteraciones hasta aproximarnos a la

solución del sistema hasta llegar al error porcentual que nos indique el problema.

𝒙 𝒏 − 𝒙𝒏−𝟏
|∈𝒂 | = | | × 𝟏𝟎𝟎%
𝒙𝒏

Ejemplo 1): Utilizar el método de Gauss-Seidel hasta que |∈𝒂 | < 𝟏% para aproximar la
solución del siguiente sistema de ecuaciones:

𝟏𝟎𝒙𝟏 + 𝟎𝒙𝟐 − 𝒙𝟑 = −𝟏
𝟒𝒙𝟏 + 𝟏𝟐𝒙 𝟐 − 𝟒𝒙 𝟑 = 𝟖
𝟒𝒙𝟏 + 𝟒𝒙𝟐 + 𝟏𝟎𝒙𝟑 = 𝟒

Ejemplo 2): Utilizar el método de Gauss-Seidel hasta que |∈𝒂 | < 𝟏% para aproximar la
solución del siguiente sistema de ecuaciones:

𝟏𝟎𝒙𝟏 − 𝒙𝟐 =𝟗
−𝒙𝟏 +𝟏𝟎𝒙𝟐 − 𝟐𝒙𝟑 = 𝟕
− 𝟐𝒙𝟐 + 𝟏𝟎𝒙𝟑 = 𝟔

SOLUCIÓN
𝟗 + 𝒙𝟐 + 𝟎
𝒙𝟏 =
𝟏𝟎
𝟕 + 𝒙𝟏 + 𝟐𝒙𝟑
𝒙𝟐 =
𝟏𝟎
𝟔 + 𝟐𝒙𝟐 + 𝟎
𝒙𝟑 =
𝟏𝟎
Iteración Nº 1:

Sea: 𝒙𝟐 = 𝟎; 𝒙𝟑 = 𝟎

(𝟏) 𝟗+𝟎 𝟗
𝑿𝟏 = = = 𝟎. 𝟗
𝟏𝟎 𝟏𝟎

(𝟏) 𝟕 + 𝟎. 𝟗 + 𝟐(𝟎)
𝑿𝟐 = = 𝟎. 𝟕𝟗
𝟏𝟎

MG. ARTURO H. GAMARRA MORENO

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA MECÁNICA
(𝟏) 𝟔 + 𝟐(𝟎. 𝟕𝟗) + 𝟎
𝑿𝟑 = = 𝟎. 𝟕𝟓𝟖
𝟏𝟎

Iteración Nº 2:

(𝟐) 𝟗 + 𝟎. 𝟕𝟗
𝑿𝟏 = = 𝟎. 𝟗𝟕𝟗
𝟏𝟎

(𝟐) 𝟕 + 𝟎. 𝟗𝟕𝟗 + 𝟐(𝟎. 𝟕𝟓𝟖)

𝑿𝟐 = = 𝟎. 𝟗𝟒𝟗𝟓
𝟏𝟎

(𝟐) 𝟔 + 𝟐(𝟎. 𝟗𝟒𝟗𝟓) + 𝟎

𝑿𝟑 = = 𝟎. 𝟕𝟖𝟗𝟗
𝟏𝟎
Verificando error:

Para 𝒙𝟏 :
𝟎. 𝟗𝟕𝟗 − 𝟎. 𝟗
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟖. 𝟎𝟔𝟗 %
𝟎. 𝟗𝟕𝟗

Para 𝒙𝟐 :
𝟎. 𝟗𝟒𝟗𝟓 − 𝟎. 𝟕𝟗
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟏𝟔. 𝟕𝟗𝟖 %
𝟎. 𝟗𝟒𝟗𝟓

Para 𝒙𝟑 :
𝟎. 𝟕𝟖𝟗𝟗 − 𝟎. 𝟕𝟓𝟖
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟒. 𝟎𝟑𝟖 %
𝟎. 𝟕𝟖𝟗𝟗

Iteración Nº 3:

(𝟑) 𝟗 + 𝟎. 𝟗𝟒𝟗𝟓
𝑿𝟏 = = 𝟎. 𝟗𝟗𝟒𝟗𝟓
𝟏𝟎

(𝟑) 𝟕 + 𝟎. 𝟗𝟗𝟒𝟗𝟓 + 𝟐(𝟎. 𝟕𝟖𝟗𝟗)

𝑿𝟐 = = 𝟎. 𝟗𝟓𝟕𝟒𝟕𝟓
𝟏𝟎

(𝟑) 𝟔 + 𝟐(𝟎. 𝟗𝟓𝟕𝟒𝟕𝟓) + 𝟎

𝑿𝟑 = = 𝟎. 𝟕𝟗𝟏𝟒𝟗𝟓
𝟏𝟎

Verificando error:

Para 𝒙𝟏 :
𝟎. 𝟗𝟗𝟒𝟗𝟓 − 𝟎. 𝟗𝟕𝟗
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟏. 𝟔𝟎𝟑 %
𝟎. 𝟗𝟗𝟒𝟗𝟓

Para 𝒙𝟐 :
𝟎. 𝟗𝟓𝟕𝟒𝟕𝟓 − 𝟎. 𝟗𝟒𝟗𝟓
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟎. 𝟖𝟑𝟐 %
𝟎. 𝟗𝟓𝟕𝟒𝟕𝟓

Para 𝒙𝟑 :

MG. ARTURO H. GAMARRA MORENO

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA MECÁNICA
𝟎. 𝟕𝟗𝟏𝟒𝟗𝟓 − 𝟎. 𝟕𝟖𝟗𝟗
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟎. 𝟐𝟎𝟏 %
𝟎. 𝟕𝟗𝟏𝟒𝟗𝟓

Iteración Nº 4:

(𝟒) 𝟗 + 𝟎. 𝟗𝟓𝟕𝟒𝟕𝟓
𝑿𝟏 = = 𝟎. 𝟗𝟗𝟓𝟕𝟒𝟕𝟓
𝟏𝟎

(𝟒) 𝟕 + 𝟎. 𝟗𝟗𝟓𝟕𝟒𝟕𝟓 + 𝟐(𝟎. 𝟕𝟗𝟏𝟒𝟗𝟓)

𝑿𝟐 = = 𝟎. 𝟗𝟓𝟕𝟖𝟕𝟑𝟕𝟓
𝟏𝟎

(𝟒) 𝟔 + 𝟐(𝟎. 𝟗𝟓𝟕𝟖𝟕𝟑𝟕𝟓) + 𝟎

𝑿𝟑 = = 𝟎. 𝟕𝟗𝟏𝟓𝟕𝟒𝟕𝟓
𝟏𝟎

Verificando error:

Para 𝒙𝟏 :
𝟎. 𝟗𝟗𝟓𝟕𝟒𝟕𝟓 − 𝟎. 𝟗𝟗𝟒𝟗𝟓
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟎. 𝟎𝟖𝟎 %
𝟎. 𝟗𝟗𝟓𝟕𝟒𝟕𝟓

Para 𝒙𝟐 :
𝟎. 𝟗𝟓𝟕𝟖𝟕𝟑𝟕𝟓 − 𝟎. 𝟗𝟓𝟕𝟒𝟕𝟓
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟎. 𝟎𝟒𝟏 %
𝟎. 𝟗𝟓𝟕𝟖𝟕𝟑𝟕𝟓

Para 𝒙𝟑 :
𝟎. 𝟕𝟗𝟏𝟓𝟕𝟒𝟕𝟓 − 𝟎. 𝟕𝟗𝟏𝟒𝟗𝟓
|∈𝒂 | = | | ∗ 𝟏𝟎𝟎 = 𝟎. 𝟎𝟏𝟎 %
𝟎. 𝟕𝟗𝟏𝟓𝟕𝟒𝟕𝟓

Solución:

𝒙𝟏 = 𝟎. 𝟗𝟗𝟓𝟕𝟒𝟕𝟓

𝒙𝟐 = 𝟎. 𝟗𝟓𝟕𝟖𝟕𝟑𝟕𝟓

𝒙𝟑 = 𝟎. 𝟕𝟗𝟏𝟓𝟕𝟒𝟕𝟓

CASO 1:
Suponga que un grupo de tres objetos está unido por una cuerda ligera mientras va en
caída libre a una velocidad de 5 m/s (figura 1). Calcule la tensión en cada sección de la
cuerda y la aceleración del grupo, dados los siguientes datos:

Objeto Masa Coeficiente de arrastre

Kg Kg/s
1 70 10
2 60 14
3 40 17

MG. ARTURO H. GAMARRA MORENO

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA MECÁNICA

MG. ARTURO H. GAMARRA MORENO