适应度函数及遗传操作

本文转载自:http://blog.csdn.net/b2b160/article/details/4680853/


这篇文章所举的例子非常给力,摘自《遗传算法原理与应用》,至少我是看过之后就明白了很多概念~


首先,什么是遗传算法?

遗传算法是为那些难以找到传统数学模型的难题找出一个解决方法。
遗传算法是仿真生物遗传学和自然选择机理,通过人工方式所构造的一类搜索算法,从某种程度上说遗传算法是对生物进化过程进行的数学方式仿真。


遗传算法基本概念

种群:初始给定的多个解的集合,它是问题解空间的一个子集。

个体:种群中的单个元素,通常由一个用于描述其基本遗传结构的数据结构来表示,如用0,1组成的长度为l的串来表示个体。
染色体:对个体进行编码后得到的编码串。染色体中的每一位称为基因,若干基因构成的有效信息段称为基因组。
适应度函数:用来对种群中个体的适应性进行度量的函数。


遗传算法执行步骤:




下面,举一个例子来简单模拟遗传计算的各个主要执行步骤

例:求下述二元函数的最大值:


(1) 个体编码
           遗传算法的运算对象是表示个体的符号串,所以必须把变量 x1, x2 编码为一种
       符号串。本题中,用无符号二进制整数来表示。
           因 x1, x2 为 0 ~ 7之间的整数,所以分别用3位无符号二进制整数来表示,将它
       们连接在一起所组成的6位无符号二进制数就形成了个体的基因型,表示一个可
       行解。
           例如,基因型 X=101110 所对应的表现型是:x=[ 5,6 ]。
           个体的表现型x和基因型X之间可通过编码和解码程序相互转换。


(2) 初始群体的产生
          遗传算法是对群体进行的进化操作,需要给其淮备一些表示起始搜索点的初始
      群体数据。
         本例中,群体规模的大小取为4,即群体由4个个体组成,每个个体可通过随机
     方法产生。
          如:011101,101011,011100,111001
         
 (3) 适应度汁算
          遗传算法中以个体适应度的大小来评定各个个体的优劣程度,从而决定其遗传
       机会的大小。
          本例中,目标函数总取非负值,并且是以求函数最大值为优化目标,故可直接
       利用目标函数值作为个体的适应度。


 (4)  选择运算

          选择运算(或称为复制运算)把当前群体中适应度较高的个体按某种规则或模型遗传到下一代群体中。一般要求适应度较高的个体将有更多的机会遗传到下一代
      群体中。                   
本例中,我们采用与适应度成正比的概率来确定各个个体复制到下一代群体中
     的数量。其具体操作过程是:
         •  先计算出群体中所有个体的适应度的总和  fi  ( i=1.2,…,M );
         •  其次计算出每个个体的相对适应度的大小 fi / fi ,它即为每个个体被遗传
             到下一代群体中的概率,
         •  每个概率值组成一个区域,全部概率值之和为1;
         •  最后再产生一个0到1之间的随机数,依据该随机数出现在上述哪一个概率区
             域内来确定各个个体被选中的次数。


(5)  交叉运算
        交叉运算是遗传算法中产生新个体的主要操作过程,它以某一概率相互交换某
    两个个体之间的部分染色体。
       本例采用单点交叉的方法,其具体操作过程是:
       • 先对群体进行随机配对;
       • 其次随机设置交叉点位置;
       • 最后再相互交换配对染色体之间的部分基因。



(6)  变异运算
         变异运算是对个体的某一个或某一些基因座上的基因值按某一较小的概率进
     行改变,它也是产生新个体的一种操作方法。
        本例中,我们采用基本位变异的方法来进行变异运算,其具体操作过程是:
        • 首先确定出各个个体的基因变异位置,下表所示为随机产生的变异点位置,
          其中的数字表示变异点设置在该基因座处;
        • 然后依照某一概率将变异点的原有基因值取反。


对群体P(t)进行一轮选择、交叉、变异运算之后可得到新一代的群体p(t+1)。


从上表中可以看出,群体经过一代进化之后,其适应度的最大值、平均值都得到了明显的改进。事实上,这里已经找到了最佳个体“111111”。 

      
[注意]需要说明的是,表中有些栏的数据是随机产生的。这里为了更好地说明问题,我们特意选择了一些较好的数值以便能够得到较好的结果,而在实际运算过程中,有可能需要一定的循环次数才能达到这个最优结果。

  • 47
    点赞
  • 230
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
好的,下面是使用鲸鱼算法(Whale Optimization Algorithm, WOA)优化支持向量机(SVM)的代码示例,其中包含最佳适应度和平均适应度: ``` import numpy as np from sklearn import datasets from sklearn.svm import SVC # 鲸鱼算法参数 max_iter = 100 # 最大迭代次数 pop_size = 30 # 种群大小 a = 2 # 收缩因子 c = 1 # 比例因子 # 数据加载 iris = datasets.load_iris() X = iris.data y = iris.target # SVM模型定义 svm = SVC(kernel='linear') # 目标函数 def fitness(X, y, svm, w): svm.coef_ = w.reshape(1, -1) svm.intercept_ = np.array([0]) svm.fit(X, y) return 1 - svm.score(X, y) # 初始化种群 pop = np.random.uniform(-1, 1, size=(pop_size, X.shape[1])) # 迭代优化 for t in range(max_iter): # 计算适应度 fit = [fitness(X, y, svm, w) for w in pop] best_fit = np.min(fit) avg_fit = np.mean(fit) # 更新最佳适应度和平均适应度 if t == 0: best_fit_all = best_fit avg_fit_all = avg_fit else: best_fit_all = min(best_fit_all, best_fit) avg_fit_all = (avg_fit_all * t + avg_fit) / (t + 1) # 更新最优解 best_idx = np.argmin(fit) best = pop[best_idx] # 更新鲸鱼位置 for i in range(pop_size): r = np.random.rand() # 随机数 A = 2 * a * r - a C = 2 * r # 控制参数 p = np.random.rand() # 随机概率 if p < 0.5: if np.abs(A) < 1: D = np.abs(C * best - pop[i]) new_pos = best - A * D else: rand_idx = np.random.choice(pop_size) D = np.abs(C * pop[rand_idx] - pop[i]) new_pos = pop[i] - A * D else: new_pos = best - C * np.random.uniform(-1, 1, X.shape[1]) # 边界处理 new_pos = np.clip(new_pos, -1, 1) # 更新种群 if fitness(X, y, svm, new_pos) < fitness(X, y, svm, pop[i]): pop[i] = new_pos # 输出结果 print('最佳适应度:', best_fit_all) print('平均适应度:', avg_fit_all) ``` 在上述代码中,我们使用鲸鱼算法来优化SVM的目标函数,即分类错误率。首先定义了鲸鱼算法的相关参数和SVM模型,然后定义了目标函数fitness,用于计算每个个体的适应度。接着初始化种群,并在每次迭代中更新最佳适应度和平均适应度。随后,我们对每个个体进行鲸鱼位置的更新,同时进行边界处理,最后根据新的位置更新种群。最终输出最佳适应度和平均适应度。 需要注意的是,本示例中的鲸鱼算法并不是专门为SVM设计的优化算法,其效果可能不如其他专门为SVM设计的优化算法(如SMO算法等)。
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值