function [ S,E,Tk ] = annealing( S,Wij,Tk,c )%ANNEALING p288隨機模擬退火算法%輸入參數：%         S：網絡節點的狀態序列，每個節點可取值+1或-1%       Wij：節點間的權值矩陣%       Tk：溫度%        c：T(k+1)=c*T(k)n=size(S,2);%有n個節點undone=true;k=1;k_max=800;while undone    for cyc=1:n*4 %所有節點輪詢四次        i=round(rand*n);        if i==0            i=1;        end        Ea=-0.5*S(i)*Wij(i,:)*S';        Eb=-Ea;        if Eb<Ea            S(i)=-S(i);        elseif exp(-(Eb-Ea)/Tk)>rand            S(i)=-S(i);        end    end    if k==k_max        undone=false;    else        k=k+1;        Tk=c*Tk;    endendfor i=1:n    Wij(i,:)=S(i)*Wij(i,:);endE=-0.5*sum(Wij*S');