Perché questa affermazione ricorsiva è sbagliata?
Domanda
Questa è una simulazione bancaria che tiene conto di 20 diverse linee di servizio con una singola coda, i clienti arrivano seguendo un tasso esponenziale e vengono serviti durante un periodo che segue una normale distribuzione di probabilità con media 40 e deviazione standard 20.
Le cose funzionavano bene fino a quando ho deciso di escludere i valori negativi dati dalla distribuzione normale usando questo metodo:
def getNormal(self):
normal = normalvariate(40,20)
if (normal>=1):
return normal
else:
getNormal(self)
Sto rovinando la chiamata ricorsiva? Non capisco perché non funzionerebbe. Ho cambiato il metodo getNormal () in:
def getNormal(self):
normal = normalvariate(40,20)
while (normal <=1):
normal = normalvariate (40,20)
return normal
Ma sono curioso di sapere perché la precedente dichiarazione ricorsiva viene sballata.
Questo è il codice sorgente completo, nel caso tu sia interessato.
""" bank21: One counter with impatient customers """
from SimPy.SimulationTrace import *
from random import *
## Model components ------------------------
class Source(Process):
""" Source generates customers randomly """
def generate(self,number):
for i in range(number):
c = Customer(name = "Customer%02d"%(i,))
activate(c,c.visit(tiempoDeUso=15.0))
validateTime=now()
if validateTime<=600:
interval = getLambda(self)
t = expovariate(interval)
yield hold,self,t #esta es la rata de generación
else:
detenerGeneracion=999
yield hold,self,detenerGeneracion
class Customer(Process):
""" Customer arrives, is served and leaves """
def visit(self,tiempoDeUso=0):
arrive = now() # arrival time
print "%8.3f %s: Here I am "%(now(),self.name)
yield (request,self,counter),(hold,self,maxWaitTime)
wait = now()-arrive # waiting time
if self.acquired(counter):
print "%8.3f %s: Waited %6.3f"%(now(),self.name,wait)
tiempoDeUso=getNormal(self)
yield hold,self,tiempoDeUso
yield release,self,counter
print "%8.3f %s: Completed"%(now(),self.name)
else:
print "%8.3f %s: Waited %6.3f. I am off"%(now(),self.name,wait)
## Experiment data -------------------------
maxTime = 60*10.5 # minutes
maxWaitTime = 12.0 # minutes. maximum time to wait
## Model ----------------------------------
def model():
global counter
#seed(98989)
counter = Resource(name="Las maquinas",capacity=20)
initialize()
source = Source('Source')
firstArrival= expovariate(20.0/60.0) #chequear el expovariate
activate(source,
source.generate(number=99999),at=firstArrival)
simulate(until=maxTime)
def getNormal(self):
normal = normalvariate(40,20)
if (normal>=1):
return normal
else:
getNormal(self)
def getLambda (self):
actualTime=now()
if (actualTime <=60):
return 20.0/60.0
if (actualTime>60)and (actualTime<=120):
return 25.0/60.0
if (actualTime>120)and (actualTime<=180):
return 40.0/60.0
if (actualTime>180)and (actualTime<=240):
return 30.0/60.0
if (actualTime>240)and (actualTime<=300):
return 35.0/60.0
if (actualTime>300)and (actualTime<=360):
return 42.0/60.0
if (actualTime>360)and (actualTime<=420):
return 50.0/60.0
if (actualTime>420)and (actualTime<=480):
return 55.0/60.0
if (actualTime>480)and (actualTime<=540):
return 45.0/60.0
if (actualTime>540)and (actualTime<=600):
return 10.0/60.0
## Experiment ----------------------------------
model()
Soluzione
Penso che tu voglia
return getnormal(self)
anziché
getnormal(self)
Se la funzione termina senza premere un'istruzione return, restituisce il valore speciale None, che è un oggetto NoneType - ecco perché Python si lamenta di un 'NoneType'. La funzione abs () vuole un numero e non sa cosa fare con un Nessuno.
Inoltre, potresti evitare la ricorsione (e il costo di creazione di un nuovo stack frame) usando
def getNormal(self):
normal = 0
while normal < 1:
normal = normalvariate(40,20)
return normal
Altri suggerimenti
Non ne sono del tutto sicuro, ma penso che tu debba cambiare il tuo metodo nel modo seguente:
def getNormal(self):
normal = normalvariate(40,20)
if (normal>=1):
return normal
else:
return getNormal(self)
Devi avere:
return getNormal(self)
anziché
getNormal(self)
Davvero, non c'è bisogno di ricorsione:
def getNormal(self):
normal = 0
while normal < 1:
normal = normalvariate(40,20)
return normal