import numpy as np
import matplotlib.pyplot as plt
# Número de clusters fractais
N = 50
# Profundidade fractal
K = 6
# Tempos de coerência locais (grande estabilidade)
T = np.linspace(50, 200, N)
# Fases iniciais aleatórias
phi0 = np.random.uniform(0, 2*np.pi, N)
# Coeficientes fractais
alpha = np.random.uniform(0.05, 0.2, (N, K))
# Distâncias fractais entre clusters
d = np.random.uniform(0.5, 2.0, (N, N))
# Tempo de simulação
t = np.linspace(0, 50, 2000)
def fase_fractal(i, t):
"""Fase fractal temporal de um cluster."""
return phi0[i] + np.sum([
alpha[i, k] * np.sin((2**k) * t)
for k in range(K)
], axis=0)
def coerencia_global(t):
"""Coerência fractal distribuída."""
C = np.zeros_like(t)
for i in range(N):
phi = fase_fractal(i, t)
C += np.exp(-t / T[i]) * np.cos(phi)
return C / N
def interacao_fractal(t):
"""Interação fractal total entre clusters."""
I = np.zeros_like(t)
for i in range(N):
phi_i = fase_fractal(i, t)
for j in range(i+1, N):
phi_j = fase_fractal(j, t)
I += (np.cos(phi_i - phi_j)) / (d[i, j]**2)
return I
# Calcular coerência e interação
C = coerencia_global(t)
I = interacao_fractal(t)
# Plot
plt.figure(figsize=(12, 6))
plt.plot(t, C, label="Coerência Fractal")
plt.plot(t, I / np.max(I), label="Interação Fractal (normalizada)")
plt.title("Decoerência Fractal Quântica")
plt.xlabel("Tempo")
plt.ylabel("Amplitude")
plt.legend()
plt.grid(True)
plt.show()







No comments:
Post a Comment