Informe Parcial ME

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

Parcial 1 - Matematicas especiales

2021-2
Integrantes: Yuranny Yeraldin Morillo Santacruz - Sara Hoyos Guisao - José Alberto Gómez
Echeverri
Problema [Tomado de Barrera.R & Díaz.L (2020). Métodos numéricos en Excel y Matlab Con
aplicaciones en ingeniería (1ra ed). Universidad De Antioquia.]
Como práctica académica, usted opta por la modalidad de trabajo social o servicio comunitario, y su
primera tarea es ayudar a diseñar un tanque de almacenamiento de agua (esférico) para un pequeño
poblado en las afueras de cierta cabecera municipal. Dadas las características del lugar, es imperativo
3
que el tanque esférico contenga 30 𝑚 de agua, pero las condiciones semiprecarias (limitada
tecnología y recursos disponibles) solo le permiten medir o controlar con certeza la altura del líquido
(nivel) dentro del tanque. Afortunadamente, de los primeros cursos en la universidad, usted sabe bien
que la siguiente ecuación le será muy útil:

(1)

Con V = volumen del tanque esférico (m3), h = altura o nivel del líquido en el tanque (m), r = radio
(m) y π = 3.1416. Utilice alguno de los métodos numéricos vistos en este capítulo para determinar
(redondeado a cuatro cifras decimales) cuál debe ser la altura o nivel del líquido que garantice el
3
volumen de 30𝑚 para r=2.08m.
Observación: para este ejercicio, el problema numérico presenta tres raíces reales, pero solo una se
ajusta a la física del problema que exige h ≤ 2r.
Resultados.
A continuación en la Tabla 1, se encuentran los resultados obtenidos utilizando cada método:
Tabla 1. Resultados

Discusión.

Para comenzar a evaluar todos lo métodos, se tomó la ecuación (1) dada por el problema, en esta se
2 3𝑟−ℎ
hizo un despeje teniendo asi la ecuacion πℎ ( 3
)-V=0 (2) esto por facilidad en los cálculos,

seguido de esto y con los datos entregados se hace el reemplazo en (2) y se tiene que
2 3(2.08)−ℎ
3. 1416ℎ ( 3
)-30=0 (3) pero podemos decir que esta ecuación es la función f que depende de
2 3(2.08)−ℎ
h, f(h)=3. 1416ℎ ( 3
)-30 (4).

Se procede a crear la tabla en excel y a seguir los pasos del método de Bipartición o Bisección
explicado en Solución de ecuación no lineal ☺ Método de Bisección/Bipartición ☺. Ejemplo en
Excel(2019, 1 diciembre)YouTube (Enlace: Solución de ecuación no lineal ☺ Método de
Bisección/Bipartición ☺ Ejemplo en Excel).
Para el ejercicio se toma valores iniciales 𝑎0=2 y 𝑏0 =4, al realizar 20 iteraciones, se observa que F(𝑎0

) es negativo, F(ℎ0) es positivo y F(𝑏0) es positivo, es decir, el cambio de signo (la raíz) se presenta

dentro del subintervalo [𝑎0,𝑥0] ([2, 3]) y, por lo tanto, el subintervalo [ℎ0, 𝑏0] ([3, 4]) puede

descartarse. De este modo, el nuevo intervalo de búsqueda [𝑎1, 𝑏1] se puede definir como 𝑎1 = 2 ,

comenzando una nueva etapa de cálculo para este nuevo intervalo [𝑎1, 𝑏1] = [2, 3]. Sucesivamente se

siguen realizando los cálculos hasta que el intervalo de búsqueda sea pequeño lo que permite que se
encuentre la raíz con una aproximación razonable, siendo en la iteración 20 el valor de la altura de
2,951636 metros lo que cumple con las condiciones del problema y al graficar la función con ayuda
de un recurso informático como es symbolab se puede observar que los valores son semejantes.
Figura 1. Gráfica de la función del problema 1 en función de la altura. [tomada de:Calculadora
gráfica - Symbolab]
Para hallar la raíz que soluciona el problema por el método de falsa posición se procede a crear la
tabla en excel y seguir los pasos del método explicado en Solución de ecuación no lineal ☺ Método
Falsa posición/Regla falsa(Regula falsi) ☺ Ejemplo en Excel(2019, 1 diciembre).YouTube. (Enlace:
Solución de ecuación no lineal ☺ Método Falsa posición/Regla falsa(Regula falsi) ☺ Ejemplo en
Excel ).
Para hallar la solución el problema se toma valores iniciales 𝑎0=-0,1 m y 𝑏0 =3,5 m, luego de digitar

los datos y las operaciones correspondientes se decide realizar 20 iteraciones para ver el
comportamiento. Se observa que en la sexta iteración la función evaluada en h (altura) se vuelve cero,
lo que significa que se encontró una raíz del problema es decir la solución ya que cumple con las
condiciones que presenta el problema 1. A partir de la séptima iteración los valores tienden a ser los
mismos en algunos casos y en otros cambia su valor a 0.

Para hallar la raíz que soluciona el problema se emplea la ecuación (4) y por el método de
sustituciones sucesivas o de punto fijo el cual suele ser de rápida convergencia, se procede a crear la
tabla en excel seguir los pasos del método explicado en Solución de ecuación no lineal ☺ Método
Sustituciones sucesivas/Punto fijo ☺ Ejemplo en Excel(2019, 2 diciembre).YouTube. (Enlace:
Solución de ecuación no lineal ☺ Método Sustituciones sucesivas/Punto fijo ☺ Ejemplo en Excel ).
Para hallar la solución el problema se toma como valor inicial ℎ0=-1 el cual no corresponde al

intervalo donde se encuentra la raíz. Al despejar la variable h se obtiene la siguiente función:


Se procede a asignar un valor inicial a la variable (ℎ0 = 1𝑚) y evaluar la función, después de varias

iteraciones realizadas al problema se nota que la diferencia de un valor al anterior no es significativo


lo que se puede decir que converge. Este valor semejante es una nueva raíz de la función
(h=5,16693671 metros en la iteración 24), pero presenta un problema ya que no puede ser la solución
del problema porque no cumple con la condición h ≤ 2r. A partir de la iteración 24 los valores tienden
a ser los mismos. Luego de dar distintos valores iniciales y aumentar el número de iteraciones, los
resultados son los mismos, lo que quiere decir que el método no da una solución válida al problema.

Se procede a crear la tabla en excel y a seguir los pasos del método de secante explicado en Solución
de ecuación no lineal ☺ Método Secante ☺ Ejemplo en Excel(2019, 2 diciembre)YouTube (Enlace:
Solución de ecuación no lineal ☺ Método Secante ☺ Ejemplo en Excel).
Para el ejercicio se toma valores iniciales 6 metros y 7 metros, al realizar 8 iteraciones el valor
evaluado en función de la altura es cero y su ℎ8 = 5, 16687828 𝑚𝑒𝑡𝑟𝑜𝑠 el cual es un valor que

corresponde a una raíz de la función, lo que quiere decir que es una solución para el problema pero
presenta el mismo problema con el anterior método y es que no puede ser la solución del problema
porque no cumple con la condición h ≤ 2r. Si se realizan más iteraciones aparece un error (#¡DIV/0!)
el cual quiere decir que la función al evaluarla será indeterminada.

Método Newton: para el desarrollo del método, se tomó 3 metros como valor inicial o valor semilla,
valor que se obtuvo reemplazando el radio escogido por el grupo de 2.08 m y el volumen requerido de
30 m3 para la construcción del tanque, este valor fue verificado con la gráfica obtenida en una
calculadora virtual, se encontró convergencia a las 3 iteraciones, se usó un criterio de convergencia
εf =1x10-6 y el resultado obtenido para la altura (h) fue de 2.9516 m, valor que cumple con el
parámetro h<2r definido en el problema propuesto.

Homotopias. Para la implementación de estos métodos se siguió el procedimiento explicado en los


videos propuestos para los ejercicios de clase, se tomó 3 m como valor inicial, M= 40 (subintervalos)
y se usó un criterio de convergencia εf =1x10-6. El resultado obtenido para la altura h fue de 2.9516 m
en los tres métodos, se encontró convergencia en 89, 111 y 120 iteraciones para Punto fijo, Afín y
Newton respectivamente.
Ahora, se empieza a crear la tabla y a seguir los pasos del método explicado en Solución de ecuación
no lineal con múltiples raíces☺Método de Deflación con Newton-Raphson en Excel. (2019, 9
diciembre).YouTube.
(https://www.youtube.com/watch?v=qQ-ZTM2uRAM&ab_channel=LUISALBERTODIAZMONTES )
Como resultados se obtuvo que en la iteración 17 se encontró la primera raíz la cual fue 5.16693671
metros pero esta no puede ser la solución del problema ya que no cumple con la condición h ≤ 2r. De
[𝑘]
esta se puede observar que el contador k en la iteración 18 comienza nuevamente en 0 y que ℎ
retoma su valor inicial 0.1. La segunda raíz se encuentra en la iteración 24 con el valor de -1.8784873
metros, esta raíz se descarta porque no tiene sentido físico por ser negativa; en esta ocurre lo mismo
[𝑘]
con la k y ℎ que con la primera raíz encontrada. Por último se encontró la tercer raíz en la iteración
27 la cual arrojó un valor de 2.95155028 metros la cual sí tiene sentido físico por ser positiva y
3
cumple con h ≤ 2r por tanto la respuesta al ejercicio es que para tener un volumen de 30𝑚 es
necesario tener una altura de 2.95155028m.
Como observación, para tener seguridad de que no se encontraran más raíces se hizo un arrastre hasta
la iteración 76 donde se obtuvo un error #¡NUM! donde se consultó en las ayudas de excel y se
descubre que el programa no puede arrojar más valores ya que no son soportados por ser demasiado
pequeños.
Para el método de deflación con buscar objetivo nuevamente usamos la ecuación (4) y basado en las
instrucciones de Solución de ecuación no lineal con múltiples raíces☺Método Deflación con Buscar
Objetivo en Excel. (2019, 9 diciembre). YouTube.(
https://www.youtube.com/watch?v=A7iDHzrFnpo&ab_channel=LUISALBERTODIAZMONTES) se
llena la tabla donde se encontraron las mismas 3 raíces a las halladas con el método anterior, donde
nuevamente podemos afirmar que la solución al problema es 2.95155028m.
Además, con este método se logró comprobar que el problema solo tenía 3 raíces, ya que cuando se
usó el método nuevamente se obtiene la raíz número 3 (5.16693671m) concluyendo que el software
ya no tenía más soluciones por encontrar.
Por último para el uso del complemento homotosolver se basó en Solución de ecuación no lineal con
múltiples raíces☺Uso de complemento Homotosolver en Excel. (2019, 12
diciembre).YouTube(https://www.youtube.com/watch?v=jXNMUZBi_YM&t=225s&ab_channel=LU
ISALBERTODIAZMONTES), se hallan las mismas raíces, además este es el último comprobante del
número de soluciones del ejercicio ya que excel indica el número de raíces que logró encontrar.

También podría gustarte