机器学习——遗传算法、粒子群算法、蒙特卡罗算法

目录

一、遗传算法——收敛速度慢

二、粒子群算法(鸟群、蚁群寻食算法)——收敛速度快

三、蒙特卡罗算法原理也是利用了穷举所有解之后做比较的思路:


几天没写总结了,说下我理解的遗传算法和粒子群算法,后面说下蒙特卡罗算法!两个算法都是用来求最优解的,如:某函数f(x)在(a,b)上的最大或小值、连通图中从A到B的最短路径、工程或活动安排(AOE关键路径)的最佳方案,同时可以结合(非)线性回归、支持向量机、神经网络分类算法用之来求解最优参数w(原本三种算法都采用梯度下降法慢慢求解最优参数)。总之,对于要求解的变量(或方案)的解空间中存在许多解(称之为误差解空间),就是存在误差解空间的模型都可以使用上面两种方法找到最优解。

其实最暴力的方法就是把所有解都列举出之后带入原函数关系式再对比结果就能得到最优解。遗传算法和粒子群算法,就是利用这一简单的思路出发的。但是他们只要使用者给列出(>2)种解,算法会自动由给出的一组的解为基础穷举出所有解情况,并且他们会朝着最优解的方向以特定步长快速穷举出中间部分解直到找到最优解。两种算法都是在如何使最初解在穷举的过程中可以快速向最优解的方向靠近提出了自己的想法。

       假如在区间(30,455)上求f(x)最大值,初始列出x1,x2,x3,x4=(31,37,34,38)四个解——种群。

一、遗传算法——收敛速度慢

  1. 将种群x代入f(x),得到f(31)<f(34)<f(37)<f(38)。选择淘汰31,令x1=x4=38,此时,x1,x2,x3,x4=(38,37,34,38)。
  2. X用二进制表示之后交叉、变异,原种群本来就把最差解给淘汰了,在经过交叉变异种群x的值整体向最优解方向变成更优的种群(可以了解一下二进制编码高位不变低位互换会使最大值变得更大)。
  3. (2)中的新种群解码带入函数比较结果大小重复(1),重复(2),直到种群x来回在某个值两边变化时得到最优解。

 

二、粒子群算法(鸟群、蚁群寻食算法)——收敛速度快

如何改变x的值使它向着最优值方向不断穷举呢?计算一个矢量速度V,使得x朝着某个方向以一定的步长不断的改变优化。鸟群x1,x2,x3,x4=(31,37,34,38)寻找食物的过程中每只鸟(1)根据自己以往经验有个判断方向(2)综合观察队友的寻找经历可排除优化方向。因此,V的方向取决于(1),(2)。

当前V的大小=上一步V大小+(个体xi历史最优值-上一步xi的值)+(整体历史最优值-上一步xi的值)。

                                            (1)

  1. 初始化粒子群体(群体规模为>2),初始化x的值和x编号的初始速度
  2. 将x代入f(x)评价得出每个粒子的适应度
  3. Find the pbest对每个粒子,将其当前适应值和个体历史最佳适应值(pbest)做比较,用较大者更个体新历史最佳x。
  4. Find the gbest对每个粒子,将其当前适应值和群体历史最佳适应值(gbest)做比较,用较大者更新群体历史最佳x。
  5. 用公式(1)计算更新每个粒子的v,Vx更新种群x向着最优值方向迁移。
  6. 不断重复(1)——(5)直到找到最佳适应度值。

三、蒙特卡罗算法原理也是利用了穷举所有解之后做比较的思路:

主要利用积分中值定理求复杂积分值或一些复杂几何的面积、体积

  1. 计算机可以多次随机自动的产生一组解空间的随机数,相对于在穷举。再计算出一个概率值。
  2. 直角坐标系中有个知道面积的长方形,现在用一条曲线分成两部分,要求各部分的面积。我们可以向长方形中随机扔n颗豆子,然数出各部分有多少豆子,就能得出面积。蒙特卡罗算法充分利用了计算机可随机生成解空间中的的未变量的特性。
  • 2
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
本框架提供了有关粒子群算法(PSO)和遗传算法(GA)的完整实现,以及一套关于改进、应用、测试、结果输出的完整框架。 本框架对粒子群算法遗传算法进行逻辑解耦,对其中的改进点予以封装,进行模块化,使用者可以采取自己对该模块的改进替换默认实现组成新的改进算法与已有算法进行对比试验。试验结果基于Excel文件输出,并可通过设定不同的迭代结束方式选择试验数据的输出方式,包括: 1. 输出随迭代次数变化的平均达优率数据(设定终止条件区间大于0)。 2. 输出随迭代次数变化的平均最优值数据(设定终止条件区间等于0)。 本框架了包含了常用基准函数的实现以及遗传算法粒子群算法对其的求解方案实现和对比,如TSP,01背包,Banana函数,Griewank函数等。并提供大量工具方法,如KMeans,随机序列生成与无效序列修补方法等等。 对遗传算法的二进制编码,整数编码,实数编码,整数序列编码(用于求解TSP等),粒子群算法的各种拓扑结构,以及两种算法的参数各种更新方式均有实现,并提供接口供使用者实现新的改进方式并整合入框架进行试验。 其中还包括对PSO进行离散化的支持接口,和自己的设计一种离散PSO方法及其用以求解01背包问题的实现样例。 欢迎参考并提出宝贵意见,特别欢迎愿意协同更新修补代码的朋友(邮箱starffly@foxmail.com)。 代码已作为lakeast项目托管在Google Code: http://code.google.com/p/lakeast http://code.google.com/p/lakeast/downloads/list 某些类的功能说明: org.lakest.common中: BoundaryType定义了一个枚举,表示变量超出约束范围时为恢复到约束范围所采用的处理方式,分别是NONE(不处理),WRAP(加减若干整数个区间长度),BOUNCE(超出部分向区间内部折叠),STICK(取超出方向的最大限定值)。 Constraint定义了一个代表变量约束范围的类。 Functions定义了一系列基准函数的具体实现以供其他类统一调用。 InitializeException定义了一个代表程序初始化出现错误的异常类。 Randoms类的各个静态方法用以产生各种类型的随机数以及随机序列的快速产生。 Range类的实现了用以判断变量是否超出约束范围以及将超出约束范围的变量根据一定原则修补到约束范围的方法。 ToStringBuffer是一个将数组转换为其字符串表示的类。 org.lakeast.ga.skeleton中: AbstractChromosome定义了染色体的公共方法。 AbstractDomain是定义问题域有关的计算与参数的抽象类。 AbstractFactorGenerator定义产生交叉概率和变异概率的共同方法。 BinaryChromosome是采用二进制编码的染色体的具体实现类。 ConstantFactorGenerator是一个把交叉概率和变异概率定义为常量的参数产生器。 ConstraintSet用于在计算过程中保存和获取应用问题的各个维度的约束。 Domain是遗传算法求解中所有问题域必须实现的接口。 EncodingType是一个表明染色体编码类型的枚举,包括BINARY(二进制),REAL(实数),INTEGER(整型)。 Factor是交叉概率和变异概率的封装。 IFactorGenerator参数产生器的公共接口。 Population定义了染色体种群的行为,包括种群的迭代,轮盘赌选择和交叉以及最优个体的保存。 org.lakeast.ga.chromosome中: BinaryChromosome二进制编码染色体实现。 IntegerChromosome整数编码染色体实现。 RealChromosome实数编码染色体实现。 SequenceIntegerChromosome整数序列染色体实现。 org.lakeast.pso.skeleton中: AbstractDomain提供一个接口,将粒子的位置向量解释到离散空间,同时不干扰粒子的更新方式。 AbstractFactorGenerator是PSO中参数产生器的公共抽象类。 AbstractParticle定义了PSO种群中粒子的基本行为,最主要是实现了如何根据现有位置计算得到下一代粒子的位置的合法值。 ConstraintSet用于在粒子迭代过程中保存和获取应用问题的各个维度的约束。 AbstractSwarm.java各种拓扑结构的PSO种群的抽象父类,主要实现了种群迭代过程中计算流程的定义以及中间数据被如何输出到测试工具类。 Domain是PSO算法求解中所有问题域必须实现的接口。 DynamicFatorGenerator若种群在迭代过程中,w,c1,c2随迭代次数发生变化,那么它们的产生器需要继承这个抽象类。 Factor封装了w,c1,c2三个参数的字面值。 Location用于保存和获取迭代中粒子的位置和速度向量的数值。 NeighborhoodBestParticle定义了采用邻域版本的PSO算法的具体实现。主要是实现了如何根据邻域版本的PSO算法计算下一迭代中的粒子速度。 RingTopoSwarm定义环拓扑结构的具体实现,主要是定义了如何获取粒子的邻域粒子的方法。 StaticTopoSwarm静态拓扑结构的PSO算法的抽象父类。 org.lakeast.pso.swarm中包含粒子群拓扑结构的各种实现,基本见名知意。 对各种问题的求解样例位于org.lakeast.main包中,以...TaskTest结尾,基本见名知意。 以ShafferF6DomainTaskTes对ShafferF6函数进行求解(采用的是PSO,遗传算法样例参见TSPValueTaskTest)为例说明求解过程如下: 1. 入口函数位于org.lakeast.main.ShafferF6DomainTaskTest中,go函数执行。 2. 在go函数中,首先指定迭代次数(numberOfIterations),测试多少轮(testCount,多次运行以得到平均达优值),种群大小(popSize),邻域大小(neighborhoodSize),迭代结束条件(exitCondition,由于制定了迭代次数,所以设定为[0,0],也就是只有达到指定迭代次数才退出)。 3. 以testCount,numberOfIterations以及迭代结束条件exitCondition为参数构建TestBatch类的实例batch。这个类用来进行管理参与测试的各种具体算法,且把数据结果按指定的格式输出为Excel文件。 4. 指定PSO中的因子产生方法,采用ExponentFactorGenerator和ConstrictFactorGenerator两种方式(实现位于org.lakeast.pso.gen包)。 5. Y表示参与测试的算法数目。 6. Testable是所有可以被TestBatch测试的类需要实现的接口,以提供TestBatch生成结果Excel文件所需要的数据。 7. Domain接口是所有可以被算法解决的问题所需要实现的接口,比如说明该问题所需要的粒子位置约束范围,速度约束范围,以及适值评估的公司等。这里的Domain被实例化为ShafferF6Domain,也就是按照ShafferF6函数评估适值。 8. RingTopoSwarm是用来封装环拓扑邻域结构的类,NeighboordBestParticle是配合该类来实现按邻域最优更新速度而不是全局最优来更新。 9. 各个测试算法都被加入到TestBatch以后,batch.run()开始执行算法比较过程并输出结果Excel文件到C盘根目录(输出路径可在Testable接口中配置,除了生成Excel文件外,还可以通过修改log4j.properties在制定的位置产生运行结果日志)。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值