dinámica de la familia de métodos iterativos de tchebychev · pdf filedin amica...
TRANSCRIPT
Dinamica de la familia de metodos iterativos deTchebychev-Halley
Ddays,Universitat Jaume I24-26 Octubre 2012
Alicia Cordero, Juan R. Torregrosa, Pura Vindel
30 de octubre de 2012
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 1 / 26
Indice
1 La familia Tchebychev-HalleyOperadorEstudio de los puntos fijosPuntos crıticos
2 Espacio de parametros
3 La cabeza del gatoLa frontera
4 El cuerpo del gatoLa frontera
5 Dentro del collar
6 En la frontera del collar
7 Fuera del conjunto gato
8 Conclusiones
9 Bibliografa
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 2 / 26
La familia Tchebychev-Halley
El operador correspondiente a la familia de metodos de Tchebychev-Halleyes:
G (z) = z −(
1 +1
2
Lf (z)
1− αLf (z)
)f (z)
f ′ (z)
donde:
Lf (z) =f (z) f ′′ (z)
(f ′ (z))2
α = 0, Tchebychev
α = 1/2, Halley
α = 1, Super-Halley
α = ±∞, Newton.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 3 / 26
La familia Tchebychev-Halley
Si p(z) = z2 + c ,
Gp (z) =z4 (−3 + 2α) + 6cz2 + c2 (1− 2α)
4z (z2 (−2 + α) + αc)
c puede eliminarse mediante la aplicacion:
h (z) =z − i
√c
z + i√
c
cuyas propiedades son h (∞) = 1, h(i√
c)
= 0, h(−i√
c)
=∞.
Op (z) =(h ◦ Gp ◦ h−1
)(z) = z3 z − 2 (α− 1)
1− 2 (α− 1) z
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 4 / 26
Puntos fijos
Op (z) = z ⇒ z = 0, z = 1, z =−3 + 2α±
√5− 12α + 4α2
2= s1, s2
Op(∞) =∞, luego ∞ tambien es un punto fijo, superatractor.
Ademas: s1 = 1s2
, s1 = s2 = ±1 ⇒ α = 12 , α = 5
2Por tanto:
α = 12 ⇒ z = −3+2α±
√5−12α+4α2
2 = −1 con multiplicidad 2 ⇒Puntos fijos: 0,∞, 1,−1
α = 52 ⇒ z = −3+2α±
√5−12α+4α2
2 = 1, entonces z = 1 conmultiplicidad 3 ⇒ Puntos fijos: 0,∞, 1α 6= 1
2 ,52 ⇒ Puntos fijos: 0,∞, 1, s1, s2
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 5 / 26
Estabilidad de los puntos fijos
O ′p (z) = 2z2 3 (1− α) + 2z(3− 4α + 2α2
)+ 3z2 (1− α)
(1− 2 (α− 1) z)2
Luego 0 e ∞ son siempre superatractores.
Lema
Estabilidad de z = 1:
Si∣∣α− 13
6
∣∣ < 13 entonces z = 1 es un atractor. Para α = 2 entonces
es superatractor.
Si∣∣α− 13
6
∣∣ = 13 , z = 1 es un punto parabolico.
Finalmente, z = 1 es un punto fijo repulsor para cualquier otro valorde α.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 6 / 26
Estabilidad de z = s1, s2
Lema
Estabilidad de los puntos extranos
s1,2 =−3 + 2α±
√5− 12α + 4α2
2
Si |α− 3| < 12 entonces s1 y s2 son puntos atractores; en particular, si
α = 3 s1,2 = 12
(3±√
5)
son superatractores
Si |α− 3| = 12 entonces s1 y s2 son parabolicos, en particular, si
α = 52 ⇒ s1 = s2 = 1.
Para cualquier otro valor de α ∈ C, s1 y s2 son puntos fijos repulsores.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 7 / 26
Puntos crıticos
O ′p (z) = 0
⇓z = 0,∞, z = 3−4α+2α2±
√−6α+19α2−16α3+4α4
3(α−1) = c1, c2
Ademas, c1 = 1c2
, por lo tanto
c1 = c2 = ±1⇔ −6α + 19α2 − 16α3 + 4α4 = 0⇒ α = 0,1
2,
3
2, 2.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 8 / 26
Puntos crıticos
El numero de punto crıticos depende de α:
Si α = 0, c1 = c2 = −1, que es un punto pre-periodico.
Si α = 12 , c1 = c2 = −1 = s1 = s2, es repulsor.
Si α = 32 , c1 = c2 = 1 es repulsor.
Si α = 2, c1 = c2 = 1 es superatractor.
Para cualquier otro valor de α ∈ C hay dos puntos crıticos distintos.
Ademas
Si α→ 1, entonces c1 → 0 y c2 →∞ ⇒ Op (z) = z4.
Si α→ ±∞ entonces c1 → 0 y c2 → ±∞ ⇒ Op (z) = z2.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 9 / 26
Espacio de parmetros
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 10 / 26
La cabeza
La cabeza del gato corresponde a valores de α donde z = 1 es atractor, esdecir
∣∣α− 136
∣∣ < 13 .
Cuando α = 2 los puntos crıticos c1 = c2 = 1, por lo que z = 1 es unpunto fijo superatractor.
−2 −1.5 −1 −0.5 0 0.5 1 1.5 2 2.5 3
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2−2 −1.5 −1 −0.5 0 0.5 1 1.5 2 2.5 3
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2 −1.5 −1 −0.5 0 0.5 1 1.5 2 2.5 3
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2
Planos de dinamicos para α = 2, α = 1,9 y α = 2,2.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 11 / 26
La frontera de la cabeza
La frontera es una circunferencia∣∣α− 13
6
∣∣ = 13 ⇒ α = 13
6 + 13 e iθ.
O ′p (1) =2e iθ + 1
2 + e iθ,∣∣O ′p (1)
∣∣ = 1
Para cada θ, α esta en la interseccion con los distintos bulbos ⇒ aparecenciclos periodicos
−1.5 −1 −0.5 0 0.5 1 1.5 2 2.5 3
−1.5
−1
−0.5
0
0.5
1
1.5
2
Plano dinamico para α = 136 + ( 1
3 + 0,01)e2iπ
3
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 12 / 26
La frontera de la cabeza
θ = π, α = 116
−1.5 −1 −0.5 0 0.5 1 1.5 2 2.5
−1.5
−1
−0.5
0
0.5
1
1.5
z = 1 tiene multiplicidad 1 ⇒ frontera de dos cuencas parabolicas.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 13 / 26
La frontera de la cabeza
Plano dinamico para α = 116 − 0,05.
z=0.95514+i0.29616
−1.5 −1 −0.5 0 0.5 1 1.5 2 2.5
−1.5
−1
−0.5
0
0.5
1
1.5
En el bulbo de la izquierda aparece un ciclo atractor de periodo 2.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 14 / 26
La frontera de la cabeza
θ = 0, α = 52 .
z=0.87689+i0.0049956
−2 −1 0 1 2 3 4
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2
z = 1 tiene multiplicidad 3 ⇒ frontera de dos cuencas atractoras.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 15 / 26
El cuerpo
El cuerpo corresponde |α− 3| < 12 .
z = 1 se convierte en repulsor y s1 y s2 en atractores ⇒ Dos cuencas deatraccion, un crıtico en cada cuenca.
−2 −1 0 1 2 3 4
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2−2 −1 0 1 2 3 4 5
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2
Planos dinamicos para α = 2,6 y α = 3,2 + 0,2iAlicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 16 / 26
El cuerpo
Si α = 3, s1 = c1 y s2 = c2 superatractores.
−2 −1 0 1 2 3 4
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 17 / 26
La frontera del cuerpo
|α− 3| = 12 , α = 3 + 1
2 e iθ ⇒ s1 y s2 son parabolicos
O ′ (s1) = O ′ (s2) = −e iθ = e i(π+θ)
para cada valor de θ ⇒ puntos de interseccion con los bulbos dondeaparecen pares de ciclos atractores.
z=3.4241+i−0.98312
−2 −1 0 1 2 3 4 5
−2
−1.5
−1
−0.5
0
0.5
1
1.5
Ciclos de periodo 3: α = 3 + 0,51eπi3 , α = 3 + 0,51e
5πi3 .
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 18 / 26
La frontera del cuerpo
θ = 0, α = 72 , s1 y s2 son parabolicos ⇒ estan en la interseccion con el
bulbo que contiene ciclos de periodo 2.
z=3.6472+i0.0013048
−2 −1 0 1 2 3 4 5 6
−2
−1.5
−1
−0.5
0
0.5
1
1.5
z=4.0623+i−4.656e−006
−2 −1 0 1 2 3 4 5 6
−2
−1.5
−1
−0.5
0
0.5
1
1.5
Planos dinamicos para α = 3,5 y α = 3,55.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 19 / 26
Dentro del collar
Lema
Si |α− 1| ≤ 12 entonces, el plano dinamico es el mismo que zn.
Si α = 1 entonces Op (z) = z4,Si α = 1
2 entonces Op (z) = z3,Si α = 3
2 entonces Op (z) = −z3.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 20 / 26
Dentro del collar
−1.5 −1 −0.5 0 0.5 1 1.5
−1
−0.5
0
0.5
1
1.5
−1.5 −1 −0.5 0 0.5 1 1.5
−1
−0.5
0
0.5
1
1.5
Planos dinamicos para α = 1 + 0,35e2πi
3 y α = 1 + 0,45eπi6 .
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 21 / 26
Dentro del collar
−2.5 −2 −1.5 −1 −0.5 0 0.5 1 1.5
−1.5
−1
−0.5
0
0.5
1
1.5
−2 −1.5 −1 −0.5 0 0.5 1 1.5
−2
−1.5
−1
−0.5
0
0.5
1
1.5
Planos dinamicos para α = 0,2 + 0,1i α = 0,4− 0,7i .
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 22 / 26
Sobre el collar
0 ≤ α < 12 y α ∈ R ⇒ antena izquierda.
Lema
El plano dinamico para valores α ∈ R y 0 ≤ α < 12 consta de dos cuencas
de atraccion, A(0) and A(∞), con infinitas pre-imagenes.
−2.5 −2 −1.5 −1 −0.5 0 0.5 1 1.5
−1.5
−1
−0.5
0
0.5
1
1.5
−1.5 −1 −0.5 0 0.5 1 1.5 2
−1.5
−1
−0.5
0
0.5
1
1.5
Planos dinamicos para α = 0,2 and α = 1,6.Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 23 / 26
Fuera del gato
El conjunto de Julia es disconexo
−4 −3 −2 −1 0 1 2
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2
−2 −1 0 1 2 3 4 5 6 7 8
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2
Planos dinamicos para α = −0,3 and α = 4,5.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 24 / 26
Conclusiones
El comportamiento dinamico de la familia Tchebychev-Halley es muyrico, como muestra su plano de parametros.
Numericamente, son especialmente interesantes:
los valores de α situados dentro del collar (especiamente α = 32 ) cuyo
conjunto de Julia es conexo, oaquellos que estan fuera del conjunto gato (Julia disconexo).
El analisis dinamico de las familias de metodos iterativos es unabuena herramienta para el analisis de su estabilidad.
Lıneas abiertas:
Profundizar en la comprension dinamico del conjunto gato.Completar el analisis con funciones no polinomicas.Extender el estudio a otras familias parametricas de metodos iterativos.
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 25 / 26
Bibliografıa
S.Amat, S.Busquier, J.M. Gutierrez. Geometric constructions ofiterative functions to solve nonlinear equations. J. of Computationaland Applied Mathematics 157 (2003) 197205.
P. Blanchard. The Dynamics of Newton’s Method. Proc. of Symposiain Applied Math. 49 (1994) 139-154.
P. Blanchard. Complex Analytic Dynamics on the Riemann Sphere.Bull. of the AMS. 11 no. 1 (1984) 85-141.
A. Douady and J.H.Hubbard. On the dynamics of polynomials-likemappings. Ann. Sci. Ec. Norm. Sup. (Paris) 18 (1985) 287-343.
N. Fagella. Invariants in dinamica complexa. Butlletı de la Soc. Cat.de Matematiques, 23(1) (2008) 29-51.
J. Milnor. Dynamics in one complex variable. Stony Brook IMSprepint (1990).
Alicia Cordero, Juan R. Torregrosa, Pura Vindel ()Dinamica de la familia de metodos iterativos de Tchebychev-Halley30 de octubre de 2012 26 / 26