ejercicio 4.6 (colas)

174 days ago by roschaefer

Usando R, simular los fallos de la máquina del problema 4.5 y mostrar dentro de la simulación, el componente que falla. La suposición es que después de una reparación, la máquina tiene igual probablidad a la máquina antes de su fracaso. problema 4.5) existen dos tipos comunes de faltas de un elemento crítico de algunas máquinas electrónicas: ya sea del componente A o del componente B pueden fallar. Si cualquiera de los componentes falla, la máquina se cae. El componente A falla de acuerdo con un proceso de Poisson con tasa media de 3,3 fallos por dia. (La compañía opera 24 / 7 con turnos de ocho horas.) El componente B falla de acuerdo con un proceso de Poisson con una tasa media de 1,2 fallos por dia. a.-Supongamos que el tiempo de reparar la máquina es insignificante en comparación con el momento en el que fracasa. 
       
T=100 t=0 i=0 j=0 k=0 A<-c() B<-c() TiempoT<-c(0) exp.A<-rexp(1,3.3) exp.B<-rexp(1,1.2) while(t<T) { if(exp.A<exp.B) { A [i]<- (t + exp.A) exp.A<-rexp(1,3.3) exp.B<-rexp(1,1.2) t = t + exp.A i = i + 1 } else { B [j]<- (t + exp.B) exp.A<-rexp(1,3.3) exp.B<-rexp(1,1.2) t = t + exp.B j = j + 1 } TiempoT [k] <- t k = k + 1 } cat("las fallas causadas por la componente A es de",i," veces","\n") cat("las fallas causadas por la componente B es de ",j," veces","\n") cat ("La maquina se daño",k,"veces.","\n") 
       
las fallas causadas por la componente A es de 204  veces 
las fallas causadas por la componente B es de  63  veces 
La maquina se daño 267 veces. 
las fallas causadas por la componente A es de 204  veces 
las fallas causadas por la componente B es de  63  veces 
La maquina se daño 267 veces. 
b.- Supongamos que el tiempo para reparar la máquina se distribuye uniformemente entre 5 minutos y 60 minutos 
       
T=100 t=0 i=0 j=0 k=0 A<-c() B<-c() TiempoT<-c(0) exp.A<-rexp(1,3.3) exp.B<-rexp(1,1.2) tarda<-runif(1,5/(60*24),1/24) while(t<T) { if(exp.A<exp.B) { A [i]<- (t + exp.A) tarda<-runif(1,5/(60*24),1/24) exp.A<-rexp(1,3.3) exp.B<-rexp(1,1.2) t = t + exp.A + tarda i = i + 1 } else { B [j]<- (t + exp.B) exp.A<-rexp(1,3.3) exp.B<-rexp(1,1.2) tarda<-runif(1,5/(60*24),1/24) t = t + exp.B + tarda j = j + 1 } TiempoT [k] <- t k = k + 1 } cat("las fallas causadas por la componente A es de",i," veces","\n") cat("las fallas causadas por la componente B es de ",j," veces","\n") cat ("La maquina se daño",k,"veces.","\n") 
       
las fallas causadas por la componente A es de 144  veces 
las fallas causadas por la componente B es de  63  veces 
La maquina se daño 207 veces. 
las fallas causadas por la componente A es de 144  veces 
las fallas causadas por la componente B es de  63  veces 
La maquina se daño 207 veces.