[ Punto fijo ] [ Ejercicio ] [ Analítico ] [ Algoritmo ] [ función ]
..
1. Ejercicio
Referencia: Burden 2.2 p41, Chapra 6.1 p143, Rodríguez 3.2 p44
Encontrar la solución a la ecuación, usando el método del punto fijo, considerando tolerancia de 0.001
f(x): e^{-x} - x = 0[ Punto fijo ] [ Ejercicio ] [ Analítico ] [ Algoritmo ] [ función ]
2. Desarrollo Analítico
Al igual que los métodos anteriores, es conveniente determinar el intervalo [a,b] donde es posible evaluar f(x). Se revisa si hay cambio de signo en el intervalo para buscar una raíz. Para el ejemplo, el rango de observación será [0,1], pues f(0)=1 es positivo y f(1)=-0.63 es negativo.
Para el punto fijo, se reordena la ecuación para para tener una ecuación con la variable independiente separada. Se obtiene por un lado la recta identidad y=x, por otro se tiene la función g(x).
f(x): e^{-x} - x = 0 x = e^{-x}Se buscará la intersección entre las dos expresiones .
Se puede iniciar la búsqueda por uno de los extremos del rango [a,b].
iteración 1
– iniciando desde el extremo izquierdo x=a,
x_0 = 0– se determina el valor de b = g(x) ,
x_1 = g(0) = e^{-0} = 1– se determina la diferencia de la aproximación o error = |b-a|
tramo =|1-0|= 1– se proyecta en la recta identidad, como el nuevo punto de evaluación
x=b.
– se repite el proceso para el nuevo valor de x, hasta que el error sea menor al tolerado. En caso que el proceso no converge, se utiliza un contador de iteraciones máximo, para evitar tener un lazo infinito.
iteración 2
x_1 = 1 x_2 = g(1) = e^{-1} = 0.3678 tramo =|0.3678 - 1|= 0.6322iteración 3
x_2 = 0.3678 x_3 = e^{-0.3678} = 0.6922 tramo =|0.6922-0.3678|= 0.3244La tabla resume los valores de las iteraciones
iteración | x | xnuevo = g(x) | |error| |
1 | 0 | 1.0 | 1 |
2 | 1.0 | 0.3678 | 0.6322 |
3 | 0.3678 | 0.6922 | 0.3244 |
4 | 0.6922 | … | |
5 |
El proceso realizado en la tabla se muestra en la gráfica, para una función que converge.
[ Punto fijo ] [ Ejercicio ] [ Analítico ] [ Algoritmo ] [ función ]
..
3. Algoritmo en Python
El algoritmo en Python del video se adjunta para seguimiento. Se requiere la función gx, el punto inicial y la tolerancia, la variable de número de iteraciones máxima, iteramax
, permite controlar la convergencia de la función con el método.
El resultado del algoritmo se presenta como:
raíz en: 0.5669089119214953 errado : 0.0006477254067881466 itera : 14
Instrucciones en Python
# Algoritmo de punto fijo # [a,b] son seleccionados desde la gráfica # error = tolera import numpy as np # INGRESO fx = lambda x: np.exp(-x) - x gx = lambda x: np.exp(-x) a = 0 b = 1 tolera = 0.001 iteramax = 15 # PROCEDIMIENTO i = 1 # iteración b = gx(a) tramo = abs(b-a) while not(tramo<=tolera or i>iteramax): a = b b = gx(a) tramo = abs(b-a) i = i + 1 respuesta = b # Valida convergencia if (i>=iteramax): respuesta = np.nan # SALIDA print('raíz en: ', respuesta) print('errado : ', tramo) print('itera : ', i)
Para obtener la gráfica básica se determinan los puntos para cada función fx y gx. Como los valore de a y b cambiaron durante el algoritmo, para la siguiente sección es necesario volver a escribirlos u obtener una copia de los valores en otra variable para reutilizarlos en ésta sección.
Se añaden las siguiente instrucciones al algoritmo anterior:
# GRAFICA import matplotlib.pyplot as plt # calcula los puntos para fx y gx a = 0 b = 1 muestras = 21 xi = np.linspace(a,b,muestras) fi = fx(xi) gi = gx(xi) yi = xi plt.plot(xi,fi, label='f(x)', linestyle='dashed') plt.plot(xi,gi, label='g(x)') plt.plot(xi,yi, label='y=x') plt.axvline(respuesta, color='magenta', linestyle='dotted') plt.axhline(0) plt.xlabel('x') plt.ylabel('f(x)') plt.title('Punto Fijo') plt.grid() plt.legend() plt.show()
[ Punto fijo ] [ Ejercicio ] [ Analítico ] [ Algoritmo ] [ función ]
..
4. Función en Python
El resultado del algoritmo como una función y mostrando la tabla de resultados por iteraciones es:
método del Punto Fijo i ['xi', 'gi', 'tramo'] 0 [0. 1. 1.] 1 [1. 0.367879 0.632121] 2 [0.367879 0.692201 0.324321] 3 [0.692201 0.500474 0.191727] 4 [0.500474 0.606244 0.10577 ] 5 [0.606244 0.545396 0.060848] 6 [0.545396 0.579612 0.034217] 7 [0.579612 0.560115 0.019497] 8 [0.560115 0.571143 0.011028] 9 [0.571143 0.564879 0.006264] 10 [0.564879 0.568429 0.003549] 11 [0.568429 0.566415 0.002014] 12 [0.566415 0.567557 0.001142] 13 [0.567557 0.566909 0.000648] raíz en: 0.5669089119214953
Instrucciones en Python
# Algoritmo de punto fijo # [a,b] son seleccionados desde # la gráfica de la función # error = tolera import numpy as np def puntofijo(gx,a,tolera, iteramax=20,vertabla=True, precision=6): """ g(x) se obtiene al despejar una x de f(x) máximo de iteraciones predeterminado: iteramax si no converge hasta iteramax iteraciones la respuesta es NaN (Not a Number) """ itera = 0 b = gx(a) tramo = abs(b-a) if vertabla==True: print('método del Punto Fijo') print('i', ['xi','gi','tramo']) np.set_printoptions(precision) print(itera,np.array([a,b,tramo])) while(tramo>=tolera and itera<=iteramax): a = b b = gx(a) tramo = abs(b-a) itera = itera + 1 if vertabla==True: print(itera,np.array([a,b,tramo])) respuesta = b # Valida respuesta if itera>=iteramax: respuesta = np.nan print('itera: ',itera, 'No converge,se alcanzó el máximo de iteraciones') return(respuesta) # PROGRAMA ---------------------- # INGRESO fx = lambda x: np.exp(-x) - x gx = lambda x: np.exp(-x) a = 0 b = 1 tolera = 0.001 iteramax = 15 # PROCEDIMIENTO respuesta = puntofijo(gx,a,tolera,iteramax,vertabla=True) # SALIDA print('raíz en: ', respuesta)
De añadirse la parte gráfica, el bloque no requiere volver a escribir los valores de a y b, al ser usados como variable internas a la función y se mantienen en la parte principal del programa.
Scipy.optimize.fixed_point
El método del punto fijo se encuentra implementado en Scipy, que también puede ser usado de la forma:
>>> import scipy.optimize as opt >>> opt.fixed_point(gx,a,xtol=0.001,maxiter=15) array(0.5671432948307147)
el valor predeterminado de iteraciones si no se escribe es 500 y la tolerancia predeterminada es xtol=1e-08
https://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.fixed_point.html
Tarea
- Revisar lo que sucede cuando el valor inicial a esta a la derecha de la raíz.
- Validar que la función converge, revisando que |g'(x)|<1, o que la función g(x) tenga pendiente menor a pendiente de la recta identidad.
[ Punto fijo ] [ Ejercicio ] [ Analítico ] [ Algoritmo ] [ función ]