UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO.
FACULTAD DE INGENIERÍA.
ANÁLISIS NUMÉRICO.
GRUPO:1
ALUMNOS:
LÓPEZ DUEÑAS ZYANYA IVONNE
RODRIGUEZ HERNANDEZ DERIAN DAVID
PROFESOR: VÍCTOR DAMIÁN PINILLA MORÁN
PROYECTO TEMA 2:
PROBLEMA DE INGENIERÍA RESUELTO A TRAVÉS DE
MÉTODOS ABIERTOS Y CERRADOS.
FECHA DE ENTREGA: 12 DE JUNIO 2021
SEMESTRE
2021-2
1
ÍNDICE
INTRODUCCIÓN.............................................................................................................................2
DESARROLLO DEL PROBLEMA.................................................................................................4
PROCEDIMIENTO Y DESARROLLO DE CÁLCULOS.............................................................5
MÉTODO CERRADO.....................................................................................................................5
RESOLUCIÓN DEL PRONLEMA POR MÉTODO DE BISECCIÓN........................................5
MÉTODO ABIERTO.......................................................................................................................8
RESOLUCIÓN DEL PROBLEMA POR MÉTODO DE NEWTON – RAPHSON...................10
CÓDIGOS.......................................................................................................................................12
CONCLUSIONES..........................................................................................................................13
REFERENCIAS.............................................................................................................................14
INTRODUCCIÓN.
MÉTODOS CERRADOS.
Se les llama Métodos cerrados a todos aquellos que requieren de un intervalo de
valores de la variable independiente [a, b] para una función f(x) que posee raíces
reales, tal que f(a) y f(b) son de signos contrarios, por lo que se cumple que
f(a)* f(b) < 0. En general, este cambio de signos en la función valuada en a y b
implica que en este intervalo existe al menos una raíz.
1. Método de bisección
El método de bisección se aplica a funciones algebraicas o trascendentes y
proporciona únicamente raíces reales. Tiene su origen en un popular
algoritmo de búsqueda de datos en arreglos vectoriales denominado
búsqueda binaria. Es un método cerrado, es decir, requiere de un intervalo
en el cual esté atrapada una raíz. Básicamente, consiste en cortar el
intervalo en dos justo por la mitad (bisectar) considerando a este punto
como una aproximación de la raíz de la función. Posteriormente, debe
determinarse si la raíz verdadera se encuentra a la derecha o a la izquierda
de la aproximación y, según corresponda, cerrar el intervalo con la
aproximación y el límite derecho o izquierdo, pero siempre manteniendo a
la raíz verdadera en el intervalo. Esta operación se repite hasta que la
diferencia entre las dos últimas aproximaciones sea menor que una
tolerancia preestablecida. La bisección es un método robusto, aunque
resulta lento en su proceso por lo oneroso de los cálculos que deben
realizarse; por otra parte, su convergencia puede en ocasiones ser
inestable.
2. Método de la Interpolación lineal.
2
Por su construcción geométrica, a este método también se le conoce como
de las cuerdas. Una cuerda, geométricamente hablando, es el segmento de
una recta que une dos puntos de un arco (internamente). Como se verá
más adelante, la geometría que se forma en torno a esta cuerda, asociada
a la búsqueda de la raíz de la función, permite establecer un método
numérico considerado más eficiente que Bisección y que, con las debidas
precauciones, suele ser más rápido.
Se aplica a funciones algebraicas y trascendentes y proporciona sólo raíces
reales.
Este método es de tipo cerrado, es decir, requiere de un intervalo [a, b] que
pertenece al dominio de la función y para el cual f(a) * f(b) < 0, lo que
implica que en el intervalo [a, b] existe al menos una raíz. [1]
MÉTODOS ABIERTOS.
Los métodos numéricos se encargan de obtener respuestas a problemas en
donde la solución analítica es complicada. En este caso, se obtendrán raíces de
ecuaciones algebraicas o trascendentes desde de una aproximación a su raíz,
obtenida a partir de la inspección de su gráfica o de su expresión analítica; a
diferencia de los métodos cerrados que requieren de un intervalo que atrape a
dicha raíz.
1. Método de Aproximaciones sucesivas.
El método de Aproximaciones sucesivas representa la esencia de los
procesos iterativos ya que permite definir una ecuación de recurrencia que,
en apariencia, no tiene sentido desde el punto de vista algebraico, pero que
resulta muy atinada iterativamente hablando ya que toma un valor inicial
que se mejora a través de las iteraciones. Sin embargo, el método como tal
no es ciento por ciento aplicable para cualquier ecuación algebraica o
trascendente, debe vigilarse estrictamente su criterio de convergencia; no
obstante, se utiliza como base para completar otros métodos abiertos.
2. Método de Newton-Raphson.
El método Newton Raphson (N-R) es, junto con la Bisección, uno de los
más utilizados. Su preferencia radica en su robustez y velocidad para
encontrar la raíz cuando la aproximación cumple con su criterio de
convergencia. Se aplica a ecuaciones algebraicas y trascendentes y
proporciona raíces reales y complejas. [2]
DESARROLLO DEL PROBLEMA.
3
Para el siguiente proyecto se resolverá el siguiente problema:
En un parque de diversiones se ensambla una montaña rusa, la cual en un
intervalo de el recorrido tiene un túnel subterráneo (imagen 1.1). Se desea saber
qué longitud tiene del punto de entrada al punto de salida; También se desea
saber cuál es el volumen del orificio que contendrá al túnel. Se sabe que el ancho
de las vías es de 0.12 [dam] la cual tiene que tener una separación de la pared de
0.11 [dam].
Nota: El problema tiene una tolerancia de 0.00001.
El intervalo se da por el siguiente gráfico dado en decámetros:
imagen (1.1)
El cual corresponde a la siguiente función:
f ( x )=xsen ( x ) +2
Dicho problema será desarrollado y resuelto por medio de los métodos cerrados y
abiertos que presentaremos a continuación.
4
PROCEDIMIENTO Y DESARROLLO DE CÁLCULOS.
MÉTODO CERRADO
Método de Bisección.
Definición del método.
A partir de una función algebraica o trascendente y de un intervalo [a, b]
que pertenece al dominio de la función y para el cual f(a) * f(b) < 0, lo que
implica que en el intervalo [a, b] existe al menos una raíz. El método
consiste en bisectar el intervalo [a, b]:
a+b
x 0=
2
Se obtiene una aproximación a la raíz x0; la función se valúa en este nuevo
valor y de acuerdo al signo de la función valuada en este punto, deberá
sustituirse uno de los extremos del intervalo de búsqueda, de tal forma que
se conserve que f(a) * f(b) < 0. De acuerdo a la geometría de la figura, la
sustitución de los intervalos deberá hacerse de la siguiente forma:
Sea a tal que f(a) < 0 y b tal que f(b) > 0:
Si f(x0) < 0, entonces x0 sustituye a a
Si f(x0) > 0, entonces x0 sustituye a b
En cada iteración debería sustituirse alguno de los límites del intervalo que
contiene a la raíz.
Repitiendo este proceso, el intervalo se reduce paulatinamente hasta que
alguna de las aproximaciones coincide razonablemente con la raíz de la
función. El proceso se detiene cuando entre la aproximación x i y la
aproximación anterior x i−1 se satisface un nivel de error (absoluto o relativo)
preestablecido (tolerancia). [3]
RESOLUCIÓN DEL PRONLEMA POR MÉTODO DE BISECCIÓN.
Para saber la longitud que hay del punto de entrada al punto de salida es
necesario obtener las coordenadas de dichos puntos, es decir las raíces, las
cuales las obtuvimos con un algoritmo que se desarrolló en Matlab.
Punto de entrada (e):
Se tomó como rango del -7 al -5 para el método de bisección.
5
e=−5. 939760208129883 con un 1.90734 x 10−6 % de error
Punto de salida (s):
Se tomó como rango del -5 al -3 para el método de bisección.
6
s=−3.710802078247070con un 1.90734 x 10−6 % de error
L= longitud = x 2−x 1
L= s−e
L= (−3.71080207 ) −(−5.93976020 )
L=2.22895813[dam ]
Por lo tanto, la longitud que hay del punto de entrada al punto de salida es de
2.22895812[dam]
El volumen se obtiene con la siguiente fórmula:
V =LWH … … . (1 )
Donde:
L = Longitud
W = Ancho
H = Alto
Sabemos que:
L = 2.22895812 [dam]
W = 0.34 [dam]
H = 4.91318 [dam]
7
Sustituimos valores en ec. (1)
V =¿V =3.72343265[dam3 ]
MÉTODO ABIERTO.
Método Newton - Raphson.
Definición del método.
El nombre original del método N-R es de las tangentes. La tangente es una
recta que intersecta a una curva en un solo punto; en consecuencia, es
perpendicular a su radio.
A partir de la imagen, se plantea que en un valor x 0 que representa una
aproximación a la raíz de la ecuación, se traza una tangente en el punto
f ¿)
Interpretación geométrica del método Newton - Raphson.
Esta recta tangente deberá cortar al eje horizontal y el punto donde esto
ocurra será la nueva aproximación x 1, de tal forma que en el punto f ( x 1) se
trace una nueva tangente. Este proceso se repetirá hasta que el corte de la
tangente en el eje horizontal coincida con la raíz de la ecuación, o bien,
cuando la diferencia entre dos aproximaciones sucesivas cumpla con una
tolerancia preestablecida. De nuevo a partir de la siguiente figura, con base
8
en las dos primeras iteraciones, se define la siguiente relación entre el
triángulo formado por la recta tangente y el ángulo theta:
f ( x 0)
tan ( θ )=
x0 −x1
Por otra parte, se conoce que:
f ´ ( x 0 ) =tan (θ)
Obtención de la ecuación de la recurrencia.
Sustituyendo ecuaciones:
f ( x0 )
f ´ ( x0)=
x 0− x1
En esta ecuación la incógnita es representada por la iteración siguiente x 1.
Despejándola y expresándola en forma iterativa para cualquier iteración:
f ( xi )
x i+1=x i−
f ´ ( xi )
Este último resultado representa la ecuación de recurrencia del método de
N-R. [4]
9
RESOLUCIÓN DEL PROBLEMA POR MÉTODO DE NEWTON – RAPHSON.
Para saber la longitud que hay del punto de entrada al punto de salida es
necesario obtener las coordenadas de dichos puntos, es decir las raíces, las
cuales las obtuvimos con un algoritmo que se desarrolló en Matlab.
Punto de entrada (e):
Se tomó como primer valor x=-6.5 para el método de Newton-Raphson:
e=−5. 939760466135034 con un 5.93976 x 10−8 % de error
Punto de salida (s):
Se tomó como primer valor x=-4 para el método de Newton-Raphson
10
s=−3.7 10802830307117 conun 2.34741 x 10−6 % de error
L= longitud = x 2−x 1
L= s−e
L= (−3.71080283)−(−5.93976046)
L=2.22895763[dam]
Por lo tanto, la longitud que hay del punto de entrada al punto de salida es de
2.22895763[dam]
El volumen se obtiene con la siguiente fórmula:
V =LWH … … ( 1 )
Donde:
L = Longitud
W = Ancho
H = Alto
Sabemos que:
L = 2.22895763 [dam]
W = 0.34 [dam]
H = 4.91318 [dam]
Sustituimos valores en ec. (1)
V =¿V =3.72343182[dam 3 ]
CÓDIGOS.
11
BISECCIÓN:
NEWTON – RAPHSON:
12
CONCLUSIONES.
El método de la Bisección converge muy lentamente lo que da pie a la
propagación del error por la cantidad de iteraciones necesarias para dicha
convergencia.
Es necesario garantizar la continuidad en el intervalo.
En el método de Newton Raphson, en casos donde el punto de inflexión
ocurre en la vecindad de una raíz las iteraciones divergen progresivamente
generando problemas en el uso de este método.
En el método de Newton Raphson, cuando se tiene un máximo o mínimo
local la tendencia del método será oscilar alrededor del máximo o mínimo
local, de enfocándose de la raíz y perdiendo el efecto de búsqueda de
raíces del método.
El método de Newton Raphson se vuelve lento cuando la derivada se va
haciendo pequeña.
Resolver el problema de una multiplicidad impar no es difícil, podría
hacerse por bisección, por ejemplo, mientras que el problema de una
multiplicidad no par.
Para el método de búsquedas incrementales es importante elegir
adecuadamente el incremento, pues de este depende que se halle
correctamente las raíces de una función dentro de un intervalo.
En ocasiones el método de regla falsa no es el mejor de los métodos
cerrados, pues no converge rápidamente, para estos casos es más
eficiente el método de bisección.
El método de punto fijo busca hallar las raíces en funciones de la forma
f(x)=0, a través de aproximaciones sucesivas que convergen a la solución
de la ecuación
13
Los métodos cerrados requieren de un intervalo que contenga la raíz o se
encuentre en sus proximidades.
Los métodos abiertos se basan en fórmulas que requieren solo de un valor
inicial que no necesariamente debe incluir una raíz.
Los métodos abiertos generalmente convergen más rápidamente que los
métodos cerrados.
REFERENCIAS.
[1]
Borras, H., Duran, R., y Iriarte, R. (1984). Apuntes de métodos numéricos (F.
de Ingeniería UNAM, Ed.).
[2] Burden, R., y Faires, D. (2011). Análisis numérico (C. Learning, Ed.). Chapra, S., y
Canale, R. (2015).
[3] Métodos numéricos para ingenieros (M. Hill, Ed.). García B., S. (2017).
[4] Análisis numérico (Alfaomega, Ed.). Luthe, R., Olivera, A., y Schutz, F. (1985).
Métodos numéricos. Olivera Salazar, A. (s.f.). Métodos numéricos (Limusa, Ed.).
Sandoval, H. (2017). Métodos numéricos.
14
