Skip to main content
4 de 4
etiquetas editadas
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

Obtener y gráficar triangulo de Sierpinski

Necesito ayuda con la generación de una secuencia a partir de

P = np.array([[0,0],
              [1,0],
              [0.5,np.sqrt(3)/2]])

P consta de tres entradas, cada una es un par ordenado que indica el vértice de un triángulo unitario

Con este triángulo necesito generar tres sub-triángulos, los vértices están en la lista T por eso obtengo una lista con 3^2 = 9 array. La cuestión es que quiero generar los siguientes sub-triángulos, lo que implica aplicar las mismas transformaciones a la nueva lista T para generar los nuevos 3^3 =27 vértices.

Este es el código que tengo actualmente:

P = np.array([[0,0],[1,0],[0.5,np.sqrt(3)/2]])

m = []
for p in P:
    m.append(p*0.5)

n=[]
for p in P:
    n.append(p*0.5 + np.array([0.5,0]))

k=[]
for p in P:
    k.append(p*0.5 + np.array([0.25,np.sqrt(3)/4]))

T = np.array([m,n,k])

donde T es de la forma

T = [ 
    array([[0.,  0.], [0.5, 0.], [0.25, 0.4330127]]), 
    array([[0.5,  0.], [1, 0.], [0.75, 0.4330127]]), 
    array([[0.25, 0.4330127], [0.75, 0.4330127],[0.5, 0.8660254]]) ]

El objetivo es obtener el triángulo de Sierpinski y graficarlo utilizando la función de matplolib siguiente:

mpatches.Polygon(T[i], fc="y")

Saludos!

Gaby Du
  • 273
  • 1
  • 3
  • 10