人工智能遗传算法试验报告

1、人工智能实验报告学号：姓名：实验名称：遗传算法实验日期：2016.1.5实验名称遗传算法【实验目的】掌握遗传算法的基本原理，熟悉遗传算法的运行机制，学会用遗传算法来求解问题。【实验原理】遗传算法(GeneticAlgorithm)是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型，是一种通过模拟自然进化过程搜索最优解的方法。遗传算法是从代表问题可能潜在的解集的一个种群开始的，而一个种群则由经过基因编码的一定数目的个体组成。每个个体实际上是染色体带有特征的实体。在一开始需要实现从表现型到基因型的映射即编码工作。由于仿照基因编码的工作很复杂，我们往往进行简化，如二进制编码，初代种

2、群产生之后，按照适者生存和优胜劣汰的原理，逐代演化产生出越来越好的近似解，在每一代，根据问题域中个体的适应度大小选择个体，并借助于自然遗传学的遗传算子进行组合交叉和变异，产生出代表新的解集的种群。这个过程将导致种群像自然进化一样的后生代种群比前代更加适应于环境，末代种群中的最优个体经过解码，可以作为问题近似最优解。遗传算法程度流程图为：【实验内容】题目：已知f(x)=x*sin(x)+1,xw0,2力求f(x)的最大值和最小值。数据结构：structpoptype(doublegenelength;/染色体doublerealnumber;/对应的实数xdoublefitness;/适应度do

3、ublerfitness;/相对适应度doublecfitness;/累计适应度;structpoptypepopulationpopsize+1;/最后一位存放max/minstructpoptypenewpopulationpopsize+1;/染色体编码:xW0,2n,变量长度为2兀，取小数点后6位，由于222<2n*106<223;因此，染色体由23位字节的二进制矢量表示，则X与二进制串(<b22b21b°>)2之间的映射如下：/4产i%2n(b22b21b0)bi*2=x'x=x'*3-Iiz0J102-1适应度函数：由于要求f(x)

4、的最值，所以适应度函数即可为f(x)。但为了确保在轮赌法选择过中，每个个体都有被选中的可能性，因此需要将所有适应度调整为大于0的值。因此，设计求最大值的适应度函数如下：evalmax=f(x)5=xsinx6;将最小问题转化为求-f(x)的最大值，同理，设计最小值的适应度函数如下：evalmin=-f(x)5=-xsinx4;种群大小：本实验默认为50,再进行种群初始化。实验参数：主要有迭代数，交叉概率，变异概率这三个参数。一般交叉概率在0.6-0.9范围内,变异概率在0.01-0.1范围内。可以通过手动输入进行调试。种群初始化主要代码如下:voidinitialize。"(sran

5、d(time(NULL);inti,j;for(i=0;i<popsize;i+)for(j=0;j<23;j+)populationi.genej=rand()%2;voidtransform。/染色体转化为实数xinti,j;for(i=0;i<=popsize+1;i+)populationi.realnumber=0;for(j=0;j<23;j+)populationi.realnumber+=populationi.genej*pow(2,j);populationi.realnumber=populationi.realnumber*2*pi/(pow(2

6、,23)-1);voidcal_fitness()/计算适应度inti;for(i=0;i<popsize;i+)populationi.fitness=populationi.realnumber*sin(populationi.realnumber)+6;voidselect()/选择操作intmem,i,j,k;doublesum=0;doublep;for(mem=0;mem<popsize;mem+)sum+=populationmem.fitness;for(mem=0;mem<popsize;mem+)populationmem.rfitness=populat

7、ionmem.fitness/sum;population0.cfitness=population0.rfitness;for(mem=1;mem<popsize;mem+)populationmem.cfitness=populationmem-1.cfitness+populationmem.rfitness;for(i=0;i<popsize;i+)/轮赌法选择机制p=rand()%1000/1000.0;if(p<population0.cfitness)newpopulationi=population0;elsefor(j=0;j<popsize;j+)i

8、f(p>=populationj.cfitness&&p<populationj+1.cfitness)newpopulationi=populationj+1;for(i=0;i<popsize;i+)/复制给下一代populationi=newpopulationi;)voidcross()/交叉操作inti,mem,one;intfirst=0;doublex;for(mem=0;mem<popsize;mem+)x=rand()%1000/1000.0;if(x<pcross)+first;if(first%2=0)Xover(one,me

9、m);/个体间染色体进行交叉函数elseone=mem;)voidmutate()/变异操作inti,j,t;doublex;for(i=0;i<popsize;i+)for(j=0;j<length;j+)x=rand()%1000/1000.0;if(x<pvariation)if(populationi.genej)populationi.genej=0;elsepopulationi.genej=1;)voidcal_max()/计算最大值inti;doublemax,sum=0;intmax_m;max=population0.fitness;for(i=0;i&l

10、t;popsize-1;i+)if(populationi.fitness>populationi+1.fitness)if(populationi.fitness>=max)max=populationi.fitness;max_m=i;)elseif(populationi+1.fitness>=max)(max=populationi+1.fitness;max_m=i+1;)if(max>populationpopsize.fitness)(iteration=0;for(i=0;i<length;i+)populationpopsize.genei=po

11、pulationmax_m.genei;populationpopsize.fitness=populationmax_m.fitness;)for(i=0;i<length;i+)sum=populationpopsize.genei-populationmax_m.genei;if(sum=0)iteration+;transform();printf("%f,%f,%f,%fn",populationpopsize.fitness,populationpopsize+1.fitness,populationpopsize.realnumber,populatio

12、npopsize+1.realnumber);)【实验结果】第9s俏:最大值旷：2.81907S-x=：2.007108010001110001101110010呼号61弋二最大值y=：2.819075,x-：2J007108染色林：01010001110001101110010第97代工最大值819610.x=：2J020339第最大值产二：281y610Ax=：2,020339fefl；：01010010010100001110011第99件二最大值*二2.&1”1呢x-：2-B20339染色林：01010010010100001110011S1眄代:最大值9=二普.81961团，x=：2.020339求取大值/取小值”射：S95代:最小值y=：-3814440，x=：491314011001000001011011110011R96代:最小值y=：-3814470x=：4913140=11001000001011011110011第971弋:最'值y=：-3814470，x=：4913140染色饪工10010