伊辛模型MATLAB源代码(共3页).docx

精选优质文档-倾情为你奉上伊辛(Ising)模型MATLAB源程序代码clear all;close all;clcN=16;Te=0:0.2:5;imax=;MT=zeros(1,length(Te);Cv=zeros(1,length(Te);xx=zeros(1,length(Te);L=1;k=1;for T=TeMt=zeros(1,imax);step=;spin=round(rand(N).*2-1;subplot(3,2,1)gspin=pcolor(spin);colormap(gray(2);title(自旋分布情况);Mt=zeros(1,imax);subplot(3,2,2)gMt=plot(1:imax,Mt);axis(0 imax -1.1 1.1);title(固定温度下的磁矩);%xlabel(时间);%ylabel(Mt);%initial energy idx=repmat(: