我正在尝试使用Python编写Ising模型
我想,我已经正确编码了,但是我在动画或绘图方面有问题。我似乎为每个配置绘制了一个新的映像,而不是更新现有的映像,从而生成了许多我不需要的已保存映像。我只想要一个正在更新的绘图,如果可能的话
我知道,我在循环中绘制,但我不记得这是一个问题,当我想绘制每个迭代时。Seaborn的热图会有问题吗
我已附上我的代码:
import numpy as np
import numpy.random as npr
import matplotlib.pyplot as plt
import seaborn as sns
#Constants
J = 1
h = 1
kbT = 1
beta = 1
#Grid
L = 20 #Dimensions
N = L**2 #Total number of grid points
#Initial configuration
spins = 2*np.random.randint(2, size = (L,L))-1
E = []
i = 0
plt.figure()
while i < 100000:
for i in range(1,N):
i += 1
s = tuple(npr.randint(0, L, 2)) # Random initial coordinate
# x and y coordinate
(sx, sy) = s
# Periodic boundary condition
sl = (sx-1, sy)
sr = ((sx+1)%L, sy)
sb = (sx, sy-1)
st = (sx, (sy+1)%L)
# Energy
E = spins[s] * ( spins[sl] + spins[sr] + spins[sb] + spins[st] )
if E <= 0 : # If negative, flip
spins[s] *= -1
else:
x = np.exp(-E/kbT) # If positve, check condition
q = npr.rand()
if x > q:
spins[s] *= -1
# Plot (heatmap)
sns.heatmap(spins, cmap = 'magma')
plt.pause(10e-10)
plt.draw()
plt.show()
我认为函数
ion
和clf
可以做到这一点使用函数
ion
可以使绘图具有交互性,因此需要:Here函数的引用
clf
的参考是here相关问题 更多 >
编程相关推荐