Comme l'a souligné Thomas Klimpel dans les commentaires, une certaine probabilité d'acceptation est souvent utilisée, qui est égale à . Voici une méthode itérative simple pour trouver une température initiale appropriée, proposée par Ben-Ameur en 2004 [1]. Dans ce qui suit, est une transition strictement positive, et sont les états après et avant la transition, la différence de coût et la probabilité de générer une transition lorsque les états d'énergie sont distribués conformément à la distribution stationnaire0.8tmin t δ t E max t - E min t π min t 1maxtmintδtEmaxt−Eminttπmint1|N(mint)|t
N(i)i
πi=|N(i)|exp(−Ei/T)∑j|N(j)|exp(−Ej/T)
, où désigne l'ensemble des voisins de .
N(i)i
Enfin, est la probabilité d'accepter une transition positive . Maintenant, nous pouvons avoir une estimation de la probabilité d'acceptation basée sur un ensemble "aléatoire" de transitions positives:t χ χ ( T ) Sexp ( - δt/ T)tχ^χ ( T)S
χ^( T)==∑t ∈ Sπmint1| N( mint) |exp ( - δt/ T)∑t ∈ Sπmint1| N( mint) |∑t ∈ Sexp ( - Emaxt/ T)∑t ∈ Sexp ( - Emint/ T).
Nous voulons trouver une température telle que , où est la probabilité d'acceptation que nous désirons. χ ( T 0 ) = χ 0 χ 0 ∈ ] 0 , 1 [T0χ ( T0) = χ0χ0∈]0,1[
S E max t E min t S T 1 T 0T0 est calculé par une méthode itérative. Certains états et un voisin pour chaque état sont générés. Cela nous donne un ensemble de transitions . Les énergies et correspondant aux états du sous-ensemble sont stockées. Ensuite, une valeur pour est choisie, qui peut être n'importe quelle valeur positive. est alors trouvé avec la formule récursiveSEmaxtEmintST1T0
Tn+1=Tnln(χ^(Tn))ln(χ0)1/p
, où est un nombre réel .
p≥1
Lorsque se rapproche de nous pouvons arrêter. est maintenant une bonne approximation de la température initiale souhaitée . Pour plus d'explications, de preuves et de discussions, veuillez consulter la première section de l'article original [1].χ^(Tn)χ0TnT0
[1] Ben-Ameur, Walid. "Calcul de la température initiale du recuit simulé." Optimisation informatique et applications 29, no. 3 (2004): 369-385.