% 結(jié)束條件為兩次差小于一個(gè)特定量

% MarkoveLength 馬爾科夫鏈長(zhǎng)度

% DecayScale 溫度衰減參數(shù)

MarkoveLength=1000;

DecayScale=0.9;

Temperature=1000; 

% initial temperaturet = 1;

% final temperature%指派矩陣，每個(gè)人做每件事所需要的費(fèi)用

assingnMatrix=[17,1,9,15,16;8,15,17,9,11;5,8,4,18,12;20,6,14,9,7;17,14,2,10,1];

% 初始解x=[1,2,3,4,5];

totalCost = 0;

for i=1:5   

totalCost=totalCost+assingnMatrix(i,x(i));

end

BestCost=totalCost;

BestAssign=x;while (Temperature > t || abs(deltaCost)<2)    

    for i = 1:MarkoveLength       

       r = randperm(5,2); 

%產(chǎn)生兩個(gè)隨機(jī)數(shù)，用來交換x中的任務(wù)分配順序        

        r1=r(1);    

    r2=r(2);    

    temp=x(r1);    

    x(r1)=x(r2);    

    x(r2)=temp;        

        %  計(jì)算        

        totalCost=0;        

        for k=1:4            

            totalCost=totalCost+assingnMatrix(k,x(k));        

        end        % 判斷費(fèi)用是否優(yōu)化       

       deltaCost=totalCost-BestCost;       

        if deltaCost <= 0            

            BestCost = totalCost; 

            % 若費(fèi)用減少則無條件接受           

            BestAssign=x;        

        else            

            if (rand()<exp(-deltaCost/Temperature))               

                BestCost = totalCost;  %否則在一定概率下接受              

                BestAssign=x;           

            end       

        end   

    end  

    Temperature=Temperature*DecayScale;

enddisp(BestCost);

disp(assingnMatrix);

disp(BestAssign);