C#遗传算法学习笔记

• 本文介绍C#遗传算法学习笔记，通过运行程序，你会发现通过不断的进化，种群的总的适应环境的能力在逐步提高。

• 以下代码实现了C#遗传算法一个简单的花朵进化的模拟过程。

  花朵的种群数量是10，共进化了50代。通过运行程序，你会发现通过不断的进化，种群的总的适应环境的能力在逐步提高（fitness的值下降）。

  C#遗传算法实现代码：

    
    
  using System;  
  using System.Collections.Generic;  
  using System.Text;  
  namespace GA  
  {  
  class Program  
  {  
  static void Main(string[] args)  
  {  
  World world = new World();  
  world.Init();  
  for (int i = 0; i < 50; i++)  
  {  
  world.Evolve();  
  Console.WriteLine(i);  
  world.Show();  
  }  
  }  
  }  
   
  class World  
  {  
  int kMaxFlowers = 11;  
  Random Rnd = new Random();  
  public int[] temperature;  
  public int[] water;  
  public int[] sunlight;  
  public int[] nutrient;  
  public int[] beneficialInsect;  
  public int[] harmfulInsect;  
  public int currentTemperature;  
  public int currentWater;  
  public int currentSunlight;  
  public int currentNutrient;  
  public int currentBeneficialInsect;  
  public int currentHarmfulInsect;  
  public World()  
  {  
  temperature = new int[kMaxFlowers];  
  water = new int[kMaxFlowers];  
  sunlight = new int[kMaxFlowers];  
  nutrient = new int[kMaxFlowers];  
  beneficialInsect = new int[kMaxFlowers];  
  harmfulInsect = new int[kMaxFlowers];  
  }  
  /** 
   
  /// 初始化第一代花朵的基因结构  
  /// 
   
  public void Init()  
  {  
  for (int i = 1; i < kMaxFlowers; i++)  
  {  
  temperature[i] = Rnd.Next(1, 75);  
  water[i] = Rnd.Next(1, 75);  
  sunlight[i] = Rnd.Next(1, 75);  
  nutrient[i] = Rnd.Next(1, 75);  
  beneficialInsect[i] = Rnd.Next(1, 75);  
  harmfulInsect[i] = Rnd.Next(1, 75);  
  }  
  currentTemperature = Rnd.Next(1, 75);  
  currentWater = Rnd.Next(1, 75);  
  currentSunlight = Rnd.Next(1, 75);  
  currentNutrient = Rnd.Next(1, 75);  
  currentBeneficialInsect = Rnd.Next(1, 75);  
  currentHarmfulInsect = Rnd.Next(1, 75);  
  }  
  /** 
   
  /// 越大说明花朵的适应环境的能力差，小说明适应环境的能力强  
  /// 
   
  ///  name="flower"> 
  /// 
       
        
  private int Fitness(int flower)  
  {  
  int theFitness = 0;  
  theFitness = Math.Abs(temperature[flower] - currentTemperature);  
  theFitnesstheFitness = theFitness + Math.Abs(water[flower] - currentWater);  
  theFitnesstheFitness = theFitness + Math.Abs(sunlight[flower] -  
  currentSunlight);  
  theFitnesstheFitness = theFitness + Math.Abs(nutrient[flower] -  
  currentNutrient);  
  theFitnesstheFitness = theFitness + Math.Abs(beneficialInsect[flower] -  
  currentBeneficialInsect);  
  theFitnesstheFitness = theFitness + Math.Abs(harmfulInsect[flower] -  
  currentHarmfulInsect);  
  return (theFitness);  
  }  
  /** 
   
  /// 排除适应能力差的花朵，让适应能力强的花朵杂交繁殖，产生下一代。同时有一定的概率变异。  
  /// 
   
  public void Evolve()  
  {  
  int[] fitTemperature = new int[kMaxFlowers];  
  int[] fitWater = new int[kMaxFlowers];  
  int[] fitSunlight = new int[kMaxFlowers];  
  int[] fitNutrient = new int[kMaxFlowers];  
  int[] fitBeneficialInsect = new int[kMaxFlowers];  
  int[] fitHarmfulInsect = new int[kMaxFlowers];  
  int[] fitness = new int[kMaxFlowers];  
  int i;  
  int leastFit = 0;  
  int leastFitIndex = 1;  
  for (i = 1; i < kMaxFlowers; i++)  
  if (Fitness(i) > leastFit)  
  {  
  leastFit = Fitness(i);  
  leastFitIndex = i;  
  }  
  temperature[leastFitIndex] = temperature[Rnd.Next(1, 10)];  
  water[leastFitIndex] = water[Rnd.Next(1, 10)];  
  sunlight[leastFitIndex] = sunlight[Rnd.Next(1, 10)];  
  nutrient[leastFitIndex] = nutrient[Rnd.Next(1, 10)];  
  beneficialInsect[leastFitIndex] = beneficialInsect[Rnd.Next(1, 10)];  
  harmfulInsect[leastFitIndex] = harmfulInsect[Rnd.Next(1, 10)];  
  for (i = 1; i < kMaxFlowers; i++)  
  {  
  fitTemperature[i] = temperature[Rnd.Next(1, 10)];  
  fitWater[i] = water[Rnd.Next(1, 10)];  
  fitSunlight[i] = sunlight[Rnd.Next(1, 10)];  
  fitNutrient[i] = nutrient[Rnd.Next(1, 10)];  
  fitBeneficialInsect[i] = beneficialInsect[Rnd.Next(1, 10)];  
  fitHarmfulInsect[i] = harmfulInsect[Rnd.Next(1, 10)];  
  }  
  for (i = 1; i < kMaxFlowers; i++)  
  {  
  temperature[i] = fitTemperature[i];  
  water[i] = fitWater[i];  
  sunlight[i] = fitSunlight[i];  
  nutrient[i] = fitNutrient[i];  
  beneficialInsect[i] = fitBeneficialInsect[i];  
  harmfulInsect[i] = fitHarmfulInsect[i];  
  }  
  for (i = 1; i < kMaxFlowers; i++)  
  {  
  if (Rnd.Next(1, 100) == 1)  
  temperature[i] = Rnd.Next(1, 75);  
  if (Rnd.Next(1, 100) == 1)  
  water[i] = Rnd.Next(1, 75);  
  if (Rnd.Next(1, 100) == 1)  
  sunlight[i] = Rnd.Next(1, 75);  
  if (Rnd.Next(1, 100) == 1)  
  nutrient[i] = Rnd.Next(1, 75);  
  if (Rnd.Next(1, 100) == 1)  
  beneficialInsect[i] = Rnd.Next(1, 75);  
  if (Rnd.Next(1, 100) == 1)  
  harmfulInsect[i] = Rnd.Next(1, 75);  
  }  
  }  
  /** 
   
  /// 显示种群中个体对环境的适应能力，还有所有个体对环境的适应能力之和。  
  /// 
   
  public void Show()  
  {  
  int sum = 0;  
  for (int i = 1; i < kMaxFlowers; i++)  
  {  
  int fitness = Fitness(i);  
  sum += fitness;  
  Console.WriteLine("No." + i + "'s fitness is " + fitness);  
  }  
  Console.WriteLine("fitness sum is " + sum);  
  }  
  }  
  } 

转载 吾搜 空间