人工智能实验报告
学 号:姓 名:实验名称:实验日期: 遗传算法 2016.1.5
【实验名称】 遗传算法 【实验目的】
掌握遗传算法的基本原理,熟悉遗传算法的运行机制,学会用遗传算法来求解问题。
【实验原理】
遗传算法(Genetic Algorithm)是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型,是一种通过模拟自然进化过程搜索最优解的方法。
遗传算法是从代表问题可能潜在的解集的一个种群开始的,而一个种群则由经过基因编码的一定数目的个体组成。每个个体实际上是染色体带有特征的实体。在一开始需要实现从表现型到基因型的映射即编码工作。由于仿照基因编码的工作很复杂,我们往往进行简化,如二进制编码,初代种群产生之后,按照适者生存和优胜劣汰的原理,逐代演化产生出越来越好的近似解,在每一代,根据问题域中个体的适应度大小选择个体,并借助于自然遗传学的遗传算子进行组合交叉和变异,产生出代表新的解集的种群。这个过程将导致种群像自然进化一样的后生代种群比前代更加适应于环境,末代种群中的最优个体经过解码,可以作为问题近似最优解。
遗传算法程度流程图为:
【实验内容】
题目:已知f(x)=x*sin(x)+1,x?[0,2?],求f(x)的最大值和最小值。
数据结构:
struct poptype {
double gene[length];//染色体 double realnumber;//对应的实数x double fitness;//适应度
double rfitness;//相对适应度 double cfitness;//累计适应度 };
struct poptype population[popsize+1];//最后一位存放max/min struct poptype newpopulation[popsize+1];//
染色体编码:
x?[0,2?],变量长度为2 π,取小数点后6位,由于222?2?*106?223;
因此,染色体由23位字节的二进制矢量表示,则X与二进制串(
?b2221b......b0?22??22? ???bi?2i??x';x?x'g232?1?i?0?10 适应度函数:
由于要求f(x)的最值,所以适应度函数即可为f(x)。但为了确保在轮赌法选择过中,每个个体都有被选中的可能性,因此需要将所有适应度调整为大于0的值。因此,设计求最大值的适应度函数如下:
evalmax?f(x)?5?xsinx?6;
将最小问题转化为求-f(x)的最大值,同理,设计最小值的适应度函数如下:
evalmin??f(x)?5??xsinx?4;
种群大小:
本实验默认为50,再进行种群初始化。
实验参数:
主要有迭代数,交叉概率,变异概率这三个参数。一般交叉概率在0.6-0.9范围内,变异概率在0.01-0.1范围内。可以通过手动输入进行调试。
主要代码如下:
void initialize()//种群初始化 {
srand(time(NULL));
int i,j;
for(i=0;i for(i=0;i<=popsize+1;i++) { } } void cal_fitness()//计算适应度 { int i; for(i=0;i void select()//选择操作 { int mem,i,j,k; double sum=0; double p; for (mem=0;mem sum+=population[mem].fitness; population[mem].rfitness=population[mem].fitness/sum; for (mem=0;mem population[0].cfitness=population[0].rfitness; for (mem=1;mem population[mem].cfitness=population[mem-1].cfitness+population[mem].rfitness; for (i=0;i p=rand()00/1000.0; if (p newpopulation[i]=population[0]; for (j=0;j if (p>=population[j].cfitness&&p newpopulation[i]=population[j+1]; else{ population[i].fitness=population[i].realnumber*sin(population[i].realnumber)+6; population[i].realnumber=0; for(j=0;j<23;j++) population[i].realnumber+=population[i].gene[j]*pow(2,j); for(j=0;j<23;j++) population[i].gene[j]=rand()%2; void transform()//染色体转化为实数x population[i].realnumber=population[i].realnumber*2*pi/(pow(2,23)-1); } for (i=0;i void cross()//交叉操作 { int i, mem, one; int first = 0; double x; for(mem=0;mem x = rand()00/1000.0; if (x ++first; if (first%2==0) Xover(one,mem);//个体间染色体进行交叉函数 else one=mem; } } } void mutate()//变异操作 { int i, j,t; double x; for (i=0;i for(j=0;j x=rand()00/1000.0; if (x if(population[i].gene[j]) population[i].gene[j]=0; else population[i].gene[j]=1; } } } void cal_max()//计算最大值 { int i; double max,sum=0; int max_m; max=population[0].fitness; for(i=0;i if(population[i].fitness>population[i+1].fitness) if(population[i].fitness>=max)