Fractales: de reglas sencillas a la complejidad de la naturaleza

¿Puede un triángulo perder casi toda su área sin desaparecer? ¿Puede una curva crecer indefinidamente dentro de una región pequeña? ¿Unas pocas instrucciones pueden generar un paisaje complejo? Los fractales permiten investigar estas preguntas con ejemplos sorprendentes y comprobables.

Empezamos con lápiz y papel, introducimos potencias y logaritmos y llegamos a las imágenes digitales. No necesitas números complejos para entender las primeras construcciones: la parte avanzada se explica por separado.

Abrir el laboratorio de fractales

1. ¿Qué es un fractal?

Un cuadrado ampliado sigue pareciendo un cuadrado, pero mirar un borde no revela otro cuadrado completo una y otra vez. En muchas construcciones fractales, partes cada vez menores reproducen la organización del conjunto: es la autosimilitud. Es exacta en el triángulo de Sierpiński; en un árbol real puede ser aproximada o estadística.

Benoît Mandelbrot introdujo el término fractal en 1975, a partir del latín fractus, roto. Ser irregular no basta para ser fractal, y no todos los fractales son copias perfectas de sí mismos ni tienen dimensión no entera. Estudiaremos modelos precisos con reglas explícitas. (B. Mandelbrot · Fractal Roughness (WPI)).

Distingamos el objeto matemático límite de su imagen. El primero se define mediante un proceso potencialmente infinito; la segunda termina tras un número finito de pasos y píxeles. Añadir detalle no significa alcanzar el infinito.

2. Sierpiński: menos área, más triángulos

Dibuja un triángulo equilátero lleno. Une los puntos medios de sus lados: aparecen cuatro triángulos iguales. Elimina el central y repite en los otros tres. El triángulo inicial es n = 0; cada nivel nuevo triplica el número de triángulos y reduce su lado a la mitad.

¿Por qué el factor del área es 3/4? Quedan tres triángulos, cada uno con lado mitad y, por tanto, un cuarto del área anterior: 3 × 1/4 = 3/4. Tras tres niveles hay 27 triángulos con lado 1/8 del inicial y área conjunta 27/64 del área inicial.

Para convertir la razón en porcentaje multiplicamos explícitamente por 100%: (27/64) × 100% = 42,1875%. En el nivel 5 quedan 243 triángulos, pero solo (3/4)⁵ × 100% ≈ 23,7305% del área inicial. Aumenta el número de piezas mientras disminuye el área. En el límite esta tiende a cero sin que el conjunto quede vacío.

Nₙ = 3ⁿ
lₙ / l₀ = (1/2)ⁿ
Aₙ / A₀ = 3ⁿ × (1/2)²ⁿ = (3/4)ⁿ
A% = (3/4)ⁿ × 100%

3. Cantor: quitar sin vaciarlo todo

Parte de un segmento de longitud 1. Elimina el tercio central: quedan dos segmentos de longitud 1/3. Repite en ambos: cuatro segmentos de longitud 1/9. Tras n pasos quedan 2ⁿ intervalos, cada uno de longitud 3⁻ⁿ.

La longitud total restante es 2ⁿ × 3⁻ⁿ = (2/3)ⁿ. En el nivel 4 vale 16/81; en porcentaje, (16/81) × 100% ≈ 19,7531%. Las filas del laboratorio representan momentos distintos: no sumes las longitudes de filas diferentes.

En el límite la longitud total es cero, pero quedan puntos: los extremos 0 y 1 nunca se eliminan. Longitud nula no equivale a ausencia de puntos. La dimensión del modelo límite es aproximadamente 0,63093: más fino que un segmento, no simplemente un segmento corto.

Nₙ = 2ⁿ
lₙ / l₀ = (1/3)ⁿ
Lₙ / L₀ = (2/3)ⁿ
L% = (2/3)ⁿ × 100%

4. Koch: longitud infinita, área finita

Divide un segmento en tres partes iguales. Sustituye el tercio central por los otros dos lados de un triángulo equilátero orientado hacia fuera. Un segmento se convierte en cuatro de longitud un tercio. Repite en todos los nuevos segmentos.

El factor de longitud es 4 × 1/3 = 4/3. Partiendo de 9 cm, tras tres pasos obtenemos 9 × (4/3)³ = 64/3 cm ≈ 21,3333 cm. Cada segmento se acorta, pero su número hace que la longitud total crezca sin límite.

Al partir de los tres lados de un triángulo equilátero aparece el copo de Koch. Con lado inicial de 10 cm, su perímetro inicial es 30 cm; el nivel 3 tiene 192 segmentos y perímetro 30 × (4/3)³ ≈ 71,1111 cm. El área encerrada tiende, en cambio, a 8/5 del área del triángulo inicial. Perímetro infinito y área finita describen el límite, no el dibujo finito.

Para el área sumamos los incrementos: el primero es A₀/3 y cada siguiente es 4/9 del anterior. La serie geométrica da Aₙ/A₀ = 8/5 − (3/5)(4/9)ⁿ. La dimensión fractal corresponde al borde del copo; su región interior llena tiene dimensión 2.

Nₙ = 4ⁿ; Lₙ = L₀ × (4/3)ⁿ

Curva y copo de Koch
Pₙ = P₀ × (4/3)ⁿ
Aₙ / A₀ = 8/5 − (3/5)(4/9)ⁿ
A∞ / A₀ = 8/5 = 1,6

5. Dimensión: ¿cómo ocupa espacio una estructura?

Dos copias a mitad de escala reconstruyen un segmento: 2 × (1/2)ᴰ = 1 da D = 1. Un cuadrado necesita cuatro copias de lado mitad: 4 × (1/2)ᴰ = 1 da D = 2. Un cubo necesita ocho y D = 3.

En modelos autosimilares con copias adecuadamente separadas se usa N × rᴰ = 1, donde N es el número de copias y r su razón de longitudes. Los logaritmos dan D = log(N)/log(1/r). La base no importa si es la misma en numerador y denominador.

Para Sierpiński, N = 3 y r = 1/2 dan D ≈ 1,58496. Para Koch, N = 4 y r = 1/3 dan D ≈ 1,26186. No apliques automáticamente la fórmula a una foto o a cualquier árbol: hay que comprobar los solapamientos y las hipótesis. El laboratorio muestra la dimensión del límite, no la de los triángulos llenos de un nivel finito.

D = log(N) / log(1/r)

6. Más ejemplos: alfombras y esponjas

La alfombra de Sierpiński empieza con un cuadrado dividido en nueve pequeños cuadrados. Elimina el central y repite en los ocho restantes. Tras n pasos quedan 8ⁿ cuadrados con área total (8/9)ⁿ de la inicial. Dimensión límite: log(8)/log(3) ≈ 1,89279.

La esponja de Menger es tridimensional. Divide un cubo en 27 cubitos y retira el central y los seis centrales de las caras: quedan 20. Tras n pasos, el volumen es (20/27)ⁿ del inicial y la dimensión límite log(20)/log(3) ≈ 2,72683.

Compara las reglas: 3 copias a media escala para el triángulo, 8 a un tercio para la alfombra y 20 a un tercio para la esponja. La complejidad nace de repetir una regla breve. Estos últimos dos modelos son ejemplos del artículo, no modos adicionales del laboratorio.

7. Construir un árbol con recursión

Dibuja un tronco y dos ramas más cortas e inclinadas en su extremo. Repite en las nuevas puntas. Con r = 0,65 y tronco de longitud 1, las generaciones tienen longitudes 1; 0,65; 0,4225; 0,274625; etc.

El nivel 0 contiene un segmento. En el nivel 4 hay 1 + 2 + 4 + 8 + 16 = 31 segmentos, incluidos 16 terminales. Aquí se conservan también las generaciones anteriores: no es el recuento de los triángulos finales de Sierpiński.

Al cambiar ángulo y razón en el laboratorio aparecen árboles estrechos, anchos o con ramas solapadas. La semejanza con un árbol natural es un modelo geométrico, no una simulación botánica. No asignamos una dimensión mediante log(2)/log(1/r), porque los solapamientos pueden alterar esa interpretación.

Nₙ = 1 + 2 + 4 + … + 2ⁿ = 2ⁿ⁺¹ − 1
Nₜ = 2ⁿ
lₙ / l₀ = rⁿ

8. Iterar una fórmula: los números complejos

Iterar es aplicar la misma regla al resultado anterior. Por ejemplo z → z² + 1, partiendo de 0, produce 0, 1, 2, 5, 26. Cada resultado prepara el paso siguiente.

Mandelbrot y Julia usan el plano complejo. El número z = x + iy corresponde al punto (x,y), con i² = −1. El parámetro c = a + ib es otro punto. Separando las partes real e imaginaria, la fórmula se traduce en dos operaciones con números reales, sin una biblioteca especial.

Una órbita es la sucesión z₀, z₁, z₂, … : la historia del cálculo, no el movimiento de un planeta. Algunas órbitas se mantienen acotadas, otras crecen sin límite. Esta diferencia organiza las imágenes.

zₙ₊₁ = zₙ² + c
(x+iy)² + (a+ib) = (x²−y²+a) + i(2xy+b)
xₙ₊₁ = xₙ² − yₙ² + a
yₙ₊₁ = 2xₙyₙ + b

9. Mandelbrot: un mapa de parámetros

En Mandelbrot fijamos z₀ = 0 y cambiamos c de un píxel a otro. Pertenecen al conjunto los parámetros cuya órbita permanece acotada. Cada píxel pregunta: ¿qué le ocurre a 0 con este c?

Ejemplos reales: c = 0 mantiene la órbita en 0; c = −1 alterna 0 y −1; c = −2 da 0, −2, 2, 2, … . Con c = 1 obtenemos 0, 1, 2, 5, … : en el paso 3 el módulo supera 2 y el escape queda confirmado. La desigualdad es estricta: alcanzar exactamente 2 no basta.

El ordenador solo prueba un número finito de iteraciones. No escapar en 160 pasos no demuestra que nunca escapará. En el laboratorio el negro significa ausencia de escape observado dentro del límite; los colores distinguen los tiempos de escape. Aumentar las iteraciones y ampliar explora detalle numérico, no demuestra la pertenencia de cada punto. (R. Devaney · Mandelbrot Set Explorer (Boston University)).

10. Julia: cambiar el punto inicial

En Julia intercambiamos lo fijo y lo variable: elegimos un único parámetro c y probamos muchos puntos iniciales z₀. El conjunto de Julia lleno contiene los puntos iniciales con órbitas acotadas; el conjunto de Julia propiamente dicho es su borde.

Con c = 0 elevamos z al cuadrado repetidamente. Si |z₀| < 1, el módulo tiende a 0; si |z₀| > 1, crece; si vale 1, permanece en 1. El conjunto lleno es el disco unidad y su borde una circunferencia. No todo parámetro produce un borde irregular.

Prueba c = −1 o c = −0,8 + 0,156i: las formas cambian mucho. El criterio de escape es |z| > max(2, |c|), seguro incluso si |c| supera 2. El test es finito y los puntos oscuros siguen siendo candidatos no escapados salvo prueba adicional, como en c = 0. (Boston University · Filled Julia sets).

11. ¿La naturaleza es realmente fractal?

El romanesco repite pequeños conos semejantes a la inflorescencia completa; los folíolos de un helecho recuerdan su fronde. Árboles y redes fluviales presentan jerarquías de ramificaciones. Muestran organización a varias escalas, no copias geométricas exactas.

Una costa medida con una regla más corta parece más recortada: entran en el cálculo pequeñas bahías antes ignoradas. El resultado depende de la escala y del método de medida. Los modelos estadísticos describen esta dependencia dentro de un intervalo de escalas. (B. Mandelbrot · Coastlines and fractal dimension (PNAS, 1975)).

Nubes, perfiles montañosos y rayos pueden mostrar complejidad multiescala; vasos sanguíneos y bronquios tienen ramificaciones. Una foto llamativa no demuestra una ley fractal: hacen falta datos, escalas declaradas y un modelo comprobado. Células, granos, moléculas y procesos físicos interrumpen la repetición.

12. Aplicaciones: para qué sirve esta geometría

La gráfica procedimental combina detalles a varias escalas para generar terrenos, nubes y texturas. Unos pocos parámetros controlan la imagen en vez de dibujar cada detalle. El ruido fractal no es el triángulo de Sierpiński, pero comparte la organización multiescala. (K. Perlin · Procedural Textures (NYU)).

Algunas antenas usan geometrías fractales para obtener estructuras compactas o comportamiento multibanda. La forma sola no garantiza mejores prestaciones: también importan materiales, frecuencias y diseño electromagnético. Es una aplicación concreta, no una propiedad mágica de cualquier dibujo fractal. (Fractal Antenna Systems · Technical overview).

El análisis de imágenes permite comparar complejidad y ramificaciones; arte y diseño aprovechan repeticiones y cambios de escala. En el recuento de cajas se superponen cuadrículas de paso ε y se cuentan las celdas que tocan la figura. La pendiente de log N(ε) frente a log(1/ε) puede estimar la dimensión solo a escalas adecuadas. Pocos píxeles o cuadrículas no producen un valor universal.

13. Experimentos guiados en el laboratorio

1. Sierpiński: parte del nivel 0. Antes de avanzar predice el número de triángulos y el área. En el nivel 5 compara 243 triángulos con aproximadamente 23,73% de área restante.

2. Cantor: cuenta los 16 intervalos del nivel 4 y verifica su longitud conjunta 16/81. 3. Koch: cambia de curva a copo y compara el perímetro creciente con la razón de áreas que se aproxima a 1,6.

4. Árbol: elige nivel 4 y razón 0,65; cambia después el ángulo. El total sigue siendo 31 aunque cambia la forma. 5. Mandelbrot: elige c = 1 y después c = −1; compara las órbitas. 6. Julia: elige c = 0, observa el disco y prueba módulos iniciales menores y mayores que 1.

Los cálculos se realizan en el navegador sin enviar parámetros al servidor. Los niveles máximos permiten usarlo en teléfonos; el dibujo complejo puede interrumpirse. La animación avanza por niveles y las tablas permiten comparar. Una condición de parada es esencial: sin ella la recursión seguiría llamándose a sí misma.

Abrir el laboratorio de fractales

14. Ponte a prueba

Haz primero una predicción o un cálculo. Las soluciones están ocultas: abre solo la que quieras comprobar.

1. ¿Qué área queda en Sierpiński en el nivel 4?

2. Un segmento de Koch empieza con 27 cm. ¿Qué longitud tiene en el nivel 3?

3. ¿Por qué Cantor puede tener longitud cero sin estar vacío?

4. ¿El árbol del nivel 6 tiene 64 o 127 segmentos?

5. ¿Un píxel negro tras 400 iteraciones pertenece con certeza a Mandelbrot?

6. En Julia con c = 0, ¿qué ocurre con z₀ = 1/2 y z₀ = 2?

Para profundizar y fuentes

• B. Mandelbrot · Fractal Roughness (WPI)
• R. Devaney · Mandelbrot Set Explorer (Boston University)
• Boston University · Filled Julia sets
• B. Mandelbrot · Coastlines and fractal dimension (PNAS, 1975)
• K. Perlin · Procedural Textures (NYU)
• Fractal Antenna Systems · Technical overview