遗传算法-01背包

遗传算法

算法思想
遗传算法(Genetic Algorithm, GA)是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型,是一种通过模拟自然进化过程搜索最优解的方法。
其主要特点是直接对结构对象进行操作,不存在求导和函数连续性的限定;具有内在的隐并行性和更好的全局寻优能力;采用概率化的寻优方法,不需要确定的规则就能自动获取和指导优化的搜索空间,自适应地调整搜索方向。
遗传算法以一种群体中的所有个体为对象,并利用随机化技术指导对一个被编码的参数空间进行高效搜索。其中,选择、交叉和变异构成了遗传算法的遗传操作;参数编码、初始群体的设定、适应度函数的设计、遗传操作设计、控制参数设定五个要素组成了遗传算法的核心内容。

设计思路

遗传算法是从代表问题可能潜在的解集的一个种群开始的,一个种群由基因编码的一定数量的个体组成。每个个体实际上是染色体带有特征的实体。
初代种群产生之后,按照适者生存和优胜劣汰的原理,逐代(generation)演化产生出越来越好的近似解,在每一代,根据问题域中个体的适应度(fitness)大小选择(selection)个体,并借助于自然遗传学的遗传算子(genetic operators)进行组合交叉(crossover)和变异(mutation),产生出代表新的解集的种群。

这个过程将导致种群像自然进化一样的后生代种群比前代更加适应于环境,末代种群中的最优个体经过解码(decoding),可以作为问题近似最优解。

遗传算法应用步骤:
    1)确定决策变量及各种约束条件,即个体的表现型X和问题的解空间;
    2)建立优化模型 (目标函数最大OR 最小) 数学描述形式 量化方法;
    3)染色体编码方法;
    4)解码方法;
    5)个体适应度的量化评价方法 F(x)
    6)设计遗传算子;
    7)确定有关运行参数。

程序代码(使用遗传算法解决01背包问题)

#include <windows.h>
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<time.h>#define  S				5			//种群的规模
#define  Pc				0.8			//交叉概率
#define  Pm				0.05			//突变概率
#define  KW             1000			//背包最大载重1000
#define  N              50			//物体总数
#define	 T				800			//最大换代数
#define	 ALIKE	        0.05			//判定相似度
int		 stop=0;						//初始化函数中初始化为所有价值之和
int		 t=0;						//目前的代数
int value[]={
220,208,198,192,180,180,165,162,160,158,155,130,125,122,120,118,115,110,105,101,100,100,98,96,95,90,88,82,80,77,75,73,72,70,69,66,65,63,60,58,56,50,30,20,15,10,8,5,3,1};
int weight[]={
80,82,85,70,72,70,66,50,55,25,50,55,40,48,50,32,22,60,30,32,40,38,35,32,25,28,30,22,25,30,45,30,60,50,20,65,20,25,30,10,20,25,15,10,10,10,4,4,2,1};struct individual				//个体结构体
{bool chromsome[N];		//染色体编码double fitness;				//适应度//即本问题中的个体所得价值double weight;			//总重量
};
int best=0;
int same=0;
individual X[S],Y[S],bestindividual;/************************************************************************/
int  comp(individual bestindividual,individual temp);	//比较函数
void checkalike(void);							//检查相似度函数
void GenerateInitialPopulation(void); 				//初始种群
void CalculateFitnessValue(void);					//适应度
void SelectionOperator(void);						//选择
void CrossoverOperator(void);					//交叉
void MutationOperator(void);					//变异
void FindBestandWorstIndividual(void);				//寻找最优解
void srand(unsigned int seed);					//随机生成
/************************************************************************/int comp(individual bestindividual,individual temp)//比较函数
{int fit=0,w=0;//第一种不变:操作后不满足重量函数,第二种:操作后适应度小于操作前for(int i=0;i<N;i++){fit+=temp.chromsome[i]*value[i];w+=temp.chromsome[i]*weight[i];}if(w>KW)return -1;return (bestindividual.fitness>fit?-1:1);//如果小于原来值或者不满足重量函数,则返回-1
}/************************************************************************/
void Checkalike(void)
{int i=0,j=0;for( i=0;i<S;i++)//相似度校验{for(j=0;j<N;j++){bool temp=X[i].chromsome[j];for(int k=1;k<S;k++){if(temp!=X[k].chromsome[j])break;}}if(j==N)same++;}if(same>N*ALIKE)//大于ALIKE作为判定为早熟{int minindex=0;for(int n=0;n<S;n++)if(X[n].fitness<X[minindex].fitness)minindex=n;//确定最小for (j=0; j<N;j++)//重新生成{bool m=(rand()%10<5)?0:1;X[minindex].chromsome[j]=m;X[minindex].weight+=m*weight[j];//个体的总重量X[minindex].fitness+=m*value[j]; //个体的总价值}}
}/************************************************************************/
void GenerateInitialPopulation(void)//初始种群,保证每个值都在符合条件的解
{int i=0,j=0; bool k;for(i=0;i<N;i++)stop+=value[i];//设置理论最优值for (i=0; i<S; i++){int w=0,v=0;for (j=0; j<N;j++){k=(rand()%10<5)?0:1;X[i].chromsome[j]=k;w+=k*weight[j];//个体的总重量v+=k*value[j]; //个体的总价值}if(w>KW) i--;	   //如果不是解,重新生成else{X[i].fitness=v;X[i].weight=w;if(v==stop){bestindividual=X[i];return;}//这种情况一般不会发生}}
}
/************************************************************************/void CalculateFitnessValue()
{int i=0,j=0;for (i=0; i<S; i++){int w=0,v=0;for (j=0; j<N;j++){w+=X[i].chromsome[j]*weight[j];//个体的总重量v+=X[i].chromsome[j]*value[j]; //个体的总价值}X[i].fitness=v;X[i].weight=w;if(v==stop){bestindividual=X[i];return;}//符合条件情况下最优解这种情况一般不会发生if(w>KW) X[i]=bestindividual;//如果不是解,找最好的一个解代之}
}
/************************************************************************/void SelectionOperator(void)
{int i, index;double p, sum=0.0;double cfitness[S];//选择、累积概率individual newX[S];for (i=0;i<S;i++)sum+=X[i].fitness;//适应度求和for (i=0;i<S; i++)cfitness[i]=X[i].fitness/sum; //选择概率for (i=1;i<S; i++)cfitness[i]=cfitness[i-1]+cfitness[i];//累积概率for (i=0;i<S;i++){p=(rand()%1001)/1000.0;//产生一个[0,1]之间的随机数index=0;while(p>cfitness[index])//轮盘赌进行选择{index++;}newX[i]=X[index];}for (i=0; i<S; i++)X[i]=newX[i];//新的种群
}/************************************************************************/
void CrossoverOperator(void)//交叉操作
{int i=0, j=0,k=0;individual temp;for(i=0; i<S; i++){int p=0,q=0;do{p=rand()%S;//产生两个[0,S]的随机数q=rand()%S;}while(p==q);int w=1+rand()%N;//[1,N]表示交换的位数double r=(rand()%1001)/1000.0;//[0,1]if(r<=Pc){for(j=0;j<w;j++){temp.chromsome[j]=X[p].chromsome[j];//将要交换的位先放入临时空间X[p].chromsome[j]=X[q].chromsome[j];X[q].chromsome[j]=temp.chromsome[j];}}if(p==best)if(-1==comp(bestindividual,X[p]))//如果变异后适应度变小X[p]=bestindividual;if(q==best)if(-1==comp(bestindividual,X[q]))//如果变异后适应度变小X[q]=bestindividual;}
}
/************************************************************************/void MutationOperator(void)
{int i=0, j=0,k=0,q=0;double p=0;for (i=0; i<S; i++){for (j=0; j<N; j++){p=(rand()%1001)/1000.0;if (p<Pm)//对每一位都要考虑{if(X[i].chromsome[j]==1)X[i].chromsome[j]=0;else	X[i].chromsome[j]=1;}}if(i==best)if(-1==comp(bestindividual,X[i]))//如果变异后适应度变小X[i]=bestindividual;}
}
/************************************************************************/void FindBestandWorstIndividual(void)
{int i;bestindividual=X[0];for (i=1;i<S; i++){if (X[i].fitness>bestindividual.fitness){bestindividual=X[i];best=i;}}
}/*主函数*****************************************************************/
int main()
{t=0;GenerateInitialPopulation(); //初始群体包括产生个体和计算个体的初始值while (t<=T){FindBestandWorstIndividual();	//保存当前最优解SelectionOperator();			//选择CrossoverOperator();			//交叉MutationOperator();			   //变异Checkalike();				  //检查相似度CalculateFitnessValue();	 //计算新种群适应度t++;}FindBestandWorstIndividual();			//找到最优解printf(" 物品价值:");for(int k=0;k<N;k++){printf(" %d ",value[k]);}printf("\n");printf(" 物品重量:");for(int k=0;k<N;k++){printf(" %d  ",weight[k]);}printf("\n");printf("背包容量 %d\n",1000);   	//输出最优值printf("-----------------------------\n");printf("最优值 %f\n",bestindividual.fitness);   	//输出最优值printf("对应重量 %f\n",bestindividual.weight);       //对应重量printf("线性解:");for(int k=0;k<N;k++){if(bestindividual.chromsome[k]==1){  //输出最优解printf(" %d ",1);}else{printf(" %d ",0);}}printf("\n");printf("\n");return 0;
}
/*结束***********************************************************************/

测试例

物品价值:

 220  208  198  192  180  180  165  162  160  158  155  130  125  122  120  118  115  110  105  101  100  100  98  96  95  90  88  82  80  77  75  73  72  70  69  66  65  63  60  58  56  50  30  20  15  10  8  5  3  1

物品重量:

 80   82   85   70   72   70   66   50   55   25   50   55   40   48   50   32   22   60   30   32   40   38   35   32   25   28   30   22   25   30   45   30   60   50   20   65   20   25   30   10   20   25   15   10   10   10   4   4   2   1

运行结果

最优值 3078.000000
对应重量 1000.000000
线性解: 1  0  0  1  1  0  1  1  1  1  1  0  1  1  1  1  1  0  1  1  0  0  0  1  1  1  1  1  1  1  0  0  0  0  1  0  1  1  0  1  0  0  0  0  0  0  1  1  1  1

分析
试验中用到的物品的重量和价值分别存储于以下两个数组之中。

int value[]={
220,208,198,192,180,180,165,162,160,158,155,130,125,122,120,118,115,110,105,101,100,100,98,96,95,90,88,82,80,77,75,73,72,70,69,66,65,63,60,58,56,50,30,20,15,10,8,5,3,1};
int weight[]={
80,82,85,70,72,70,66,50,55,25,50,55,40,48,50,32,22,60,30,32,40,38,35,32,25,28,30,22,25,30,45,30,60,50,20,65,20,25,30,10,20,25,15,10,10,10,4,4,2,1};

遗传算法有许多不足,算法对新空间的探索能力是有限的,也容易收敛到局部最优解。涉及到大量个体的计算,当问题复杂时,计算时间很长。
但是他也有很多优点:

  1. 与问题领域无关切快速随机的搜索能力。
  2. 搜索从群体出发,具有潜在的并行性,可以进行多个个体的同时比较,robust.
  3. 搜索使用评价函数启发,过程简单
  4. 使用概率机制进行迭代,具有随机性。
  5. 具有可扩展性,容易与其他算法结合。

参考:
短短的路走走停停被抢注啦
鹤鹤有明
wangqiuyun

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/536106.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Angular Web App部署Linux Nginx Https

Angular Web App部署Linux Nginx Https 提示:这篇文章是基于内网的 互联网就开始将 WEB 服务从 HTTP 迁移到 HTTPS,而现在为了更快的推进 HTTPS 的普及,Chrome 将从 2018 年 7 月起标记所有的 HTTP 网站为不安全链接。 HTTPS 会逐渐成为 WEB 服务的标配,最最重要的是,它能…

SOLO算法简读

论文链接&#xff1a;https://arxiv.org/abs/1912.04488 代码链接&#xff1a;https://github.com/WXinlong/SOLO 摘要 提出一种新的实例分割方法。与语义分割等其他密集预测任务相比&#xff0c;实例分割的难度要大得多。为了预测每个实例的掩码&#xff0c;主流方法要么遵…

Rxjs的flatMap使用

Rxjs的flatMap使用 flatMap是Rxjs比较绕的一个概念&#xff0c;这里我们只是讲解如何使用。在Rxjs 4.0版本时叫flatMap,在Rxjs 5.0时被更名为margeMap,现在flatMap作为margeMap的别名使用&#xff0c;这是考虑向下兼容。 官方flatMap的定义&#xff1a; Projects each sourc…

关于Loss的简单总结

Dice Loss 参考&#xff1a;https://blog.csdn.net/l7H9JA4/article/details/108162188 Dice系数&#xff1a; 是一种集合相似度度量函数&#xff0c;通常用于计算两个样本的相似度&#xff0c;取值范围为[0,1]。 s2∣X∩Y∣∣X∣∣Y∣s \frac{2|X ∩ Y|}{|X||Y|} s∣X∣∣Y…

Angular_PWA使用+Demo

Angular_PWA使用+Demo 什么是PWA PWA(Progressive Web App)利用TLS,webapp manifests和service workers使应用程序能够安装并离线使用。 换句话说,PWA就像手机上的原生应用程序,但它是使用诸如HTML5,JavaScript和CSS3之类的网络技术构建的。 如果构建正确,PWA与原生应…

SOLOv2论文简读

论文&#xff1a;SOLOv2: Dynamic, Faster and Stronger 代码&#xff1a;https://github.com/WXinlong/SOLO 摘要 主要提出了作者在SOLOv2中实现的优秀的实例分割方法&#xff0c;旨在创建一个简单、直接、快速的实例分割框架&#xff1a; 通过提出动态学习对象分割器的mas…

Angular6_PWA

Angular6_PWA Angular正式发布了V6.0,我们已经可以利用对应的@angular/cli V6.0来直接开发PWA应用了。 第一步:安装@angular/cli V6.0 如果你机器上有老版本,请先卸载。 打开你的终端,执行: npm install -g @angular/cli 或 cnpm install -g @angular/cli 安装成功…

Ubuntu18.04 关于使用vnc的踩坑

由于种种原因&#xff0c;手上多了一台可使用的桌面版Ubuntu&#xff0c;正好用来测试代码&#xff0c;方便调试。因为只能远程&#xff0c;所以需要配置远程连接。因此就打算使用vnc进行远程连接&#xff0c;谁料一路坎坷&#xff0c;特此记录。 安装 设置桌面共享 需要注意…

App_Shell模型

App_Shell模型 App Shell 架构是构建 Progressive Web App 的一种方式,这种应用能可靠且即时地加载到您的用户屏幕上,与本机应用相似。 App shell是支持用户界面所需的最小的 HTML、CSS 和 JavaScript,如果离线缓存,可确保在用户重复访问时提供即时、可靠的良好性能。这意…

Angular6_服务端渲染SSR

Angular6_服务端渲染 在使用服务端渲染之前,需要安装最新版本的Angular。 npm install -g @angular/cli 或 cnpm install -g @angular/cli github项目 创建项目 ng new PWCat --routing 为项目添加universalng g universal --client-project=PWCat 或

Jenkins自定义主题教程

Jenkins自定义主题 由于Jenkins自带的样式比较丑陋&#xff0c;所以有很多第三方的样式库&#xff0c;这里针对jenkins-material-theme样式库做一个安装教程。 下载样式库 下载连接 Select your color 选择一个你喜欢的主题颜色。Choose your company logo 上传你自定义的…

IndexedDB_Web 离线数据库

IndexedDB_Web 离线数据库 本文会从头剖析一下 IndexedDB 在前端里面的应用的发展。 indexedDB 目前在前端慢慢得到普及和应用。它正朝着前端离线数据库技术的步伐前进。以前一开始是 manifest、localStorage、cookie 再到 webSQL&#xff0c;现在 indexedDB 逐渐被各大浏览器认…

Angular 单元测试讲解

Angular_单元测试 测试分类 按开发阶段划分按是否运行划分按是否查看源代码划分其他ATDD,TDD,BDD,DDD ATDDTDDBDDDDDAngular单元测试 Karma的介绍jasmine介绍单元测试的好处使用jasmine和karma创建一个Angular项目Karma配置Test.ts文件测试体验测试Form测试服务service常用断言…

基于 Docker 的微服务架构

基于 Docker 的微服务架构-分布式企业级实践前言Microservice 和 Docker服务发现模式客户端发现模式Netflix-Eureka 服务端发现模式ConsulEtcdZookeeper 服务注册自注册模式 Self-registration pattern第三方注册模式 Third party registration pattern小结一 服务间的 IPC 机制…

funcode游戏实训,java及C/C++,网上整理

软件&#xff0c;常见错误都有。 所有资源可到公众号获取(源码也是)&#xff0c;不再直接分享

Docker 容器部署 Consul 集群

Docker 容器部署 Consul 集群 Consul 介绍 Consul 提供了分布式系统的服务发现和配置的解决方案。基于go语言实现。并且在git上开放了源码consul-git。consul还包括了分布式一致协议的实现&#xff0c;健康检查和管理UI。Consul和zk相比较起来&#xff0c;更加轻量级&#xf…

swing皮肤包 substance

分享一下swing皮肤包substance 资源可到公众号获取

基于Android的聊天软件,Socket即时通信,实现用户在线聊天

基于Android的聊天软件&#xff0c;Socket即时通信&#xff0c;单聊&#xff0c;聊天室&#xff0c;可自行扩展功能&#xff0c;完善细节。 【实例功能】 1.运行程序&#xff0c;登录界面, 注册账号功能 2.进入主界面&#xff0c;有通讯录, 个人信息。 3.点击好友会话框&#…

用Docker搭建Elasticsearch集群

用Docker搭建Elasticsearch集群 对于用Docker搭建分布式Elasticsearhc集群的一个介绍&#xff0c;以及一些实施中遇到问题的总结 搜索服务简述 结合业务的场景&#xff0c;在目前的商品体系需要构建搜索服务&#xff0c;主要是为了提供用户更丰富的检索场景以及高速&#xf…

Go实现简单的RESTful_API

Go实现简单的RESTful_API 何为RESTful API A RESTful API is an application program interface (API) that uses HTTP requests to GET, PUT, POST and DELETE data. A RESTful API – also referred to as a RESTful web service – is based on representational state t…