示例#1
0
nosh.config.samples = 3000  # 2900
nosh.config.elites = 300  # 200
nosh.config.smoother = 0.5
nosh.config.iters = 51

mymeth = ("DP/LBH", "EMSR/NV", "EMSR/SL", "HCR/OSA", "HAR/OSA")  #'EMSR/HCR', 'EMSR/HAR',


def scenas(scen, ratios):
    f = scen.f
    rat = f[2] / float(f[1])
    for r in ratios:
        newscen = copy(scen)
        f1 = f[0] * r
        newscen.f = (f[0], f1, f1 * rat)
        newscen.nid = r
        yield newscen


pickle = "fareratio3.pkl"
gld = [s for s in scenas(sina, (i / 10.0 for i in range(3, 8)))]
# nosh.enuSim(gld, 10000, pickle, mymeth)
DISP = "ORfareratio3"
xlab = "fare ratio"
nosh.drawSubFigs(DISP, xlab, *nosh.loadResults(pickle), no_percentiles=True)
# nosh.drawFigs(DISP, xlab, *nosh.loadResults(pickle))
# nosh.drawPolicies(DISP,xlab,*nosh.loadPolicies(pickle))
if DISP is None:
    pylab.show()
示例#2
0
nosh.config.LBH=True
nosh.config.BINO=True
nosh.config.svl =0.001

nosh.config.samples = 3000
nosh.config.elites  = 300 
nosh.config.smoother = 0.7

#mymeth = ('OBSA/CR','EMSR/CR','EMSR/NV', 'CRSA/NV', 'DP/LBH', 'EMSR/SL')
mymeth = ('DP/LBH', 'EMSR/NV', 'EMSR/SL', 'HCR/OSA', 'HAR/OSA')

def scenas(scen, ps):
	for C in ps:
		newscen = copy(scen)
		newscen.C = newscen.nid = C 
		print "## LOAD Factor:", newscen.demandFactor()
		yield newscen

from ORinstance3 import sina as osin
sina = nosh.NormalScena(osin)

pickle = 'df3normal.pkl'
gld = [s for s in scenas(sina, range(100, 201, 10))]
nosh.enuSim(gld, 10000, pickle, mymeth)
DISP = 'ORdf3normal'
xlab = 'total number of seats'
nosh.drawSubFigs(DISP, xlab, *nosh.loadResults(pickle))
nosh.drawFigs(DISP, xlab, *nosh.loadResults(pickle))
nosh.drawPolicies(DISP,xlab,*nosh.loadPolicies(pickle))
if DISP is None: pylab.show()