Taller 3
Taller 3
función g(x), si iniciamos desde el valor x = 5 y evaluamos la función, será tomar la primera flecha vertical
que parte desde (5,0) a un punto de la gráfica de g, después obtenemos la flecha horizontal que parte del
punto anterior a un punto de la recta y = x, esto es tomar el valor 0,0067, guardarlo y después evaluarlo
1
en la función g(x), obteniendo la segunda flecha vertical. Este proceso se continua, hasta iterar 9 veces.
Observe que las flechas cada vez son más cerca al punto de interseccion de la recta y la gráfica de g(x).
Definición 1. El número real r es un punto fijo de la función g(x) si g(r) = r.
Por ejemplo la función anterior tiene un punto fijo aproximado a 0,5602. De la anterior definición se
concluye que debemos resolver la ecuación g(x) = x o lo que es mismo la ecuación f (x) = g(x) − x = 0.
Ejemplo 2. Los siguientes ejemplos consisten en tomar la ecuación g(x) = x y solucionarla de manera
algebraica, se deja de ejercicio resolver estas ecuaciones.
1. Para g(x) = x3 , la función tiene tres puntos fijo: x = 0, x = 1 y x = −1, los cuales son soluciones de
la ecuación x3 − x = 0 o x3 = x.
2
2. La función g(x) = xx−1+1
tiene un único punto fijo x = −1.
El objetivo para solucionar una ecuación f (x) = 0 es llevarla a una ecuación de la forma g(x) = x,
donde la función auxiliar g(x) es la que debemos de calcular un punto fijo. Esto es, el punto fijo de la
función g(x) es la raı́z de la ecuación f (x) = 0. Ahora como aseguramos que la función auxiliar tiene un
único punto fijo en algún intervalo [a, b], para garantizar este hecho debemos de suponer que la función
g(x) no varı́a demasiado. Este resultado se sigue del siguiente teorema.
Teorema 1. Sea g(x) una función continua en [a, b] y diferenciable en (a, b), donde |g 0 (x)| < 1, para todo
x ∈ (a, b), entonces la función g(x) tiene un único punto fijo r en [a, b].
Para obtener un punto fijo de la función auxiliar g(x), utilizamos el siguiente algoritmo de punto fijo
en Matlab véase figura 3.
Ejemplo 3. Queremos resolver la ecuación f (x) = x2 − x − 1 = 0, para esto buscaremos varias funciones
auxiliares g(x) que permitan calcular una de las raı́ces de la ecuación cuadrática.
√
1. De la ecuación despejemos x2 , esto es x2 = x + 1, De aquı́ x = x + 1, obteniendo la primer
función auxiliar g(x). Se utiliza el programa con los datos colocados en la ventana comando windows
de Matlab, observe la manera de ejecutar el programa como la función r = f pi(g, x0, k) (véase la
gráfica derecha de la figura 3) obteniendo el valor aproximado de r = 1,6180, el cual es un valor
aproximado de la ecuación f (x) = 0. Se deja de ejercicio verificar que |g 0 (x)| ≤ 1 para todo x > −3
4
.
Un gráfico detallado de la iteración empezando del valor x0 = 0, puede verse a continuación (véase
el gráfico de la izquierda)
2
obtenemos el mismo valor r = 1,6180. También observe que |g 0 (x)| < 1. Esta convergencia puede
verse a partir del gráfico de la mitad de la figura anterior.
3. ¿Que sucede si se trabaja con la función auxiliar g(x) = x2 − 1 tomando el valor inicial de x0 = 1?
Observe del siguiente gráfico que nunca converge la iteración de punto fijo, modificando un poco el
programa de la figura 1 se puede concluir que esta convergencia es repetitiva entre los valores −1 y
0. Además, se puede concluir que |g 0 (x)| < 1 para − 12 < x < 12 .
En consecuencia, si |g 0 (x)| < 1, entonces los errores disminuyen en cada iteración. Si |g 0 (x)| > 1,
entonces los errores crecen. Este tipo de problemas dan origen a la siguiente definición.
Definición 2. Un punto fijo r de una función g(x) es estable (o atractor) si |g 0 (x)| < 1 o inestable (o
repulsor) si |g 0 (x)| > 1.
La figura 4 presenta la clasificación de los puntos fijos, el primero es inestable (gráfica izquierda)
y el segundo es estable (gráfica derecha), observe que no importa el valor inicial que se tome x0 , el
comportamiento de la derivada es de suma importancia para que pueda atraerse la iteración de punto fijo.
3
2. Ejercicios 3 2 3
1. Muestre que 1, 2 y 3 son puntos fijos de la siguiente g(x). (a) x6x−10
+x−6
(b) 6+6x11−x
2. Use el hecho de que |g 0 (x)| ≤ 1 para determinar si la iteración de punto fijo de g(x) converge
localmente al punto fijo dado r. (a) g(x) = 2x−1
x2
, r = 1 (b) g(x) = e2x − 1, r = 0 (c) g(x) = senx − x,
r = 0.
3. Encuentre cada punto fijo y decida si la iteración de punto fijo converge localmente a éste. (a)
g(x) = x2 − 32 x + 32 (b) g(x) = x2 − 14 x + 38 .
4. Exprese cada ecuación como un problema de punto fijo x = g(x) en tres formas diferentes. (a)
x3 − x + ex = 0 (b) 3x−2 + 9x3 = x2 .
q
15 2 5 3
5. Considere la iteración de punto fijo aplicado a g(x) = 1 − 5x + 2 x − 2 x (a) Muestre que 1 − 35
q
y 1 + 35 son puntos fijos de g(x). (b) Muestre que ninguno de los tres puntos fijos es localmente
convergente, esto es |g 0 (r)| ≥ 1
6. Utilice la iteración de punto fijo para las siguientes funciones con los valores iniciales indicados
√
(a)g(x) = 2x3 − 11,7x2 + 17,7x − 5 x0 = 3 (b) g(x) = sen( x) − x, x0 = 0,5
7. Aplique la iteración de punto fijo para encontrar la solución de cada ecuación hasta ocho posiciones
decimales correctas, (a)x5 + x = 1 (b) senx = 6x + 5 (c) ln|x| + x2 = 3.
8. Calcule las raı́ces cuadradas de los números siguientes hasta ocho cifras decimales correctas usando la
iteración de punto fijo: (a) 3 (b) 5. Establezca su estimación inicial y el número de pasos necesarios.
9. Calcule las raı́ces cúbicas de los números siguientes hasta ocho posiciones decimales correctas usando
2x+ A2
la iteración de punto fijo con g(x) = 3
x
donde A es (a) 2 (b) 3 (c) 5. Establezca su estimación
inicial y el número de pasos necesarios.
10. Deduzca tres diferentes g(x) para encontrar las raı́ces hasta seis cifras decimales correctas de las
siguientes f (x) = 0 mediante la iteración de punto fijo. Ejecute la iteración de punto fijo (IPF) para
cada g(x), reporte los resultados y si hay convergencia o divergencia. Cada ecuación f (x) = 0 tiene
tres raı́ces. Si es necesario deduzca más g(x) hasta encontrar todas las raı́ces mediante la IPF. Para
cada serie convergente. determine el valor de S a partir de los errores ei+1ei
. (a) f (x) = 2x3 − 6x − 1
(b) f (x) = ex−2 + x3 − x (c) f (x) = 1 + 5x − 6x3 − e2x
11. Encuentre el error hacia adelante y hacia atrás para las siguientes funciones, donde la raı́z es 3/4
y la raı́z aproximada es xa = 0,74: (a)f (x) = 4x − 3, (b)f (x) = (4x − 3)2 , (c)f (x) = (4x − 3)3 ,
1
(d)f (x) = (4x − 3) 3 .
12. Encuentre el error hacia adelante y hacia atrás para las siguientes funciones, donde la raı́z es 1/3
y la raı́z aproximada es xa = 0,3333: (a)f (x) = 3x − 1, (b)f (x) = (3x − 1)2 , (c)f (x) = (3x − 1)3 ,
1
(d)f (x) = (3x − 1) 3 .
13. (a) Encuentre la multiplicidad de la raı́z r = 0 de f (x) = 1 − cosx. (b) Encuentre los errores hacia
adelante y hacia atrás de la raı́z aproximada xa = 0,0001.
4
14. (a) Encuentre la multiplicidad de la raı́z r = 0 de f (x) = x2 sen2 x. (b) Encuentre los errores hacia
adelante y hacia atrás de la raı́z aproximada xa = 0,01.