首页 >> 数学 >> 生物 >> 计算机科学 >> 文章

100个透明三角形拼成的Firefox

这是个真实的故事。

从前在海岸边有一群扇贝在悠哉游哉地生活繁衍着。它们自然是衣食不愁,连房子也有了着落。它们担忧的只有一件事:每隔一段时间,总有一个人来挖走它们之中的一部分。当然啦,挖回去干什么这大家都知道。但扇贝们不知道的是,这人的家族图腾是Firefox的图标,所以他总是选择那些贝壳花纹长得比较不像Firefox图标的扇贝。

这种状况持续了好几十万代。大家应该也猜到扇贝们身上发生什么事情了:它们的贝壳上都印着很像Firefox图标的图案。

可能有些读者会说:你这不是糊弄我们么,Firefox才有多少年历史,你就搞了个几十万代的扇贝?

确有其事,但是这些扇贝不是真实的,它们在我电脑的内存里边生活。它们是一个遗传算法程序的一部分,这个程序的目的就是用100个半透明三角形把Firefox的图标尽可能像地画出来。

什么是遗传算法呢?

简单地说,遗传算法是一种解决问题的方法。它模拟大自然中种群在选择压力下的演化,从而得到问题的一个近似解。

在二十世纪五十年代,生物学家已经知道基因在自然演化过程中的作用了,而他们也希望能在新出现的计算机上模拟这个过程,用以尝试定量研究基因与进化之间的关系。这就是遗传算法的滥觞。后来,有人将其用于解决优化问题,于是就产生了遗传算法。

那么,具体来说,在计算机里边是怎么模拟进化过程的呢?

我们还是以开头提到的程序为例。

首先,我们知道,生物个体长什么样子很大程度上是由染色体上的基因决定的。同样,如果我们把100个半透明三角形组成的东西看成一个生物个体的话(为了说话方便,称为扇贝吧),我们也可以说它的样子是由这些三角形的具体位置和颜色决定的。所以,我们可以把一个一个的半透明三角形看作是这些扇贝的“基因”。而组成扇贝的这100个基因就组成了每个扇贝个体的“染色体”(chromosome)。

从下面的图可以大约看出来这些基因是怎么决定扇贝的样子的(为了观看方便,我们只画出其中五个三角形):

基因组成的个体

然后,扇贝们除了生活,当然还要繁衍后代。生物界中的繁衍无非就是父母的基因组合产生新的个体,而在这个程序里边我们当然也这么办:选择两个原有的扇贝,然后从这两个扇贝的染色体中随机选取一共100个基因组成新个体的染色体。如图所示:(仍然是将扇贝看成是五个三角形组成的)

为了产生新的基因,使基因的种类更多样化,在组合的时候,新的扇贝的基因有一定的概率发生变异。也就是说,其中的一个透明三角形的位置或者颜色会随机改变,如图(仍然是五个三角形……我偷工减料……):

扇贝基因的变异

其次,为了使扇贝的样子向Firefox图标靠近,我们要给它们加上一点选择压力,就是文章开头故事中提到的那个人的行动:在每一代把最不像Firefox的扇贝淘汰出去,同时也给新的个体留下生存的空间。怎么评价这个扇贝像不像Firefox呢?最直接的方法就是一个一个像素比较,颜色相差得越多就越不像。这个评价的函数叫做“适应函数”,它负责评价一个个体到底有多适应我们的要求。

在淘汰的过程中,为了便于编程,我们通常会在淘汰旧个体和产生新个体的数目上进行适当的调整,使种群的大小保持不变。淘汰的作用就是使适应我们要求的个体存在的时间更长,从而达到选择的目的。

最后,在自然界中,种群的演化是一个无休止的过程,但程序总要停下来给出一个结果。那么,什么时候终止演化输出结果呢?这就要订立一个终止条件,满足这个条件的话程序就输出当前最好的结果并停止。最简单的终止条件就是,如果种群经过了很多代(这里的“很多”是一个需要设定的参数)之后仍然没有显著改变适应性的变异的话,我们就停止并输出结果。我们就用这个终止条件。

好了,现在是万事俱备只欠东风了。定义好基因,写好繁衍、变异、评价适应性、淘汰和终止的代码之后,只需要随机产生一个适当大小的种群,然后让它这样一代代的繁衍、变异和淘汰下去,到最后终止我们就会获得一个惊喜的结果:(这回是完整的了,图片下的数字表示这个扇贝是第几代中最好的)

怎么样?虽说细节上很欠缺,但是也算是不错了。要不,你来试试用100个透明三角形画一个更像的?就是这样的看上去很简单的模拟演化过程也能解决一些我们这些有智慧的人类也感到棘手的问题。

实际上,在生活和生产中,很多时候并不需要得到一个完美的答案;而很多问题如果要得到完美的答案的话,需要很大量的计算。所以,因为遗传算法能在相对较短的时间内给出一个足够好能凑合的答案,它从问世伊始就越来越受到大家的重视,对它的研究也是方兴未艾。当然,它也有缺点,比如说早期的优势基因可能会很快通过交换基因的途径散播到整个种群中,这样有可能导致早熟(premature),也就是说整个种群的基因过早同一化,得不到足够好的结果。这个问题是难以完全避免的。

其实,通过微调参数和繁衍、变异、淘汰、终止的代码,我们有可能得到更有效的算法。遗传算法只是一个框架,里边具体内容可以根据实际问题进行调整,这也是它能在许多问题上派上用场的一个原因。像这样可以适应很多问题的算法还有模拟退火算法、粒子群算法、蚁群算法、禁忌搜索等等,统称为元启发式算法(Meta-heuristic algorithms)。

另外,基于自然演化过程的算法除了在这里说到的遗传算法以外,还有更广泛的群体遗传算法和遗传编程等,它们能解决很多棘手的问题。这也从一个侧面说明,我们不一定需要一个智能才能得到一个构造精巧的系统。

无论如何,如果我们要将遗传算法的发明归功于一个人的话,我会将它归功于达尔文,进化论的奠基人。如果我们不知道自然演化的过程,我们也不可能在电脑中模拟它,更不用说将它应用于实际了。

向达尔文致敬!

100个半透明三角形组成的达尔文肖像

Bonus:如果你看明白了这篇文章,而且你懂英语的话,你自然明白下面这个冷笑话:(来自http://xkcd.com/534/

Just to make sure you don't have it maximize instead of minimize.

Just to make sure you don't have it maximize instead of minimize.

相关文章
  1. 关于交配的问题,如果只是将每个三角形作为固定的基因保存下来,然后彼此更换,会不会导致基因库不够丰富而快速陷入进化陷阱呢?

    如果将三角形的三个点颜色和透明度,作为5个基因保存下来并且相互更换呢?

    我自己做了个程序玩,但是效果不理想~苦恼啊~

  2. 好神奇,可能我的算法有问题,但是为何最开始生成的随机三角形,就和文中给出的图不太一样?而且图中的图案从第10代开始就有收敛的倾向,颜色还这么均匀。确认算法里面没有作弊么?

  3. Sorry for not being able to type Chinese. Genetic algorithm is nothing mysterious, simply an improvement to brutal force search with a biology-inspired heuristics, while the correctness of this heuristics is hard to prove. So I would like to think GA as a purely technical solution --- it can sometimes lead you to a solution much faster than brutal force search, but you may find the heuristics is arbitrary. This does not definitely mean a problem is really solved.

  4. Pingback: AC » Archive » genetic programming 多边形近似

  5. 以前不太了解遗传算法,刚看了看,感觉没什么,就是模拟自然界罢了,其实模拟自然界倒是容易,更困难的是如何根据电脑的优势,拥有更好的表现。

  6. 我论文做的就是人工免疫算法 克隆小生境算法

  7. 推荐《复杂》这本书,里面有讲到基因算法的由来。豆瓣网址http://book.douban.com/subject/1030898/

  8. Pingback: 科学松鼠会 » 人算不如天算·那一阵风告诉我们的

  9. 遗传算法是很好的算法,在经济管理中相当有用,
    也就是对最优化问题是一个很好的解决方法!
    现在统计中常用的实数编码遗传算法相当好!

  10. 另外,我一直弄不明白模拟退火算法的原理,
    为什么模拟退火算法可以得到全局最优呢?
    希望知道的朋友说一说!非常感谢

    • 模拟退火中有一定概率“跳”到其它位置,以一定概率覆盖住了所有解空间

      • 不知道这么想对不对:在一定意义上来说,这类模糊算法都能找到全局最优解,只不过由于变异率的设置导致一部分算法很难跳出局部最优解(遗传算法如果变异率太高就没有意义)

  11. 你这个文章也就学过编程的能看懂,外行看了肯定懵……

    遗传算法太难写了~

  12. 笑话是说终结者里面的那个“天网”。是说不要让进化算法最终演变成为人类的敌人……

  13. 能不能把源代码公开一下,自己算算玩玩,自己编了一个算的太慢

  14. 按照这种算法臆测。上帝是不是按照某以参照物,通过这种算法模式制造了我们还有这个世界?

  15. 唉唉唉,我写的程序效果不好呢?到底该用什么选择、杂交、变异算子呢?

  16. 唉唉唉,我写的程序效果不好呢?到底该用什么选择、杂交、变异算子呢?