Prácticas de Matemáticas con Mathematica . Fundamentos de Matemáticas III . Grado en Ingeniería Civil. Práctica nº 2 (continuación). Métodos de resolución aproximada de ecuaciones. Departamento de Matemática Aplicada. E.P.S. de Zamora Universidad de Salamanca Ejemplo 6: La función f H xL = x 5 3 - 2 x 3 Sin[x] -20Sin 2 [x] tiene una raiz en el intervalo (2,6). Deter- mine el intervalo donde hay que elegir el punto de arranque x 0 en el método de Newton para obtener conver- gencia a esa raiz. f@x_D := Hx^ H1 3LL ^5 - 2x^3Sin@xD - 20 Sin@xD ^2 El Mathematica no trabaja con el valor real para las raices de índice impar y radicando negativo, así que definiremos nuestra propia función para hallar la raiz cúbica. Plot@x^ H1 3L, 8x, - 4, 4<DH* obsérvese que no dibuja la parte de la gráfica correspondiente a valores negativos de la variable independiente*L -4 -2 2 4 0.5 1.0 1.5
Se utiliza el programa Mathematica para determinar parejas de puntos que conducen a un bucle en la aplicación del método de Newton para resolver ecuaciones de forma aproximada.
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Prácticas de Matemáticas con Mathematica .
Fundamentos de Matemáticas III . Grado en Ingeniería Civil.
Práctica nº 2 (continuación). Métodos de resolución aproximada de ecuaciones.
Departamento de Matemática Aplicada.E.P.S. de Zamora
Universidad de Salamanca
Ejemplo 6: La función f HxL = x53
- 2 x3 Sin[x] -20Sin 2[x] tiene una raiz en el intervalo (2,6). Deter-
mine el intervalo donde hay que elegir el punto de arranque x0 en el método de Newton para obtener conver-gencia a esa raiz.
El Mathematica no trabaja con el valor real para las raices de índice impar y radicando negativo, así quedefiniremos nuestra propia función para hallar la raiz cúbica.
y ahora sí que se dibuja la gráfica completa*LPlot@raiz3@xD, 8x, -4, 4<D
-4 -2 2 4
-1.5
-1.0
-0.5
0.5
1.0
1.5
Utilizando esta raiz cúbica la función f(x), que denotaremos ahora por fm(x) se escribe:
fm@x_D := raiz3@xD^5 - 2 x^3 Sin@xD - 20 Sin@xD^2
Plot@fm@xD, 8x, -5, 5.2<D
-4 -2 2 4
50
100
150
200
250
Para hallar los valores a,b que den lugar a un posible bucle en el método de Newton hemos de resolver elsistema siguiente, tomando como puntos de partida valores en el intervalo (2,6):
b ® 5.0227263182941152467076575641193896601247901061503<La función de iteración del método de Newton se escribe:
2 3 FM III (26-9-2012) Bucles en el método de Newton.nb
g@x_D = x - f@xD � f'@xD
x - Ix5�3- 2 x3 Sin@xD - 20 Sin@xD2M �
5 x2�3
3- 2 x3 Cos@xD - 6 x2 Sin@xD - 40 Cos@xD Sin@xD
Pero como nos vuelven a aparecer raices cúbicas, utilizamos la funcion raiz cúbica definida por nosotros, yqueda la función de iteración modificada gm(x) siguiente:
Cuando iteremos y la función de iteración se aplique a valores positivos podremos utilizar la función g(x) ogm(x) pues es lo mismo, pero cuando aparezcan valores negativos deberemos usar la gm(x):
Si partimos del valor inicial x0 = a =2.1654089718798745299479799727 o del valorx0 = b =5.0227263182941152467076575641 obtendremos un bucle:
5.0227263182941152467076575641193896633476051304794<Para valores de arranque dentro del intervalo (a,b) el procedimiento será convergente a la raiz en el intervalo
dado:
4 3 FM III (26-9-2012) Bucles en el método de Newton.nb
3.0317352187764541326941980514675378799438476562500<Que no haya convergencia a la raiz buscada no significa que no pueda haber convergencia hacia otra raiz. El
valor x0 = 5.1 no está en el intervalo (a,b), y sin embargo sí que da lugar a un procedimiento convergente hacia laraiz 0.0001250000013671876699337948135593
3 FM III (26-9-2012) Bucles en el método de Newton.nb 5
0.00012500000136718766993379481355930238350993022322655<Una gráfica detallada en un intervalo que contenga la raiz anterior nos muestra que en efecto hay una raiz en
ese punto, la cual en una gráfica con una escala mayor resulta inapreciable.
6 3 FM III (26-9-2012) Bucles en el método de Newton.nb
La función dada tiene infinitas raíces debido a su carácter oscilatorio, y para hallar cualquiera de sus raíces seprocedería de manera parecida a como se ha hecho anteriormente.
8 3 FM III (26-9-2012) Bucles en el método de Newton.nb
La función dada tiene infinitas raíces debido a su carácter oscilatorio, y para hallar cualquiera de sus raíces seprocedería de manera parecida a como se ha hecho anteriormente.
Plot@fm@xD, 8x, -20, 20<D
-20 -10 10 20
-5000
5000
10 000
Ejemplo 7: La función f HxL = Exp@xD - 4 + 10 Cos@xD tiene tres raíces en el intervalo (-3,3). Determine sihay bucles en el método de Newton para hallar dichas raíces.
Vamos a intentar localizar la existencia de bucles en el cálculo de las tres raíces más cercanas al origen de lafunción f[x_] = Exp[x] - 4 + 10 Cos[x] mediante el método de Newton
f@x_D := Exp@xD - 4 + 10 Cos@xD;
Plot@f@xD, 8x, -3, 3<DH* en esta gráfica observamos
la situación de cada una de las tres raíces *L
-3 -2 -1 1 2 3
-10
-5
5
g@x_D = x - f@xD � f'@xD
x --4 + ãx + 10 Cos@xD
ãx - 10 Sin@xDPara la raiz más a la izquierda encontramos que hay un bucle, de manera que cuando aplicamos el método deNewton tomando como punto de arranque uno de los puntos del bucle no encontramos la raiz
En el caso de la segunda raiz, también es posible localizar un bucle, de manera que tomando como punto dearranque uno de los puntos del bucle, no encontramos la raiz buscada:
En el caso de la raiz más a la derecha no es posible encontrar un bucle, pues en un entorno de esta raiz lafunción es cóncava hacia arriba. El intento de localizar un bucle conduce a obtener la raiz buscada.
10 3 FM III (26-9-2012) Bucles en el método de Newton.nb
A continuación aparece un dibujo donde se muestra la dinámica del método de Newton. Si r es una raiz def(x) entonces se tiene que g(r)=r-f(r)/f’(r)=r, es decir, las raíces de f(x) son puntos fijos de g(x), que son lospuntos donde g(x) corta a la recta diagonal y=x. Los puntos en rojo y verde señalan los extremos de los buclespara cada una de las raíces.
3 FM III (26-9-2012) Bucles en el método de Newton.nb 11
En los puntos donde se anula el denominador de g(x) habrá asíntotas verticales.A continuación aparece undibujo donde se muestra el denominador de g(x) indicando que hay infinitas asíntotas verticales hacia laizquierda.
12 3 FM III (26-9-2012) Bucles en el método de Newton.nb
Plot@ãx
- 10 Sin@xD, 8x, -20, 4<D
-20 -15 -10 -5
-10
10
20
30
3 FM III (26-9-2012) Bucles en el método de Newton.nb 13