Python 需要使用SimPy对列车系统建模的帮助吗

我需要帮助在SimPy中对火车系统（如地铁系统）进行建模，问题是我的“模拟”老师希望我使用Python+SimPy，而我不知道如何使用它，在阅读了一周SimPy之后，我设法理解了银行示例，这很简单，但是现在我仍然不知道如何模拟老师给我的问题

这个问题是大而广泛的，但我只需要一点帮助，而不是整个项目的完成，所以，如果有人可以帮助我，我会非常感激

帮助我入门的简化问题如下：

假设你有3个车站（A、B、C），相隔一定距离（比如100米），我有两列火车，一列在A上，一列在C上（1列从A到C，2列从C到A），每列火车都有一个最大速度（比如50米/秒）和加速度（比如5米/秒^2和-5米/秒^2）他们所要做的就是在每一站停留几秒钟（比如说24秒），然后继续到下一站，当他们到达终点时，他们再等待20秒（换轨），然后重新开始

车站的乘客数量是有限的（不需要模拟）。唯一的问题是每个车站都有随机数量的乘客，当火车到达时，一些乘客上车，另一些乘客下车

所以，基本上我需要一个小模型的火车和车站与这3个车站。。。如果有人能帮助我，请

---

真正的问题有22个车站，2个额外的换道站，随机等待时间，42列不同最高速度和容量的列车，根据车站和一天中的时间产生不同的乘客数量，等等。。。所有这些我以后都能处理，但火车站系统的建模逻辑我似乎无法理解。。。谢谢你的帮助

我绝对不是SimPy的专家，但这应该让你开始：

"""
Simulation of a train network
"""
from SimPy.Simulation import *
from math import sqrt
from random import randint
from itertools import cycle

def timeTo(A, maxV, d):
    """
    Given a trapezoidal velocity envelope defined by
    A       constant acceleration, m/s/s
    maxV    maximumum velocity, m/s
    return time in seconds required to travel
    d       distance, m
    """
    tA = float(maxV)/A          # time to accelerate to maxV
    dA = A*tA*tA                # distance traveled during acceleration from 0 to maxV and back to 0
    if (d < dA):                # train never reaches full speed?
        return sqrt(4.0*d/A)        # time needed to accelerate to half-way point then decelerate to destination
    else:
        return 2*tA + (d-dA)/maxV   # time to accelerate to maxV plus travel at maxV plus decelerate to destination

class Train(Process):
    def __init__(self, name, sim, accel=1.0, maxV=50.0, passengers=0, maxPassengers=400):
        super(Train, self).__init__(name, sim)
        self.accel = accel
        self.maxV  = maxV
        self.p     = passengers
        self.maxP  = maxPassengers

    def roll(self, route):
        here = route.next()     # starting location
        for dest in route:
            # travel to next station
            print "{:.1f}s: {} leaving {} for {}".format(self.sim.now(), self.name, here, dest)
            yield hold, self, timeTo(self.accel, self.maxV, here.distanceTo[dest])
            # arrive at next station
            here = dest
            print "{:.1f}s: {} at {}".format(self.sim.now(), self.name, here)
            yield hold, self, here.arrive(self)

    def getOff(self, num):
        if self.p >= num:
            print "  {} passengers got off".format(num)
            self.p -= num
        else:
            num = self.p
            print "  train is empty - only {} passengers got off".format(num)
            self.p = 0

    def getOn(self, num):
        if (self.maxP is None) or (self.p + num <= self.maxP):
            print "  {} passengers got on".format(num)
            self.p += num
        else:
            num = self.maxp - self.p
            print "  train is full - only {} passengers got on".format(num)
            self.p = self.maxp

class TrackNode(object):
    def __init__(self, name, delay=5.0):
        self.name = name
        self.delay = delay
        self.distanceTo = {}
    def arrive(self, train):
        pass
    def __str__(self):
        return self.name

class Station(TrackNode):
    def arrive(self, train):
        train.getOff(randint(1,15))
        train.getOn(randint(1,15))
        return self.delay

class Switch(TrackNode):
    def arrive(self, train):
        print("  switching tracks")
        return self.delay

class SampleRailroad(Simulation):
    def run(self, maxTime=100.0):
        self.initialize()
        # create places
        x = Switch("switch x", 20.0)
        A = Station("Station A", 24.0)
        B = Station("Station B", 27.0)
        C = Station("Station C", 25.0)
        y = Switch("switch y", 18.0)
        # distances between places
        x.distanceTo[A] = 50.0
        A.distanceTo[B] = 5000.0
        B.distanceTo[C] = 2000.0
        C.distanceTo[y] = 80.0
        y.distanceTo[C] = 80.0
        C.distanceTo[B] = 2000.0
        B.distanceTo[A] = 5000.0
        A.distanceTo[x] = 50.0
        # set up first train
        sf = Train("Santa Fe 219", self)
        self.activate(sf, sf.roll(cycle([A,B,C,y,C,B,A,x])), at=0.0)
        # set up second train
        cn = Train("Canadian National 41", self, maxPassengers=200)
        self.activate(cn, cn.roll(cycle([C,B,A,x,A,B,C,y])), at=5.0)
        # start simulating!
        self.simulate(until=maxTime)

def main():
    rr = SampleRailroad()
    rr.run(800.0)

if __name__=="__main__":
    main()

“”“
列车网络仿真
"""
从SimPy.Simulation导入*
从数学导入sqrt
从随机导入randint
从itertools导入周期
def timeTo（A、maxV、d）：
"""
给定梯形速度包络线，由
恒定加速度，m/s/s
最大速度最大速度，m/s
行驶所需的返回时间（秒）
d距离，m
"""
tA=浮动（最大值）/A#加速到最大值的时间
dA=A*tA*tA#从0加速到maxV再加速到0时行驶的距离
如果（d=num：
打印“{}乘客下车”。格式（num）
self.p-=num
其他：
num=self.p
打印“列车空车-只有{}名乘客下车”。格式（num）
self.p=0
def getOn（self，num）：
如果（self.maxP为None）或（self.p+num），我创建了一个框架，用于使用simpy简化这些类型的模拟。请参阅：。