José Ezequiel Arvizu Dzul Tarea 1 MMAIP

Universidad Iberoamericana Puebla

Trimestre 24I
Primer parcial de Métodos Matemáticos
Nombre del profesor: Dr. Roel Hernández Rodríguez

Nombre del alumno: Arvizu Dzul José Ezequiel

Trimestre que cursa: 24-I
Problema 1 Problema 2 Problema 3
Calificación

Instrucciones: Leer con cuidado cada una de las instrucciones en cada problema. Para cada problema se
debe reportar la solución encontrada y los errores relativos porcentuales, ambos para cada iteración con
cada método utilizado.

Problema 1 (30 puntos)

1. ¿Qué son los métodos numéricos?

Son un conjunto de sucesiones de operaciones matemáticas (algoritmo) ordenadas o lógicas para
encontrar una solución aproximada a un problema determinado.
2. ¿En el contexto de Matlab, que instrucciones permiten tomar decisiones dentro un algoritmo?
Explicar.
Las expresiones condicionales, si la expresión es verdadera, el grupo o bloque de instrucciones se
ejecutan. Si la expresión es falsa, no se ejecuta o se salta el grupo de instrucciones.
3. ¿Qué es una sentencia condicional?
Es una instrucción que permiten tomar decisiones sobre si se ejecuta un tipo de comandos o
instrucciones que cumplen una condición o, por el contrario, omitirlos.
Existen 3 tipos de sentencias condicionales:
1) if – end
Si la expresión es verdad se ejecutan las sentencias y sigue a end. Si la
condición es falsa, el problema salta inmediatamente a la línea después del
comando end.

1
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
2) if – else – end

Permite ejecutar un conjunto de instrucciones u operaciones

si la condición es verdadera y un conjunto diferente de
operaciones si la condición es falsa.

3) if – elseif – else -end

Esta estructura o sentencia incluye dos sentencias
condicionales, lo que hace posible ejecutar 1 de entre 3
grupos de instrucciones diferentes.

Ejemplos:

Ejemplo de las instrucciones 1), 2), 3)

4. En el contexto de Matlab, ¿qué instrucciones permiten realizar operaciones repetitivas?

Ciclos o loops. Repite tantas veces como sea necesario los pasos del algoritmo. Permite repetir o no un
número N de veces un conjunto de operaciones ordenadas.
Existen 2 tipos de ciclos: 1) Ciclo controlado por un controlador y 2) Ciclo terminado por un condicional.
Ejemplo del 1):
for: incrementos, se usa como si fuera un vector, usando un contador: i=1:1:N, estos
incrementos siempre deben ser enteros.

5. ¿Qué significa encontrar la raíz de una ecuación algebraica?

Es encontrar uno o varios puntos en donde la función pase por cero, esto significa que encontramos una
solución a esa ecuación.

2
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
6. ¿Qué es una serie de Taylor y para qué sirve?

Es una expresión matemática que permite aproximar una función 𝑓(𝑥) mediante la combinación de
derivadas de la función y una ley de potencia de la forma (𝑥 − 𝑥0)

Teorema de Taylor: si una función 𝑓(𝑥) y sus 𝑛 + 1 derivadas son continuas en un intervalo que
contiene los puntos 𝑥 y 𝑥0 , el valor de la función en algún valor de 𝑥 se puede aproximar con
respecto o alrededor del punto 𝑥0 como:
∞
𝑑𝑛𝑓 (𝑥 − 𝑥0 )𝑛
𝑓( 𝑥 ) = ∑ | =⋯
𝑑𝑥 𝑥 0
2
𝑛!
𝑛=0

7. ¿Qué es una subrutina?

Es un programa adyacente al programa principal que tiene una forma diferente de iniciarse además de
tener tareas muy particulares (como escribir una función para no tener que repetirla las veces que la
necesitemos, en lugar de esto sólo la llamamos), tiene la siguiente estructura:
Function[Datossalida]=nombre_de_subrutina(Datosentrada)

Llamada a subrutina. Puede hacerse de dos maneras:

1) [Datossalida]=nombre_de_subrutina(Datosentrada)
2) [Datossalida]=feval(‘nombre_de_subrutina’,Datosentrada
8. Define que es exactitud y que es precisión en métodos numéricos.
Exactitud: es la diferencia entre un valor aproximado y un valor real asociado.
Precisión: que tan cerca están las aproximaciones unas de otras. Indica la cantidad de cifras
significativas en una estimación.
9. Como se define el error absoluto y el error relativo.

Error Absoluto: Sea 𝑥 el valor real de un número, y 𝑥∗ el valor aproximado del mismo, se denota al
error absoluto (mide la exactitud) como:

𝐸𝑥∗ = |𝑥 − 𝑥∗|

Error Relativo: Sea 𝑥 el valor real de un número, y 𝑥∗ el valor aproximado del mismo, se denota al error
relativo como:

|𝑥 − 𝑥∗|
𝜀𝑥∗ = %𝜀𝑥∗ = 𝜀𝑥∗ ∗ 100
|𝑥 |

3
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
Problema 2 (35 Puntos)
El coeficiente de fricción p para el flujo turbulento en un tubo está dado por la siguiente ecuación, la cual
es una ecuación algebraica no lineal
1
1.14  2
log  9.35  (Correlación de Colebrook)
 
e 
10  0 Re
p D p


donde Re es el número de Reynolds, e es la rugosidad de la superficie del tubo y D es el diámetro del

tubo. Usando el método gráfico, y el método de bisección, determine el valor del factor de fricción p
para los dos siguientes casos

i. D  0.1m , e  0.0025 , Re  3x104

ii. D  0.1m , e  0.0001, Re  5x106
En el segundo método el de error relativo porcentual debe ser menor que 0.1% . Como resultados de este
ejercicio se espera tener

 El programa en MATLAB que resuelve el problema mediante el método gráfico

 La gráfica de cada caso donde se señale la solución mediante el método gráfico
 El porcentaje de error relativo como función de las iteraciones para el segundo y tercer método
 El valor de p final, el valor del error relativo porcentual final, el número de iteraciones finales
necesarias para que el error sea menor que la tolerancia.
NOTA: Para presentar los resultados de los métodos iterativos puede utilizar una tabla, como las usadas
en clase.
Problema 2 inciso i)
Línea de código:

clear all, close all, clc

%inciso 1

p=0.001:0.00001:0.1;
D=0.1; %metros (Diámetro del tubo)
e=0.0025; %adimensional (Rugosidad)
Re=30000; %adimensional (No. de Reynolds)

fp=(1./sqrt(p))-1.14+2*log10((e./D)+(9.35./Re.*sqrt(p))); %Correlación de Colebrook

figure(1)
plot(p,fp,'g'), grid on
xlabel('coeficiente de ficción (p)'), ylabel('f(p)'), title('Probelma 2');

4
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
Con el método gráfico se obtiene lo siguiente:

Gráfica 1. Método gráfico problema 2 inciso i)

El valor aproximado con el método gráfico es 𝒑 ≈ 𝟎. 𝟎𝟓𝟑𝟎𝟏

Con el método de bisección:
Como se observa en la gráfica 1 se sabe que existe una raíz en el intervalo [0.05, 0.06]

Iteración % error
ai bi mi f(ai) f(bi) f(mi)
i relativo
1 0.05 0.06 0.0550 0.1304 -0.259 -0.07757
2 0.05 0.0550 0.0525 0.1304 -0.07757 0.02272 4.76
3 0.0525 0.0550 0.0537 0.02272 -0.07757 -0.0263 2.23
4 0.0525 0.0537 0.0531 0.02272 -0.0263 -0.001998 1.13
5 0.0525 0.0531 0.0528 0.02272 -0.001998 0.01031 0.57
6 0.0528 0.0531 0.05295 0.01031 -0.001998 0.00414097 0.28

Para la sexta iteración se tiene que 𝒑 ≈ 𝟎. 𝟎𝟓𝟐𝟗𝟓 con error relativo de 0.28%

5
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
Problema 2 inciso ii)
Línea de código:

clear all, close all, clc

%inciso 2

p=0.001:0.0001:0.1;
D=0.1; %metros (Diámetro del tubo)
e=0.0001; %adimensional (Rugosidad)
Re=5000000; %adimensional (No. de Reynolds)

fp=(1./sqrt(p))-1.14+2*log10((e./D)+(9.35./Re.*sqrt(p))); %Correlación de Colebrook

figure(1)
plot(p,fp,'cyan'), grid on
xlabel('coeficiente de ficción (p)'), ylabel('f(p)'), title('Probelma 2');

Con el método gráfico se obtiene lo siguiente:

Gráfica 2. Método gráfico problema 2 inciso ii)

El valor aproximado con el método gráfico es 𝒑 ≈ 𝟎. 𝟎𝟏𝟗𝟔

Con el método de bisección:
Como se observa en el la gráfica 2 se sabe que existe una raíz en el intervalo [0.01, 0.03]

6
José Ezequiel Arvizu Dzul Tarea 1 MMAIP

% error
Iteración i ai bi mi f(ai) f(bi) f(mi)
relativo
1 0.01 0.03 0.02 2.86016 -1.36622 -0.0687025
2 0.01 0.02 0.015 2.86016 -0.0687025 1.02516 33.33
3 0.015 0.02 0.0175 1.02516 -0.0687025 0.419504 14.29
4 0.0175 0.02 0.01875 0.419504 -0.0687025 0.16319 6.67
5 0.01875 0.02 0.019375 0.16319 -0.0687025 0.0444381 3.23
6 0.019375 0.02 0.0196875 0.0444381 -0.0687025 0.01280057 1.60

Para la sexta iteración se tiene que 𝒑 ≈ 𝟎. 𝟎𝟏𝟗𝟔𝟖𝟕𝟓 con error relativo de 1.60%

Problema 3 (35 Puntos)

La ecuación de Ergun es una ecuación frecuentemente utilizada para describir el flujo de un fluido en
medios porosos o lechos empacados, la cual está dada por
P Dp 
3

 150 1 1.75

Go2 L 1   DpGo


En esta ecuación, P es la diferencia de presión,  es la densidad, Go es la velocidad másica (el cociente
del flujo de masa dividido entre el área de la sección transversal de flujo), D es el diámetro de las
p

partículas dentro del cauce, es la viscosidad dinámica del fluido, L es la longitud del medio poroso
Go
o del lecho empacado y  es la porosidad del medio poroso o lecho empacado. Dados los siguientes
valores para los siguientes grupos de parámetros en la ecuación,
P Dp
 1000
2
Go L
DpGo
 10


encuentre la porosidad del lecho empacado usando el método gráfico y el método de bisección. En el
segundo método, el de error relativo porcentual debe ser menor que 0.1% %. Como resultados de este
ejercicio se espera tener

 El programa en MATLAB que resuelve el problema mediante el método gráfico

 La grafica donde señale la solución
 El porcentaje de error relativo como función de las iteraciones usando el método de bisección.

7
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
 El valor de  final, el valor del error relativo porcentual final, el número de iteraciones final para
alcanzar ese valor.

8
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
Problema 3
La ecuación de Ergun igualada a cero, con los parámetros sustituidos es la siguiente:

𝜀3
(1000) + 15𝜀 − 16.75 = 0
1−
𝜀
Debido a que era complicado leer el valor de 𝜀 con el método gráfico, decidí hacer una gráfica del tipo
Log(y) vs x para poder leer mejor los valores, a continuación, se muestra la línea de código,
posteriormente se muestra la gráfica 3 en donde se muestra el resultado del método gráfico.
Subrutina:

function[y]=ec_Ergun(epsilon)
y=(1000.*(epsilon.^3./(1-epsilon)))+(15.*epsilon)-16.75;

Programa principal:

clear all, close all, clc

epsilon=0.0001:0.0001:0.9;

[y]=ec_Ergun(epsilon);

figure(1)
semilogy (epsilon,y,'m'), grid on
xlabel('Porosidad (epsilon)'), ylabel('Log[f(epsilon)]'), title('Problema 3');

Gráfica 3. Método gráfico del problema 3

9
José Ezequiel Arvizu Dzul Tarea 1 MMAIP
Como se observa en el la gráfica 2, existe una raíz en el intervalo [0.0, 0.3]

% error
i ai bi mi f(ai) f(bi) f(mi)
relativo
1 0 0.3 0.15 -16.75 26.3214 -10.5294
2 0.15 0.3 0.225 -10.5294 26.3214 1.32258 33.33
3 0.15 0.225 0.1875 -10.5294 1.32258 -5.82452 20.00
4 0.1875 0.225 0.20625 -5.82452 1.32258 -2.60279 9.10
5 0.20625 0.225 0.215625 -2.60279 1.32258 -0.73438 4.35
6 0.215625 0.225 0.2203125 -0.73438 1.32258 0.26972 2.13
7 0.215625 0.2203125 0.21796875 -0.73438 0.26972 -0.238316005 1.08
8 0.21796875 0.2203125 0.219140 -0.238316005 0.26972 0.014055536 0.53
9 0.21796875 0.219140 0.218854375 -0.238316005 0.014055536 -0.047765982 0.13
10 0.218854375 0.219140 0.2189971875 -0.047765982 0.014055536 -0.01687764 0.06

Para la decima iteración se tiene que 𝜺 ≈ 𝟎. 𝟐𝟏𝟖𝟗𝟗𝟕𝟏𝟖𝟕𝟓 con error relativo de 0.06%

1
0