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

上传人:H****o 文档编号:33670635 上传时间:2022-08-12 格式:PDF 页数:3 大小:53.55KB
返回 下载 相关 举报
2022年遗传算法-求函数零点-C语言代码 .pdf_第1页
第1页 / 共3页
2022年遗传算法-求函数零点-C语言代码 .pdf_第2页
第2页 / 共3页
点击查看更多>>
资源描述

《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 页 - - - - - - - - -

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 技术资料 > 技术总结

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁