1 / 61

Faculdade de Engenharia - Campus de Guaratinguetá

Faculdade de Engenharia - Campus de Guaratinguetá. Pesquisa Operacional Livro: Introdução à Pesquisa Operacional Capítulo 6 – Teoria de Filas Fernando Marins – fmarins@feg.unesp.br Departamento de Produção. Sumário. Introdução Estrutura de um Sistema de Filas Características Básicas

iniko
Download Presentation

Faculdade de Engenharia - Campus de Guaratinguetá

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Faculdade de Engenharia - Campus de Guaratinguetá Pesquisa Operacional Livro: Introdução à Pesquisa Operacional Capítulo 6 – Teoria de Filas Fernando Marins – fmarins@feg.unesp.br Departamento de Produção

  2. Sumário Introdução Estrutura de um Sistema de Filas Características Básicas Notação de Kendall Medidas de Desempenho de Sistemas de Filas Tipos de Sistemas de Filas Estatística e Sistemas de Filas Problemas de Decisão em Sistemas de Filas Processos de Nascimento e Morte (P-N-M) Modelos Markovianos Exercícios

  3. Introdução Objetivo: Desenvolvimento de modelos matemáticos para permitir prever o comportamento de Sistemas de Prestação de Serviços (Sistemas de Filas). Motivação: Possibilidade de testar possíveis modificações em Sistemas de Filas que contribuam para melhorar seu rendimento, e obter subsídios para escolha da melhor alternativa de ação.

  4. Exemplos de Sistemas de Prestação de Serviços (Filas) • Congestionamentos de sistemas telefônicos (Erlang); • Escolha do tipo de sinalização para intersecções de vias urbanas (pare, semáforo); • Otimização de ciclo de um semáforo; • Dimensionamento de equipes de manutenção; • Análise de congestionamentos em aeroportos.

  5. Estrutura de um Sistema de Filas Capacidade do Sistema Processo de Chegadas Processo de Atendimento Tamanho da População Disciplina de Atendimento

  6. Características Básicas Processo de Chegadas: Quantos clientes solicitaram serviço em [0, t]? • Caracterizado pela distribuição de probabilidades dos intervalos entre chegadas consecutivas (ou pela distribuição de probabilidades do número de chegadas em intervalos disjuntos). • Chegadas podem ser individuais ou em grupos. • Caso importante: Número de Chegadas segue um Processo de Poisson com média de chegadas .

  7. Características Básicas Processo de Atendimento: Quanto tempo demora o atendimento? • Caracterizado pela distribuição de probabilidades da duração do atendimento (ou distribuição de probabilidades do número de atendimentos em intervalos disjuntos). . • A estação de serviços pode ser formada por um ou mais servidores. • Atendimento pode ser individual ou em grupos. • Caso importante: duração do atendimento é uma variável aleatória com distribuição de probabilidades dada por uma exponencial negativa com média 1/.

  8. Características Básicas Tamanho da População: Finita ou Infinita Capacidade do sistema: Quantos clientes podem estar no sistema (fila + posições de atendimento) ao mesmo tempo. Disciplina de Atendimento: - Forma como os clientes são selecionados da fila para serem atendidos. Tipos: FIFO, LIFO, FEFO, com prioridades, aleatória. • Modelos sofisticados: • Troca de filas • Desistências • Redes de filas • Atendimento/chegada dependente do tamanho da fila

  9. Notação de Kendall   Fila Campo I/Campo II/Campo III/Campo IV/Campo V/Campo VI com Campo I - informação sobre a distribuição de probabilidades dos intervalos entre chegadas (D, M, EK ,G); Campo II - informação sobre a distribuição de probabilidades dos tempos de atendimentos (D, M, EK ,G); Campo III - informação sobre o número de atendentes (em paralelo) para a fila; Campo IV - informação sobre a capacidade do sistema; Campo V - informação sobre o tamanho da população; Campo VI - informação sobre a disciplina de atendimento. Exemplo: Fila M/M/3/20/300/LIFO com = 3/h e  = 5/h. Obs: Quando capacidade do sistema e o tamanho da população puderem ser admitidas como sendo , e a disciplina for a FIFO os Campos IV, V , VI podem ser omitidos.

  10. Medidas de Desempenho (KPIs) de Sistemas de Filas (Queues) • Tempo que o cliente fica na fila - Wq • Tempo que o cliente fica no sistema – Ws (ou W) • Número de clientes na fila - Lq • Número de clientes no sistema – Ls (ou L) • Ociosidade dos servidores - P0 • Observe-se que todas estas medidas são variáveis aleatórias. Convencionou-se analisar seus valores médios medidos numa situação de estado de equilíbrio de funcionamento do sistema.

  11. Tipos de Sistemas de Filas Markovianos - para caracterizar o comportamento futuro do sistema basta conhecer o seu estado atual (por exemplo: o número de clientes no sistema no momento). Não-Markovianos - para caracterizar o comportamento futuro do sistema é necessário conhecer o seu estado atual e se ter informações sobre o passado (por exemplo: o tempo que o cliente que está sendo atendido já demandou de serviço).

  12. Estatística e Sistemas de Filas Inferência Estatística: Qual modelo se adapta a uma dada situação? Necessário aplicar técnicas estatísticas de estimação de parâmetros e testes de aderência para a escolha da distribuição de probabilidades adequada para cada caso. Importante: as medidas de eficiência do sistema de filas são dependentes das taxas de chegadas () e de atendimentos ().

  13. Problemas de Decisão em Sistemas de Filas Projeto - Sistema deve ter algumas características desejadas  escolher o número de servidores, qual o tamanho máximo admissível para a fila, qual a disciplina de atendimento adequada... Controle - Como e quando alterar as características básicas do sistema de forma a otimizar algum critério econômico de interesse.

  14. Problemas de Decisão em Sistemas de Filas Variáveis comumente observadas - tamanho da fila ou a quantidade de trabalho acumulado por executar pelos servidores. Custos comumente envolvidos - custo de espera (cliente), custo do atendimento (servidor), custo de perda de clientes em potencial, custo de ociosidade do servidor. Decisões possíveis - fechar o sistema para novos clientes e/ou acelerar o atendimento. Importante: na impossibilidade de se desenvolver modelos matemáticos para situações complexas recomenda-se o uso da técnica de Simulação de Sistemas.

  15. Processos de Nascimento e Morte (P-N-M): Definições Estado do sistema = número de clientes no sistema = n Nascimentorepresenta entrada de cliente no sistema. Morte representa saída de cliente do sistema. n = taxa média de chegadas quando há n clientes no sistema. n = taxa média de atendimentos quando há n clientes no sistema. P-N-M: chegadas e atendimentos aleatórios, e n e n dependem apenas do estado do sistema.

  16. Processos de Nascimento e Morte (P-N-M): • Postulados (Funçãoo(t) é um nulo quandot tende a zero) • Admita o sistema no estado n no instante t • Nascimento - Probabilidade de ocorrer exatamente 1 nascimento no intervalo de tempo t é dada por n T + o(t). • (II)Morte - Probabilidade de ocorrer exatamente 1 morte no intervalo de tempo t é dada por nT + o(t). • (III)Salto múltiplo - Probabilidade de ocorrer número de nascimentos e mortes superior a 1 é dado por o(t).

  17. Processos de Nascimento e Morte (P-N-M): Consequência: o estado n do sistema em t + t pode ter sido alcançado a partir de 4 situações mutuamente exclusivas do sistema no instante t: n, n - 1, n + 1, outra. Seja pn(t) = probabilidade do sistema estar no estado n no instante t = probabilidade de existir n clientes no sistema no instante t.

  18. Processos de Nascimento e Morte (P-N-M): Rearranjando a expressão pn(t + t), dividindo por t e passando ao limite quando t  0, tem-se: = n-1pn-1(t) + n+1pn+1(t) - (n + n)pn(t). Analisando o sistema em regime estacionário, tem-se = 0 Que resulta em:

  19. Processos de Nascimento e Morte (P-N-M): Da expressão (A) tem-se: Da expressão (B) tem-se: Aplicando conjuntamente ( C) e (D) tem-se:

  20. Processos de Nascimento e Morte (P-N-M): Como sabe-se que, considerando todos os estados possíveis do sistema, tem-se: Substituindo a expressão de Pn+1 vem: P0 = proporção do tempo que o sistema fica vazio (servidores ociosos)

  21. Modelos Markovianos 1. Fila M/M/1 2. Fila M/M/S  3. Fila M/M/1/N* 4. Fila M/M/S/N* 5. Fila M/M/1 com população finita p 6. Fila M/M/S com população finita p Obs: No livro texto há 9 modelos descritos.

  22. Modelo M/M/1 Como não há limitação para a formação de fila, e a população é infinita, tem-se: n=  e n=  para todo n = 0, 1, 2, 3... Assim:

  23. Modelo M/M/1 Defina  = / como o fator de utilização ou congestionamento do Sistema de Filas. Tem-se então: Pn= n(1- ) para n = 0, 1, 2, 3,... Cálculo do número médio de clientes no sistema: Cálculo do número médio de clientes na fila:

  24. Modelo M/M/1 Utilizando-se a Fórmula de Little (1961) – LS, Q = WS, Q, tem-se: Tempo médio de um cliente no sistema – Tempo médio de um cliente na fila –

  25. Modelo M/M/S Neste caso o modelo de filas tem s servidores em paralelo para atendimento dos clientes: Onde  = taxa média de atendimento/servidor n =  para n = 0, 1, 2, 3...

  26. Modelo M/M/S Para este sistema tem-se: com  = /s < 1 Obs: ver Diagrama 1 com os valores de P0 para S= 1, 2, 3, 4, 5, 7,10, 15, 20, 25. PN = -1

  27. Diagrama 1 - Valores de P0 para a Fila M/M/S 7

  28. Modelo M/M/S 7 Fórmula C de Erlang WQ= LQ/ WS= WQ + 1/ LS = LQ + / Obs: ver Diagrama 2 com valores de LS para S = 1, 2, 3, 4, 5, 7, 10, 15, 20, 25.

  29. Diagrama 2 - Valores de LS para a Fila M/M/S

  30. FILA M/M/1/N* Processo de Chegadas: Poisson, taxa  Tempo de Serviço: Exponencial, taxa  Número de Atendentes: 1 Capacidade do Sistema: N*

  31. FILA M/M/1/N* • A restrição de capacidade pode surgir ou por uma restrição física do sistema ou por iniciativa do cliente. • Como existe essa restrição de capacidade, neste caso tem-se: • e para

  32. FILA M/M/1/N* c) Não será preciso a restrição  = /s < 1 Seguindo a mesma seqüência da fila M/M/1, tem-se:

  33. FILA M/M/1/N* para

  34. Para a Fila M/M/1/N*, o tempo médio que um usuário passa no sistema pode ser de dois tipos: tempo médio, levando em conta todos os usuários, atribuindo-se tempo zero aos usuários que foram recusados; ele pode ser calculado diretamente por: tempo médio, considerando-se apenas os usuários que se juntam ao sistema; deve-se tomar algum cuidado no seu cálculo, pois, neste caso, o processo de chegadas não possui mais taxa , pois toda vez que um consumidor encontra N* usuários no sistema, ele é perdido. Pode-se mostrar que a taxa de chegadas é agora igual a e portanto:

  35. FILA M/M/S/N* Observação: para 35

  36. FILA M/M/S/N* 36

  37. FILA M/M/S/N* 37

  38. FILA M/M/1 , com população finita p Em alguns casos o número de clientes em potencial para a estação de serviços é pequeno. Se este valor for tão pequeno que a chegada de um cliente para ser atendido ou um atendimento afeta a probabilidade de futuras chegadas, não será mais válido o pressuposto de uma população infinita. Assim, tem-se: onde p = número de usuários na população

  39. FILA M/M/1 , com população finita p Pode-se obter as seguintes expressões: Obs: ver Diagrama 3 com os valores de P0 para p= 1, 2, 3, 4, 5, 6, 7, 8, 9,10

  40. M/M/1, com população finita p p=1 p=2 p=3 p=4 p=5 p=6 p=7 p=8 p=9 p=10

  41. FILA M/M/S , com população finita p Aqui tem-se: para

  42. FILA M/M/S , com população finita p Pode-se obter:

  43. FILA M/M/S , com população finita p

  44. Exercícios Uma copiadora para uso de escritório é usada e operada por pessoas deste mesmo escritório que precisam fazer cópias, principalmente secretárias. Como o trabalho a ser copiado varia de tamanho (no de páginas do original) e quanto ao no de cópias, a taxa de atendimento é distribuída aleatoriamente, mas se aproxima de um processo de Poisson tendo uma taxa de atendimento médio de 10 trabalhos por hora. Geralmente as necessidades de uso são aleatórias durante as 8 horas de trabalho diário, mas chegam a uma taxa de 5 por hora.

  45. Exercícios Várias pessoas notaram que surge uma linha de espera às vezes, e elas têm questionado a política de se manter apenas uma unidade. Se o tempo de uma secretária custa $3,50 por hora, faça uma análise para determinar: (a)  A utilização da copiadora; (b)  O percentual de vezes que uma secretária tem de esperar para usar a copiadora; (c)  O tempo médio que uma secretária permanece no sistema de filas da copiadora; (d)  O custo total médio de espera e operação da copiadora pelas secretárias num dia. L Po

  46. Exercícios • A situação apresentada no exercício (1) foi questionada e resolveu-se considerar a possibilidade de se instalar duas máquinas ou de se alugar uma máquina maior, sendo que os dados adicionais são fornecidos a seguir. • Recalcule todos os itens de (a) a (d). Qual é a melhor opção entre as três analisadas: 1 máquina pequena, 1 máquina grande ou 2 máquinas pequenas?

  47. Exercícios • Um empresário tem plano para abrir um serviço automático de lavagem de carros numa determinada região de uma cidade, para tanto realizou uma pesquisa que proporcionou os seguintes dados: • Número de Clientes em potencial deverá seguir uma distribuição de Poisson, com uma chegada cada 5 minutos, desde que haja lugar na área de estacionamento do sistema de lavagem de carros. • Tempo para lavar um carro  deverá seguir uma distribuição Exponencial Negativa, com média de 4 minutos.

  48. Exercícios Para ajudar na decisão de onde abrir o negócio, o empresário identificou 3 locais disponíveis, L1, L2 e L3. Estes locais têm capacidades diferentes, com respeito à área para estacionamento, conforme abaixo: (a)  L1 - não tem espaço para estacionamento, só cabe a máquina para lavar os carros; (b)  L2 – há espaço para 2 carros estacionarem, além daquele que está sendo atendido na máquina; (c)  L3 - há espaço para 4 carros estacionarem, além daquele que está sendo atendido na máquina.

  49. Exercícios Evidentemente, os valores dos aluguéis de cada local são diferentes, sendo o do L1 o mais barato e o de L3 o mais caro. Para ter mais uma informação para a sua tomada de decisão, além dos valores dos aluguéis, o empresário deseja comparar porcentagem de fregueses perdidos por não haver espaço no estacionamento do local escolhido (não incluindo o carro sendo lavado).

  50. Exercícios • Caminhões chegam num armazém para descarga segundo uma distribuição de Poisson no ritmo de 3 caminhões por hora. A distribuição do tempo de atendimento é aproximadamente uma exponencial com média de 15 minutos. Calcular: • (a)  o número de caminhões na fila • (b)  o número de caminhões no sistema • (c)  o tempo médio de espera na fila • (d)  o tempo médio de espera no sistema • (e)  a probabilidade de 6 caminhões estarem no sistema • (f)  o fator de utilização. L Po

More Related