UNIERSIDAD POLITECNICA SALESIANA

SEDE QUITO

CARRERA:
INGENIERIA MECATRONICA

METODOS NUMERICOS Y ELEMENTOS FINITOS

TEMA:
SOLUCIÓN NUMÉRICA DE SISTEMAS DE ECUACIONES
LINEALES

REALIZADO POR:
JUAN DAVID BENAVIDES AVILA
AGUILAR SILVA SANTIAGO ANDRES
PABLO DAVID MIRANDA NUÑEZ

TUTOR:
ING. LORENA ALEJANDRA FERNANDEZ YANEZ

GRUPO 1
NIVEL 5

Quito, 20 de noviembre del 2020

 TAREA – UNIDAD 1

TEMA: SOLUCIÓN NUMÉRICA DE SISTEMAS DE

ECUACIONES LINEALES

TRABAJO GRUPAL

INDICACIONES
 Elaborar el encabezado, colocando todos los datos de los integrantes
del grupo
 Desarrollar cada una de las preguntas planteadas con su respectivo
proceso
 Aquellos ejercicios donde requiera generar programas debe constar lo
siguiente

 El código del programa copiado desde el Matlab (Octave). Si utiliza

uno por método coloque al inicio de todos los ejercicios y explique su
funcionamiento. Si utiliza un código diferente por ejercicio, cada
código con su explicación debe estar con cada ejercicio.
 El código debe contener los comentarios necesarios y suficientes para
entender su funcionamiento.
 En cada ejercicio debe constar el procedimiento para obtener las
ecuaciones en caso de que no estén definidas explícitamente.
 En cada ejercicio debe presentarse una captura de pantalla con la
corrida del programa y el resultado.

 Una vez finalizado, debe escanear y cargar al AVAC

 La valoración es de 10 puntos.
PARTE 1. Cifras Significativas y Errores
1.-Convierta los números que se presentan en base 2 a base 10.

a) 101101

b) 110.011

c) 0.01101

2. Realice el siguiente ejercicio de cifras significativas

 a) Evalúe el polinomio 𝑦 = 𝑥3 − 5𝑥2 + 6𝑥 + 0.55 en x=1.37. Utilice tres cifras
significativas en el resultado. Evalúe el error relativo porcentual.

b) Repita el inciso a) pero exprese y como 𝑦 = ((𝑥 − 5) 𝑥 + 6) 𝑥 + 0.55. Evalúe el error

y compare con el del inciso a).

3. Ejercicio 4.1 literal b) Libro de Chapra.

b) Use la serie de Taylor para estimar 𝒇(𝒙) = 𝒆−𝒙 en xi+1 = 1 para xi = 0.2. Emplee
versiones de cero, primero, segundo y tercer orden, y calcule |εt| para cada caso.

4. Ejercicio 4.5 Libro de Chapra.

Emplee la expansión de la serie de Taylor de cero hasta tercer orden para predecir f(3)
si f(x) = 25x3 – 6x2 + 7x – 88

5. Ejercicio 4.6 Libro de Chapra.

Use la expansión de la serie de Taylor de cero al cuarto orden para estimar f(2.5) si f(x)
= ln x utilizando x = 1 como punto base. Calcule el error relativo porcentual verdadero
et para cada aproximación. Analice los resultados.
PARTE 2. Sumatorias
6. Genere programas para resolver las siguientes sumatorias, con k=1 donde
corresponda.
Código
% PARTE 2 SUMATORIAS Ejercicio 6
clear all;
clc;
syms X
f=input('Valor Final: ');
n=input('Valor Inicial: ');
A=input('Ingrese la sumatoria: ');
B=0;
% Calculo de la sumatoria
for t=n:f
X=t; %Reemplaza el valor de X de uno en uno hasta llegar a n
B=B+subs(A,X); %Acumula el valor de las sumatorias en una variable B
end
fprintf('La sumatoria es: %f',B);

Breve descripción del funcionamiento de la programación:

Se utilizó un solo programa para la resolución de todos los literales.

El usuario debe ingresar el valor final de las sumatorias, el valor inicial y la sumatoria que se desea realizar.
Dentro de un bucle for se realizarán las sumatorias ingresadas por el usuario, por cada repetición el valor de
X (variable simbólica) tomara un nuevo valor, dentro del bucle en una variable definida como B se
almacenarán los resultados, por cada repetición la variable B suma el valor anterior al valor actual

EJEMPLO DE EL PROGRAMA CORRIDO:

 SUMATORIA RESULTADOS

A ∑ 2𝑘
𝑘=1

B ∑(3𝑘 + 1)
𝑘=1

C ∑(6𝑘 − 5)
𝑘=1

D ∑ 𝑗(𝑗 + 1)
𝑗=1

E ∑(2𝑗 − 1)2
𝑗=1

F ∑ 2𝑗(2𝑗 + 2)
𝑗=1

G ∑(2𝑖 − 1)3
𝑖=1

H ∑ 𝑖(𝑖 + 1)(𝑖 + 2)
𝑖=1

𝑛
I
∑ 𝑖 2 (2𝑖 + 3)
𝑖=1
7.- Dados los valores de 𝑥1 = 2; 𝑥2 = 5; 𝑥3 = 9, 𝑥4 = −9, 𝑥5 = 0, 𝑥6 = −5. Realice un programa para
calcula el valor de numérico de:

Código:
%% PARTE 2 SUMATORIAS Ejercicio 7

clear all;
clc;
syms X
Y=[];
f=input('Valor Final: ');
n=input('Valor Inicial: ');

%Valore para X1, X2, X3

for t=n:f
Y(t)=input('Ingrese el valor de X: ');
end

A=input('Ingrese la Función: ');

B=0;

% Cálculo de la sumatoria

for t=n:f
X=Y(t); %Reemplaza el valor de X de uno en uno hasta llegar a n
B=B+subs(A,X) %Acumula el valor de las sumatorias en una variable B
end
fprintf('La sumatoria es: %f',B);

Breve descripción del funcionamiento de la programación:

Se utilizó un solo programa para la resolución de todos los literales.

El usuario debe ingresar el valor final de las sumatorias, el valor inicial de la sumatoria.

Dentro de un bucle for se generara un arreglo donde se ingresaran los valores de X en este caso el número de
valores para X es el número de repeticiones de la sumatoria

Posteriormente se debe ingresar la sumatoria a realizar

Dentro de un bucle for se realizarán las sumatorias ingresadas por el usuario, por cada repetición el valor de
X (variable simbólica) tomara un nuevo valor (Valores de X ingresados por el usuario), dentro del bucle en
una variable definida como B se almacenarán los resultados, por cada repetición la variable B suma el valor
anterior al valor actual.
EJEMPLO DE LA CORRIDA DEL PROGRAMA:

SUMATORIA RESULTADOS

A ∑(𝑥𝑖 − 2)
𝑖=1

La segunda repetición de la sumatoria es indefinida debido a que el

denominador es igual a cero
𝑥2 =5
5+2
0

3
𝑥𝑖 + 2
∑
B (𝑥𝑖 − 5)2
𝑖=1
 La quinta repetición de la sumatoria es indefinida debido a que el
denominador es igual a cero.
𝑥2 =0
02 − 4
0

C 𝑥𝑖2 − 4
∑
𝑥𝑖
𝑖=1

PARTE 3. Métodos de solución de Sistemas de Ecuaciones Lineales.

Halle la solución a los siguientes ejercicios usando
a) Método de Eliminación de Gauss
Breve descripción del funcionamiento de la programación:
Para la resolución de los sistemas de ecuaciones se utilizó el método de Gauss simple y Gauss Seidel a través
de un programa generado en Matlab.
Para Gauss simple se ingresa una matriz A donde se almacenan los coeficientes y una matriz B que almacena
los términos independientes.
El programa analiza si la matriz A es cuadrada caso de no serlo envía un mensaje indicando que no lo es, en
caso de que A es cuadrada el programa realiza una eliminación hacia adelante y el remplazo de valores hacia
atrás.
Para Gauss Seidel se ingresa una matriz A donde se almacenan los coeficientes y una matriz B que almacena
los términos independientes, una matriz donde se almacenan valores iniciales, una tolerancia de error y un
número de iteraciones.
El programa analiza si la matriz A es cuadrada caso de no serlo envía un mensaje indicando que no lo es, en
caso de que A es cuadrada el programa realiza el despeje de las variables de la fila dominante e ira remplazando
los valores iniciales hasta conseguir un valor aproximado a la respuesta exacta (depende de la tolerancia de
error).
En caso de que el programa supere un número de iteraciones ingresado por el usuario saldrá un mensaje de
error expresando que se ha sobrepasado el número máximo de iteraciones
%% PARTE 3. Métodos de solución de Sistemas de Ecuaciones Lineales
%Gauss simple

clear all;
clc;
A=input('MATRIZ A: ');
b=input('MATRIZ B:');
[C,F]=size(A);
[n,~]=size(A);
x=zeros(n,1);
MC=['X';'Y';'Z';'T';'W';'U'];

% Valores para imprimir caracteres X Y Z T W U

ES=[' '];
MA=[];

% Verifica que la matriz A sea cuadrada

if C~=F
disp('Matriz A no es cuadrada')
end
if A(1,1)==0
disp('Ingrese la matriz de tal forma que el valor A(1,1) sea diferente de
cero');
end
if C==F
for i=1:C
MA(i,1)=MC(i,1);
ES(i,1)=ES(1,i);
end
for i=1:n-1
m= A(i+1:n,i)/A(i,i);
A(i+1:n,:) = A(i+1:n,:) - m*A(i,:);
b(i+1:n,:) = b(i+1:n,:) - m*b(i,:);
end
x(n,:) = b(n,:)/A(n,n);
for i=n-1:-1:1
x(i,:) = (b(i,:)-A(i,i+1:n)*x(i+1:n,:))/A(i,i);
end
x;
x=num2str(x);
x=[MA,ES,x]
end

Capturas de pantalla de los resultados del código:

MATRIZ CAPTURA DE PANTALLA

2𝑥 − 𝑦 + 3𝑧 = 5
A { 2𝑥 + 2𝑦 + 3𝑧 = 7
−2𝑥 + 3𝑦 = −3
 𝑥 + 2𝑦 + 3𝑧 = 5
B { −3𝑥 − 2𝑦 − 𝑧 = 7
4𝑥 + 4𝑦 + 4𝑧 = −3

C 3𝑥 − 𝑦 + 7𝑧 = 1
{
5𝑥 +𝑧 =2

𝑥 − 𝑦 + 5𝑧 = √2
√5𝑥 + 𝑧 = √3
D 2 5
𝑥 + 3𝑦 + 2𝑧 =
{5 2

𝑥 + 2𝑦 − 3𝑧 − 𝑡 = 0
−3𝑦 + 2𝑧 + 6𝑡 = −8
E −3𝑥 − 𝑦 + 3𝑧 + 𝑡 = 0
{2𝑥 + 3𝑦 + 2𝑧 − 𝑡 = −8
F (1 + 𝑖 )𝑥 − 𝑖𝑦 = 3𝑖
{
2𝑥 + 𝑖𝑦 = −𝑖

𝑥 + 2𝑦 + 3𝑧 = 5
−3𝑥 − 2𝑦 − 𝑧 = 7
4𝑥 + 4𝑦 + 4𝑧 = −3
 b) Método de Gauss Seidel

%% PARTE 3. Método de Gauss Seidel

clear all
clc
%Método de resolución de Gauss Seidel
fprintf('Sistema de eliminación de Gauss Seidel \n')
A=input('Matriz A:');%Matriz de coeficientes
B=input('Matriz B:');%Matriz de términos independientes
X0=input('Ingrese los valores iniciales:');%Valores iniciales para remplazar
tol=input('Ingrese la tolerancia:');%Tolerancia
maxiter=input('Ingrese el valor máximo de iteraciones:');%Limite para las
iteraciones
N=length(B);

%incr es el valor criterio de parada, que se compara con la tolerancia para

obtener el resultado más aproximado

iter=0;
incr=tol+1;
v=[]; %Para ir guardando los incrementos por iteración

%Verifica que A sea cuadrado

[C,F]=size(A);
if C~=F
disp('Matriz A no es cuadrada')
end

%En caso de que A sea cuadrada

if C==F
while iter <= maxiter && incr > tol
for j=1:N
if j==1
X(1)=(B(1)-A(1,2:N)*X0(2:N))/A(1,1);
elseif j==N
X(N)=(B(N)-A(N,1:N-1)*X(1:N-1)')/A(N,N);
else
%X tiene la aproximación k-enésima
%X0 tiene la aproximación k-1-enesima
X(j)=(B(j)-A(j,1:j-1)*X(1:j-1)'-A(j,j+1:N)*X0(j+1:N))/A(j,j);
end
end
incr=norm(X'-X0,inf)/norm(X',inf);
v=[v,incr];
iter=iter+1;
X0=X';
end
A
%X es la matriz solución aproximada generada por el método iterativo de Gauss
Seidel
X=X'
fprintf('El número de iteraciones es:%f',iter);
if iter > maxiter
disp('Error, usted sobrepasó el número máximo de iteraciones y el error es
mayor que la tolerancia');
end
end
Capturas de pantalla de los resultados del código:
MATRIZ CAPTURA DE PANTALLA

2𝑥 − 𝑦 + 3𝑧 = 5
A { 2𝑥 + 2𝑦 + 3𝑧 = 7
−2𝑥 + 3𝑦 = −3

C { 3𝑥 − 𝑦 + 7𝑧 = 1
5𝑥 +𝑧 =2
 𝑥 − 𝑦 + 5𝑧 = √2
√5𝑥 + 𝑧 = √3
D 2 5
𝑥 + 3𝑦 + 2𝑧 =
{5 2

𝑥 + 2𝑦 − 3𝑧 − 𝑡 = 0
−3𝑦 + 2𝑧 + 6𝑡 = −8
−3𝑥 − 𝑦 + 3𝑧 + 𝑡 = 0
E {2𝑥 + 3𝑦 + 2𝑧 − 𝑡 = −8

𝑥 + 2𝑦 + 3𝑧 = 5
{ −3𝑥 − 2𝑦 − 𝑧 = 7
4𝑥 + 4𝑦 + 4𝑧 = −3
14. En una columna de cinco platos, se requiere absorber benceno que está contenido en una corriente
de gas V, con aceite L que circula a contracorriente del gas. Considérese que el benceno transferido no
altera sustancialmente el número de moles de V y L fluyendo a contracorriente, que la relación de
Equilibrio está dada por la ley de Henry (𝛾 = 𝑚𝑥) y que la columna opera a régimen permanente. Calcule
la composición de benceno en cada plato.
Datos:
V= 100 moles/min
L= 500 moles /min
𝛾𝑜 = 0.09 fracción molar de benceno en V
𝑥𝑜 = 0.00 fracción molar de benceno en L
𝑚 = 0.12 Los balances de materia para el benceno en cada plato son

ECUACION DESPEJE

500(0 − 𝑥5 ) + 100(0.12𝑥4 − 0.12𝑥5 ) = 0

𝐿 (𝑥0 − 𝑥5 ) + 𝑉 (𝑦4 − 𝑦5 ) = 0 −500𝑥5 + 12𝑥4 − 12𝑥5 = 0
−512𝑥5 + 12𝑥4 = 0

500(𝑥5 − 𝑥4 ) + 100(0.12𝑥3 − 0.12𝑥4 ) = 0

𝐿 (𝑥5 − 𝑥4 ) + 𝑉 (𝑦3 − 𝑦4 ) = 0 500𝑥5 − 500𝑥4 + 12𝑥3 − 12𝑥4 = 0
500𝑥5 − 512𝑥4 + 12𝑥3 = 0

500(𝑥4 − 𝑥3 ) + 100(0.12𝑥2 − 0.12𝑥3 ) = 0

𝐿(𝑥4 − 𝑥3 ) + 𝑉(𝑦2 − 𝑦3 ) = 0 500𝑥4 − 500𝑥3 + 12𝑥2 − 12𝑥3 = 0
500𝑥4 − 512𝑥3 + 12𝑥2 = 0
 500(𝑥3 − 𝑥2 ) + 100(0.12𝑥1 − 0.12𝑥2 ) = 0
𝐿(𝑥3 − 𝑥2 ) + 𝑉(𝑦1 − 𝑦2 ) = 0 500𝑥3 − 500𝑥2 + 12𝑥1 − 12𝑥2 = 0
500𝑥3 − 512𝑥2 + 12𝑥1 = 0

500(𝑥2 − 𝑥1 ) + 100(0.09 − 0.12𝑥1 ) = 0

𝐿(𝑥2 − 𝑥1 ) + 𝑉(𝑦0 − 𝑦1 ) = 0 500𝑥2 − 500𝑥1 + 9 − 12𝑥1 = 0
500𝑥2 − 512𝑥1 = −9

−512𝑥5 + 12𝑥4 + 0 + 0 + 0 = 0
500𝑥5 − 512𝑥4 + 12𝑥3 + 0 + 0 = 0
ECUCIONES 0 + 500𝑥4 − 512𝑥3 + 12𝑥2 + 0 = 0
0 + 0 + 500𝑥3 − 512𝑥2 + 12𝑥1 = 0
0 + 0 + 0 + 500𝑥2 − 512𝑥1 = −9

ECUACION PARA REEMPLAZAR EN EL CODIGO

% Matriz A: [-512 12 0 0 0; 500 -512 12 0 0; 0 500 -512 12 0; 0 0 500 -

512 12; 0 0 0 500 -512]
% Matriz B:[0;0;0;0;-9]

CAPTURA RESPUESTAS EQUIVALENCIA

𝑥1 = 𝑊 = 0.018
𝑥2 = 𝑇 = 0.00043199
𝑥3 = 𝑍 = 1.0356−5
𝑥4 = 𝑌 = 2.373−7
𝑥5 = 𝑋 = −5. .56 17−9

NOTA: El código utilizado para resolver este ejercicio fue el mismo utilizado en la PARTE 3. Métodos de
solución de Sistemas de Ecuaciones Lineales, por el método de eliminación de Gauss Simple.
15. Un bloque de 10Kg está suspendido en un sistema de cables y anillos, tal como
muestra la figura. Determine la tensión en cada cable del sistema
 ECUACIONES
FCF FCB FCA FBE FAB FAD
3/5 24/25 -1 - - - 0
4/5 -7/25 - - - - 0
-3/5 - - 5/13 - - 0
-4/5 - - -12/13 1 - 0
- -24/25 - -5/13 - 15/17 0
- 7/25 - 12/13 - 8/17 98.1

MATRIZ PARA REEMPLAZAR EN EL CODIGO

% Matriz A: [3/5 24/25 -1 0 0 0;4/5 -7/25 0 0 0 0;-3/5 0 0

5/13 0 0; -4/5 0 0 -12/13 1 0;0 -24/25 0 -5/13 0 15/17;0 7/25
0 12/13 0 8/17]
% Matriz B: [0; 0; 0; 0; 0; 98.1]

CAPTURA DE PANTALLA DE LAS

RESPUESTAS

NOTA: El código utilizado para resolver este ejercicio fue el mismo utilizado en la PARTE 3. Métodos de
solución de Sistemas de Ecuaciones Lineales, por el método de eliminación de Gauss Simple.

16. Determine las concentraciones molares de una mezcla de cinco componentes en

solución a partir de los siguientes datos espectrofotométricos.
Asuma que la longitud de la trayectoria óptica es unitaria y que el solvente no absorbe a estas longitudes de
onda Considere la ley de Beer para generar las ecuaciones del sistema.
Se considera que se cumple la ley de Beer, entonces a una longitud de onda, i.
5

𝐴𝑟𝑜𝑟𝑖 = ∑ 𝜀𝑖𝑗 𝐶𝑗
𝑗=1

Donde:
𝐴𝑟𝑜𝑟𝑖 Es la absorbancia total observada a la longitud de onda i.
𝜀𝑖𝑗 Es la absorbida mola del componente j a la longitud de onda i.

𝐶𝑗 Es la concentración molar del componente j en la mezcla.

C1 C2 C3 C4 C5
1 98 9 2 1 0.5 0.1100
2 11 118 9 4 0.88 0.2335
3 27 27 85 8 2 0.2800
4 1 3 17 142 25 0.3000
5 2 4 7 17 118 0.1400

ECUACIONES CAPTURA DE PANTALLA

98𝐶1 + 9𝐶2 + 2𝐶3 + 𝐶4 + 0.5𝐶5 = 0.1100

11𝐶1 + 118𝐶2 + 9𝐶3 + 4𝐶4 + 0.88𝐶5 = 0.2335
27𝐶1 + 27𝐶2 + 85𝐶3 + 8𝐶4 + 2𝐶5 = 0.2800
𝐶1 + 3𝐶2 + 17𝐶3 + 142𝐶4 + 25𝐶5 = 0.3000
2𝐶1 + 4𝐶2 + 7𝐶3 + 17𝐶4 + 118𝐶5 = 0.1400

MATRIZ PARA REEMPLAZAR EN EL CODIGO

Matriz A: % [98 9 2 1 0.5; 11 118 9 4 0.88; 27 27 85 8 2; 1 3 17 142 25; 2 4
7 17 118]
Matriz B: % [0.11; 0.2335; 0.28; 0.3; 0.14]

NOTA: El código utilizado para resolver este ejercicio fue el mismo utilizado en la PARTE 3. Métodos de
solución de Sistemas de Ecuaciones Lineales, por el método de eliminación de Gauss Simple.
17. Determine la intensidad de corriente en cada rama del circuito mostrado

Los valores de tensiones, resistencias internas de fuente y resistencias del circuito se muestran en la siguiente
tabla.
 ECUACIONES
I1 I2 I3 I4 I5 I6
61 - - - - 20 -6
- - 53 - - 9 9.8
20 - 9 - - 49 -6
- - 1 - -1 1 0
1 1 -1 - - - 0
-1 - - -1 - -1 0

MATRIZ ORDENADA CAPTURA DE PANTALLA

ECUACION PARA REEMPLAZAR EN EL CODIGO.

% MATRIZ A: [61 0 0 0 0 20;1 1 -1 0 0 0;20 0 9 0 0 49;-1 0 0 -1 0 -1;0 0 1 0 -1 1;0 0

53 0 0 9]
% MATRIZ B:[-6;0;-6;0;0;9.8]

NOTA: El código utilizado para resolver este ejercicio fue el mismo utilizado en la PARTE 3. Métodos de
solución de Sistemas de Ecuaciones Lineales, por el método de eliminación de Gauss Simple.