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];  
}  
/**//// <summary> 
/// 初始化***代花朵的基因结构  
/// </summary> 
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);  
}  
/**//// <summary> 
/// 越大说明花朵的适应环境的能力差，小说明适应环境的能力强  
/// </summary> 
/// <param name="flower"></param> 
/// <returns></returns> 
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);  
}  
/**//// <summary> 
/// 排除适应能力差的花朵，让适应能力强的花朵杂交繁殖，产生下一代。同时有一定的概率变异。  
/// </summary> 
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);  
}  
}  
/**//// <summary> 
/// 显示种群中个体对环境的适应能力，还有所有个体对环境的适应能力之和。  
/// </summary> 
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);  
}  
}  
} 
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
61.
62.
63.
64.
65.
66.
67.
68.
69.
70.
71.
72.
73.
74.
75.
76.
77.
78.
79.
80.
81.
82.
83.
84.
85.
86.
87.
88.
89.
90.
91.
92.
93.
94.
95.
96.
97.
98.
99.
100.
101.
102.
103.
104.
105.
106.
107.
108.
109.
110.
111.
112.
113.
114.
115.
116.
117.
118.
119.
120.
121.
122.
123.
124.
125.
126.
127.
128.
129.
130.
131.
132.
133.
134.
135.
136.
137.
138.
139.
140.
141.
142.
143.
144.
145.
146.
147.
148.
149.
150.
151.
152.
153.
154.
155.
156.
157.
158.
159.
160.
161.
162.
163.