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

遗传算法求函数的零点 本文适合初学者 刚刚学会一点遗传算法，并成功用代码实现，希望与大家分享。本文主要以求函数零点来说明遗传算法，没有大段的理论描述，代码是主要的！ 例： 已知函数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]; } void out() { for(int i=1;i=size;i++) printf(x=%lf adapt=%lf up=%lf:\n,pop[i],adapt(pop[i]),pop_adapt[i]); } int main() { unsigned st,end; srand(time(0)); st=clock(); create_pop(); pop_adapt[0]=0; for(int j=0;j170;j++) { for(int i=1;i=size;i++) { pop_adapt[i]=pop_adapt[i-1]+adapt(pop[i]); if(adapt(pop[i])0.999) {out();end=clock();printf(yes time=%u\n j=%d\n,end-st,j-1); return 0;} } for(int i=1;i=si