模拟 Ising 模型的 Matlab 代码，用于模拟一个二维正方形晶格上自旋的随机热力学行为

这是一个模拟 Ising 模型的 Matlab 代码，用于模拟一个二维正方形晶格上自旋的随机热力学行为，其中自旋是取值为+1或-1的变量。该模型在固体物理和统计物理中得到了广泛的应用，它的主要特点是磁性自组装、相变以及相变的临界行为。该代码使用 Metropolis 算法模拟 Ising 模型的热力学演化，每次迭代随机选取一个晶格点，计算翻转该点自旋带来的能量变化，并根据 Metropolis 准则以一定概率接受或拒绝该翻转操作。迭代次数越多，晶格达到的热力学平衡态越稳定，所显示的图像也越接近于该温度下的平衡态自旋构型。初始时，晶格中的自旋是随机排列的。在迭代中，用函数 randi() 随机选取一个晶格点，然后计算翻转该点自旋带来的能量变化 $\Delta E$，并根据 Metropolis 准则接受或拒绝该翻转操作。最后通过 imagesc() 函数将晶格的自旋状态可视化。