伊辛(ising)模型MATLAB源程序

clear all；
close all；
clc
N=16，
Te=[0：0.2：5]；
imax=500000；
MT=zeros(1， length(Te) ) ；
Cv=zeros(1， length(Te) ) ；
xx=zeros(1， length(Te) ) ；
L=1k=1
forT=Te
Mt=zeros(1， imax)
step= 100000
spin=round(rand(N) ) .*2-1；
subplot(3.2.1)
gsp in=p color(spin) ：
colormap(gray(2) )
title(自旋分布情况) ；
Mt=zeros(1， imax) ；
subplot(3， 2.2)
gMt=plot([1：imax] .Mt) ：
22.a0.00
axis([O imax-1.11.1] )
title('固定温度下的磁矩) ；
%xlabel('时间) ；
%ylabel(‘Mt”) ；
%initial energy
idx=repmat((：) ， ndims(spin) ， 1) ；
idx(1) =[N1：N-1] ；
E total=spin(idx() )