function [ S,E,Tk ] = c_annealing( S,Wij,Tk,c )%C_ANNEALING p291確定性模擬退火算法%輸入參數：%         S：網絡節點的狀態序列，每個節點可取值+1或-1%       Wij：節點間的權值矩陣%       Tk：溫度%        c：T(k+1)=c*T(k)undone=true;n=size(S,2);%有n個節點k=1;k_max=800;while undone    i=round(rand*n);    if i==0        i=1;    end    li=Wij(i,:)*S';    liT=li/Tk;    %S(i)=round((1-exp(-liT))/(1+exp(-liT)));    S(i)=1/(1+exp(-liT));    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');