【随机模拟退火算法】matlab实现

《模式分类》第二版课本上机练习p314 第二题

T=10，c=0.9

%初始化六个节点的状态
s=rand(1,6);
s(s>0.5)=1;s(s<0.5)=-1;
s=s';
disp(['初始构型s为：',num2str(s')]);%num2str函数讲数值转化为字符串


T=10;%初始温度
w=[0 5 -3 4 4 1;5 0 -1 2 -3 1;-3 -1 0 2 2 0;4 2 2 0 3 -3;4 -3 2 3 0 5;1 1 0 -3 5 0];
%权值矩阵
c=0.9;%温度下降速度
e=10^(-5);%阈值
k=1;
kmax=100;%最大迭代次数
Ea= inf;
Eb= inf;
E=0;
while k<=kmax
    for t=1:1000
        i=randi([1,6]);
        Ea=(-0.5)*s(i)*(w(i,:)*s);%计算si不变的系统总能量
        Eb=-Ea;%改变i的状态计算总能量
        if Eb<Ea
            s(i)=-s(i);%以1的概率接受改变
        elseif exp(-(Eb-Ea)/T)>rand(1)
            s(i)=-s(i);%以p的概率接受改变
        else
            s(i)=s(i);
        end
    end
%计算总能量
for it=1:6
    for jt=1:6
        E=E+w(it,jt)*s(it)*s(jt);
    end
end
E=(-0.5)*E;
    k=k+1;
    T=c*T;%温度更新
end
disp(['随机模拟退火算法最终能量为：',num2str(E)]);
disp(['随机模拟退火算法最终构型为：',num2str(s')]);

第一次随机结果

 第二次随机结果