Click here to Skip to main content
11,922,559 members (32,490 online)
Click here to Skip to main content
Add your own
alternative version


35 bookmarked

Simulated Annealing Example in C#

, 30 Jun 2006 CPOL
Rate this:
Please Sign up or sign in to vote.
Artificial intelligence algorithm: simulated annealing


Simulated annealing (SA) is an AI algorithm that starts with some solution that is totally random, and changes it to another solution that is “similar” to the previous one. It makes slight changes to the result until it reaches a result close to the optimal.

Simulated annealing is a stochastic algorithm, meaning that it uses random numbers in its execution. So every time you run the program, you might come up with a different result. It produces a sequence of solutions, each one derived by slightly altering the previous one, or by rejecting a new solution and falling back to the previous one without any change.

When SA starts, it alters the previous solution even if it is worse than the previous one. However, the probability with which it will accept a worse solution decreases with time,(cooling process) and with the “distance” the new (worse) solution is from the old one. It always accepts a new solution if it is better than the previous one.

The probability used is derived from The Maxwell-Boltzmann distribution which is the classical distribution function for distribution of an amount of energy between identical but distinguishable particles. It's value is:


Besides the presumption of distinguishability, classical statistical physics postulates further that:

  • There is no restriction on the number of particles which can occupy a given state.
  • At thermal equilibrium, the distribution of particles among the available energy states will take the most probable distribution consistent with the total available energy and total number of particles.
  • Every specific state of the system has equal probability.

The name “simulated annealing” is derived from the physical heating of a material like steel. This material is subjected to high temperature and then gradually cooled. The gradual cooling allows the material to cool to a state in which there are few weak points. It achieves a kind of “global optimum” wherein the entire object achieves a minimum energy crystalline structure.

If the material is rapidly cooled, some parts of the object, the object is easily broken (areas of high energy structure). The object has achieved some local areas of optimal strength, but is not strong throughout, with rapid cooling.

In my program, I took the example of the travelling salesman problem: file tsp.txt.The matrix designates the total distance from one city to another (nb: diagonal is 0 since the distance of a city to itself is 0).

As for the program, I tried developing it as simple as possible to be understandable.

You could change the starting temperature, decrease or increase epsilon (the amount of temperature that is cooling off) and alter alpha to observe the algorithm's performance.

The program calculates the minimum distance to reach all cities(TSP). The best minimal distance I got so far using that algorithm was 17. Can you calculate a better distance?

The Code

public string StartAnnealing()
    ArrayList list = new ArrayList();
    //primary configuration of cities
    int [] current={0,1,2,3,4,5,6,7,8,9,10,11,12,13,14};
    //the next configuration of cities to be tested
    int []next=new int[15];
    int iteration =-1;
    //the probability
    double proba;
    double alpha =0.999;
    double temperature = 400.0;
    double epsilon = 0.001;
    double delta;
    double distance = TspDataReader.computeDistance(current);

    //while the temperature did not reach epsilon
    while(temperature > epsilon)
        //get the next random permutation of distances 
        //compute the distance of the new permuted configuration
        delta = TspDataReader.computeDistance(next)-distance;
        //if the new distance is better accept it and assign it
            distance = delta+distance;
            proba = rnd.Next();
            //if the new distance is worse accept 
            //it but with a probability level
            //if the probability is less than 
            //E to the power -delta/temperature.
            //otherwise the old value is kept
            if(proba< Math.Exp(-delta/temperature))
                distance = delta+distance;
        //cooling process on every iteration
        temperature *=alpha;
        //print every 400 iterations
        if (iteration%400==0)
        return "best distance is"+distance;
        return "error";

/// <summary>
/// compute a new next configuration
/// and save the old next as current
/// </summary>
/// <param name="c">current configuration</param>
/// <param name="n">next configuration</param>
void computeNext(int[] c, int[] n)
    for(int i=0;i<c.Length;i++)
    int i1 = (int)(rnd.Next(14))+1;
    int i2 = (int)(rnd.Next(14))+1;
    int aux = n[i1];

Make sure the debug window is opened to observe the algorithm's behavior through iterations.

Happy programming!


This article, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)


About the Author

Assaad Chalhoub
Lebanon Lebanon
No Biography provided

You may also be interested in...

Comments and Discussions

Generalpoor explanation Pin
Omar Gameel Salem15-Sep-14 15:24
memberOmar Gameel Salem15-Sep-14 15:24 
GeneralThanks - and some information Pin
mike griggs22-Nov-09 11:33
membermike griggs22-Nov-09 11:33 
GeneralRe: Thanks - and some information Pin
peterkjgaard27-Apr-12 2:57
memberpeterkjgaard27-Apr-12 2:57 
QuestionIs this a Bug? Pin
Member 897296-Nov-09 12:18
memberMember 897296-Nov-09 12:18 
AnswerRe: Is this a Bug? Pin
MikeDC22-May-12 15:54
memberMikeDC22-May-12 15:54 
Generaldelta problem PinPopular
alitozan16-Nov-06 1:14
memberalitozan16-Nov-06 1:14 
GeneralNice summary and concise explanations PinPopular
Hal Angseesing30-Jun-06 2:13
memberHal Angseesing30-Jun-06 2:13 
GeneralRe: Nice summary and concise explanations Pin
Assaad Chalhoub2-Jul-06 21:48
memberAssaad Chalhoub2-Jul-06 21:48 

General General    News News    Suggestion Suggestion    Question Question    Bug Bug    Answer Answer    Joke Joke    Praise Praise    Rant Rant    Admin Admin   

Use Ctrl+Left/Right to switch messages, Ctrl+Up/Down to switch threads, Ctrl+Shift+Left/Right to switch pages.

| Advertise | Privacy | Terms of Use | Mobile
Web02 | 2.8.151120.1 | Last Updated 30 Jun 2006
Article Copyright 2006 by Assaad Chalhoub
Everything else Copyright © CodeProject, 1999-2015
Layout: fixed | fluid