solvebeta.m

用Monte-Carlo方法解决二维Ising模型

%估计M*M格点的临界beta
function beta = solvebeta(M)
%调用格式: beta=solvebata(M)
%自变量:   M - 二维Ising模型的边长
%输出:    beta - 临界温度值


%请选择合适的策略,否则程序运行非常缓慢!
t=20; %重复次数
p=50; %beta个数
N=500;N1=500;N2=3;number=1; %组策略



P=zeros(t,1);
for k=1:t
    [U,C]=Ising(M,p,N,N1,N2,number,2,0);
    P(k)=(maxN(C)-1)/p;
end
beta=drawimage(P);
end

function P=maxN(C)
n=length(C);
P=1;
for j=2:n
    if C(j)>C(P);
        P=j;
    end
end
end