73.1. Del ejemplo a la estructura común

Los capítulos anteriores construyeron dos fractales clásicos mediante reglas aparentemente diferentes. En el conjunto de Cantor se conservaban dos copias reducidas del intervalo. En la curva de Koch un segmento se reemplazaba por cuatro copias reducidas, rotadas y trasladadas. En ambos casos aparecía una misma idea, según la cual el objeto completo podía reconstruirse como unión de imágenes contraídas de sí mismo.

Los sistemas de funciones iteradas, o IFS por iterated function systems, formalizan esa arquitectura.

La abstracción es útil porque separa tres niveles que conviene no confundir. El primero son las transformaciones individuales, que actúan sobre puntos. El segundo es el operador inducido, que actúa sobre conjuntos. El tercero es el conjunto límite, que aparece como punto fijo de ese operador.

En los ejemplos clásicos, cada transformación individual es contractiva. Si se itera una sola de ellas sobre un punto, la trayectoria converge a un punto fijo ordinario. El fractal aparece porque se aplican simultáneamente varias contracciones al conjunto completo. La unión de sus imágenes puede conservar una estructura extensa y no trivial aunque cada mapa por separado reduzca distancias.

Este capítulo desarrollará esa afirmación con precisión. Definiremos contracciones en espacios métricos completos, introduciremos la métrica de Hausdorff sobre compactos no vacíos y demostraremos que el operador de Hutchinson es contractivo. El teorema del punto fijo de Banach garantiza entonces la existencia y unicidad del atractor del IFS.

La construcción permitirá reinterpretar Cantor y Koch, generar el triángulo de Sierpiński y comprender por qué diferentes condiciones iniciales producen el mismo conjunto límite.

73.2. Transformaciones contractivas

Sea (X,d)(X,d) un espacio métrico. Una aplicación

f:X→Xf:X\to X

es contractiva si existe una constante

0≤s<10\leq s<1

tal que

d(f(x),f(y))≤s d(x,y)d(f(x),f(y)) \leq s\,d(x,y)

para todo x,y∈Xx,y\in X.

El menor valor admisible de ss, cuando existe, puede interpretarse como factor de Lipschitz de la transformación. La propiedad que importa es la desigualdad estricta, s<1.s<1.

Una aplicación con s=1s=1 puede preservar distancias y no tiene por qué poseer un punto fijo atractivo. Una rotación del círculo es un ejemplo inmediato.

En Rm\mathbb R^m, una transformación afín

f(x)=Ax+bf(x)=Ax+b

es contractiva respecto de la norma euclídea si

∥A∥2<1,\|A\|_2<1,

donde ∥A∥2\|A\|_2 es la norma operatorial inducida. En una similitud pura,

f(x)=rRx+b,f(x)=rRx+b,

con RR ortogonal y

0≤r<1,0\leq r<1,

se tiene exactamente

∥f(x)−f(y)∥=r∥x−y∥.\|f(x)-f(y)\| = r\|x-y\|.

Las transformaciones de Cantor y Sierpiński pertenecen a esta clase.

73.3. El teorema del punto fijo de Banach

La herramienta fundamental es el principio de contracción de Banach.

Sea XX un espacio métrico completo y

f:X→Xf:X\to X

una contracción con factor s<1s<1. Entonces existe un único punto

x∗∈Xx^\ast\in X

tal que

f(x∗)=x∗.f(x^\ast)=x^\ast.

Además, para cualquier condición inicial x0x_0, xn+1=f(xn)x_{n+1}=f(x_n) converge a x∗x^\ast.

La demostración contiene la tasa de convergencia. Como d(xn+1,xn)≤snd(x1,x0),d(x_{n+1},x_n) \leq s^n d(x_1,x_0), la distancia entre iterados separados puede acotarse por una serie geométrica. La completitud garantiza que la sucesión de Cauchy tiene límite en XX. La continuidad de una contracción permite pasar al límite y obtener el punto fijo.

La unicidad también es inmediata. Si pp y qq fueran puntos fijos distintos,

d(p,q)=d(f(p),f(q))≤s d(p,q).d(p,q) = d(f(p),f(q)) \leq s\,d(p,q).

Como s<1s<1, esto obliga a d(p,q)=0.d(p,q)=0.

Por tanto, p=q.p=q.

El teorema se refiere todavía a puntos. Para producir fractales necesitamos aplicar la misma lógica a un espacio cuyos «puntos» sean conjuntos compactos.

73.4. El espacio de compactos

Sea XX un espacio métrico completo. Denotemos por

K(X)\mathcal K(X)

el conjunto de todos los subconjuntos compactos no vacíos de XX.

Queremos medir la distancia entre dos conjuntos

A,B∈K(X).A,B\in\mathcal K(X).

La distancia ordinaria entre puntos no basta. Una posibilidad natural es preguntar cuánto debe expandirse uno de los conjuntos para cubrir al otro.

Para un punto xx y un conjunto BB, definimos

d(x,B)=inf⁡b∈Bd(x,b).d(x,B) = \inf_{b\in B}d(x,b).

La distancia dirigida de AA a BB es

h(A,B)=sup⁡a∈Ad(a,B).h(A,B) = \sup_{a\in A}d(a,B).

Esta cantidad no es simétrica. Puede ocurrir que todos los puntos de AA estén cerca de BB, mientras BB tiene regiones lejanas de AA.

La distancia de Hausdorff corrige ese problema,

dH(A,B)=max⁡{h(A,B),h(B,A)}.d_H(A,B) = \max\{h(A,B),h(B,A)\}.

Para compactos no vacíos, dHd_H es una métrica.

Un resultado central es que si XX es completo, entonces K(X)\mathcal K(X), equipado con dHd_H, también es completo.

Ese hecho permite aplicar Banach a transformaciones que actúan sobre conjuntos.

73.5. Un IFS

Un sistema de funciones iteradas consiste en una familia finita de aplicaciones contractivas

F={f1,…,fN}\mathcal F= \{f_1,\ldots,f_N\}

sobre un espacio métrico completo XX.

Supongamos que cada fif_i tiene factor de contracción sis_i,

d(fi(x),fi(y))≤sid(x,y),0≤si<1.d(f_i(x),f_i(y)) \leq s_i d(x,y), \qquad 0\leq s_i<1.

Definimos

s=max⁡isi.s=\max_i s_i.

Entonces s<1.s<1.

El sistema completo no se itera escogiendo necesariamente una sola función. Primero construiremos un operador que actúa sobre conjuntos y aplica todas las transformaciones simultáneamente.

73.6. El operador de Hutchinson

Para

A∈K(X),A\in\mathcal K(X),

definimos

H(A)=⋃i=1Nfi(A).\mathcal H(A) = \bigcup_{i=1}^N f_i(A).

Esta aplicación

H:K(X)→K(X)\mathcal H:\mathcal K(X)\to\mathcal K(X)

se denomina operador de Hutchinson.

La unión sigue siendo compacta porque cada imagen fi(A)f_i(A) es compacta y existe un número finito de ellas.

Cantor aparece inmediatamente. Si

f0(x)=x3,f1(x)=x+23,f_0(x)=\frac{x}{3}, \qquad f_1(x)=\frac{x+2}{3},

entonces

H([0,1])=[0,13]∪[23,1].\mathcal H([0,1]) = \left[0,\frac13\right] \cup \left[\frac23,1\right].

Es decir,

H(C0)=C1.\mathcal H(C_0)=C_1.

Repitiendo,

Hn(C0)=Cn.\mathcal H^n(C_0)=C_n.

La construcción sustractiva de C71 puede reescribirse como iteración positiva de un operador sobre conjuntos.

73.7. Por qué el operador también es contractivo

El paso principal es demostrar

dH(H(A),H(B))≤s dH(A,B).d_H(\mathcal H(A),\mathcal H(B)) \leq s\,d_H(A,B).

Comencemos con una sola función fif_i. Para cualquier a∈Aa\in A, existe un punto de BB tan cercano como queramos a la distancia d(a,B)d(a,B). Por contractividad, d(fi(a),fi(B))≤sid(a,B).d(f_i(a),f_i(B)) \leq s_i d(a,B).

Al tomar supremos, h(fi(A),fi(B))≤sih(A,B).h(f_i(A),f_i(B)) \leq s_i h(A,B).

Repitiendo en la dirección inversa, h(fi(B),fi(A))≤sih(B,A).h(f_i(B),f_i(A)) \leq s_i h(B,A).

Entonces

dH(fi(A),fi(B))≤sidH(A,B).d_H(f_i(A),f_i(B)) \leq s_i d_H(A,B).

Para la unión de todas las imágenes, cada punto de H(A)\mathcal H(A) pertenece a alguna rama fi(A)f_i(A). Puede compararse con la rama correspondiente fi(B)f_i(B), contenida en H(B)\mathcal H(B). Por ello,

h(H(A),H(B))≤s dH(A,B).h(\mathcal H(A),\mathcal H(B)) \leq s\,d_H(A,B).

La misma desigualdad vale en la otra dirección. En consecuencia,

dH(H(A),H(B))≤s dH(A,B).d_H(\mathcal H(A),\mathcal H(B)) \leq s\,d_H(A,B).

El operador de Hutchinson es una contracción sobre el espacio completo K(X)\mathcal K(X).

73.8. Existencia y unicidad del atractor del IFS

Ahora puede aplicarse el teorema de Banach.

Existe un único conjunto compacto no vacío

A∗∈K(X)A^\ast\in\mathcal K(X)

tal que

H(A∗)=A∗.\mathcal H(A^\ast)=A^\ast.

Es decir,

A∗=⋃i=1Nfi(A∗).A^\ast = \bigcup_{i=1}^N f_i(A^\ast).

Este conjunto se denomina habitualmente atractor del IFS.

Además, para cualquier conjunto compacto no vacío inicial A0A_0,

An+1=H(An)A_{n+1} = \mathcal H(A_n)

converge a A∗A^\ast en métrica de Hausdorff.

La tasa satisface

dH(An,A∗)≤sndH(A0,A∗).d_H(A_n,A^\ast) \leq s^n d_H(A_0,A^\ast).

La conclusión es más fuerte que «la imagen parece estabilizarse». La existencia del límite, su unicidad y la independencia respecto del conjunto inicial se deducen de la contractividad.

Esta independencia debe interpretarse correctamente. Las aproximaciones intermedias dependen de A0A_0. El límite no.

73.9. Un atractor de conjuntos no es un punto fijo individual

La palabra «atractor» puede inducir una confusión.

Cada mapa contractivo fif_i posee su propio punto fijo ordinario. Sin embargo, el atractor del IFS no es, en general, ninguno de esos puntos.

En el IFS de Cantor, f0(x)=x/3f_0(x)=x/3 tiene punto fijo 00, mientras f1(x)=(x+2)/3f_1(x)=(x+2)/3 tiene punto fijo 11.

El atractor del sistema completo es el conjunto de Cantor,

C=f0(C)∪f1(C).C=f_0(C)\cup f_1(C).

Por tanto,

punto fijo de fi≠punto fijo de H.\text{punto fijo de }f_i \neq \text{punto fijo de }\mathcal H.

El primero es un punto de XX. El segundo es un elemento de K(X)\mathcal K(X), es decir, un conjunto compacto.

Esta distinción es central. El fractal surge porque elevamos el problema desde el espacio de estados al espacio de conjuntos.

73.10. Cantor como IFS contractivo

Volvamos al conjunto de Cantor con

f0(x)=x3,f1(x)=x+23.f_0(x)=\frac{x}{3}, \qquad f_1(x)=\frac{x+2}{3}.

Ambas transformaciones poseen razón de contracción

s0=s1=13.s_0=s_1=\frac13.

El operador de Hutchinson satisface

dH(H(A),H(B))≤13dH(A,B).d_H(\mathcal H(A),\mathcal H(B)) \leq \frac13 d_H(A,B).

El atractor único es el conjunto compacto CC que satisface

C=13C∪(23+13C).C = \frac13C \cup \left( \frac23+\frac13C \right).

Si iniciamos con A0=[0,1],A_0=[0,1], obtenemos exactamente las aproximaciones de C71, An=Cn.A_n=C_n.

En este caso puede calcularse incluso la distancia exacta al límite,

dH(Cn,C)=12 3n+1.d_H(C_n,C) = \frac{1}{2\,3^{n+1}}.
Convergencia en métrica de Hausdorff de las aproximaciones de Cantor hacia el atractor del IFS.
Figura 73.3. Convergencia en métrica de Hausdorff de las aproximaciones de Cantor hacia el atractor del IFS. Ampliar figura ↗

La figura 73.3 muestra esta convergencia geométrica. Cada iteración reduce la distancia por un factor

13,\frac13,

coherente con el factor de contracción del IFS.

73.11. Koch como sistema de cuatro contracciones

La curva de Koch también puede describirse mediante cuatro similitudes del plano. Cada una reduce longitudes por un factor

13.\frac13.

Las cuatro ramas incluyen traslaciones y rotaciones para situar las copias sobre los cuatro segmentos del generador.

Si las denotamos

f1,f2,f3,f4,f_1,f_2,f_3,f_4,

la curva límite KK satisface

K=f1(K)∪f2(K)∪f3(K)∪f4(K).K = f_1(K)\cup f_2(K)\cup f_3(K)\cup f_4(K).

Todas las transformaciones son contractivas con factor 1/31/3. Por el teorema del IFS, existe un único compacto que satisface esa ecuación.

Esto aporta una lectura distinta de C72. Allí construimos la curva sustituyendo segmentos. Ahora podemos afirmar que la curva de Koch es el punto fijo de una contracción que actúa sobre conjuntos compactos.

La equivalencia entre ambas construcciones no depende de la apariencia de las aproximaciones. Se deriva de la misma ecuación de invariancia.

73.12. El triángulo de Sierpiński como ejemplo nuevo

Consideremos los vértices de un triángulo equilátero,

v1=(0,0),v_1=(0,0),
v2=(1,0),v_2=(1,0),
v3=(12,32).v_3= \left( \frac12,\frac{\sqrt3}{2} \right).

Definimos

fi(x)=12x+12vi,i=1,2,3.f_i(x) = \frac12 x+\frac12v_i, \qquad i=1,2,3.

Cada mapa lleva el triángulo completo a uno de sus tres subtriángulos de esquina y reduce todas las distancias por un factor

s=12.s=\frac12.

El atractor SS satisface

S=f1(S)∪f2(S)∪f3(S).S = f_1(S)\cup f_2(S)\cup f_3(S).

Ese conjunto es el triángulo de Sierpiński.

Si iniciamos con un solo punto, el operador produce primero tres puntos, luego nueve, después veintisiete. Tras nn iteraciones aparecen 3n3^n puntos si no se eliminan duplicados.

Octava iteración determinista del operador de Hutchinson para el triángulo de Sierpiński.
Figura 73.1. Octava iteración determinista del operador de Hutchinson para el triángulo de Sierpiński. Ampliar figura ↗

La figura 73.1 muestra la octava iteración. No es el atractor exacto, sino una aproximación finita al punto fijo en K(R2)\mathcal K(\mathbb R^2).

73.13. Convergencia desde diferentes conjuntos iniciales

Una consecuencia fuerte del teorema es que no necesitamos comenzar con el triángulo correcto.

Podemos partir de un punto, un círculo, una nube aleatoria o cualquier compacto no vacío apropiado. La iteración

An+1=H(An)A_{n+1} = \mathcal H(A_n)

convergerá al mismo atractor.

La geometría temprana puede ser muy distinta. Si A0A_0 es grande, las primeras etapas contienen copias grandes de su forma. Si A0A_0 es un único punto, las aproximaciones son nubes finitas.

Aun así, la distancia al atractor satisface

dH(An,A∗)≤sndH(A0,A∗).d_H(A_n,A^\ast) \leq s^n d_H(A_0,A^\ast).

Esto permite separar la dependencia transitoria de la unicidad asintótica.

Es el análogo, en el espacio de conjuntos, de la convergencia hacia un punto fijo estable en un sistema dinámico ordinario.

73.14. Qué papel desempeña la completitud

El teorema de Banach requiere completitud.

Una sucesión puede ser de Cauchy dentro de un espacio y, sin embargo, converger a un punto que no pertenece a ese espacio si el espacio es incompleto.

Un ejemplo elemental es el conjunto de números racionales con la distancia usual. Existen sucesiones racionales de Cauchy que convergen a números irracionales.

En el problema del IFS queremos que la sucesión de compactos

A0,A1,A2,…A_0,A_1,A_2,\ldots

tenga su límite dentro del mismo espacio sobre el que actúa H\mathcal H.

Por eso no basta con escribir informalmente «las imágenes se hacen cada vez más parecidas». La completitud es la hipótesis que permite transformar el comportamiento de Cauchy en existencia de un límite.

Para espacios euclídeos, esta condición está satisfecha y K(Rm)\mathcal K(\mathbb R^m) con la métrica de Hausdorff es completo.

73.15. El juego del caos, o iterar una rama por vez

El operador de Hutchinson aplica todas las ramas simultáneamente. Existe otra forma de aproximar el mismo atractor.

Partimos de un punto x0x_0. En cada paso elegimos un índice

In∈{1,…,N}I_n\in\{1,\ldots,N\}

y calculamos xn+1=fIn(xn).x_{n+1} = f_{I_n}(x_n).

Si todas las ramas relevantes se eligen repetidamente bajo condiciones apropiadas, la nube de puntos termina visitando el atractor.

Este procedimiento es conocido como juego del caos.

Para Sierpiński, si elegimos los tres mapas con probabilidad 1/31/3, la nube resultante reproduce visualmente el triángulo fractal después de descartar un transitorio breve.

La explicación intuitiva es que cada composición

fin∘⋯∘fi1f_{i_n}\circ\cdots\circ f_{i_1}

selecciona una región de escala aproximadamente sns^n. Secuencias largas de símbolos recorren regiones cada vez más finas del atractor.

Triángulo de Sierpiński aproximado mediante 100 000 iteraciones del juego del caos.
Figura 73.2. Triángulo de Sierpiński aproximado mediante 100 000 iteraciones del juego del caos. Ampliar figura ↗

La figura 73.2 utiliza 100 000100\,000 iteraciones aleatorias.

73.16. El juego del caos no cambia la ecuación del soporte

Supongamos que asignamos probabilidades

p1,…,pN,pi>0,∑ipi=1.p_1,\ldots,p_N, \qquad p_i>0, \qquad \sum_i p_i=1.

Las transformaciones geométricas siguen siendo las mismas. El compacto invariante definido por

A∗=⋃ifi(A∗)A^\ast = \bigcup_i f_i(A^\ast)

no incorpora esas probabilidades.

Lo que cambia es la frecuencia con que la órbita aleatoria visita distintas regiones.

Con probabilidades iguales, el triángulo de Sierpiński clásico recibe una distribución simétrica de visitas. Con probabilidades desiguales, algunas ramas se visitan con mayor frecuencia.

Por tanto, debemos separar el

soporte geomeˊtrico\text{soporte geométrico}

de la

medida sobre el soporte.\text{medida sobre el soporte}.

Esta diferencia será desarrollada con mayor detalle en C78. Aquí basta advertir que una imagen generada por juego del caos puede parecer «más oscura» en ciertas regiones por la medida de muestreo, aunque el atractor geométrico sea el mismo.

73.17. El código simbólico reaparece

Cada punto generado por una composición de mapas puede asociarse a una palabra

i1i2…in.i_1i_2\ldots i_n.

La palabra indica qué ramas se eligieron.

Para una secuencia infinita

ω=(i1,i2,…),\omega=(i_1,i_2,\ldots),

consideremos

fi1∘fi2∘⋯∘fin(x0).f_{i_1}\circ f_{i_2}\circ\cdots\circ f_{i_n}(x_0).

Debido a la contractividad, el efecto de la condición inicial desaparece cuando n→∞n\to\infty. La composición converge entonces a un punto del atractor bajo las hipótesis estándar.

Esto conecta C73 con C67 y C71. El código simbólico excede la condición de etiqueta, ya que describe una secuencia de contracciones que localiza progresivamente un punto dentro de regiones cada vez menores.

Cuando las imágenes de las ramas se solapan, diferentes códigos pueden corresponder al mismo punto. La codificación no tiene por qué ser inyectiva.

Por tanto,

secuencia simboˊlica uˊnica\text{secuencia simbólica única}

no es una propiedad automática de cualquier IFS.

73.18. El papel del solapamiento

Cantor y Sierpiński tienen ramas con separación clara. En otros IFS, las imágenes

fi(A∗)f_i(A^\ast)

pueden solaparse.

La existencia y unicidad del atractor no requieren que las ramas estén separadas. La contractividad basta para el teorema básico de Hutchinson.

Sin embargo, el solapamiento sí afecta otras propiedades, entre ellas la unicidad de la codificación simbólica, el cálculo simple de dimensión, la distribución de medidas y la estructura local del conjunto.

Esto explica por qué la ecuación de similitud

∑iris=1\sum_i r_i^s=1

no puede utilizarse sin más para cualquier sistema de similitudes solapadas. En C71 y C72 las condiciones geométricas eran especialmente favorables.

El teorema del IFS proporciona el conjunto invariante. La dimensión exige hipótesis adicionales.

73.19. Transformaciones contractivas que no son similitudes

No toda contracción reduce todas las direcciones por el mismo factor.

Una transformación afín

f(x)=Ax+bf(x)=Ax+b

puede contraer más en una dirección que en otra. Por ejemplo,

A=(0.4000.7)A= \begin{pmatrix} 0.4 & 0\\ 0 & 0.7 \end{pmatrix}

contrae ambos ejes de manera anisotrópica.

El IFS sigue siendo contractivo si

∥A∥<1.\|A\|<1.

El atractor existe y es único por el mismo argumento.

Sin embargo, deja de ser estrictamente autosimilar en el sentido de copias obtenidas por una única razón escalar. Puede aparecer autoafinidad, que será tratada en C76.

Esta generalización es importante. El poder de los IFS no depende de que cada rama sea una miniatura rígida del conjunto completo.

73.20. Un contraejemplo con una familia no contractiva

Considérese

f1(x)=x,f2(x)=x+12f_1(x)=x, \qquad f_2(x)=\frac{x+1}{2}

sobre [0,1][0,1].

La primera aplicación tiene factor de Lipschitz s1=1.s_1=1.

Por tanto, la familia no satisface la hipótesis contractiva estricta con s<1.s<1.

El teorema de Hutchinson–Banach, en la forma usada aquí, ya no puede invocarse.

Esto no significa que todo IFS con una rama no contractiva carezca de conjuntos invariantes. Significa algo más preciso, a saber, que la existencia y unicidad global obtenidas por el argumento de contracción dejan de estar garantizadas.

Esta es la forma correcta de leer una hipótesis matemática. Cuando falla una condición suficiente, no se concluye automáticamente el resultado contrario; se pierde el derecho a usar ese teorema.

73.21. Un segundo contraejemplo con una sola rama

Supongamos que tomamos el IFS de Sierpiński, aunque en cada paso aplicamos siempre f1f_1.

Entonces xn+1=f1(xn)x_{n+1}=f_1(x_n) converge al vértice v1v_1.

No aparece el triángulo de Sierpiński.

El atractor fractal pertenece al sistema completo de ramas. Para aproximarlo mediante una sola trayectoria aleatoria necesitamos una regla de selección que explore las ramas.

Esto aclara una confusión frecuente entre dos procedimientos,

An+1=⋃ifi(An)A_{n+1} = \bigcup_i f_i(A_n)

y xn+1=fIn(xn).x_{n+1} = f_{I_n}(x_n).

El primero es una iteración determinista sobre conjuntos. El segundo es una iteración sobre puntos guiada por símbolos o azar.

Pueden aproximar el mismo soporte, aunque no son el mismo algoritmo.

73.22. Cotas a posteriori para estimar el error sin conocer el atractor

La desigualdad

dH(An,A∗)≤sndH(A0,A∗)d_H(A_n,A^\ast) \leq s^n d_H(A_0,A^\ast)

es conceptualmente clara, aunque contiene una cantidad incómoda, la distancia inicial al atractor A∗A^\ast, que es el objeto que a menudo no conocemos.

El principio de contracción permite obtener una cota a posteriori basada solo en dos iterados consecutivos.

Como

A∗=H(A∗),A^\ast=\mathcal H(A^\ast),

tenemos

dH(An,A∗)≤dH(An,An+1)+dH(An+1,A∗).d_H(A_n,A^\ast) \leq d_H(A_n,A_{n+1}) + d_H(A_{n+1},A^\ast).

Por contractividad,

dH(An+1,A∗)=dH(H(An),H(A∗))≤s dH(An,A∗).d_H(A_{n+1},A^\ast) = d_H(\mathcal H(A_n),\mathcal H(A^\ast)) \leq s\,d_H(A_n,A^\ast).

Por tanto,

(1−s)dH(An,A∗)≤dH(An,An+1),(1-s)d_H(A_n,A^\ast) \leq d_H(A_n,A_{n+1}),

y finalmente

dH(An,A∗)≤dH(An,An+1)1−s.d_H(A_n,A^\ast) \leq \frac{d_H(A_n,A_{n+1})}{1-s}.

Esta desigualdad es especialmente útil en cómputo. Si conocemos el factor contractivo y podemos estimar cuánto cambió el conjunto en la última iteración, obtenemos una cota sobre la distancia al límite sin disponer de una representación exacta del atractor.

También puede escribirse una cota desde el comienzo,

dH(An,A∗)≤sn1−sdH(A0,A1).d_H(A_n,A^\ast) \leq \frac{s^n}{1-s} d_H(A_0,A_1).

La demostración utiliza la serie geométrica formada por las diferencias entre iterados sucesivos.

En el IFS de Cantor,

s=13.s=\frac13.

Además,

dH(Cn,Cn+1)=12 3n+1.d_H(C_n,C_{n+1}) = \frac{1}{2\,3^{n+1}}.

Entonces la cota a posteriori produce

dH(Cn,C)≤11−1/312 3n+1=14 3n.d_H(C_n,C) \leq \frac{1}{1-1/3} \frac{1}{2\,3^{n+1}} = \frac{1}{4\,3^n}.

La distancia exacta es aún menor,

dH(Cn,C)=12 3n+1,d_H(C_n,C) = \frac{1}{2\,3^{n+1}},

de modo que la cota es conservadora y correcta.

Esto ilustra una diferencia metodológica que reaparecerá muchas veces. Una cota demostrada no tiene que ser numéricamente ajustada para ser útil, porque su función primaria es garantizar que el error no excede cierto valor.

73.23. El teorema del collage y el problema inverso

Hasta ahora partimos de un IFS y preguntamos cuál es su atractor. Puede plantearse el problema inverso. Dado un conjunto objetivo BB, ¿podemos encontrar un IFS cuyo atractor se parezca a BB?

Supongamos que H\mathcal H es contractivo con factor s<1s<1 y que conseguimos que

dH(B,H(B))≤ε.d_H(B,\mathcal H(B)) \leq \varepsilon.

Es decir, el conjunto objetivo se parece bastante a la unión de sus imágenes bajo las transformaciones elegidas.

Sea A∗A^\ast el atractor del IFS. Entonces

dH(B,A∗)≤dH(B,H(B))+dH(H(B),H(A∗)).d_H(B,A^\ast) \leq d_H(B,\mathcal H(B)) + d_H(\mathcal H(B),\mathcal H(A^\ast)).

El segundo término satisface

dH(H(B),H(A∗))≤s dH(B,A∗).d_H(\mathcal H(B),\mathcal H(A^\ast)) \leq s\,d_H(B,A^\ast).

Por tanto,

dH(B,A∗)≤ε+s dH(B,A∗),d_H(B,A^\ast) \leq \varepsilon + s\,d_H(B,A^\ast),

y

dH(B,A∗)≤ε1−s.d_H(B,A^\ast) \leq \frac{\varepsilon}{1-s}.

Esta desigualdad se conoce habitualmente como teorema del collage en la literatura de IFS.

La idea es poderosa, porque no necesitamos calcular primero el atractor exacto para saber que estará cerca del objetivo. Basta conseguir que el objetivo sea casi invariante bajo el operador elegido.

El nombre «collage» refleja el procedimiento geométrico, que consiste en intentar cubrir el objeto objetivo mediante copias transformadas y contraídas de sí mismo.

Sin embargo, esta estrategia tiene límites importantes.

Primero, una aproximación geométrica buena no implica que el IFS represente el mecanismo que generó el objeto observado.

Segundo, diferentes sistemas de transformaciones pueden producir atractores parecidos dentro de la resolución disponible. El problema inverso puede ser no identificable.

Tercero, la cota depende de

11−s.\frac{1}{1-s}.

Si s≈1,s\approx1, una pequeña discrepancia de collage puede amplificarse en una cota poco informativa. Las contracciones débiles convergen lentamente y ofrecen garantías geométricas menos ajustadas.

Por tanto, el teorema convierte el ajuste visual en una desigualdad, sin resolver automáticamente la identificación causal ni la selección entre modelos rivales.

73.24. Contracción fuerte, convergencia rápida y pérdida de memoria

La constante ss cumple varios papeles simultáneos.

Controla cuánto se contraen las distancias en una iteración,

dH(H(A),H(B))≤s dH(A,B).d_H(\mathcal H(A),\mathcal H(B)) \leq s\,d_H(A,B).

Controla la tasa de pérdida de memoria del conjunto inicial,

dH(Hn(A),Hn(B))≤sndH(A,B).d_H(\mathcal H^n(A),\mathcal H^n(B)) \leq s^n d_H(A,B).

También controla las cotas de error a posteriori mediante el factor

11−s.\frac{1}{1-s}.

Si ss es pequeño, dos inicializaciones distintas se vuelven cercanas rápidamente. Si ss está próximo a uno, la memoria transitoria puede persistir durante muchas iteraciones.

Esto permite cuantificar una frase que de otro modo sería vaga, la de que «el atractor no depende de la condición inicial». Matemáticamente, la independencia es asintótica y la velocidad con que se vuelve visible está gobernada por ss.

Por ejemplo, con

s=12,s=\frac12,

una diferencia inicial entre dos compactos se reduce, como máximo, por un factor 2−n2^{-n} después de nn iteraciones del operador.

Con s=0.99,s=0.99, la cota correspondiente es 0.99n,0.99^n, mucho más lenta.

En una simulación finita, ambos sistemas poseen un atractor único, mientras que el primero puede parecer convergido en pocas etapas y el segundo conserva durante largo tiempo rasgos del conjunto inicial.

La distinción entre unicidad asintótica y velocidad de convergencia es otra razón para no inferir propiedades del límite solo a partir de una figura obtenida con un número arbitrario de iteraciones.

73.25. Demostración del teorema del IFS

Podemos resumir la cadena lógica.

Sea (X,d)(X,d) completo y

F={f1,…,fN}\mathcal F=\{f_1,\ldots,f_N\}

una familia de contracciones con factores si<1.s_i<1.

Definimos

s=max⁡isi<1.s=\max_i s_i<1.

Sea

H(A)=⋃i=1Nfi(A).\mathcal H(A) = \bigcup_{i=1}^N f_i(A).

El espacio adecuado para trabajar es

K(X),\mathcal K(X),

que es completo con la métrica de Hausdorff. Cada rama contrae conjuntos, según

dH(fi(A),fi(B))≤sidH(A,B).d_H(f_i(A),f_i(B)) \leq s_i d_H(A,B).

La unión conserva la cota máxima, de modo que

dH(H(A),H(B))≤sdH(A,B),d_H(\mathcal H(A),\mathcal H(B)) \leq s d_H(A,B),

y por tanto H\mathcal H es una contracción. Aplicando Banach se obtiene que existe un único

A∗∈K(X)A^\ast\in\mathcal K(X)

tal que

H(A∗)=A∗.\mathcal H(A^\ast)=A^\ast.

Para la convergencia, tomemos cualquier A0∈K(X)A_0\in\mathcal K(X). Entonces

An=Hn(A0)A_n=\mathcal H^n(A_0)

satisface

dH(An,A∗)→0.d_H(A_n,A^\ast)\to0.

Además,

dH(An,A∗)≤sndH(A0,A∗).d_H(A_n,A^\ast) \leq s^n d_H(A_0,A^\ast).

La existencia del fractal como límite queda así reducida a un teorema de punto fijo sobre un espacio de conjuntos.

73.26. Verificación computacional

El paquete reproducible implementa tres pruebas distintas.

Primero, reconstruye Cantor mediante el operador de Hutchinson y verifica la cantidad y escala de los intervalos. Para ese caso se registra además la distancia exacta

dH(Cn,C)=12 3n+1.d_H(C_n,C) = \frac{1}{2\,3^{n+1}}.

Segundo, verifica numéricamente que cada mapa de Sierpiński reduce distancias exactamente por un factor 1/21/2.

Tercero, genera pares de nubes finitas A,BA,B y comprueba que la aproximación numérica a la distancia de Hausdorff satisface

dH(H(A),H(B))≤12dH(A,B).d_H(\mathcal H(A),\mathcal H(B)) \leq \frac12 d_H(A,B).

También se genera el triángulo mediante iteración determinista del conjunto y mediante juego del caos.

Se ejecutaron 93 comprobaciones automáticas y todas resultaron satisfactorias.

Los ensayos sobre nubes no sustituyen la demostración general del teorema. Verifican que la implementación concreta respeta la cota que ya fue derivada analíticamente.

73.27. Qué aporta la formalización

Cantor y Koch podían estudiarse sin hablar de espacios de compactos. Sin embargo, el lenguaje de IFS revela que ambos son manifestaciones del mismo mecanismo.

La formalización aporta cuatro cosas. La existencia del conjunto límite deja de depender de la confianza en una sucesión de dibujos. La unicidad garantiza que, bajo contractividad, existe un único compacto invariante. La convergencia global asegura que cualquier compacto inicial converge al mismo atractor. La tasa establece que la métrica de Hausdorff decrece al menos de acuerdo con el factor de contracción.

Estas propiedades convierten una receta gráfica en un problema bien planteado de análisis.

73.28. Qué no aporta

El teorema básico de IFS no responde todas las preguntas sobre el atractor.

No determina automáticamente su dimensión, si distintas ramas se solapan, si la codificación es única, qué medida probabilística es natural sobre el conjunto, cómo estimarlo desde datos empíricos ni si una estructura observada en la naturaleza fue generada efectivamente por un IFS.

Tampoco afirma que todo fractal sea atractor de un IFS finito contractivo.

Un modelo generativo puede reproducir una geometría sin ser el mecanismo causal de un objeto real. Esta distinción entre representación matemática y explicación empírica será importante cuando el proyecto pase de fractales ideales a datos.

73.29. De IFS a dimensión

La formalización de este capítulo deja preparado el siguiente paso.

En un IFS de similitudes con razones

r1,…,rN,r_1,\ldots,r_N,

aparece naturalmente la ecuación

∑i=1Nris=1.\sum_{i=1}^N r_i^s=1.

Para Cantor, 2(1/3)s=1.2(1/3)^s=1.

Para Koch, 4(1/3)s=1.4(1/3)^s=1.

Para Sierpiński, 3(1/2)s=1.3(1/2)^s=1.

En ejemplos con separación apropiada, el exponente ss coincide con dimensiones geométricas importantes.

C73 termina deliberadamente antes de convertir esta ecuación en una definición universal. C74 introducirá la dimensión por conteo de cajas desde una perspectiva operacional. C75 mostrará por qué dimensión de similitud, dimensión de cajas y dimensión de Hausdorff pueden coincidir en estos modelos canónicos y divergir en situaciones más generales.

Referencias

Banach, S. (1922). Sur les opérations dans les ensembles abstraits et leur application aux équations intégrales. Fundamenta Mathematicae, 3, 133–181. https://doi.org/10.4064/fm-3-1-133-181

Hutchinson, J. E. (1981). Fractals and self-similarity. Indiana University Mathematics Journal, 30(5), 713–747. https://doi.org/10.1512/iumj.1981.30.30055

Barnsley, M. F. (1993). Fractals Everywhere (2nd ed.). Academic Press Professional. ISBN 978-0-12-079061-6.

Falconer, K. (2003). Fractal Geometry: Mathematical Foundations and Applications (2nd ed.). Wiley. https://doi.org/10.1002/0470013850

Recursos reproducibles

Descargar recursos del capítulo (ZIP, 324 kB).

2 archivos CSV · 3 figuras · 1 programa Python.

El programa ilustra cálculos del capítulo. Las figuras y los CSV están incluidos, pero el script no vuelve a guardar todos esos archivos.

Ver instrucciones y contenido del paquete.