一个关于遗传算法优化的简单例子

在课程上学了一些关于遗传算法的思想的,想用这个思想来写一个简单的小例子。

先来说遗传算法的思想:遗传算法是模拟生物的遗传、变异、选择、进化来对问题的解进行优化,可以理解为将一组初始解看成是“基因”,在求解的开始设置一个过滤器,对“基因”进行筛选,通过如果目前生成的“基因”暂不满足上述条件,那么“基因”就要开始“变异”,在迭代过程中通过产生的随机数,对“基因”进行更改,达到“变异”的目的,也就是”遗传“给了下一代。而因为变异的随机性,计算机能够相当大的范围内对问题的解进行搜索,直至随着迭代的代数继续增加而解几乎不再变化为止。这时,我们可以说,我们的得到了进化后的最优解。

例:用遗传算法思想对函数 f(x)=x^2+sin(x) 进行优化,找出该函数的极值。

(使用的Anaconda3的spyder编写的)

 1 import random as R
 2 import math
 3 import matplotlib.pyplot as plt #通常遗传算法会优图像演示需求
 4 f=lambda x:x*x+math.sin(x) #定义lambda函数,即我们要求最优解得函数
 5 x1=x0=float(input('enter a value of x:'))
 6 y1=y0=f(x0)
 7 count=0 #记录迭代代数
 8 delta=0.001 #每次变异,对x的修改量,可以使之尽量小,从而避免在遗传变异中错过最优值
 9 xi=[] #通过下x,y来记录迭代的点坐标,方便画图
10 yi=[]
11 while(count<100000): #迭代代数,可以自己设置
12     sign=R.random()
13     if sign<=0.5:
14         x0+=delta
15         if abs(y0-f(x0))<1e-6: #进化完成条件,极值处的导数为0,此时,我们将精度定位1e-6
16             break
17         count+=1
18         if y0>f(x0):
19             y0=f(x0)
20     else:
21         x1-=delta
22         if abs(y0-f(x1)<1e-6):
23             break
24         count+=1
25         if y0>f(x1):
26             y0=f(x1)  
27     if count%25==0:
28         xi.append(count)
29         yi.append(y0)
30 plt.plot(xi,yi,'r-')
31 print('极值结果为:{:.6f},迭代代数为{}代'.format(y0,count))

以上便是关于遗传算法的介绍以及一个Python小程序,虽然简单,但整个算法的基本步骤基本都包含了,适合复习时浏览一下。

上一篇:git的简单用法


下一篇:窗口拖动