Categoría: Investigación de Operaciones
Indice
Fernando Marrero Delgado. Doctor en Ciencias Técnicas. Máster en Informática Aplicada. Ingeniero Industrial. Profesor Auxiliar. Departamento de Ingeniería Industrial. Universidad Central de Las Villas. Santa Clara. Cuba. Email: fmarrero@fce.uclv.edu.cu
Javier Asencio García. Doctor en Ciencias Técnicas. Ingeniero Industrial. Profesor Titular. Departamento de Ingeniería Industrial. Universidad Central de Las Villas. Santa Clara. Cuba. Email: asencio@fce.uclv.edu.cu
René Abreu Ledón. Máster en Ingeniería Industrial. Ingeniero Industrial. Profesor Asistente Departamento de Ingeniería Industrial. Universidad Central de Las Villas. Santa Clara. Cuba. Email: rabreu@fce.uclv.edu.cu
René Orozco Sánchez. Ingeniero Industrial. Aspirante a Máster del Departamento de Ingeniería Industrial. Universidad Central de Las Villas. Santa Clara. Cuba. Email: fmarrero@fce.uclv.edu.cu
Hugo R. Granela Martín. Doctor en Ciencias Técnicas. Ingeniero Industrial. Profesor Auxiliar. Departamento de Ingeniería Industrial. Universidad Central de Las Villas. Santa Clara. Cuba. Email: hugran@fce.uclv.edu.cu
Teoría de Colas, Investigación de Operaciones, Optimización
Muchas industrias de productos y de servicios tienen un sistema de colas en el que los ¨productos¨ (o clientes) llegan a una ¨estación ¨ esperan en una ¨fila¨ (o cola), obtienen algún ¨servicio¨ y luego salen del sistema. Considere los siguientes ejemplos:
Los problemas administrativos relacionados con tales sistemas de colas se clasifican en dos grupos básicos:
Basándose en estas preguntas, los gerentes tomarán decisiones como emplear o no más gente, agregar una estación de trabajo adicional para mejorar el nivel de servicio, o si es necesario o no aumentar el tamaño del área de espera.
Estas decisiones de diseño se toman mediante la evaluación de los méritos de las diferentes alternativas, respondiendo a las preguntas de análisis del grupo 1 y luego seleccionando la alternativa que cumpla con los objetivos administrativos. En el presente capítulo se proporcionan las técnicas para analizar un sistema de colas dado. Sin embargo, las técnicas matemáticas específicas dependen de la clase de sistema a la cual pertenece su modelo de colas.
Características de un sistema de colas
Para analizar un sistema de colas, es mejor primero identificar las características importantes que aparecen en la siguiente sección características claves, y que se ilustran en la figura 1.1
Características claves
Las siguientes características se aplican a los sistemas de colas:
Una población de clientes, que es el conjunto de los clientes posibles.
Un proceso de llegada, que es la forma en que llegan los clientes de esa población
Un proceso de colas, que está conformado por (a) la manera que los clientes esperan para ser atendidos y (b) la disciplina de colas, que es la forma en que son elegidos para proporcionarles el servicio.
Un proceso de servicios, que es la forma y la rapidez con la que es atendido el cliente
Proceso de salida, que son de los siguientes dos tipos:
a. Los elementos abandonan completamente el sistema después de ser atendidos, lo que tiene como resultado un sistema de colas de un paso. Por ejemplo como se muestra en la figura 1.2 (a) los clientes de un banco esperan en una sola fila, son atendidos por uno de los tres cajeros y, después que son atendidos abandonan el sistema.
b. Los productos, ya que son procesados en una estación de trabajo, son trasladados a alguna otra parte para someterlos a otro tipo de proceso, lo que tiene como resultado una red de colas. Por ejemplo, los productos que se muestran en la figura 1.2 (b) primero son procesadas en la estación de trabajo A y después son enviadas a la estación de trabajo B o C. Los productos terminados en ambas estaciones, B y C, luego son procesados en la estación D, antes de abandonar el sistema.
En el presente capítulo solamente se considerarán sistemas de un paso Se necesitan diferentes análisis matemáticos para cada uno de estos dos tipos de procesos de salida.
La población de clientes
Al tomar en cuenta la base de clientes, la principal preocupación es el tamaño de la población. Para problemas como los de un banco o un supermercado, en donde el número de clientes potenciales es bastante grande (cientos de miles), el tamaño de la población se considera, para fines prácticos, como si fuera infinita.
Al contrario, considere una fábrica que tiene cuatro máquinas, que a menudo se descomponen y requieren servicio de reparación en un taller especializado. En este caso, es de solamente cuatro. El análisis de poblaciones finitas (es decir de tamaño limitado ) es más complicado que el análisis en donde la base de población se considera infinita.
El proceso de llegada
El proceso de llegada es la forma en que los clientes llegan a solicitar un servicio. La característica más importante del proceso es el tiempo entre llegadas, que es la cantidad de tiempo entre dos llegadas sucesivas. Este lapso es importante porque mientras menor sea el intervalo de tiempo, con más frecuencia llegan los clientes, lo que aumenta la demanda de servidores disponibles.
Características claves
Existen dos clases básicas de tiempo entre llegadas:
Determinístico, en el cual clientes sucesivos llegan en un mismo intervalo de tiempo, fijo y conocido. Un ejemplo clásico es el de una línea de ensamble, en donde los artículos llegan a una estación en intervalos invariables de tiempo (conocido como ciclos de tiempo)
Probabilístico, en el cual el tiempo entre llegadas sucesivas es incierto y variable. Los tiempos entre llegadas probabilísticos se describen mediante una distribución de probabilidad.
En el caso probabilístico, la determinación de la distribución real, a menudo, resulta difícil. Sin embargo, una distribución , la distribución exponencial, ha probado ser confiable en muchos de los problemas prácticos. La función de densidad, para una distribución exponencial depende de un parámetro, digamos l (letra griega lambda), y está dada por:
f(t)=(1/ l )e- l t
en donde l (lambda) es el número promedio de llegadas en una unidad de tiempo.
Con una cantidad, T, de tiempo usted puede hacer uso de la función de densidad para calcular la probabilidad de que el siguiente cliente llegue dentro de las siguientes T unidades a partir de la llegada anterior, de la manera siguiente:
P(tiempo entre llegadas <=T)=1-e- l t
El proceso de colas
Parte del proceso de colas tiene que ver con la forma en que los clientes esperan para ser atendidos. Los clientes pueden esperar en una sola fila, como en un banco, observe la figura 1.3(a) éste es un sistema de colas de una sola línea. Al contrario, los clientes pueden elegir una de varias filas en las que deben esperar para ser atendidos, como en las cajas cobradoras de un supermercado observe la figura 1.3(b), éste es un sistema de colas de líneas múltiples
Otra característica del proceso de colas es el número de espacios de espera en cada fila, es decir, el número de clientes que pueden esperar (o que esperarán ) para ser atendidos en cada línea. En algunos casos, como en un banco, ese número es bastante grande y no significa ningún problema práctico, pues para cuestiones de análisis la cantidad de espacio de espera se considera infinita. En contraste, un sistema telefónico puede mantener un número finito de llamadas (es decir limitado) de llamadas, después del cual las llamadas subsecuentes no tienen acceso al sistema. Las condiciones de espacio de espera infinito y finito requieren análisis matemáticos diferentes
Características claves
Otra característica del proceso de colas es la disciplina de cola, es decir la forma en que los clientes que esperan son seleccionados para ser atendidos. A continuación presentamos algunas de las formas más comunes.
Primero en entrar, primero en salir (PEPS). Los clientes son atendidos en el orden en que van llegando a la fila. Los clientes de un banco y de un supermercado son atendidos de esa manera.
Último en entrar, primero en salir (VEPS). El cliente que ha llegado más recientemente es el primero en ser atendido. Un ejemplo de esta disciplina se da en un proceso de producción en el que los productos llegan a una estación de trabajo y son apilados uno encima del otro. El trabajador elige, para su procesamiento, el producto que esta encima de la pila, que fue el último que llegó para ser procesado o para brindarle un servicio.
Selección de prioridad. A cada cliente que llega se le da una prioridad y se le elige según ésta para brindarle el servicio. Un ejemplo de esta disciplina son los pacientes que llegan a la sala de urgencias de un hospital. Mientras más severo sea el caso, mayor será la prioridad del ¨cliente ¨
El proceso de servicio
El proceso de servicio define cómo son atendidos los clientes. En algunos casos, puede existir más de una estación en el sistema en el cual se proporcione el servicio requerido. Los bancos y los supermercados, de nuevo, son buenos ejemplos de lo anterior. Cada ventanilla y cada registradora son estaciones que proporcionan el mismo servicio. A tales estructuras se les conoce como sistemas de colas de canal múltiple. En dichos sistemas, los servidores pueden ser idénticos, en el sentido en que proporcionan la misma clase de servicio con igual rapidez, o pueden no ser idénticos. Por ejemplo, si todos los cajeros de un banco tienen la misma experiencia, pueden considerarse como idénticos. En este capítulo, se tomarán en cuenta servidores idénticos.
Al contrario de un sistema de canal múltiple, considere un proceso de producción con una estación de trabajo que proporciona el servicio requerido. Todos los productos deben pasar por esa estación de trabajo; en este caso se trata de un sistema de colas de canal sencillo. Es importante hacer notar que incluso en un sistema de canal sencillo pueden existir muchos servidores que, juntos, llevan a cabo la tarea necesaria. Por ejemplo, un negocio de lavado a mano de automóviles, que es una sola estación, puede tener dos empleados que trabajan en un auto de manera simultánea
Otra característica del proceso de servicio es el número de clientes atendidos al mismo tiempo en una estación. En los bancos y en los supermercados (sistema de canal sencillo), solamente un cliente es atendido a la vez. Por el contrario, los pasajeros que esperan en una parada de autobús son atendidos en grupo, según la capacidad del autobús que llegue. En este capítulo se verá el servicio de uno a la vez.
Otra característica más de un proceso de servicio es si se permite o no la prioridad, esto es ¿puede un servidor detener el proceso con el cliente que está atendiendo para dar lugar a un cliente que acaba de llegar? Por ejemplo, en una sala de urgencia, la prioridad se presenta cuando un médico, que está atendiendo un caso que no es crítico es llamado a atender un caso más crítico. Cualquiera que sea el proceso de servicio, es necesario tener una idea de cuánto tiempo se requiere para llevar a cabo el servicio. Esta cantidad es importante debido a que cuanto más dure el servicio, más tendrán que esperar los clientes que llegan. Como en el caso del proceso de llegada, este tiempo pude ser determinístico o probabilístico . Con un tiempo de servicio determinístico, cada cliente requiere precisamente de la misma cantidad conocida de tiempo para ser atendido. Con un tiempo de servicio probabilístico, cada cliente requiere una cantidad distinta e incierta de tiempo de servicio. Los tiempos de servicio probabilísticos se describen matemáticamente mediante una distribución de probabilidad. En la práctica resulta difícil determinar cuál es la distribución real, sin embargo, una distribución que ha resultado confiable en muchas aplicaciones , es la distribución exponencial .En este caso, su función de densidad depende de un parámetro, digamos (la letra griega my) y esta dada por
s(t)=(1/ m )e-m t
en la que
m =número promedio de clientes atendidos por unidad de tiempo, de modo que
1/ m =tiempo promedio invertido en atender a un cliente
En general, el tiempo de servicio puede seguir cualquier distribución, pero, antes de que pueda analizar el sistema, usted necesita identificar dicha distribución.
Clasificaciones de los modelos de cola
Como se menciono al inicio del presente capitulo, para aplicar las técnicas matemáticas apropiadas, usted debe identificar las características de un sistema de colas, basado en la población de clientes y en los procesos de llegada, de colas y de servicio. El método de clasificación presentado aquí pertenece a un sistema de colas en el que el tamaño de la población de clientes es infinita, los clientes que llegan esperan en una sola fila y el espacio de espera en cada línea es efectivamente infinito.
Características claves
En este modelo, los símbolos describen las características del sistema
Cuando de espera y/o el tamaño de la población de clientes es finito, los dos siguientes símbolos adicionales se incluyen para indicar estas limitaciones:
Medidas de rendimiento para evaluar un sistema de colas
El objetivo último de la teoría de colas consiste en responder cuestiones administrativas pertenecientes al diseño y a la operación de un sistema de colas. El gerente de un banco puede querer decidir si programa tres o cuatro cajeros durante la hora de almuerzo. En una estructura de producción, el administrador puede desear evaluar el impacto de la compra de una nueva máquina que pueda procesar los productos con más rapidez.
Cualquier sistema de colas pasa por dos fases básicas. Por ejemplo, considere un día como se muestra en la figura 1.4. Cuando el banco abre en la mañana, no hay nadie en el sistema, de modo que el primer cliente es atendido de forma inmediata. Conforme van llegando más clientes, lentamente se va formando la cola y la cantidad de tiempo que tienen que esperar se empieza a aumentar. A medida que avanza el día, el sistema
Llega a una condición en la que el efecto de la falta inicial de clientes ha sido eliminado y el tiempo de espera de cada cliente ha alcanzado niveles bastante estables.
Como se indica en la figura 1.4, la fase inicial, que conserva los efectos de las condiciones iniciales, se conoce como fase transitoria .Después de que los efectos de las condiciones son eliminados, el sistema entra en un estado estable. A pesar de que las preguntas pertenecientes a ambas faces son importantes, esta sección trata solamente el comportamiento del estado estable.
Algunas medidas de rendimiento comunes
Existen muchas medidas de rendimiento diferentes que se utilizan para evaluar un sistema de colas en estado estable, algunas de las cuales se describen en la presente sección. Para diseñar y poner en operación un sistema de colas, por lo general, se los administradores se preocupan por el nivel de servicio que recibe un cliente, así como el uso apropiado de las instalaciones de servicio de la empresa. Algunas de las medidas que se utilizan para evaluar el rendimiento surgen de hacerse las siguientes preguntas:
Preguntas relacionadas con el tiempo, centradas en el cliente, como:
Preguntas cuantitativas relacionadas al número de cliente, como:
Preguntas probabilísticas que implican tanto a los clientes como a los servidores, por ejemplo:
Preguntas relacionadas con los costos, como:
El cálculo específico de estas medidas de rendimiento depende de la clase de sistema de colas. Algunas de estas medidas están relacionadas entre sí. Conocer el valor de una medida le permita encontrar el valor de una medida relacionada. Tales relaciones generales se describen primeramente en la sección 1.1.5, antes de que se presenten los métodos utilizados para calcular estas medidas de rendimiento para un sistema de colas dado
Relaciones entre medidas de rendimiento
El cálculo de muchas de las medidas de rendimiento depende de los procesos de llegadas y de servicio del sistema de colas en específico. Recuerde de la sección 1.1, que el caso probabilístico, estos procesos son descritos matemáticamente mediante distribuciones de llegada y de servicio. Incluso sin conocer la distribución especifica, las relaciones entre algunas de las medidas de rendimiento pueden obtenerse para ciertos sistemas de colas, únicamente mediante el uso de los siguientes parámetros de los procesos de llegada y de servicio.
l =número promedio de llegadas por unidad de tiempo
m = número promedio de clientes atendidos por unidad de tiempo en una sección
Suponga una población de clientes infinita y una cantidad limitada de espacio de espera en la fila. El tiempo total que un cliente invierte en el sistema es la cantidad de tiempo invertido en la fila más el tiempo durante el cual es atendido:
El tiempo promedio en el sistema y el tiempo promedio de espera están representados por las cantidades W y Wq, respectivamente. El tiempo promedio de servicio puede expresarse en términos de parámetros de &. Por ejemplo, si & es cuatro clientes por hora, entonces , en promedio, cada cliente requiere un cuarto de hora para ser atendido. En general, el tiempo de servicio es 1/&, lo cual nos conduce a la siguiente relación :
Considere ahora la relación entre el número promedio de clientes en el sistema y el tiempo promedio que cada cliente pasa en el sistema. Imagine que un cliente acaba de llegar y se espera que permanezca en el sistema un promedio de media de hora. Durante esta media hora, otros clientes siguen llegando a una tasa ¿¿¿, digamos doce por hora. Cuando el cliente en cuestión abandona el sistema, después de media hora, deja tras de sí un promedio de (1/2)*12 = 6 clientes nuevos. Es decir, en promedio, existen seis clientes en el sistema en cualquier tiempo dado. En términos de & y de las medidas de rendimiento, entonces:
Numero promedio número promedio tiempo promedio de clientes = de llegadas por * en el sistema
en el sistema unidad de tiempo de modo que:
L =l *W
Utilizando una lógica parecida se obtiene la relación entre el número promedio de clientes que esperan en la cola y el tiempo promedio de espera en la fila:
Lq =l * Wq
ANÁLISIS ECONÓMICO DE LOS SISTEMAS DE COLAS
En la sección anterior se vio la ventaja de tener más de un servidor, a saber la reducción del tiempo de espera y el número de clientes que esperan para ser atendidos. Claramente, mientras más servidores se tengan, mejor será el servicio a los clientes. Sin embargo, cada servidor implica costos de operación. ¿De que manera evalúa usted este equilibrio entre el nivel de servicio y el costo?
EJEMPLO 1.1: Problema de colas de American Weavers, Inc.
American Weavers, Inc, tiene una fabrica de manufactura en Georgia. La planta tiene un gran numero de maquinas tejedoras que con frecuencia se atascan. Estas maquinas son reparadas basándose en el procedimiento de la primera en entrar, la primera en ser revisada, por uno de los 7 miembros del personal de reparación. Durante varios recorridos, la gerente de producción ha observado que, en promedio, de 10 a 12 maquinas están fuera de operación en cualquier momento debido a que están atascadas. Ella sabe que contratar personal de reparaciones adicional bajaría el número de máquinas sin funcionar, lo cual traería como consecuencia un aumento en la producción, pero no sabe a cuantas personas más debería contratar. Se desea determinar dicho número.
MODELO Y ANALISIS DEL SISTEMA DE COLA ACTUAL.
El primer paso que se debe dar consiste en analizar las condiciones de operación actuales. Se debe reconocer que las maquinas tejedoras conforman un modelo de colas. Los clientes están constituidos por las maquinas que se atascan de vez en cuando. Existe un gran numero de tales maquinas, de modo que se podría suponer razonablemente, que la población de clientes es infinita. Se tienen 7 servidores independientes e idénticos que reparan las maquinas basándose en una estrategia de primera en entrar, primera en darle servicio. Se puede pensar en estas maquinas formando una sola fila en espera de pasar con el siguiente servidor que este disponible.
Para modelar esta operación, el siguiente paso consiste en reunir y analizar los datos correspondientes a los procesos de llegada y de servicio. Se supone lo siguiente:
Con estas observaciones, el sistema actual puede modelarse como un sistema de colas M / M / 7, con l = 25 maquinas por hora, m = 4 maquinas por hora y una población y un área de espera infinita.
TABLA 1:Medidas de rendimiento obtenidas con Queuing Analysis en el WinQSB .
| 02-08-2000 | Performance Measure | Result |
| 23:12:22 | ||
| 1 | System: M/M/7 | From Formula |
| 2 | Customer arrival rate (lambda) per hour = | 25.0000 |
| 3 | Service rate per server (mu) per hour = | 4.0000 |
| 4 | Overall system effective arrival rate per hour = | 25.0000 |
| 5 | Overall system effective service rate per hour = | 25.0000 |
| 6 | Overall system utilization = | 89.2857 % |
| 7 | Average number of customers in the system (L) = | 12.0973 |
| 8 | Average number of customers in the queue (Lq) = | 5.8473 |
| 9 | Average number of customers in the queue for a busy system (Lb) = | 8.3333 |
| 10 | Average time customer spends in the system (W) = | 0.4839 hours |
| 11 | Average time customer spends in the queue (Wq) = | 0.2339 hours |
| 12 | Average time customer spends in the queue for a busy system (Wb) = | 0.3333 hours |
| 13 | The probability that all servers are idle (Po) = | 0.1017 % |
| 14 | The probability an arriving customer waits (Pw or Pb) = | 70.1674 % |
| 15 | Average number of customers being balked per hour = | 0 |
Como se puede ver, el gerente de producción había estimado con bastante precisión el hecho de que entre 10 y 12 maquinas están atascadas, en promedio; en cualquier momento. De hecho, ese numero en el informe es de 12.09. La línea 10 del reporte indica que las maquinas atascadas están fuera de operación durante un tiempo promedio de 0.4839 horas, aproximadamente 29 minutos.
Es necesario determinar el numero de reparadores adicionales que se necesitarían contratar. Se conocen las medidas de rendimiento de un total de 7 trabajadores.
¿ De que manera cambian las medidas de rendimiento si se aumenta el personal de reparación ? Las medidas de rendimiento asociadas para un numero entre 7 y 11 reparadores se muestran en la TABLA 2
A medida que aumenta el tamaño del personal de 7 a 11, el numero promedio de maquinas fuera de operación disminuye de aproximadamente 12 a 6.333. Similarmente , la cantidad promedio de tiempo que una maquina esta fuera de operación disminuye de 0.4839 horas ( aproximadamente 29 minutos ) a 0.2533 horas ( aproximadamente 15 minutos ).
Ahora se necesita información sobre los costos para determinar cuantos reparadores adicionales, si se requieren, deben contratarse.
TABLA 2: Medidas de rendimiento con diferentes tamaños de personal de reparación.
| Numero | de | reparadores | |||
| 7 | 8 | 9 | 10 | 11 | |
| Utilización (%) | 89.2857 | 78.1250 | 69.4444 | 62.5000 | 56.8182 |
| Numero esperado en la cola | 5.8473 | 1.4936 | 0.5363 | 0.2094 | 0.0830 |
| Numero esperado en el sistema | 12.0973 | 7.7436 | 6.7863 | 6.4594 | 6.3330 |
| Probabilidad de que un cliente tenga que esperar | 0.7017 | 0.4182 | 0.2360 | 0.1257 | 0.0630 |
| Tiempo esperado en cola | 0.2339 | 0.0597 | 0.0215 | 0.0084 | 0.0033 |
| Tiempo esperado en el sistema | 0.4839 | 0.3097 | 0.2715 | 0.2584 | 0.2533 |
ANALISIS DE COSTOS DEL SISTEMA DE COLAS.
Al analizar los méritos de contratar personal de reparación adicional en American Weavers, Inc., se deben identificar dos componentes importantes:
Costo total por = Costo por hora para cada * Numero promedio la espera maquina fuera de operación máquina fuera de operación
Para seguir adelante, se necesita ahora conocer el costo por hora de cada miembro del personal de reparación ( denotado con Cs ) y el costo por hora de una maquina fuera de operación ( denotado Ce ), que es el costo de una hora de producción perdida. Suponga que el departamento de contabilidad le informa que cada reparador le cuesta a la compañía $ 50 por hora, incluyendo impuestos, prestaciones, etc. El costo de una hora de producción perdida deberá incluir costos explícitos, como la cantidad de ganancias no obtenidas, y costos implícitos, como la perdida de voluntad del cliente no se cumple con la fecha limite de entrega.
Sin embargo, suponga que el departamento de contabilidad estima que la compañía pierde $ 100 por cada hora que una maquina este fuera de operación.
Ahora se puede calcular un costo total para cada uno de los tamaños de personal. Para un personal de 7 reparadores, el numero esperado de maquinas en el sistema es 12. 0973.
Costo total = Costo del personal + Costo de la espera
Costo por hora por * Numero de + Costo por hora por * Número persona reparadores cada maquina fuera esperado de de operación máquinas fuera de operación = ( 50 * 7 ) + ( 100 * 12.0973 ) = $ 1559.73 por hora.
Realizando cálculos parecidos para cada uno de los tamaños de personal restantes se tiene como resultado los costos por hora de cada alternativa presentada en la TABLA 3.
De los resultados, se puede ver que la alternativa que tiene el menor costo por hora, $ 1128.63, es tener un total de 9 reparadores. En consecuencia, la recomendación a la gerencia de producción, es contratar a dos reparadores adicionales. Estos dos nuevos empleados tendrán un costo de $ 100 por hora, pero este costo adicional esta mas que justificado por los ahorros que se tendrán con menos maquinas fuera de operación. La recomendación reducirá el costo por hora de $1559.73 a $ 1128.63, un ahorro de aproximadamente $ 430 por hora, mayor que la cantidad que cubre sus honorarios.
TABLA 3: Costo por hora para diferentes tamaños de personal de reparación.
| Tamaño de personal | Numero esperado en el sistema | Costo por hora ($) |
| 7 | 12.0973 | ( 50 * 7 ) + ( 100 * 12.0973 ) = 1559.73 |
| 8 | 7.7436 | ( 50 * 8 ) + ( 100 * 7.7436 ) = 1174.36 |
| 9 | 6.7863 | ( 50 * 9 ) + ( 100 * 6.7863 ) = 1128.63 |
| 10 | 6.4594 | ( 50 * 10 ) + ( 100 * 6.4594 ) = 1145.94 |
| 11 | 6.3330 | ( 50 * 11 ) + ( 100 * 6.3330 ) = 1183.30 |
Características claves
En resumen, para evaluar un sistema de colas en el que usted controla el número servidores o su tasa de servicio, se necesitan las siguientes estimaciones de costos y medidas de rendimiento:
Modelos que caracterizan los diferentes procesos de servicio.
A continuación se analizarán los diferentes modelos que caracterizan los procesos de servicio.
Se estudiarán a continuación modelos para líneas de espera que combinan procesos de llegadas y salidas. Solo analizaremos líneas de espera donde los clientes son atendidos por c servidores que ofrecen servicios iguales desde el punto de vista del tiempo que requieren para prestar el mismo.
El objetivo final de analizar situaciones de espera consiste en generar medidas de desempeño para evaluar sistemas reales. Se hace necesario para analizar los sistemas de espera decidir con anterioridad si nos interesa analizar el sistema en condiciones transitorias o de estado estable. Las líneas de espera que combinan salidas y llegadas se inician en condiciones transitorias y llegan gradualmente al estado estable después de haber transcurrido un tiempo lo suficientemente grande, siempre que los parámetros permitan se alcance el estado estable.
Utilizando el modelo general de Poisson, obtenemos la ecuación general de equilibrio:
l n-1pn-1 + m n + 1pn+1 = (l n+m n)pn , n =1,2,...
Las ecuaciones de equilibrio se resuelven en forma recursiva comenzando con p, y procediendo por inducción para determinar pn. En general, podemos demostrar por inducción que la probabilidad de que el sistema alcance el estado estable es:
pn = [ ( l n-1* l n-2* ¼ l 0) / ( m n* m n-1* ¼ m 1) ] * p0 , n = 1,2
donde el valor de p0 se determina con la siguiente ecuación:
pn = 1
Las medidas de desempeño del estado estable se pueden usar para analizar la operación de las líneas de espera con el fin de hacer recomendaciones sobre el diseño del sistema. Se puede decir que existe una íntima relación entre Ls y Ws de manera que cualquier medida
Se determina automáticamente a partir de la otra. Siendo l e, f la tasa promedio efectiva de llegadas, entonces:
Ls = l e, f * Ws Y Lq = l e, f * Wq
También existe una relación directa entre Ws y Wq pues por definición:
Ws = Wq + 1/m
siendo m la tasa de servicio por servidor activo, mientras que el tiempo estimado de servicio es 1/m ..
- Modelos para líneas de espera especializadas de Poisson.
Cada modelo de los que a continuación se analizaran se describe en términos de la notación extendida por Kendall, como la deducción de pn es completamente independiente de la disciplina de la línea de espera, es apropiado usar el símbolo DG (disciplina general) en la notación de Kendall.
Este es un modelo de servidor único sin límites en la capacidad del sistema o de la fuente de llamadas, con llegadas y salidas de Poisson con tasa medias l y m ..
Definiendo r = l / m obtenemos la siguiente fórmula general para este modelo:
Pn = (1- r ) * r n , n = 0,1,2,¼
que es una distribución geométrica, donde además p0 = 1- r .
El requisito matemático de que r > 1 necesario para garantizar la convergencia de la serie geométrica (1 + r + r 2 +¼ ) , conduce a un elemento intuitivo. O sea r > 1 significa que l < m lo que establece que la tasa de llegadas debe ser estrictamente menor que la tasa de servicio en la instalación, para que el sistema alcance estabilidad.
Para este modelo las medidas de básicas de desempeño se calculan de la siguiente forma :
Ls = E{ n} = (r ) / (1-r ) Ws = Ls / l = 1 / [ m ( 1-r ) ]
Lq = Ls - l / m = r 2 /(1-r ) Wq = Lq / l = r / [ m ( 1-r ) ]
La diferencia de este modelo y el anterior radica en que el número máximo de clientes (para este modelo) permitidos en el sistema es N (longitud máxima de la línea de espera es = N-1). Esto significa cuando haya N clientes en el sistema, se impiden todas las nuevas llegadas o no se les permite unirse a la línea de espera.
En este modelo, haciendo r = l / m obtenemos que:
(1-r ) / ( 1-r N+1) ; r ¹ 1
Po = 1 / N+1 ; r = 1
entonces las formulas para pn pueden resumirse como:
[ (1-r ) / ( 1-r N+1) ] r n ; r ¹ 1
pn = n = 0,1,2, ...
1 / (N+1) ; r = 1
Para este modelo no se hace necesario que r > 1 pues el número de unidades en el sistema esta controlado por la longitud de la línea de espera (= N-1).
Usando el valor anterior de pn, encontramos que el número esperado de unidades en el sistema se calcula como sigue:
{ r [ 1-( N+1)* r N + Nr N+1] } / [ ( 1-r ) ( 1-r N+1) ] ; r ¹ 1
Ls = N / 2 ; r = 1
Las medidas Lq, Ws y Wq se pueden calcular a partir de Ls, una vez que se determina la tasa efectiva de llegadas l ef de la forma siguiente:
l e, f = l ( 1-pn)
usando Ls y l ef obtenemos las formulas para calcular, Lq, Wq y Ws:
Lq = Ls-( l e, f / m ) = Ls - [ ( 1-pN )] / m pN = Probabilidad de que una unidad no sea
capaz de unirse al sistema.
Wq = Lq / l e, f = Ls / [ l ( 1-pN ) ]
Ws = Wq +1/ m = Ls / [ l ( 1-pN ) ]
En este modelo los clientes llegan con una tasa constante l y un máximo de c unidades puede ser atendidos simultáneamente. La tasa de servicio por servidor activo es también constante e igual a m y l e, f = l ..
El efecto último de usar c servidores paralelos es acelerar la tasa de servicio al permitir servicios simultáneos. Si el numero de clientes en el sistema, n, es igual o excede a c, la tasa combinadas de salidas de la instalación es c* m .. Por otra parte, si n es menor que c, la tasa de servicio es igual a n* m . Así, en terminos del modelo generalizado, tenemos:
l n = l n ³ 0
nm , n £ c
m n = cm , n ³ c
Si hacemos r = l ¤ m ; el valor de pn y p0 se calcula de la siguiente forma:
( r n/ n!) p0 0 £ n £ c
Pn = ( r n/ cn-cc!) p0 n > c
p0 = {
( r n/ n!) + r c/ [ c!( 1-r / c) ] } -1
Los valores de las medidas de desempeño se obtienen como sigue :
Lq = [ r c+1/ (c-1)(c-r ) ] p0 = [ cr / ( c-r ) 2] pc
Ls = Lq + r
Wq = Lq/?
Ws = Wq + 1/m
Las operaciones asociadas con este modelo pueden ser tediosas. Morse (1958) da dos aproximaciones útiles para po y Lq. Para r mucho menor que 1,
P0 @ 1-r y Lq @ r c+1/ c2
Y para r / c muy próxima a 1,
P0 @ [ (c-r )(c-1)] ! / cc y Lq @ r / c-r
- Modelo (M/M/c) : (DG/N/¥ ) , c £ N.
Esta situación de espera difiere de la anterior pues se impone un límite N sobre la capacidad del sistema (es decir, tamaño máximo de la línea de espera = N-c). En términos del modelo generalizado, l n y m n para el modelo actual están dadas por :
l , 0 £ n < N
l n = 0, n ³ N
nm , 0 £ n £ c
m n = cm , c£ n £ N
Sustituyendo por l n y m n en la expresión general de pn y observando que r = l / m ; se obtiene :
( r n/ n!) p0, 0 £ n £ c
Pn = [ r n/ ( c!* cn-c) ] p0, c £ n £ N
donde :
{
( r n/ n!) + r c[ 1-( r / c) N-c+1] / c!(1-r / c) } -1, r / c ¹ 1
p0 =
[
( r n/ n!)
+ ( r c/ c!) ( N-c+1)] -1 , r / c = 1
La única diferencia entre pn en este modelo y el anterior ocurre en la expresión de p0 y el factor de uso r / c no necesita ser menor que 1. Las medidas de desempeño se calculan de la siguiente forma :
p0[ r c+1/ ( c-1)!(c-r )2] { 1-( r / c) N * c - ( N-c)(r / c) N-c[ 1- ( r / c) ] } , r / c ¹ 1
Lq = P0[ r c( N-c)(N-c+1)/ 2c! , r / c = 1
Ls = Lq + (c+
) = Lq + l e,f /
m
donde :
= número estimado de servidores inactivos =
(c-n)pn
l e, f
= l ( 1-pn) = m ( c-
)
l e, f = La tasa efectiva de llegada.
En este modelo el numero de servidores es ilimitado porque el cliente mismo es también el servidor. Este es normalmente el caso en los establecimientos de autoservicio. Una vez más en términos del modelo generalizado se tiene :
l n = l , para toda n ³ 0
m n = nm , para toda n ³ 0
La sustitución directa en la expresión de pn produce que :
Pn = ( l n / n!m n ) p0 = (r n/ n!)p0
ya que
pn = 1 se deduce que :
P0 = 1 / 1+ r + ( r 2/ 2 ) + ¼ = 1/ er = e-r
Como resultado se obtiene que :
Pn = ( e-r r n) / n! , n = 0,1,2,¼ , además pn sigue una distribución de Poisson con media
E{ n} = r .
Las medidas de desempeño se pueden obtener mediante la siguientes expresiones :
Ls = E{ n} = r .
Ws =1/m .
Lq = Wq = 0.
Nótese que Wq = 0 porque cada cliente se atiende a si mismo. Esta es la razón por la que Ws es igual al tiempo de servicio medio 1/m .
A modo de obtener una conclusión parcial de lo analizado hasta aquí, se puede decir que los resultados del modelo (M/M/¥ ) : (DG/¥ / ¥ ) se pueden utilizar para determinar aproximadamente los del (M/ M / c) : ( DG/¥ / ¥ ) esto cuando c crece "lo suficiente". La ventaja que esto ofrece es que las operaciones son mucho más sencillas en el modelo (M/M/¥ ) : (DG/¥ / ¥ ) . También vale la pena destacar que los cálculos demuestran que cuando r se hace chica (es decir), l es mucho mayor que m , el modelo (M/M/¥ ) es una aproximación bastante exacta del modelo (M/M/c) aun para c tan chica como 10.
- Modelo de servicio de máquinas (M/M/R) : (DG/K/K), R > k.
En este modelo se supone que se dispone de R técnicos en reparaciones para dar servicio a k máquinas. Como una maquina descompuesta no puede generar nuevas llamadas mientras esta en servicio, el modelo es ejemplo de una fuente de llamadas finita.
Si se define l como la tasa de descompostura por máquina, se tiene :
(k-n)l , 0 £ n £ k
l n = 0, n ³ k
nm , 0 £ n £ R
m n = Rm , R £ n £ k
0, n > k
Sustituyendo por l n y m n en la expresión de pn se obtiene :
k r np0 , 0 £ n £ R n
Pn =k (n!r n / R!Rn-R ) p0, R £ n £ k n
p0 = {
k r n +
k (n!r n / R!Rn-R
) n n
las otras medidas de desempeño están dadas por :
Lq =
(n-R)pn
Ls = Lq + (R-
) = Lq +l ef /m
Donde :
= número esperado de técnico
s = ![]()
l ef
= m (R-
) = l (k-Ls)
La segunda expresión de l ef se obtiene como sigue. Como la tasa de llegadas dadas n máquinas en el sistema es l ( k-n)(donde l es la tasa de descompostura por máquina), en condiciones de estado estable :
l ef = E{ l ( k-n) } = l ( k-Ls)
Los resultados se aplican al caso de un solo técnico haciendo sencillamente R = 1. En este caso se puede probar que :
Lq = k - [ 1+( 1/ r ) ] ( 1 - p0)
Ls = k -( 1 - p0) / r
Los modelos de líneas de espera donde los procesos de llegada y/o salida no siguen las hipótesis de Poisson, conducen a resultados complejos y tal vez poco manejables. En estos modelos el tiempo de servicio se describe por medio de una distribución general de probabilidades con media E{ t} y varianza var{ t} . El análisis de esta situación no proporciona una expresión analítica manejable para las probabilidades pn. En vez de esto, los resultados del modelo sólo proporcionan las medidas básicas de desempeño, Ls, Lq, Ws, y Wq.
Para l igual a la tasa de llegadas a una instalación con un solo servidor y dadas E{ t} y var{ t} como la media y la varianza de la distribución del tiempo de servicio, se puede demostrar usando el análisis de cadena de Markov que :
Ls = l E{ t} +[ l 2( E2{ t} +var{ t} ) ] / 2 (1-l E{ t} ) , (M/G/1) : DG/¥ / ¥ )
donde : l E{ t} < 1 , de esta fórmula se obtienen :
Ws = Ls /l
Lq = Ls-l E{ t}
Wq = Lq /l
Aquí la tasa de servicio está dada por m = 1/ E{ t} y l e, f = l para este modelo. Para el caso cuando el tiempo de servicio es aproximadamente constante, var{ t} = 0 la fórmula de Ls anterior se reduce a :
Ls = r + r 2/ [ 2 (1- r ) ] , (M/D/1) : (DG/¥ / ¥ )
donde r = l / m y m es la tasa constante de servicio.
Las otras medidas de desempeño permanecen iguales.
Si el tiempo de servicio es tipo Erlang con parámetros m y m , con E{ t} = 1/ m y var{ t} =1/(mm 2) , la fórmula de Ls sería :
Ls = r + ( 1+m) / (2m) [ r / ( 1- r ) ] , (M/Em/1) : (DG/¥ / ¥ )
En los modelos de espera con prioridad, se suponen que se forman varias líneas de espera en paralelo, incluyendo los clientes que pertenezcan a cierto orden de prioridad. Si la instalación tiene m, líneas de espera, suponiendo que la línea de espera 1 tiene más alta prioridad de servicio y la línea de espera m incluye a clientes con más baja prioridad. Las tasas de llegadas y servicio pueden variar para las diferentes filas de prioridad.
El servicio de prioridad puede seguir una de estas dos reglas:
A continuación se analizaran dos modelos que se rigen por la regla de no prioridad, que se aplican a servidores únicos y múltiples. En el caso de servidor único supone llegadas de Poisson y distribuciones de servicio arbitrarias. En el caso de servidores múltiples, las llegadas y salidas siguen la distribución de Poisson. El símbolo NPRP se utiliza con la notación de Kendall para representar la disciplina de no prioridad, Mi y Gj representan distribuciones de Poisson y arbitraria.
Sea Fi (t) la FDA de la distribución de tiempo de servicio arbitraria para la i-ésima línea de espera (i = 1,2, m) y sea E{ t} y var{ t} la media y la varianza ,además sea l i la tasa de llegada en la i-ésima fila por el tiempo unitario. Los valores de las medidas de desempeño se obtendrán de las siguientes expresiones :
Wq(k) = [ å l i (E2{ t} +var{ t} )] / 2( 1-Sk-1) ( 1-Sk)
Lq(k) = l kWq(k)
Ws(k) = Wq(k) + Ek{ t}
Ls(k) = Lq(k) + r k
donde :
r k = l kE{ t}
Sk = å r i< 1 , k = 1,2,3 , m
S0 = 0
y aquí se definen Lq(k), Wq(k), Ls(k) y Ws(k) como las medidas de desempeño para la k-ésima línea de espera. Además el tiempo de espera estimado en la línea de espera para cualquier cliente sin importar su prioridad está dado por :
Wq =
(
l k/ l ) Wq(k)
donde :
l =
l i y l k/ l es
el peso relativo de Wq(k). Un resultado similar se
aplica a Ws.
- Modelo de no prioridad con varios servidores (M/M/c) : (NPRP/¥ / ¥ ) .
Este modelo supone que todos los clientes tienen la misma distribución del tiempo de servicio, los c canales tienen una distribución de servicio exponencial idéntica con tasa de servicio m . Las llegadas de K-ésima línea de espera con prioridad ocurren según una distribución de Poisson con una tasa de llegada l K, k = 1,2, ,m. Se puede demostrar para la k-ésima línea de espera que:
Wq(k) = E{ x 0} / ( 1-Sk-1) ( 1-Sk) , k = 1,2, ,m
donde :
S0 = 0
Sk =
l i / cm < 1
E{ x 0} = 1/ cm [ r -c( c-r ) ( c-1) !( å r n/ n!+ 1 ) , r = l / m
y la espera numérica estimada en la línea de espera para todo el sistema viene dada por :
Lq = l Wq.
A continuación se analizaran líneas de espera de Poisson con estaciones de servicio dispuestas en serie, de manera que todo cliente debe pasar por las dos estaciones antes de completar el servicio.
Para este modelo que describe esta situación los tiempos de servicio en cada estación están exponencialmente distribuidos con la misma tasa de servicio m . Las llegadas ocurren según una distribución de Poisson con tasa l y no se permiten colas en ninguna de las dos estaciones.
Para plantear este modelo se requiere de identificar primeramente los estados del sistema. Los símbolos 0, 1, b, representan los estados libres, ocupado y bloqueado respectivamente. Sean i y j los estados en la estación 1 y 2. Entonces los estados del sistema se pueden representar :
{ ( i, j) } = { ( 0, 0) ; ( 1, 0) ; ( 0, 1) ; ( 1, 1) ; ( b, 1) }
Definiendo pij(t) como la probabilidad de que el sistema se halle en el estado (i, j) en el tiempo t. Las probabilidades de transición entre los tiempos t y t+h se resumen en la siguiente tabla:
| Estados t | (0, 0) | (0, 1) | (1, 0) | (1, 1) | (b, ¡) |
| (0, 0) | 1- l * h | - | l * h | - | - |
| (0, 1) | m * h( 1-l * h) | 1-m * h-l * h | - | l * h( 1-m h) | - |
| (1, 0) | - | m * h( 1-l * h) | 1 - m h | - | - |
| (1, 1) | - | m * h( 1-l * h) | m * h | ( 1-m h) ( 1-m h) | m * h |
| (b, 1) | - | m * h( 1-l |