这是一个用Matlab编写的Ising模型代码，用于模拟一个二维正方形晶格上自旋的随机热力学行为。其中，自旋的取值为+1或-1。Ising模型在固体物理和统计物理中被广泛应用，其特点为磁性自组装、相变以及相变的临界行为。该代码使用Metropolis算法模拟Ising模型的热力学演化，迭代计算晶格点的自旋变化，并根据Metropolis准则接受或拒绝该变化操作。最后，可以使用imagesc()函数将晶格的自旋状态可视化。该代码适用于物理学和计算机科学领域的研究和教学。