一维伊辛模型MATLAB(总2页).docx

一维伊辛模型MATLAB源代码clear;clc;%本次模拟中，J=1，k=1kT=0.1:0.1:5;kT_len=length(kT);num=102;Ising=sign(2*rand(1,num)-1);%初始构型N=2*105;%-初始化，预分配内存-%E_total=zeros(1,kT_len);E_squ=zeros(1,kT_len);C=zeros(1,kT_len);M_total=zeros(1,kT_len);for i=1:kT_len E_sum=0; E_squ_sum=0; M_sum=0; for turn=1:N E_old=E_calc(Ising); ri=randi(1,num,1); if ri=1 ri=num-1; end if ri=num ri=2;