Taller 3 Solución de Ecuaciones No-Lineales-2023 - 014559

Unidad 2: Solución de Ecuaciones No-Lineales
UNIVERSIDAD DEL MAGDALENA

FACULTAD DE INGENIERÍA
ANÁLISIS NUMÉRICO
TALLER 3: SOLUCIÓN DE ECUACIONES NO-LINEALES
Mg. Alvaro Espinosa Pérez
Valor: 40 Puntos
La solución de estos ejercicios debe incluir su gráfica, así como la solución realizada con las
funciones de OCTAVE en la consola. Se debe especificar el método utilizado y definir el criterio de
pare. Debe mostrar como mínimo dos iteraciones al usar los métodos para verificar la aplicación
de cada una de las fórmulas de iteración.
1. La función 2𝑥(1 − 𝑥 2 + 𝑥)𝑙𝑛 𝑥 = 𝑥 2 − 1 tiene una única raíz real en el intervalo [0, 1].
Utilice los métodos de bisección, de la regula falsi, de la secante hasta obtener la convergencia
con cuatro cifras decimales.
1
2. Investigue la naturaleza de la iteración de punto fijo cuando 𝑔(𝑥) = −4 + 4𝑥 − 𝑥 2
2
a. Resuelva 𝑔(𝑥) = 𝑥 y prueba que los puntos fijos son 2 y 4.
b. Tome como valor inicial 𝑥0 = 1.9 y calcule 𝑥1 , 𝑥2 𝑦 𝑥3 .
c. Tome como valor inicial 𝑥0 = 3.7 y calcule 𝑥1 , 𝑥2 𝑦 𝑥3 .
d. Halle los errores absolutos 𝐸𝑘 y relativos 𝐸𝑟𝑘 para cada 𝑥𝑘 obtenidos en b) y c).
e. ¿Qué conclusiones puedes obtener?
3. Aplique el método de Newton y analice el tipo de divergencia de las siguientes funciones:
a. 𝑓(𝑥) = 𝑥𝑒 −𝑥 .
b. 𝑓(𝑥) = 𝑥 3 − 𝑥 − 3.
c. 𝑓(𝑥) = arctan(𝑥).
4. Considere la función 𝑓 (𝑥) = 2 + 𝑠𝑒𝑛𝑥 − 𝑥. Utilice el método de Newton con 𝑥0 = 2 para
aproximar un cero de 𝑓 . Genera términos hasta |𝑥𝑛+1 − 𝑥𝑛 | < 0,0002.
a. Construye una Secuencia de Aitken {𝑥𝑛 }∞ 𝑛=0 para aproximar la solución. ¿Ha mejorado
la convergencia?
b. Utilice el método de Steffensen para aproximar la solución con un error de 10−5 .
5. Las raíces del polinomio 𝑃(𝑥) = 𝑥 3 − 3.9𝑥 2 + 4.79𝑥 − 1.881, son 𝑥1 = 0.9, 𝑥2 =
1.1 𝑦 𝑥3 = 1.9. Aun cuando 𝑥0 = 1 está entre 𝑥1 𝑦 𝑥2 , y cerca de ambos, el método de
Newton-Raphson no conduce a ninguna de estas raíces.
a. ¿Por qué ocurre esto?
b. ¿Si utilizamos el método de Horner se resolvería el problema? ¿Por qué?
c. ¿Para qué valor de 𝑥0 el método de Newton-Raphson aproxima a 𝑥2 ?
6. La función descrita por 𝑓(𝑥) = 𝑙𝑛(𝑥 2 + 1) − 𝑒 0.4𝑥 𝑐𝑜𝑠 𝜋𝑥 tiene un número infinito de ceros.
a. Determine, dentro de 10−6 , el único cero negativo.
b. Determine una aproximación inicial razonable para determinar, dentro de 10−6 , el
vigésimo quinto cero positivo más pequeño de 𝑓.
7. Implemente el método de Newton con los puntos iniciales adecuados para encontrar una
solución al siguiente sistema:
𝑓(𝑥, 𝑦) = 𝑥 4 − 4𝑦 2 − 4
𝑔(𝑥, 𝑦) = (𝑥 − 1)2 + 𝑦 2 − 4
8. APLICACIONES:
a. Un objeto que cae verticalmente a través del aire está sujeto a una resistencia viscosa
además de la fuerza de gravedad. Supongamos que se deja caer un objeto con masa 𝑚
desde una altura 𝑠0 y que la altura del objeto después de t segundos es:
__________________________________________________________________________________________________________
Análisis Numérico-2023 1
Unidad 2: Solución de Ecuaciones No-Lineales
𝑚𝑔 𝑚2 𝑔
𝑠(𝑡) = 𝑠0 − 𝑡 + 2 (1 − 𝑒 −𝑘𝑡/𝑚 ),
𝑘 𝑘
donde 𝑔 = 32.17 𝑝𝑖𝑒𝑠/𝑠2 y 𝑘 representa el coeficiente de resistencia del aire en lb-
s/pie. Supongamos que 𝑠0 = 300 𝑝𝑖𝑒𝑠, 𝑚 = 0,25 𝑙𝑖𝑏𝑟𝑎𝑠 y 𝑘 = 0,1 libras-s/pie.
Encuentre, con una precisión de 0,01 s, el tiempo que le toma al objeto chocar con el
suelo.
b. Considere el problema físico de hallar la porción de una esfera de radio 𝑟 que queda
sumergida al meter la esfera en el agua. Suponga que la esfera está construida de un
material que tiene una densidad de 𝜌 = 0.638 𝑔/𝑐𝑚3 y que su radio mide 𝑟 = 10 𝑐𝑚.
¿Cuánto vale la profundidad ℎ a la que está sumergido el polo sur de la esfera? (Sug: Use
el principio de Arquímedes).
c. Para medir la altura ℎ de un pozo se recurre al método de lanzar una piedra y medir el
tiempo 𝑡 hasta que se escuche el golpe contra el agua. En este proceso se observan dos
fases. En la primera se tiene en cuenta el tiempo t1 de caída libre del cuerpo que se rige
1
por la ecuación ℎ = 𝑔𝑡1 2 , con g = 9,81m/s2 la aceleración de la gravedad y en segundo,
2
el tiempo t2 que tarda el sonido en recorrer la altura ℎ = 𝑐𝑡2 con c = 343 m/s la velocidad
del sonido. Es evidente entonces que 𝑡 = 𝑡1 + 𝑡2 . De ahí que:
2ℎ ℎ
𝑡=√ +
𝑔 𝑐
A. A partir de la ecuación anterior deduzca que:
𝑐2
ℎ2 − 2 (𝑡𝑐 + ) ℎ + 𝑐 2𝑡 2 = 0
𝑔
B. Construye una función ℎ(𝑡) de la ecuación anterior y dibuja una gráfica que muestre
la altura ℎ en función del tiempo 𝑡 medido en segundos en [0, 10] y calcula una raíz.
9. PROBLEMAS DE CÓMPUTO
A continuación, realizaran una implementación sencilla de un algoritmo en OCTAVE que
permita obtener un script en el que aparecen los comandos necesarios para aproximar la raíz
de la ecuación 𝑦 = 𝑥 2 − 1 + 𝑒 𝑥 ; primero mostrando la gráfica, y segundo dando a conocer
las iteraciones incluyendo su criterio de pare usando una tolerancia de 0,001 para los
siguientes métodos:
a. Método de la Bisección.
b. Método de Newton.
c. Método de iteración de Punto Fijo.
Ejemplo: Podemos implementar el algoritmo de la Bisección en OCTAVE de la siguiente manera:
====================================
====== 𝑀𝐸𝑇𝑂𝐷𝑂 𝐷𝐸 𝐿𝐴 𝐵𝐼𝑆𝐸𝐶𝐶𝐼𝑂𝑁 =======
====================================
𝐷𝑖𝑔𝑖𝑡𝑒 𝑙𝑎 𝑓𝑢𝑛𝑐𝑖𝑜𝑛: 𝑥. ^3 + 4.∗ 𝑥. ^2 − 10
𝐷𝑖𝑔𝑖𝑡𝑒 𝑒𝑙 𝑙𝑖𝑚𝑖𝑡𝑒 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟 𝑑𝑒𝑙 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜: 1
𝐷𝑖𝑔𝑖𝑡𝑒 𝑒𝑙 𝑙𝑖𝑚𝑖𝑡𝑒 𝑠𝑢𝑝𝑒𝑟𝑖𝑜𝑟 𝑑𝑒𝑙 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜: 2
𝐷𝑖𝑔𝑖𝑡𝑒 𝑙𝑎 𝑡𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎: 0.001
𝐸𝑙 𝑛𝑢𝑚𝑒𝑟𝑜 𝑚í𝑛𝑖𝑚𝑜 𝑑𝑒 𝑖𝑡𝑒𝑟𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝑒𝑠: 10
𝑖 𝑎 𝑏 𝑐 𝑓(𝑐) 𝐸𝑎(%)
1 1.000000 2.000000 1.500000 2.375000
2 1.000000 1.500000 1.250000 − 1.796875 20.00000
3 1.250000 1.500000 1.375000 0.162109 9.09091
4 1.250000 1.375000 1.312500 − 0.848389 4.76190
5 1.312500 1.375000 1.343750 − 0.350983 2.32558
6 1.343750 1.375000 1.359375 − 0.096409 1.14943
7 1.359375 1.375000 1.367188 0.032356 0.57143
8 1.359375 1.367188 1.363281 − 0.032150 0.28653
9 1.363281 1.367188 1.365234 0.000072 0.14306
10 1.363281 1.365234 1.364258 − 0.016047 0.07158
𝐿𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑒𝑠: 1.364258.
__________________________________________________________________________________________________________
Análisis Numérico-2023 2

Taller 3 Solución de Ecuaciones No-Lineales-2023 - 014559

Cargado por

Copyright:

Formatos disponibles

Taller 3 Solución de Ecuaciones No-Lineales-2023 - 014559

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Taller 3 Solución de Ecuaciones No-Lineales-2023 - 014559

Cargado por

Copyright:

Formatos disponibles

Unidad 2: Solución de Ecuaciones No-Lineales

UNIVERSIDAD DEL MAGDALENA

𝐿𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑒𝑠: 1.364258.

También podría gustarte