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!