Eliminación Gausiana
Eliminación Gausiana
Eliminación Gausiana
SOLUCIÓN:
a) Verificando que el sistema es compatible determinado:
m=4: Número de ecuaciones
n=4: Número de incógnitas
Con el siguiente script en Matlab, hallamos los rangos de la matriz A y su ampliada, así
como la solución del sistema:
clc; A=[6 -2 2 4; 12 -8 6 10; 3 -13 9 3; -6 4 1 -18];
b=[12 34 27 -38]';
fprintf('El rango de A es: %d\n',rank(A));
fprintf('El rango de la ampliada de A es: %d\n',rank([A b]));
fprintf('La solución del sistema de ecuaciones es:\n');
disp(inv(A)*b)
𝑃𝑖𝑣𝑜𝑡𝑒 𝑒𝑠 6 ⟶ 6 −2 2 4 12
𝑓2 = 𝑓2 − 𝑓1 . 12/6 0 −4 2 2 10
[ | ]
𝑓3 = 𝑓3 − 𝑓1 . 3/6 0 −12 8 1 21
𝑓4 = 𝑓4 − 𝑓1 . −6/6 0 2 3 −14 −26
6 −2 2 4 12
𝑃𝑖𝑣𝑜𝑡𝑒 𝑒𝑠 − 4 ⟶
𝑓3 = 𝑓3 − 𝑓2 . −12/−4 [0 −4 2 2 10
| ]
𝑓4 = 𝑓4 − 𝑓2 . 2/−4 0 0 2 −5 −9
0 0 4 −13 −21
6 −2 2 4 12
𝑃𝑖𝑣𝑜𝑡𝑒 𝑒𝑠 2 ⟶ 0 −4 2 2 10
[ | ]
𝑓4 = 𝑓4 − 𝑓3 . 4/2 0 0 2 −5 −9
0 0 0 −3 −3
Aplicamos sustitución regresiva:
La matriz ampliada equivalente representa al siguiente sistema de ecuaciones equivalente:
−3𝑥4 = −3 ⟶ 𝑥4 = 1
2𝑥3 − 5𝑥4 = −9 ⟶ 𝑥3 = −2
−4𝑥2 + 2𝑥3 + 2𝑥4 = 10 ⟶ 𝑥2 = −3
6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 12 ⟶ 𝑥1 = 1
𝑥1 1
𝑥2 −3
[𝑥 ] = [ ]
3 −2
𝑥4 1
6 −2 2 4 12
12 −8 6 10 34
𝐴𝑏 = [ ]
3 −13 9 3 27
−6 4 1 −18 −38
12 −8 6 10 34
𝑓1 𝑥 𝑓2
6 −2 2 4 12
𝑃𝑖𝑣𝑜𝑡𝑒 𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜 [ | ]
3 −13 9 3 27
𝑑𝑒 6 𝑎 12
−6 4 1 −18 −38
⟶ 12 −8 6 10 34
𝑓2 = 𝑓2 − 𝑓1 . 6/12 0 2 −1 −1 −5
𝑓3 = 𝑓3 − 𝑓1 . 3/12 [ | ]
0 −11 7.5 0.5 18.5
𝑓4 = 𝑓4 − 𝑓1 . −6/12 0 0 4 −13 −21
12 −8 6 10 34
𝑓2 𝑥 𝑓3
0 −11 7.5 0.5 18.5
𝑃𝑖𝑣𝑜𝑡𝑒 𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜 [ | ]
0 2 −1 −1 −5
𝑑𝑒 2 𝑎 − 11
0 0 4 −13 −21
⟶ 12 −8 6 10 34
𝑓3 = 𝑓3 − 𝑓1 . 2/−11 [ 0 −11 7.5 0.5
|
18.5
]
𝑓4 = 𝑓4 − 𝑓1 . 0/−11 0 0 0.3636 −0.9091 −1.6364
0 0 4 −13 −21
12 −8 6 10 34
𝑓3 𝑥 𝑓4
0 −11 7.5 0.5 18.5
𝑃𝑖𝑣𝑜𝑡𝑒 𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜 [ | ]
0 0 4 −13 −21
𝑑𝑒 0.3636 𝑎 4
0 0 0.3636 −0.9091 −1.6364
12 −8 6 10 34
⟶ 0 −11 7.5 0.5 18.5
[
𝑓4 = 𝑓4 − 𝑓1 . 0.3636/6.7273 0 | ]
0 4 −13 −21
0 0 0 0.2727 0.2727
6 −2 2 4 12
12 −8 6 10 34
𝐴𝑏 = [ | ]
3 −13 9 3 27
−6 4 1 −18 −38
𝑓1 𝑥 𝑓4 −18 4 1 −6 −38
𝑐1 𝑥𝑐4 10 −8 6 12 34
[ | ]
𝑃𝑖𝑣𝑜𝑡𝑒 𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜 3 −13 9 3 27
𝑑𝑒 12 𝑎 − 18 4 −2 2 6 12
⟶ −18 4 1 −6 −38
𝑓2 = 𝑓2 − 𝑓1 . 10/−18 0 −5.7778 6.5556 8.6667 12.8889
𝑓3 = 𝑓3 − 𝑓1 . 3/−18 [ | ]
0 −12.3333 9.1667 2 20.6667
𝑓4 = 𝑓4 − 𝑓1 . 4/−18 0 −1.1111 2.2222 4.6667 3.5556
−18 4 1 −6 −38
𝑓2 𝑥 𝑓3
0 −12.3333 9.1667 2 20.6667
𝑃𝑖𝑣𝑜𝑡𝑒 𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜 [ | ]
0 −5.7778 6.5556 8.6667 12.8889
𝑑𝑒 − 5.7778 𝑎 − 12.3333
0 −1.1111 2.2222 4.6667 3.5556
⟶ −18 4 1 −6 −38
𝑓3 = 𝑓3 − 𝑓1 . −5.7778/−12.3333 [ 0 −12.3333 9.1667 2
|
20.6667
]
𝑓4 = 𝑓4 − 𝑓1 . −1.1111/−12.3333 0 0 2.2613 7.7297 3.2072
0 0 1.3964 4.4865 1.6937
𝑐3 𝑥 𝑐4 −18 4 −6 1 −38
0 −12.3333 2 9.1667 20.6667
𝑃𝑖𝑣𝑜𝑡𝑒 𝑐𝑎𝑚𝑏𝑖𝑎𝑑𝑜 [ | ]
0 0 7.7297 2.2613 3.2072
𝑑𝑒 2.2613 𝑎 7.7297
0 0 4.4865 1.3964 1.6937
−18 4 −6 1 −38
⟶ 0 −12.3333 2 9.1667 20.6667
𝑓4 = 𝑓4 − 𝑓1 . 4.4865/7.7297 [ 0 0
|
7.7297 2.2613 3.2072
]
0 0 0 0.0839 −0.1678
Considerando los siguientes cambios de columnas:
Distribución original 𝒄𝟏 𝒙𝒄𝟒 𝒄𝟑 𝒙 𝒄𝟒
𝑥1 𝑥4 𝑥4
𝑥2 𝑥2 𝑥2
𝑥3 𝑥3 𝑥1
𝑥4 𝑥1 𝑥3