用Matlab模拟二维正方形晶格上自旋的随机热力学行为

该代码使用Metropolis算法模拟Ising模型的热力学演化，每次迭代随机选取一个晶格点，计算翻转该点自旋带来的能量变化，并根据Metropolis准则以一定概率接受或拒绝该翻转操作。迭代次数越多，晶格达到的热力学平衡态越稳定，所显示的图像也越接近于该温度下的平衡态自旋构型。初始时，晶格中的自旋是随机排列的。在迭代中，用函数randi()随机选取一个晶格点，计算翻转该点自旋带来的能量变化$\Delta E$，并根据Metropolis准则接受或拒绝该翻转操作。最后通过imagesc()函数将晶格的自旋状态可视化。