可爱的 Python：SimPy 简化了复杂模型

　原文http://www.fzs8.net/python/2007-03-26/1870.html

人们常常难以理解或预知实际系统的随机行为。有时可能精确论证系统的统计特性（诸如一般的、最差和最好的性能特性）。但在其它时候，只有实际运行（或模拟）系统时，具体设计的缺陷才会显现出来。本文中，David 讨论了 SimPy，它是一种 Python 包，允许您非常方便地创建离散事件系统的模型。

　　
　　在我遇到 SimPy 包的其中一位创始人 Klaus Miller 时，从他那里知道了这个包。Miller 博士阅读过几篇提出使用 Python 2.2+ 生成器实现半协同例程和“轻便”线程的技术的可爱的 Python 专栏文章。特别是（使我很高兴的是），他发现在用 Python 实现 Simula-67 样式模拟时，这些技术很有用。
　　
　　结果表明 Tony Vignaux 和 Chang Chui 以前曾创建了另一个 Python 库，它在概念上更接近于 Simscript，而且该库使用了标准线程技术，而不是我的半协同例程技术。该小组在一起研究时，认为基于生成器的样式更有效得多，并且最近在 SourceForge 上发起了使用 GPL 的项目，称为 SimPy（请参阅参考资料，获得 SimPy 主页的链接），目前处于 beta 测试版状态。Vignaux 教授希望他在惠灵顿维多利亚大学（University of Victoria）的将来大学教学中使用统一的 SimPy 包；我相信该库也非常适合应用到各类实用问题中。
　　
　　我承认在近期的通信交流和调查研究之前，我对编程领域的模拟方面没有任何基础知识。我猜想本专栏文章的大部分读者也和我一样，对这方面的知识知之甚少。尽管有人会认为这种样式编程的方式有些新奇，但在理解资源有限的实际系统的行为时，模拟是很有用的。不管您感兴趣的是有限带宽网络、汽车交通行为、市场和商业性优化、生物／进化的交互作用还是其它“随机”系统，SimPy 对这样的建模都提供了简单的 Python 工具。
　　
　　 随机的定义
　　与“连接”相类似，它是那些最适合形容其作业的词汇之一 — 再也找不到更适合的了：
　　随机（stochastic），源自希腊语 stokhastikos（形容词）
　　1）推测的、与推测相关的或者具有推测特点的；好推测的。
　　2）在统计学上：涉及或包含一个随机变量或多个随机变量，或涉及偶然性或概率。
　　
　　在本专栏文章中，我将一直使用食品杂货店内具有多条通道的付款区域这个相当简单的示例。通过使用所演示的模拟，我们可以根据对扫描器技术、购物者习惯、人员配备需求等进行的各种更改所产生的经济上和等待时间上的含义提出问题。这个建模的优点是在您对所做的更改产生的含义有清晰的想法时，它让您能提前制定策略。很明显，大多数读者不会专门经营一家食品杂货店，但这些技术可以广泛地应用于各类系统中。
　　
　　 模拟的概念
　　SimPy 库只提供了三个抽象／父类，并且它们对应于模拟的三个基本概念。有许多其它常规函数和常量用于控制模拟的运行，但重要的概念都与这些类结合在一起。
　　
　　模拟中的核心概念是进程。一个进程只是一个对象，它完成某些任务，随后在它准备完成下一个任务之前有时会等待一会儿。在 SimPy 中，您还可以“钝化”进程，这意味着在一个进程完成一个任务后，只有当其它进程要求该进程完成其它任务时，它才会去做。把进程当作尝试完成一个目标，常常是很有用的。在编写进程时，通常把它编写成可以在其中执行多个操作的循环。在每个操作之间，可以插入 Python“yield”语句，它让模拟调度程序在返回控制之前执行每个等待进程的操作。
　　
　　进程执行的许多操作取决于资源的使用。资源只是在可用性方面受到限制。在生物学模型中，资源可能是食物供应；在网络模型中，资源可以是路由器或有限带宽通道；在我们的市场模拟中，资源是付款通道。资源执行的唯一任务是在任何给定的时间内将它的使用限于一个特定的进程上。在 SimPy 编程模型下，进程单独决定它要保留资源的时间有多长，资源本身是被动的。在实际系统中，SimPy 模型可能适合概念性方案，也可能不适合；很容易想象到资源在本质上会限制其利用率（例如，如果服务器计算机在必需的时间帧内没有获得满意的响应，则它会中断连接）。但作为编程问题，进程或资源是否是“主动”方就不是特别重要（只要确保您理解了您的意图）。
　　
　　最后一个 SimPy 类是监控程序。实际上监控程序不是很重要，只不过它很方便。监控程序所做的全部任务就是记录向它报告的事件，并保存有关这些事件的统计信息（平均值、计数、方差等）。该库提供的 Monitor 类对记录模拟措施是个有用的工具，但您也可以通过您想使用的其它任何技术来记录事件。事实上，我的示例使 Monitor 子类化，以提供某些（稍微）增强的能力。
　　
　　 设置商店：对模拟编程
　　在我所撰写的大部分文章中，我都会马上给出样本应用程序，但在本例中，我认为带您经历食品杂货店应用程序的每个步骤会更有用。如果您愿意的话，可以把每个部分剪贴在一起；SimPy 创造者们将在将来的发行版中包含我的示例。
　　
　　SimPy 模拟中的第一步是几个常规的导入（import）语句：
　　
　　清单 1. 导入 SimPy 库 #!/usr/bin/env python
　　from __future__ import generators
　　from SimPy import Simulation
　　from SimPy.Simulation import hold, request, release, now
　　from SimPy.Monitor import Monitor
　　import random
　　from math import sqrt
　　
　　有些 SimPy 附带的示例使用 import * 样式，但我更喜欢使我填充的名称空间更清晰。对于 Python 2.2（SimPy 所需的最低版本），将需要如指出的那样，导入生成器特性。对于 Python 2.3 以后的版本，不需要这样做。
　　
　　对于我的应用程序，我定义了几个运行时常量，它们描述了在特定的模拟运行期间我感兴趣的几个方案。在我更改方案时，我必须在主脚本内编辑这些常量。要是这个应用程序的内容更充实，那么我就可能用命令行选项、环境变量或配置文件来配置这些参数。但就目前而言，这个样式已经足够了：
　　
　　清单 2. 配置模拟参数 AISLES = 5　　　　 # Number of open aisles
　　ITEMTIME = 0.1　　 # Time to ring up one item
　　AVGITEMS = 20　　　# Average number of items purchased
　　CLOSING = 60*12　　# Minutes from store open to store close
　　AVGCUST = 1500　　 # Average number of daily customers
　　RUNS = 10　　　　　# Number of times to run the simulation
　　
　　我们的模拟需要完成的主要任务是定义一个或多个进程。对于模拟食品杂货店，我们感兴趣的进程是在通道处付款的顾客。
　　
　　清单 3. 定义顾客的操作 class Customer(Simulation.Process):
　　def __init__(self):
　　Simulation.Process.__init__(self)
　　# Randomly pick how many items this customer is buying
　　self.items = 1 + int(random.expovariate(1.0/AVGITEMS))
　　def checkout(self):
　　start = now()　　　　　 # Customer decides to check out
　　yield request, self, checkout_aisle
　　at_checkout = now()　　 # Customer gets to front of line
　　waittime.tally(at_checkout-start)
　　yield hold, self, self.items*ITEMTIME
　　leaving = now()　　　　 # Customer completes purchase
　　checkouttime.tally(leaving-at_checkout)
　　yield release, self, checkout_aisle
　　
　　每位顾客已经决定采购一定数量的商品。（我们的模拟不涉及从食品杂货店通道上选择商品；顾客只是推着他们的手推车到达付款处。）我不能确定这里的指数变量分布确实是一个精确的模型。在其低端处我感觉是对的，但我感到对实际购物者究竟采购了多少商品的最高极限有点失实。在任何情况下，您可以看到如果可以使用更好的模型信息，则调整我们的模拟是多么简单。
　　
　　顾客采取的操作是我们所关注的。顾客的“执行方法”就是 .checkout()。这个进程方法通常被命名为 .run() 或 .execute()，但在我的示例中，.checkout() 似乎是最可描述的。您可以对它起任何您希望的名称。Customer 对象所采取的实际操作仅仅是检查几个点上的模拟时间，并将持续时间记录到 waittime 和 checkouttime 监控程序中。但在这些操作之间是至关重要的 yield 语句。在第一种情况中，顾客请求资源（付款通道）。只有当顾客获得了所需的资源之后，他们才能做其它操作。一旦来到付款通道，顾客实际上就在付款了 — 所花时间与所购商品的数量成比例。最后，经过付款处之后，顾客就释放资源，以便其他顾客可以使用它。
　　
　　上述代码定义了 Customer 类的操作，但我们需要在运行模拟之前，创建一些实际的顾客对象。我们可以为一天中将要购物的每位顾客生成顾客对象，并为每位顾客分配相应的付款时间。但更简洁的方法是“在每位顾客到商店时”，让工厂对象生成所需的顾客对象。实际上模拟并不会同时对一天内将要购物的所有顾客感兴趣，而是只对那些要同时争用付款通道的顾客感兴趣。注意：Customer_Factory 类本身是模拟的一部分 — 它是一个进程。尽管对于这个客户工厂，您可能联想到人造的机器工人（la Fritz Lang 的 Metropolis），但还是应该只把它看作编程的便利工具；它并不直接对应已建模域中的任何事物。
　　
　　清单 4. 生成顾客流 class Customer_Factory(Simulation.Process):
　　def run(self):
　　while 1:
　　c = Customer()
　　Simulation.activate(c, c.checkout())
　　arrival = random.expovariate(float(AVGCUST)/CLOSING)
　　yield hold, self, arrival
　　
　　正如我前面提到的，我想收集一些当前 SimPy Monitor 类没有解决的统计信息。也就是，我并不仅仅对平均付款时间感兴趣，而且还对给定方案中最