[R-es] Teorias de Cola
Olivier Nuñez
onunez en iberstat.es
Jue Ene 27 13:29:26 CET 2011
Por si te siver de ayuda para empezar, en su momento, redactŽ el
siguiente c—digo para la simulaci—n de una cola MM1.
lamb=10;mu=9 #valores de lambda y mu
N=1000 #nº total de clientes que entran en la cola
T=rexp(N,lamb) #tiempos entre dos llegadas
X=rexp(N,mu) #tiempos de servicio
T[1]=0;X[1]=0
Tlleg=cumsum(T) #Tiempos de llegada de los clientes
Tsal=rep(0,N)
for(n in 2:N) Tsal[n]=max(Tlleg[n],Tsal[n-1])+X[n] #instantes de
llegadas al sistema
#instantes de salidas del sistema con politica FCFS (First Come First
Serve)
Tmax=max(Tlleg) #instante de ultima llegada
instantes=c(Tlleg,Tsal[Tsal<Tmax]) #instantes de saltos hasta la
ultima llegada
eventos=data.frame(instantes=instantes,eventos=c(rep(1,N),rep(-1,
(length(instantes)-N))))
Eventos=eventos[order(eventos$instantes),] #eventos ordenados
Nt= cumsum(Eventos$eventos) #numero de clientes en el sistema en los
instantes de saltos
salida=cbind(Instantes=Eventos$instantes[-1],NClientes=Nt[-1])
plot(salida,type="l",xlab="Tiempo",ylab="Nœmero de clientes en el
sistema",
main="Evoluci—n del nœmero de clientes en el sistema a lo largo del
tiempo",col="blue4")
grid()
Un saludo. Olivier
--
____________________________________
Olivier G. Nu–ez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es
____________________________________
El 27/01/2011, a las 13:09, daniel pacheco gomez escribi—:
>
> Hola,
>
> ÀUna preguntilla, alguien ha trabajdo con R para teorias de colas
> modelos simples o mas complejos ? M/M/1 o M/M/c. Por si existe
> alguna libreria o script que pueda servir de ayuda o algo parecido.
>
> Saludos!
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
Más información sobre la lista de distribución R-help-es