[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