Clase 5 PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 2

Dimensión Fractal

Felipe Arancibia , Carlos Muñoz


September 2019

1 El algoritmo de conteo de cajas


from random import *
import matplotlib.pyplot as plt
import numpy as np

X = []
Y = []
r = 0.3
N = 100
M = np.zeros((N,N))

for k in range(1000):
x = random()
for m in range(1000):
y = random()
if (x-0.5)**2+(y-0.5)**2<=r**2:
X.append(x)
Y.append(y)

for i in range(N):
for j in range(N):
if M[i][j]==1:
break
else:
for x in X:
for y in Y:
if 1/5 + (0.3/N)*i < x and x <= 1/5 + (0.3/N)*(i+1):
if 1/5 + (0.3/N)*j < y and y <= 1/5 + (0.3/N)*(j+1):
M[i][j] = 1
elif 1/2 + (0.3/N)*i < x and x <= 1/2 + (0.3/N)*(i+1):
if 1/2 + (0.3/N)*j < y and y <= 1/2 + (0.3/N)*(j+1):
M[i][j] = 1

1
plt.imshow(M,cmap='hot')
plt.show

Debiamos estimar la dimensión fractal en dos dimensiones de un circulo de radio


0.3 centrado en las coordenadas (1/2, 1/2). Para esto generamos dos puntos al
azar simultaneamente, los cuales debı́an cumplir que la raiz de la suma de sus
cuadrados fuera menor o igual a 0.3, los pares de datos que cumplian la condición
eran graficados mediante la función scatter de matplotlib generando el circulo
relleno solicitado.

2 AQUÍ LA EXPLICACIÓN DEL PROGRAMA


DE EL SEGUNDO FOR
Intentamos correr el programa en la plataforma de Cocalc y en el programa de
Python 3 de 64 bits instalado en Windows 10, pero en ambos casos no paraba
de calcular; incluso disminuyendo las iteraciones, por lo que no fue posible ver
el resultado del programa. Sin embargo, usted revisó nuestro programa en clase
y comentó que debı́a funcionar, por lo que asumimos que con la capacidad de
computo necesaria nuestro programa entregarı́a el resultado requerido.

También podría gustarte