Marzo2013 Prof.RománLaraC.
TeoríadeColas
EsunaramadelasMatemáticasqueseocupadel estudiodelas«colasdeespera», estudiodelas«colasdeespera»,conunalcancemás conunalcancemás generalquelossistemasdecomun generalquelossistemasdecomunicaciones icaciones.. Suobjetivoeselestudiodesistemas Suobjetivoeselestudiodesistemascompuestos,po compuestos,por r unaomásunidadesllamadasprocesadoreso unaomásunidadesllamadas procesadoreso servidores,encargadosderealizarlastareas servidores,encargadosderealizarlastareas encomendadasporotrasunidades,llamadasclientes encomendadasporotrasunidades,llamadas clientes,, conlaparticularidaddequesiduran conlaparticularidaddequesidurantealgúninterva tealgúnintervalo lo detiempolallegadadeclientessuperalacapacidad deprocesamientodelsistema,dichosclientes permanecenencola permanecenencolahastaqueseanservidos. hastaqueseanservidos.
TeoríadeColas
EsunaramadelasMatemáticasqueseocupadel estudiodelas«colasdeespera», estudiodelas«colasdeespera»,conunalcancemás conunalcancemás generalquelossistemasdecomun generalquelossistemasdecomunicaciones icaciones.. Suobjetivoeselestudiodesistemas Suobjetivoeselestudiodesistemascompuestos,po compuestos,por r unaomásunidadesllamadasprocesadoreso unaomásunidadesllamadas procesadoreso servidores,encargadosderealizarlastareas servidores,encargadosderealizarlastareas encomendadasporotrasunidades,llamadasclientes encomendadasporotrasunidades,llamadas clientes,, conlaparticularidaddequesiduran conlaparticularidaddequesidurantealgúninterva tealgúnintervalo lo detiempolallegadadeclientessuperalacapacidad deprocesamientodelsistema,dichosclientes permanecenencola permanecenencolahastaqueseanservidos. hastaqueseanservidos.
TeoríadeColas Recursos
Clientes
Fuentesdetráfico
f T r á i c
o
de
lared
Servidores
Aplicacióndelateoríadecolas
Loscamposdeaplicacióndelateoríadecolasenla IngenieríadeTelecomunicaciónsonnumerosos, pudiendoagruparseenlassiguientesgrandesáreas:
Análisisydiseñodeprotocolosdecomunicación. Controlygestiónderecursosenredesde telecomunicaciones. Dimensionadoderedesdecomunicación. Modeladodeltráficoquecirculaporlasredes:voz, videoydatos.
Lascolas… Lascolassonfrecuentesennuestra vidacotidiana: Enunbanco Enunrestaurantedecomidas rápidas Almatricularenlauniversidad Losautosenunlavacar
Lascolas…
Engeneral,anadielegustaesperar Cuandolapacienciallegaasulímite, lagentesevaaotrolugar Sinembargo,unserviciomuyrápido tendríauncostomuyelevado Esnecesarioencontrarunbalance adecuado
Modelobásico
Unsistemadecolaspuededividirseen doscomponentesprincipales: Lacola Lainstalacióndelservicio Losclientesollegadasvienenenforma individualpararecibirelservicio
Modelobásico
Losclientesollegadaspuedenser: Personas Automóviles Máquinasquerequierenreparación Documentos Entremuchosotrostiposde artículos
Modelobásico
Sicuandoelclientelleganohaynadie enlacola,pasadeunavezarecibirel servicio Sino,seunealacola Esimportanteseñalarquelacolano incluyeaquienestárecibiendoel servicio
Modelobásico Lasllegadasvanalainstalacióndel serviciodeacuerdoconladisciplinade lacola Generalmenteéstaes primeroenllegar, primeroenserservido(FIFO) Peropuedenhaberotrasreglasocolas conprioridades
Modelobásico Sistema de colas
Llegadas
Cola
Disciplina de la cola
Instalación Salidas delservicio
Estructuras:picasdesistemasde colasunalínea,unservidor Sistema de colas
Llegadas
Cola
Servidor
Salidas
Estructuras:picasdesistemasdecolasuna línea,múl@plesservidores Sistema de colas
Servidor Llegadas
Cola
Servidor Servidor
Salidas Salidas Salidas
Estructuras:picasdecolasvarias líneas,múl@plesservidores Sistema de colas
Cola Llegadas
Cola Cola
Servidor Servidor Servidor
Salidas Salidas Salidas
Estructuras:picasdecolasunalínea, servidoressecuenciales Sistema de colas Llegadas
Cola Servidor Cola Servidor
Salidas
Costosdeunsistemadecolas 1. Costodeespera:Eselcostoparael
clientealesperar
Representaelcostodeoportunidad deltiempoperdido Unsistemaconunbajocostode esperaesunafuenteimportantede competitividad
Costosdeunsistemadecolas 2. Costodeservicio:Eselcostode
operacióndelserviciobrindado
Esmásfácildeestimar Elobjetivodeunsistemadecolases encontrarelsistemadelcostototal mínimo
Costosdeunsistemadecolas
Eltiempoquetranscurreentredos llegadassucesivasenelsistemade colassellamatiempoentrellegadas Eltiempoentrellegadastiendeaser muyvariable Elnúmeroesperadodellegadaspor unidaddetiemposellamatasamedia dellegadas(λ)
SistemasdecolasLasllegadas
Eltiempoesperadoentrellegadases1/ λ
Porejemplo,silatasamediade llegadasesλ=20clientesporhora Entonceseltiempoesperadoentre llegadases1/λ=1/20=0.05horaso3 minutos
SistemasdecolasLasllegadas Ademásesnecesarioestimarla distribucióndeprobabilidaddelos tiemposentrellegadas Generalmentesesuponeuna distribuciónexponencial Estodependedelcomportamientode lasllegadas
SistemasdecolasLasllegadas– Distribuciónexponencial
Laformaalgebraicadeladistribución exponenciales:
P (tiempo de servicio ≤ t ) 1 − e =
t
− µ
Dondet representaunacantidad expresadaendetiempounidadesde tiempo(horas,minutos,etc.)
SistemasdecolasLasllegadas– Distribuciónexponencial P(t)
0
Media
Tiempo
SistemasdecolasLasllegadas– Distribuciónexponencial
Ladistribuciónexponencialsupone unamayorprobabilidadparatiempos entrellegadaspequeños Engeneral,seconsideraquelas llegadassonaleatorias Laúltimallegadanoinfluyeenla probabilidaddellegadadelasiguiente
SistemasdecolasLasllegadas-Distribuciónde Poisson
Esunadistribucióndiscretaempleada conmuchafrecuenciaparadescribirel patróndelasllegadasaunsistemade colas Paratasasmediasdellegadas pequeñasesasimétricaysehacemás simétricayseaproximaalabinomial paratasasdellegadasaltas
SistemasdecolasLasllegadas- DistribucióndePoisson
Suformaalgebraicaes: k
P (k )
λ e =
− λ
k ! Donde: P(k):probabilidaddekllegadaspor unidaddetiempo λ:tasamediadellegadas e=2,7182818…
SistemasdecolasLasllegadas- DistribucióndePoisson P
0
Llegadas por unidad de tiempo
SistemasdecolasLacola
Elnúmerodeclientesenlacolaesel númerodeclientesqueesperanel servicio Elnúmerodeclientesenelsistemaes elnúmerodeclientesqueesperanen lacolamáselnúmerodeclientesque actualmenterecibenelservicio
SistemasdecolasLacola Lacapacidaddelacolaeselnúmero máximodeclientesquepuedenestar enlacola Generalmentesesuponequelacolaes infinita Aunquetambiénlacolapuedeser finita
SistemasdecolasLacola
Ladisciplinadelacolaserefiereal ordenenqueseseleccionanlos miembrosdelacolaparacomenzarel servicio LamáscomúnesFIFO:primeroen llegar,primeroensalir Puededarse:selecciónaleatoria, prioridades,LIFO,entreotras.
SistemasdecolasElservicio
Elserviciopuedeserbrindadoporun servidoroporservidoresmúltiples Eltiempodeserviciovaríadeclientea cliente Eltiempoesperadodeservicio dependedelatasamediadeservicio (µ)
SistemasdecolasElservicio
Eltiempoesperadodeservicio equivalea1/µ Porejemplo,silatasamediade servicioesde25clientesporhora Entonceseltiempoesperadode servicioes1/µ=1/25=0.04horas,o2.4 minutos
SistemasdecolasElservicio
Esnecesarioseleccionaruna distribucióndeprobabilidadparalos tiemposdeservicio Haydosdistribucionesque representaríanpuntosextremos: Ladistribuciónexponencial (σ=media) Tiemposdeservicioconstantes(σ=0)
SistemasdecolasElservicio
Unadistribuciónintermediaesla distribuciónErlang Estadistribuciónposeeunparámetro deformak quedeterminasu desviaciónestándar: σ
1 =
k
media
SistemasdecolasElservicio
Sik =1,entoncesladistribuciónErlang esigualalaexponencial
Sik =∞,entoncesladistribución Erlangesigualaladistribución degeneradacontiemposconstantes LaformadeladistribuciónErlang varíadeacuerdoconk
SistemasdecolasElservicio P(t)
k=∞
k=8
k=2
k=1
0
Media
Tiempo
Sistemasdecolas DistribuciónErlang Distribución
Desviación estándar
Constante
0
Erlang, k = 1 Erlang, k = 2 Erlang, k = 4 Erlang, k = 8 Erlang, k = 16 Erlang, cualquier k
media 1/
2 media
1/2 media 1/
8 media
1/4 media 1/
k media
NotacióndeKendall
Utilizadaparaclasificarlossistemasespecificandolas característicasdeloselementosquelocomponen. A/ B/C/K/m/z
A:Distribucióndetiemposentrellegadas(τ) B:Distribucióndetiemposdeservicio(s) D(determinista),M(memory-less;exponencial),Ek(Erlang-k),Hk (hiperexponencialdekestados)oG(general;seaproximaráaunadelas anterioressegún↓↑ 2 (0,1,1/ko≥1)). :Númerodeservidores(canales) C
K =capacidadtotalmáximadelsistema(usuariosencola+ servidores).
m=tamañodelapoblación.
SiK=C=>noexistecola. Pordefecto,infinito. Pordefecto,infinito.
z=disciplinadelacola.
Estadodelsistemadecolas
Enprincipioelsistemaestáenun estadoinicial Sesuponequeelsistemadecolasllega aunacondicióndeestadoestable (nivelnormaldeoperación) Existenotrascondicionesanormales (horaspico,etc.) Loqueinteresaeselestadoestable
Desempeñodelsistemadecolas
Paraevaluareldesempeñosebusca conocerdosfactoresprincipales: 1. Elnúmerodeclientesqueesperan enlacola 2. Eltiempoquelosclientesesperan enlacolayenelsistema
ProcesosdePoisson
Tambiénllamadosprocesostotalmentealeatorios,modelande formaadecuadalallegadadeusuariosasistemasreales.
Otraopción=>Procesosautosimilares:característicassimilares endistintosinstantesdetiempo.
Másajustadosalarealidadperodepeormanejomatemático.
Características:
P(llegadadeusuarioent)≠funcióndellegadasanteriores.
Probabilidaddellegadaenunintervalodirectamenteproporcional alalongituddeéste. Probabilidaddemásdeunallegadaenunintervalolo suficientementepequeñoesdespreciable. Lallegadaenunintervaloesindependientedellegadaspasadaso futuras.
CasoparticulardeprocesodeMarkov:
Probabilidaddesiguienteestadosólodependedelestadoactualy nodelahistoria.
Medidasdeldesempeñodelsistemade colas
Deinterésparalaexplotacióndelsistemaconelmáximo beneficioylamenorinversión.
IntensidaddeTráfico= A(=I )[Erlangs]
Demandaderecursosrealizadaporlosusuarios. Utilizacióndelosrecursosdesplegados. InterpretaciónErlang=númerodecanalespermanentemente ocupadosnecesariosparacursartodoeltráfico. AO =tráficoofrecido A perdido =tráficonocursado AC =tráficocursado = /
Factordeutilización= ρ
Probabilidaddequeunservidorestéocupadooporcentajede tiempoenqueelservidorestáocupado. λ’=tasaefectiva(cursada)dellegada.
Parámetrosorientadosalsistema
Throughput(caudal)=Th[usuarios/segundo]
Medidadelaproductividaddelsistema. Númeromediodeusuariosservidosporunidaddetiempo.
Sinpérdidas=>Th=λ Thmax=μC
Th= µ
Volumendetráficocursadoporunservidor=V
Tiempototaldeocupacióndeeseservidorenunintervalodetiempode
referencia(T ).=>V≤T
Tambiénsepuededefinir A=V/T . Varíaalolargodeldía=>franjashorariasypromedioenvariosdías. HC=HoraCargada=horadeldíaconmayortráfico.
Apartirdeltráficoenestahoraserealizalaplanificación.
Unidades:
Devolumen:[LLR](llamadade120seg)y[CCS](llamadade100seg).
Parámetrosorientadosalsistema
MedidadelaQoSpercibidaporelusuario. Tiempomediodeesperaencola=W
Tiempomedioenelsistema=T
Esunav.a.yvaríaconeltiempo.
Númeromediodeusuariosencola=N =λW
Esunav.a.yvaríaparacadausuario.Ti= W i+Si Régimenpermanente=>procesosestocásticosestacionarios.
Númeromediodeusuariosenelsistema=N =λT
Esunav.a.yvaríaparacadausuario.
Esunav.a.yvaríaconeltiempo.
Otros:
Probabilidaddequeexistaunservidorlibre. Probabilidaddequelacolasupereciertovalor. Tiempomediodeesperaparalosqueentranencola.
Medidasdeldesempeñodelsistemadecolas fórmulasgenerales
T = W + E [ S ] = W +
1 µ
N = λ T N q
=
λ W
N = N q
λ +
µ
Medidasdeldesempeñodelsistemade colasejemplo
Supongaunaestacióndegasolinaalacual lleganenpromedio45clientesporhora Setienecapacidadparaatenderen promedioa60clientesporhora Sesabequelosclientesesperanen promedio3minutosenlacola Determinarlosparámetrostotalesdel sistema
Medidasdeldesempeñodelsistemade colasejercicio
Supongaunrestaurantdecomidas rápidasalcuallleganenpromedio100 clientesporhora Setienecapacidadparaatenderen promedioa150clientesporhora Sesabequelosclientesesperanen promedio2minutosenlacola Calculelasmedidasdedesempeñodel sistema
Probabilidadescomomedidasdel desempeño
Beneficios: Permitenevaluarescenarios Permiteestablecermetas Notación: P n:probabilidaddetenernclientes enelsistema P (T s ≤ t ) : probabilidad de que un
cliente no espere en el sistema más de t horas
ac or e u sistema
zac n e
Dadalatasamediadellegadasλyla tasamediadeservicioµ,sedefineel factordeutilizacióndelsistemaρ. Generalmenteserequierequeρ<1 Sufórmula,conunservidoryconC servidores,respectivamente,es:
ρ
λ =
µ
ρ
λ
`
=
C µ
Factordeu@lizacióndelsistema-ejemplo
Conbaseenlosdatosdelejemplo anterior,λ=0.75,µ=1 Elfactordeutilizacióndelsistemasise mantuvieraunservidores ρ=λ/µ=0.75/1=0.75=75% Condosservidores(C =2): ρ=λ/C µ=0.75/(2*1)=0.75/2=37,5%
Ejemplo
Unaparatoregistradordetráficotomamedidascada3minutos, durantelaHC,delnúmerodecircuitosocupadosenungrupo. Lasmedidasobtenidasserepresentanenelsiguientegráfico:
Ejemplo
Enunhazdecuatrocircuitos,cadaunoestá ocupadouncuartodehoradiferentedelaHora Cargada. a)¿Cuáleseltráficocursadoporcadacircuito? ¿Yporelhaz? b)¿Ysicoincidenloscuatrocuartosdehora?
Análisiseconómicodelíneasdeespera Costos
Costo total
Costo del servicio
Costo de espera Tasa óptima
Tasa de servicio
Procesosde NacimientoyMuerte
Resultaninteresantesporquelamayorpartedelossistemasde esperacontiemposdellegadaydeservicioexponencialesse puedenmodelarcomoprocesosdenacimientoymuerte. SonuncasoespecialdelosprocesosdeMarkovdondesólose realizan transicionesaestadosadyacentes. Estadodelsistema:númerodeelementosdelsistema. Laevoluciónentreestadosdelsistemasólodependedelestado actual. Transiciones Gráficamente:
Estado“n”delsisteman Eventos
n
n
Nacimiento:llegadadeunelementoalsistema Muerte:salidadeunelementodelsistema
Consideraremosquelasllegadasysalidassonindependientesentresí
Procesosde NacimientoyMuerte
Latransiciónentreestadostienelugarconunadeterminadaprobabilidad Probabilidad de que suceda un nacimiento es la probabilidad de que, estando el sistema λ n-1 en el estado n-1,pasealestadon n-1
μn
n
Probabilidad de que suceda una muerte es la probabilidad de que, estando el sistema en el estado n,pasealestadon-1
Ecuacióndeconservacióndelflujo(n>0)
Enrégimenpermanente,elflujodeentradaysalidadecadaestadocoinciden
λn-1 pn-1= μn pnè pn=λn−1 pn−1/µn
Procesosde NacimientoyMuerte λ0
λ1
λ2
pn=λn−1 pn−1/µn 0
1
2
μ 1μ 2μ 3 Paraobtener pn enfunciónde p0:
yp ?Sielsistemaesestable
p 1=λ0 /µ1 p0
∑ ≥0↑▒pn =1 p 2=λ1 /µ2 p1 =λ1λ0/µ1µ2 p0
0
SistemasconPérdidas SinReintento
SistemasconPérdidasySinReintento: Sistemasincola. Silademandasuperalacapacidaddeservicio,laspeticiones serechazanypierden. èLaspeticionesquelleganalsistemaynoencuentranun servidorlibre,serechazan.
Lasllamadasquenosepuedencursarinmediatamente,se rechazan.Sesuponequenohayreintento(¿otro sistema?). Estructuradeunsistemaconpérdidas El modelo de tráfico que estudia esta situación es el modelo de Erlang B (ErB)
SistemasconPérdidas SinReintento
HipótesisdelmodeloErlang-B: 1.Tamañodelapoblaciónesinfinito. Característica: latasadepeticionesdeservicionoseve afectada por el estado del sistema. 2.RégimendellegadasdePoisson. 3.Lav.a.tiempodeservicio(s)sigueunadistribución exponencial. 4.Númerodeservidores:c. 5.Nohayespera(Q=0). Notación Kendall: M/M/c/c/∞
Sesuponequelasunidades,una vezservidas,regresanala población
SistemasconPérdidas SinReintento
SistemasconPérdidas SinReintento
SistemasconPérdidas SinReintento
Paraconectarcuatrogruposdeterminalesdedatosaunordenadorcentralseproponen dosconfiguraciones,representadasenlasfiguras(a)y(b).
Sabiendoquecadagrupotiene22terminalesyque,portérminomedio,estánactivosel 10%deltiempo,determineelnúmerodecircuitosquesenecesitanencadacasosi probabilidaddebloqueomáximaesdel5%. (Nota:Modeleelsistemadeformaquelasllamadasbloqueadassepierden.)
SistemasconPérdidas SinReintento
SistemasconPérdidas ConReintento
SistemasconPérdidasyConReintento:
Elanálisissimplificadodeestossistemassebasaentres hipótesis:
Silapeticióndeconexiónnosecursa,lallamadavuelvealsistema comounreintento. H1)Todaslasllamadasqueelsistemarechazacuandoestá bloqueadosecursanenposterioresreintentos. H2)Eltiempoquetranscurreentreelinstanteenqueunapetición encuentrabloqueoyelreintentoesaleatorioyestadísticamente independiente. H3)Eltiempomedioentrereintentosesmayorqueeltiempo mediodeservicio.
Desdeelp.d.v.analítico,distinguiremosentrellamadade primerintentoyllamadadereintento.
SistemasconPérdidas ConReintento
SistemasconPérdidas ConReintento
SistemasconPérdidas ConReintento
SistemasconPérdidas ConReintento Determinelaprobabilidaddebloqueodeunenlacede10circuitosentreuna centralitayunacentrallocalsabiendoqueeltráficoofrecido,sisepudieracursar ensutotalidad,esde7Erlangs. Nota:Apliqueelmodelodereintentosconllamadasperdidasysuponga poblacióninfinita.
SistemasconPérdidas ConReintento
SistemasconPérdidas ConReintento
Sistemasconespera
Sistemasconespera
Sistemasconespera
Sistemasconespera
Sistemasconespera
Sistemasconespera
Sistemasconespera Unacentralitasirvedepuenteentre100 Unacentralitasirvedepuenteentre100extensionesy20 extensionesy20 líneasde salidaconlassiguientescaracterísticas: Encasodesaturación,lasextensionesesperaneltiempo queseanecesarioa queunadelaslíneasdesalidaquedelibre. Cadaextensióngenera3llamadasdurantelaHC. Lamediadeduracióndelasllamadasesde2 minutos. Calcule: PD. PD. Tiempomediodeespera.
Sistemasconespera
Sistemasconespera Unacentralitasirvedepuenteentre100extensionesy20 líneasde salidaconlassiguientescaracterísticas: Encasodesaturación,lasextensionesesperaneltiempo queseanecesarioa queunadelaslíneasdesalidaquedelibre. Cadaextensióngenera3llamadasdurantelaHC. Lamediadeduracióndelasllamadasesde2 minutos. Calcule: PD. Tiempomediodeespera.