CURSO
MÉTODOS NUMÉRICOS PARA
INGENIERÍA QUÍMICA
TEMA:
SISTEMAS DE ECUACIONES LINEALES
[Link]. Ing. Alex Willy Pilco Nuñez
Introducción
Los equipos fundamentales en una refinería de petróleo son las columnas de destilación, en donde se lleva a cabo la
separación de los diversos componentes de petróleo. La operación de dichas torres se sustenta en un balance de materia
y energía dando lugar éstos a sistemas de ecuaciones lineales, no lineales y de ecuaciones diferenciales.
Refinería de petróleo de Talara
Nota. [Link]
minem-nueva-refineria-de-talara-procesara-gasolina-y-diesel-de-
alta-calidad-para-peruanos
Sistemas de ecuaciones lineales
Muchos problemas prácticos de ingeniería se reducen a la solución de un sistema de ecuaciones lineales.
Un sistema de 𝑚 ecuaciones lineales en 𝑛 incógnitas tiene la forma general
𝑎1,1 𝑥1 + 𝑎1,2 𝑥2 + … + 𝑎1,𝑛 𝑥𝑛 = 𝑏1
𝑎2,1 𝑥1 + 𝑎2,2 𝑥2 + … + 𝑎2,𝑛 𝑥𝑛 = 𝑏2
(1)
. . . .
. . . .
. . . .
𝑎𝑚,1 𝑥1 + 𝑎𝑚,2 𝑥2 + … + 𝑎𝑚,𝑛 𝑥𝑛 = 𝑏𝑚
Con la notación matricial, se puede escribir la ecuación anterior como
𝑎1,1 𝑎1,2 . . . 𝑎1,𝑛 𝑥1 𝑏1
𝑎2,1 𝑎2,2 . . . 𝑎2,𝑛 𝑥2 𝑏2
. . . . . . . . (2)
. . . . . . . = .
. . . . . . . .
𝑎𝑚,1 𝑎𝑚,2 . . . 𝑎𝑚,𝑛 𝑥𝑛 𝑏𝑛
y concretamente como 𝐴 x = b.
Donde 𝐴 es la matriz coeficiente del sistema, x el vector incógnita y b el vector de términos independientes.
Sistemas de ecuaciones lineales
Sistema de ecuaciones lineales
𝐴x=b
Rango 𝐴 ≠ Rango 𝐵 Rango 𝐴 = Rango 𝐵
Inconsistente Consistente
Sin solución Rango 𝐴 = 𝑛 Rango 𝐴 < 𝑛
Solución única Número infinito
de soluciones
Sistemas de ecuaciones lineales
Métodos para la solución de
un sistema de ecuaciones
lineales
Métodos directos Métodos iterativos
• Eliminación de Gauss • Método de Jacobi
• Eliminación de Gauss con • Método de Gauss-Seidel
pivoteo
• Método de Thomas
• Factorización directa
• Factorización con pivoteo
• Método de Doolitle
• Factorización de matrices
simétricas
• Método de Cholesky
Método de Jacobi
Ejemplo 𝑥1 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘) Τ4
𝑥2 (𝑘+1) = 1Τ4 + 𝑥1 (𝑘) Τ4 + 𝑥3 (𝑘) Τ4
Resuelva el siguiente sistema por el método de Jacobi
𝑥3 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘) Τ4 + 𝑥4 (𝑘) Τ4
4𝑥1 − 𝑥2 =1 𝑥4 (𝑘+1) = 1Τ4 + 𝑥3 (𝑘) Τ4
−𝑥1 + 4𝑥2 − 𝑥3 =1
(1) donde 𝑘 = 0, 1, 2, . . .
−𝑥2 + 4𝑥3 − 𝑥4 = 1
−𝑥3 + 4𝑥4 = 1
El cálculo de x(1) en el método de Jacobi se obtiene
Use como vector inicial x(0) = [0 0 0 0]T reemplazando x(0) en cada una de las Ecuaciones 2.
Solución 𝑥1 (1) = 1Τ4 + 𝑥2 (0) Τ4
𝑥2 (1) = 1Τ4 + 𝑥1 (0) Τ4 + 𝑥3 (0) Τ4
Despejando 𝑥1 de la primera ecuación, 𝑥2 de la segunda
𝑥3 (1) = 1Τ4 + 𝑥2 (0) Τ4 + 𝑥4 (0) Τ4
ecuación, etc., se obtiene
𝑥4 (1) = 1Τ4 + 𝑥3 (0) Τ4
𝑥1 = 1Τ4 + 𝑥2 Τ4
𝑥2 = 1Τ4 + 𝑥1 Τ4 + 𝑥3 Τ4 Usando el vector inicial x(0) = [0 0 0 0]T
(2)
𝑥3 = 1Τ4 + 𝑥2 Τ4 + 𝑥4 Τ4
𝑥4 = 1Τ4 + 𝑥3 Τ4 𝑥1 (1) = 1Τ4 + 0Τ4 = 1Τ4
𝑥2 (1) = 1Τ4 + 0Τ4 + 0Τ4 = 1Τ4
De acuerdo al método de Jacobi, el sistema de 𝑥3 (1) = 1Τ4 + 0Τ4 + 0Τ4 = 1Τ4
ecuaciones lineales 2 en forma iterativa 𝑥4 (1) = 1Τ4 + 0Τ4 = 1Τ4
Método de Jacobi
𝑥1 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘) Τ4
𝑥2 (𝑘+1) = 1Τ4 + 𝑥1 (𝑘) Τ4 + 𝑥3 (𝑘) Τ4
𝑥3 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘) Τ4 + 𝑥4 (𝑘) Τ4
𝑥4 (𝑘+1) = 1Τ4 + 𝑥3 (𝑘) Τ4
donde 𝑘 = 0, 1, 2, . . . , 10
Método de Gauss-Seidel
Ejemplo 𝑥1 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘) Τ4
𝑥2 (𝑘+1) = 1Τ4 + 𝑥1 (𝑘+1) Τ4 + 𝑥3 (𝑘) Τ4
Resuelva el siguiente sistema por el método Gauss-
𝑥3 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘+1) Τ4 + 𝑥4 (𝑘) Τ4
Seidel
𝑥4 (𝑘+1) = 1Τ4 + 𝑥3 (𝑘+1) Τ4
4𝑥1 − 𝑥2 =1
−𝑥1 + 4𝑥2 − 𝑥3 =1 donde 𝑘 = 0, 1, 2, . . .
(1)
−𝑥2 + 4𝑥3 − 𝑥4 = 1
−𝑥3 + 4𝑥4 = 1 El cálculo de x(1) en el método de Jacobi se obtiene
reemplazando x(0) en cada una de las Ecuaciones 2.
Use como vector inicial x(0) = [0 0 0 0]T
𝑥1 (1) = 1Τ4 + 𝑥2 (0) Τ4
Solución
𝑥2 (1) = 1Τ4 + 𝑥1 (1) Τ4 + 𝑥3 (0) Τ4
𝑥3 (1) = 1Τ4 + 𝑥2 (1) Τ4 + 𝑥4 (0) Τ4
Despejando 𝑥1 de la primera ecuación, 𝑥2 de la segunda
ecuación, etc., se obtiene 𝑥4 (1) = 1Τ4 + 𝑥3 (1) Τ4
𝑥1 = 1Τ4 + 𝑥2 Τ4 Usando el vector inicial x(0) = [0 0 0 0]T
𝑥2 = 1Τ4 + 𝑥1 Τ4 + 𝑥3 Τ4 𝑥1 (1) = 1Τ4 + 0Τ4 = 1Τ4
(2)
𝑥3 = 1Τ4 + 𝑥2 Τ4 + 𝑥4 Τ4 1
𝑥4 = 1Τ4 + 𝑥3 Τ4 𝑥2 (1) = 1Τ4 + (4)Τ4 + 0Τ4 = 5Τ16
5
𝑥3 (1) = 1Τ4 + (16)Τ4 + 0Τ4 = 21Τ64
De acuerdo al método de Gauss-Seidel, el sistema de 21
ecuaciones lineales 2 en forma iterativa 𝑥4 (1) = 1Τ4 + (64)Τ4 = 85Τ256
Método de Gauss-Seidel
𝑥1 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘) Τ4
𝑥2 (𝑘+1) = 1Τ4 + 𝑥1 (𝑘+1) Τ4 + 𝑥3 (𝑘) Τ4
𝑥3 (𝑘+1) = 1Τ4 + 𝑥2 (𝑘+1) Τ4 + 𝑥4 (𝑘) Τ4
𝑥4 (𝑘+1) = 1Τ4 + 𝑥3 (𝑘+1) Τ4
donde 𝑘 = 0, 1, 2, . . . , 10
Comparación de los métodos iterativos
Método de Jacobi Método de Gauss-Seidel
Criterio de convergencia de un proceso iterativo
¿Cuándo detener el proceso iterativo?
Sean los vectores
𝑇
𝒙(𝑘) = 𝑥1 (𝑘) 𝑥2 (𝑘) . . . 𝑥𝑛 (𝑘)
𝑇
𝒙(𝑘+1) = 𝑥1 (𝑘+1) 𝑥2 (𝑘+1) . . . 𝑥𝑛 (𝑘+1)
y 𝜀 el criterio de exactitud, cuyo valor conduce a un número de iteraciones. Si el valor de 𝜀 tiende a cero, esto puede
significar un número mayor de iteraciones. Para casos prácticos su valor puede ser 0.1, 0.01 o 0.001.
Detener el proceso iterativo cuando el valor absoluto 𝒙(𝑘+1) − 𝒙(𝑘) < 𝜀, donde
𝒙(𝑘+1) − 𝒙(𝑘) = 𝑥1 (𝑘+1) − 𝑥1 (𝑘) 2 + 𝑥2 (𝑘+1) − 𝑥2 (𝑘) 2 + . . . + 𝑥𝑛 (𝑘+1) − 𝑥𝑛 (𝑘) 2
¿Cuándo se asegura la convergencia del proceso iterativo?
𝑛 𝑛
𝑎𝑖,𝑖 > 𝑎𝑖,𝑗 1≤𝑖≤𝑛 o 𝑎𝑖,𝑖 > 𝑎𝑗,𝑖 1≤𝑖≤𝑛
𝑗=1 𝑗=1
𝑗≠𝑖 𝑗≠𝑖
