遗传算法-求函数零点-C语言代码.pdfVIP

遗传算法求函数的零点 本文适合初学者 刚刚学会一点遗传算法， 并成功用代码实现， 希望与大家分享。 本文主要以求函 数零点来说明遗传算法，没有大段的理论描述，代码是主要的！ 例： 已知函数 y = 2^x + x*x - 2 ; 求它在区间 (0,1)上， y=0 时 x 的值。 解 ：函数 y = 2^x + x*x – 2 的图像为 区间 (0,1)上，函数的零点大致为 0.6 左右。 我们采用遗传算法来求方程 2^x + x*x – 2=0 在区间（ 0，1）上的根。 网上有很多关于遗传算法的理论知识， 估计这并不是初学者所关心的， 所以这里 就直接给出代码 . #includestdio.h #includestdlib.h #includemath.h #includetime.h #includestring.h #define size 40//种群大小 #define wei 14// 染色体长度 #define pc 0.9// 交叉概率 #define pm 0.03// 变异概率 //y=pow(2.0,x)+x*x-2;(0,1) 9999=10011100001111 14位 double pop[size+1],pop_adapt[size+1]; double prd()// 产生 [0,1)之间的随机数 { return (rand()%10000)/10000.0; } void create_pop()// 建立初始种群 { for(int i=1;i=size;i++) pop[i]=prd(); } double adapt(double x)// 适应度 { double temp; temp=fabs(pow(2.0,x)+x*x-2); return pow(8.0,-temp); } double chose()// 自然选择，适者生存 { double x; x=prd(); for(int i=1;i=size;i++) if(xpop_adapt[i])return pop[i]; else continue; } int variation(int x)// 变异 { int length,t; length=rand()%wei; t=1length; return x^t; } void cross(double x,double y,double *x1,double *y1)// 交叉 { int xx,yy,length; int tp1,tp2,tp3,tp4; xx=int(x*10000);yy=int(y*10000); length=rand()%13+1; tp1=(xxlength)length; tp2=tp1^xx; tp3=(yylength)length; tp4=tp3^yy; if(prd()pm){*x1=variation(tp1^tp4)/10000.0;} else *x1=(tp1^tp4)/10000.0; if(prd()pm){*y1=variation(tp2^tp3)/10000.0;} else *y1=(tp2^tp3)/10000.0; } void update_next()// 产生下代 { double temp[size+1],x,y; for(int i=1;i=size;i+=2) { x=chose();y=chose(); if(prd()pc)cross(x,y,temp[i],temp[i+1]); else {temp[i]=x;temp[i+1]=y;} } for(int i=1;i=size;i++) pop[i]=temp[i]; } v