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

遗传算法求函数的零点本文适合初学者刚刚学会一点遗传算法， 并成功用代码实现， 希望与大家分享。 本文主要以求函数零点来说明遗传算法，没有大段的理论描述，代码是主要的！例：已知函数 y = 2x + x*x - 2 ; 求它在区间 (0,1)上，y=0时 x 的值。解：函数 y = 2x + x*x 2 的图像为区间(0,1)上，函数的零点大致为0.6 左右。我们采用遗传算法来求方程2x + x*x 2=0在区间（ 0，1）上的根。网上有很多关于遗传算法的理论知识，估计这并不是初学者所关心的，所以这里就直接给出代码 . #include #include #include #include #include #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 popsize+1,pop_adaptsize+1; double prd()/ 产生0,1)之间的随机数 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 3 页 - - - - - - - - - return (rand()%10000)/10000.0; void create_pop()/建立初始种群 for(int i=1;i=size;i+) popi=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_adapti)return popi; else continue; int variation(int x)/ 变异 int length,t; length=rand()%wei; t=1length)length)length; tp4=tp3yy; if(prd()pm)*x1=variation(tp1tp4)/10000.0; else *x1=(tp1tp4)/10000.0; if(prd()pm)*y1=variation(tp2tp3)/10000.0; else *y1=(tp2tp3)/10000.0; void update_next()/ 产生下代 double tempsize+1,x,y; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 3 页 - - - - - - - - - for(int i=1;i=size;i+=2) x=chose();y=chose(); if(prd()pc)cross(x,y,&tempi,&tempi+1); else tempi=x;tempi+1=y; for(int i=1;i=size;i+) popi=tempi; void out() for(int i=1;i=size;i+) printf(x=%lf adapt=%lf up=%lf:n,popi,adapt(popi),pop_adapti); int main() unsigned st,end; srand(time(0); st=clock(); create_pop(); pop_adapt0=0; for(int j=0;j170;j+) for(int i=1;i0.999) out();end=clock();printf(yes time=%un j=%dn,end-st,j-1); return 0; for(int i=1;i=size;i+) pop_adapti/=pop_adaptsize; update_next(); / double t=0.6534;printf(%lfn,adapt(t); return 0; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 3 页 - - - - - - - - -