《2022年遗传算法-求函数零点-C语言代码 .pdf》由会员分享,可在线阅读,更多相关《2022年遗传算法-求函数零点-C语言代码 .pdf(3页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、遗传算法求函数的零点本文适合初学者刚刚学会一点遗传算法, 并成功用代码实现, 希望与大家分享。 本文主要以求函数零点来说明遗传算法,没有大段的理论描述,代码是主要的!例:已知函数 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 #i
2、nclude #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 页 - - - - - - -
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; i
4、nt 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; 名师资料总结 - - -精品资料欢迎下载
5、- - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 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,pop
6、i,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 页 - - - - - - - - -