SISTEMA DE ECUACIONES LINEALES (SEL) Y SUS SOLUCIONES

Un Sistema de Ecuaciones Lineales (SEL), es un conjunto de ecuciones cuyas variables tienen como máximo exponente la unid
El SEL contiene un número de vriables y la misma cantidad de ecuaciones. #var = # Ecuaciones, si est condición no se cump

La forma denotación general de un SEL es:

■8(𝑎_11 𝑋1+&𝑎_12 𝑋2+ …& 𝑎_1𝑛
𝑋𝑛=𝑏1@𝑎_21 𝑋1+&𝑎_22 𝑋2+ …
&𝑎_2𝑛 𝑋𝑛=𝑏2@█(…_ @𝑎_𝑛1
𝑋1+)&█(…@𝑎_𝑛2 𝑋2+ …)&█(…@𝑎_𝑛𝑛
𝑋𝑛=𝑏𝑛))
SOLUCIONES
MÉTODO DE LA DETERMINANTE (Crammer)
■8(𝑎_11 𝑋1+&𝑎_12 𝑋2+ …& 𝑎_1𝑛
Un SEL está formado con los siguientes componentes: 𝑋𝑛=𝑏1@𝑎_21 𝑋1+&𝑎_22 𝑋2+ …
&𝑎_2𝑛 𝑋𝑛=𝑏2@█(…_ @𝑎_𝑛1
𝑋1+)&█(…@𝑎_𝑛2 𝑋2+ …)&█(…@
𝑋𝑛=𝑏𝑛))
aij =

X=

b=

Sea un SEL de 3 x 3: a11 X + a12 Y + a13 Z = b1

a21 X + a22 Y + a23 Z = b2
a31 X + a32 Y + a33 Z = b3

Obteniendo a11 a12 a13

[A]= a21 a22 a23
a31 a32 a33

X
[X]= Y
Z

b1
[b] = b2
b3

Para resolver un SEL por este método, debemos obtener las Determinantes, de la Matriz A, de una nueva matriz

a11 a12 a13

matriz [A]= a21 a22 a23 hallando su Determinante, por Sarrus
a31 a32 a33
 b1 a12 a13
[ Ax ] = b2 a22 a23
b3 a32 a33

a11 b1 a13
[ Ay ] = a21 b2 a23
a31 b3 a33

a11 a12 b1
[ Az ] = a21 a22 b2
a31 a32 b3
======================================================================================
Hallando los valores de las variables:
Ejemplo: Resolver el siguiente sistema.

X = | Ax | / | A |

Y = | Ay | / | A | 4 -9
[A]= 2 -4
Z = | Az | / | A | 1 -1

X1 6.95
[X]= X2 2.5
X3 -0.15

5
[b]= 3
4
 Y SUS SOLUCIONES
omo máximo exponente la unidad, sólo la unidad, en caso contrario conformarían ecuaciones NO LINEALES
nes, si est condición no se cumple, entonces tedríamos infinitas soluciones o ninguna solución.

1+&𝑎_12 𝑋2+ …& 𝑎_1𝑛

21 𝑋1+&𝑎_22 𝑋2+ …
𝑏2@█(…_ @𝑎_𝑛1
@𝑎_𝑛2 𝑋2+ …)&█(…@𝑎_𝑛𝑛

[A]

[X]

[b]

Matriz A, de una nueva matriz en Ax, Ay, Az

a11 a12 a13

terminante, por Sarrus |A|= a21 a22 a23 =a11*a22*a33+a21*a32*a13+a31*a12*a23-a31*a22*a13-
a31 a32 a33
 b1 a12 a13
| Ax | = b2 a22 a23 =b1*a22*a33+b2*a32*a13+b3*a12*a23-b3*a22*a13-b2*a
b3 a32 a33

a11 b1 a13
| Ay | = a21 b2 a23 =a11*b2*a33+a21*b3*a13+a31*b1*a23-a31*b2*a13-a21*
a31 b3 a33

a11 a12 b1
| Az | = a21 a22 b2 =a11*a2*b3+a21*a32*b1+a31*a12*b2-a31*a22*b1-a21*a
a31 a32 b3
==================== = 3

4X1 - 9X2 + 2X3 = 5

2X1 - 4X2 + 6X3 = 3
X1 - X2 + 3X3 = 4
2 4 -9 2
6 |A|= 2 -4 6 -20
3 1 -1 3

Solución 5 -9 2
Solución | Ax | = 3 -4 6 -139
Solución 4 -1 3

4 5 2
| Ay | = 2 3 6 -50
1 4 3

4 -9 5
| Az | = 2 -4 3 3
1 -1 4
3+a31*a12*a23-a31*a22*a13-a21*a12*a33-a11*a32*a23
+b3*a12*a23-b3*a22*a13-b2*a12*a33-b1*a32*a23

+a31*b1*a23-a31*b2*a13-a21*b1*a33-a11*b3*a23

a31*a12*b2-a31*a22*b1-a21*a12*b3-a11*a32*b2
 Método Matricial
Este método utiliz definiciones de álgebra matricial
De la expresión general de un SEL, se despeja el vector de las variables y se opera una multiplicación matricialmente.

Sea un SEL 3 x 3
a11 X + a12 Y + a13 Z = b1
a21 X + a22 Y + a23 Z = b2 extraemos sus componentes:
a31 X + a32 Y + a33 Z = b3

la expresión matricial será:

[A]*[X]=[b]
despejando el vector de las variables tendremos:

[ X ] = [ A ]-1 * [ b ]

Ejemplo: resolver el SEL siguiente.

10 x1 + x2 - 5 x3 = 1
- 20 x1 + 3 x2 + 20 x3 = 2 [A]=
5 x1 + 3 x2 + 5 x3 = 6

[X]=

[b]=

F2, ctrl + shift + ENTER

[ A ] -1
cación matricialmente.

[A]
[X]
[b]

10 1 -5
-20 3 20
5 3 5

x1
x2
x3

1
2
6

-0.36 -0.16 0.28 x1 = 1

1.6 0.6 -0.8 x2 = -2
-0.6 -0.2 0.4 x3 = 1.4
 MÉTODO ITERACTIVO DE GAUSS SEIDEL
Método iteractivo que resuelve un Sistema de Ecuaciones Lineales (SEL)
-Requiere tener definido un SEL, donde estén un número de ecuaciones y de la misma cantidad de variables, todo ellos igualad
- El método iteractivo es un cálculo de los siguientes valores de cada variable mediante un ajuste y reajuste de los valores que
- Este método requiere que se obtenga otro sistema de ecuaciones lineales, donde cada ecuación es producto del despeje de

Tendremos un SEL de 3 x 3:

a11 x1 + a12 x2 + a13 x3 = b1 x1 = ( b1 - a12 x - a13 x3 ) / a11

a21 x1 + a22 x2 + a23 x3 = b2 x2 = ( b2 - a21 x1 - a23 x3 )/ a22
a31 x1 + a32 x2 + a33 x2 = b3 x3 = ( b3 - a31 x1 - a32 x2 ) / a33

Los puntos iniciales ( semilla) [ 1, 1, 1, ….] ó [ 0, 0, 0, ….. ], salvo que el problema suministre estos valores, o in
[x1o, x2o, x3o, …. ], [xo, yo, zo, ….]
Modo de hallar los siguientes valores de las variables [x11, x21, x31, …. ], [x1, y1, z1, …

a11 x + a12 y + a13 z = b1 x1 = ( b1 - a12 yo - a13 zo) / a11

a21 x + a22 y + a23 z = b2 y1 = ( b2 - a21 x1 - a23 zo) / a22
a31 x + a32 y + a33 z = b3 z1 = ( b3 - a31 x1 - a32 y1 )/a33

ALGORITMO DEL MÉTODO ITERACTIVO DE GAUSS SEIDEL

0.- Inicio
1.- Determinar el SEL
2.- Definir los puntos iniciales [x1, y1, z1, … ]
3.- Despejar de cada ecuación la variable correspondiente. x1 = ( b1 - a12 yo - a13 zo
4- Formular el nuevo SEL en función a los valores iniciales y los anteriores calculados. y1 = ( b2 - a21 x1 - a23 zo)
5.- Calcular los nuevos valores [x1, y1, z1, … ] z1 = ( b3 - a31 x1 - a32 y1
6.- calcualr la distancia vectorial𝐷𝑉=√(( 〖𝑥 1−𝑥𝑜) 〗 ^2+(𝑦1−𝑦𝑜)^2+(𝑧1−𝑧𝑜)^2+…) )
7.- Es DV < EPS ?
SI X1, Y1, Z1, … solución , fin
No Xo = X1, Yo = Y1, Zo = Z1, … continua en 5.

Ejemplo Resolver el siguiente SEL

10 x + y - 5z = 1 inicio del SEL sea [1 ,1,1, ]

-20 x+3 Y+20 z=2
5x + 3 y +5 z = 6
TABLA DE CONTROL DEL MÉTODO ITERACTIVO DE GAUSS SEIDEL
x1 = ( 1 - yo + 5 zo) /10
y1 = (2 + 20 x1 - 20 zo) / 3 iter X Y Z
z1 = (6 - 5 x1 -3y1) / 5 0 1 1 1
1 0.5 -2.66666667 2.3
2 1.51666667 -4.55555556 2.41666667
3 1.76388889 -3.68518519 1.64722222
 4 1.29212963 -1.70061728 0.92824074
5 0.7341821 -0.62705761 0.84205247
6 0.583732 -1.05546982 1.2495499
7 0.83032193 -2.12818644 1.64658993
8 1.13611361 -2.73650882 1.70579168
9 1.22654672 -2.52829972 1.49043311
10 1.09804653 -1.94924389 1.27149981
11 0.93067429 -1.60550343 1.23262776
12 0.87686422 -1.70509026 1.34618993
13 0.94360399 -2.0172396 1.46673977
14 1.03509384 -2.21097283 1.49148985
15 1.06684221 -2.16431762 1.43174836
16 1.03230594 -1.9962828 1.36546374
17 0.98236015 -1.88735725 1.3500542
18 0.96376283 -1.90860918 1.38140268
19 0.98156226 -1.99893615 1.41779943
20 1.00879333 -2.06004067 1.42723107
21 1.0196196 -2.05074313 1.41082627
22 1.01048745 -2.00225883 1.39086785
23 0.99565981 -1.9680536 1.38517236
24 0.98939154 -1.97187212 1.39373173
25 0.99405308 -1.9978577 1.40466154
26 1.00211654 -2.01696667 1.40806346
27 1.0057284 -2.01556709 1.40361186
28 1.00336264 -2.00166146 1.39763424
29 0.99898327 -1.99100649 1.39562063
30 0.99691096 -1.99139777 1.3979277
31 0.99810363 -1.99882714 1.40119266
32 1.00047905 -2.00475744 1.40237542
33 1.00166345 -2.00474644 1.40118441
34 1.00106685 -2.00078374 1.39940339
35 0.99978007 -1.99748882 1.39871322
36 0.99910549 -1.99738486 1.39932542
37 0.9994012 -1.99949484 1.4002957
38 1.00009734 -2.00132246 1.40069614
39 1.00048031 -2.00143882 1.40038298
40 1.00033537 -2.00031738 1.39985506
41 0.99995927 -1.99930527 1.39962389
42 0.99974247 -1.99920947 1.39978321
43 0.99981255 -1.99980438 1.40007008
44 1.00001548 -2.000364 1.40020293
45 1.00013786 -2.00043375 1.40012239
46 1.00010457 -2.00011879 1.39996671
47 0.99999523 -1.99980983 1.39989066
48 0.99992631 -1.99976233 1.39993108
49 0.99994177 -1.99992872 1.40001546
d de variables, todo ellos igualados a una constante.
te y reajuste de los valores que se calculan, teniendo como tope al EPS.
ión es producto del despeje de cada variable en su correspondiente ecuación.

2 x - a13 x3 ) / a11
1 x1 - a23 x3 )/ a22
1 x1 - a32 x2 ) / a33

ma suministre estos valores, o indique cómo obtenerlos.

x11, x21, x31, …. ], [x1, y1, z1, …. ]

1 = ( b1 - a12 yo - a13 zo) / a11

1 = ( b2 - a21 x1 - a23 zo) / a22
1 = ( b3 - a31 x1 - a32 y1 )/a33

CTIVO DE GAUSS SEIDEL

DV EPS

3.92230091 0.0001
2.14828384
1.18773236
 2.16286868
1.21295681
0.6101036
1.17011467
0.68342484
0.31290332
0.63225766
0.38429417
0.1603416
0.34120913
0.21567457
0.08218097
0.18390788
0.12081726
0.04219279
0.0989975
0.06755923
0.02174327
0.05322088
0.03771333
0.01127542
0.02857343
0.02101774
0.00590094
0.01531976
0.01169453
0.00312572
0.00820232
0.00649692
0.00167972
0.00438531
0.00360397
0.00091687
0.00234113
0.00199629
0.00050822
0.00124793
0.0011042
0.00028558
0.00066417
0.00060992
0.00016226
0.00035291
0.00033645
9.2948E-05
solución
 MÉTODO DE REDUCCIÓN DE GAUSS
- Este método es de cálculo director de los valores de las variables que intervienen en el SEL.
- Utiliza los componentes de un SEL, como son: la diagonal principal, los triángulos de elementos superior e inferior de un SEL
- La solución se encuentra de la siguiente forma:

La notación general de un SEL: ■8(𝑎_11 𝑋1+&𝑎_12 𝑋2+ …& 𝑎_1𝑛

𝑋𝑛=𝑏1@𝑎_21 𝑋1+&𝑎_22 𝑋2+ …
&𝑎_2𝑛 𝑋𝑛=𝑏2@█(…_ @𝑎_𝑛1
𝑋1+)&█(…@𝑎_𝑛2 𝑋2+ …)&█(…@𝑎_𝑛𝑛
𝑋𝑛=𝑏𝑛))

Sea el SEL 3x3 a11 x + a12 y + a13 z = b1

a21 x + a22 y + a23 z = b2 [A] * [X] = [b]
a31 x + a32 y + a33 z = b3

- De acuedo al método de reducción de Gauss, ubicar el triángulo inferior de elementos

La matriz A será: a11 a12 a13 b1

a21 a22 a23 Vector b es: b2
a31 a32 a33 b3

-Señalizar los elementos del triángulo inferior, para convertirlos en CERO

a11 a12 a13

a21 a22 a23
a31 a32 a33

Formular una matriz extendida conformada por la unión de la [ A ] y [ b ], que finalmente tendremos un resultado de la forma

a11 a12 a13 b1

Matriz extendida = 0 a22' a23' b2'
0 0 a33'' b3''

a11 x + a12 y + a13 z = b1

Expresar el nuevo SEL= a22'y + a23'z = b2'
a33" z = b3"

Despejando una variable abajo hacia arriba, tendremos

z = b3"/a33"
y = (b2' - a23' * (b3"/ a33"))/a22'
x = (b1 - a12 * (b2'- a23'*(b3"/ a33")/a22') - a13 *(b33" / a33"))
CONVIRTIENDO EN CERO LOS VALORES DE LOS ELEMENTOS DEL TRIÁNGULO INFERIOR
Para convertir en CERO a21, debo
L1 a11 a12 a13 b1
L2 a21 a22 a23 b2 = (-a21/a11) * L1 + L2
L3 a31 a32 a33 b3 =(-a31/a11) * L1 + L3

L1 a11 a12 a13 b1

L2' 0 a22' a23' b2' Para convertir en CERO a32', deb
L3' 0 a32' a33' b3' = (-a32'/a22') * L2' + L3'

L1 a11 a12 a13 b1

L2' 0 a22' a23' b2'
L3' 0 0 a33" b3"
GAUSS
os superior e inferior de un SEL

x
vector X = y
z

remos un resultado de la forma siguiente.

ara convertir en CERO a21, debo calcualar el FACTOR = - a21/a11 Para convertir en CERO a31, debo calcualar el FACTOR = - a

(-a21/a11) * L1 + L2 L1 = a11 a12 a13 b1

(-a31/a11) * L1 + L3 Factor = * (-a21/a11) * (-a21/a11) * (-a21/a11) * (-a21/a11)
L2 = '+ a21 + a22 + a23 + b2
L2' = 0 a22' a23' b2'

ara convertir en CERO a32', debo calcular el FACTOR = - a32'/a22'

(-a32'/a22') * L2' + L3'
L2' = a22' a23' b2'
Factor = * (-a32'/a22')* (-a32'/a22')* (-a32'/a22')
L3' = + a32' + a33' + b3'
L3" = 0 a33" b3"
debo calcualar el FACTOR = - a31/a11

L1 = a11 a12 a13 b1

Factor = * (-a31/a11) * (-a31/a11) * (-a31/a11) * (-a31/a11)
L3 = + a31 + a32 + a33 + b3
L3' = 0 a32' a33' b3'
 APLICACIÓN
PROBLEMA 1 Con los datos del diagrama siguiente (donde los porcentajes están dados en peso), en
de las corrientes M1, M2 y M3. M4 = 100 kilogramos.

M1 + M2 + M3 = M4 = 100

Etanol : 0,83*M1 + 0,55 * M3 = 58 0.83 0 0.55

Metanol: 0,61 * M2 + 0,24 * M3 = 21 0 0.61 0.24
Agua : 0,17 * M1 + 0,39 * M2 + 0,21 * M3 = 21 0.17 0.39 0.21

Convirtiendo en CERO los elementos del triángulo inferior

L1 0.83 0 0.55 58
L2 0 0.61 0.24 21
L3 0.17 0.39 0.21 21 =(-0,17/0,83) * L1 + L3

L1 0.83 0 0.55 58
L2 0 0.61 0.24 21
L3 0 0.39 0.0973493 9.120482 =(-0,39/0,61) * L2 + L3

L1 0.83 0 0.55 58
L2 0 0.61 0.24 21
L3 0 0 -0.056093 -4.305748

Formulando el nuevo SEL

0,83 M1 + 0 M2 + 0,55 M3 = 58
0 M1 + 0,61 M2 + 0,24 M3 = 21
0 M1 + 0 M2 -0,056093 M3 = -4,305748
Hallando valores de la variables.
M3 = 76.7608793
M2 = 4.22522783
M1 = 19.0138752
M4 = 99.9999823
=======================================================================================================
PROBLEMA 4 (Niev 3.1 pp 254) En una columna de cinco platos, se requiere absorber benceno que
V, con un aceite L que circula a contracorriente del gas. Considérese que el benceno transferido no
sustancialmente el número de moles de V y L fluyendo a contracorriente, que la relación de equilib
dada por la ley de Henry (y = mx) y que la columna opera a régimen permanente. Calcule la compo
del benceno en cada plato.

Datos: V = 100 moles/min;

L = 500 moles/min.
y0 = 0.09 fracción molar de benceno en V.
x0 = 0.0 fracción molar de benceno en L (el aceite entra por el domo sin benceno).
m = 0.12.

Solucionando, realizaremos balance de materia en cada plato de la columna.

5to plato L xo + V y4 - L x5 - V y5 = 0
4to. Plato L x5 + V y3 - L x4 - V y4 = 0
3er Plato L x4 + V y2 - Lx3 - V y3 = 0
2do Plato L x3 + V y1 - L x2 - V y2 = 0
1er Plato L x2 + V yo - L x1 - V y1 = 0

reemplazando valores de acuerdo a la Ley de Henry

L xo + V*0,12*x4 - L x5 - V*0,12*x5 = 0 0x1 + 0x2 + 0x3 + 12x4 - 512x5 = 0

L x5+ V*0,12*x3 - L x4 - V*0,12*x4 = 0 0x1 + 0x2 + 12x3 - 512x4 + 500x5 = 0
L x4 + V*0,12*x2 -L x3 - V*0,12*x3 = 0 0x1 + 12x2 - 512x3 + 500x4 + 0x5 = 0
L x3 + V*0,12*x1 - L X2 - V*0,12*x2 = 0 12x1 - 512x2 + 500x3 + 0x4 + 0x5 = 0
L x2 + V*0.09 - L x1 - V*0,12*x1 = 0 -512x1 + 500x2 + 0x3 + 0x4 + 0x5 = -9

Reunificar y reordenar el SEL. CALCULAR……

=======================================================================================================
Propuesto
PROBLEMA 9.- (Tema 9 LIN.1.) Balances de materia en estado estacionario para una secuencia de
separación.
El p-xileno, estireno, tolueno y benceno se separan en la siguiente secuencia
de columnas de destilación.

a) Calcula el flujo molar de las corrientes D1, D2, B1 y B2.

b) Reduce el flujo original de la corriente de alimento en un 1% y
luego en un 2% y repite el apartado (a). Explica los resultados.
c) Determina el flujo molar y las composiciones de las corrientes B y
D del apartado (a).
-1 -9.2591 -2.9014 0.1625
-1 -9.6304 -2.785 0.0408
-1 -9.9547 -2.4932 -0.07696
-1 -10.2598 -2.1633 -0.01655
están dados en peso), encuentre posibles valores

58
21
21

(-0,17/0,83) * L1 + L3

(-0,39/0,61) * L2 + L3
====================================
ere absorber benceno que está contenido en una corriente de gas
l benceno transferido no altera
que la relación de equilibrio está
anente. Calcule la composición

0 0 0
0 0 12
[A] = 0 12 -512
12 -512 500
-512 500 0

x1
x2
[x] = x3
x4
x5

0
0
[b] = 0
0
-9

-0.001952 -0.001998848 -0.0019999724

-0.001998848 -0.00204682 -0.0020479717
[A] -1 = -0.001999972 -0.002047972 -0.0020491237
-0.001999999 -0.002047999 -4.915132E-05
-0.002 -4.8E-05 -1.151984E-06
========================================================

0.07D1 + 0.18B1 + 0.15D2 + 0.24B2 = 10.5

0.04D1 + 0.24B1 + 0.10D2 + 0.65B2 = 17.5
0.54D1 + 0.42B1 + 0.54D2 + 0.10B2 = 28
0.35D1 + 0.16B1 + 0.21D2 + 0.01B2 = 14

A)
0.07 0.18 0.15 0.24
0.04 0.24 0.1 0.65
[A] = 0.54 0.42 0.54 0.1
0.35 0.16 0.21 0.01

10.5
17.5
[b] = 28
14

-2.81048387 1.221774194 -1.802419355 6.06048387097

94.3629032 -29.83064516 -36.68548387 41.1370967742
[A] -1 = -66.0443548 20.24596774 30.52822581 -36.205645161
-24.5080645 9.362903226 8.959677419 -9.9919354839

B) 0.07 0.18 0.15 0.24

0.04 0.24 0.1 0.65
0.54 0.42 0.54 0.1
0.35 0.16 0.21 0.01

10.395
17.325
27.72
13.86

-2.81048387 1.221774194 -1.802419355 6.06048387097

94.3629032 -29.83064516 -36.68548387 41.1370967742
[A] -1 = -66.0443548 20.24596774 30.52822581 -36.205645161
 -24.5080645 9.362903226 8.959677419 -9.9919354839

0.07 0.18 0.15 0.24

0.04 0.24 0.1 0.65
0.54 0.42 0.54 0.1
0.35 0.16 0.21 0.01

10.1871
16.9785
27.1656
13.5828

0.07 0.18 0.15 0.24

[A] -1 = 0.04 0.24 0.1 0.65
0.54 0.42 0.54 0.1
0.35 0.16 0.21 0.01

0.07D1 + 0.18B1 = 43.75

0.04D1 + 0.24B1 = 17.5
0.54D1 + 0.42B1 + 0.54D2 + 0.10B2 = 28
0.35D1 + 0.16B1 + 0.21D2 + 0.01B2 = 14
 12 -512
-512 500
500 0
0 0
0 0

-0.0019999993 -0.002 x1= 0.0179999999

-0.0020479993 -4.799998E-05 x2= 0.0004319999
-4.915132E-05 -1.151984E-06 x3= 1.036786E-05
-1.178969E-06 -2.763207E-08 x4= 2.486887E-07
-2.763207E-08 -6.476268E-10 x5= 5.828641E-09
 26.25
17.5
8.75
17.5

25.9875
17.325
8.6625
 17.325

11.103939
16.027704
28.659708
12.122649