ESTUDO PARA A COMPARAÇÃO
ENTRE OS MODELOS DE FILAS M/M/C
E M/EK/C APLICADA EM UMA
PANIFICADORA
IGOR MICHEL SANTOS LEITE (UNAMA)
Jorge Luiz Okabe Auad (UNAMA)
Sadraque Araujo da Silva Junior (UNAMA)
ROMULO JOSE PINHEIRO DA SILVA (UEPA)
Yvelyne Bianca Iunes Santos (UNAMA)
Este artigo buscou estudar a comparação entre dois modelos de filas, o
M/M/c, com chegada por modelo de Poisson e atendimento por
Exponencial Negativa,e o modelo M/Ek/c, com chegado por Poisson e
atendimento por Erlang de grau k, em uma filla formada em uma
panificado de pequeno porte, visando encontrar aquela cujo os
resultadosmais se adequaram a realidade. Porém conclui-se que os
dois modelos são muitos proximos quando o parâmetro k = 1, não
obtendo nenhum resultado com diferença sensível, para o parâmetro
utilizado pode-se dizer que M/M/c e M/E1/c são idênticos. Contudo em
uma situação simulada, fazendo com que os picos no tempo de
atendimento fossem mais centralizados, o modelo M/Ek/c obteve um
melhor desempenho (Erlang a grau 3) quando comparado ao M/M/c.
Palavras-chaves: Filas, Modelo deErlang, Exponencial Negativa,
Distribuição
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
2
1. Introdução
As filas estão distribuídas em todos os seguimentos da sociedade, seja em um supermercado,
ou em uma ligação telefônica, ou em ambientes produtivos. Em certas situações elas podem
ser abstratas, mas não abdicando de possuir uma grandiosa importância.
A formação de uma fila é aleatória e se não estiver bem dimensionada pode causar grandes
transtornos, tanto para os clientes como para o estabelecimento, tais como:
congestionamentos, desperdício financeiros, dentre outros.
Por uma fila ser concebida de maneira aleatória seu estudo é feito através de processos
matemáticos estatísticos (processo estocástico) buscando um equilíbrio entre o custo do
atendimento e a satisfação daqueles que fazem parte desse sistema.
Dentro do contexto apresentado, o presente trabalho tem como objetivo comparar modelos de
atendimento (Exponencial Negativa e Erlang), geradas no acesso ao balcão de atendimento de
uma padaria, identificando o modelo que mais se adéqua a realidade.
2. Conceito das filas
Segundo Tôrres (2010), a fila é formada quando um determinado posto de serviço não pode
atender prontamente um grupo de clientes, ocasionando então a espera para dispor dos
serviços oferecidos pelo posto, tendo sua formação caracterizada basicamente pelos seguintes
elementos:
O regime de chegada: onde cada unidade (cliente) ocupa um lugar de acordo com o
momento.
O regime de serviço: que depende da disponibilidade de atendimento variando de
acordo com o seu intervalo de tempo, pois em alguns lugares o serviço é sempre
disponível.
Capacidade do sistema: de acordo com o número de clientes que podem ser atendidos
simultaneamente. Observando-se a duração do tempo de serviço.
Disciplina de Atendimento: representa a ordem na qual os clientes são atendidos,
sendo normalmente empregada a disciplina FIFO (First in First out),
3. Descrição do Sistema
O local onde foram feitos o estudo e a coleta de dados trata-se de uma panificadora que dispõe
de dois atendentes (caixa).
Tanto o atendimento do pedido quanto a função de caixa é feito por ambos os atendentes,
existindo apenas uma fila para os dois funcionários, onde os clientes aguardam a sua vez de
acordo com a ordem de sua chegada.
Sendo que o sistema operacional da panificadora não é informatizado. O pagamento do
produto é feito somente em dinheiro.
4. Modelagem Analítica do Sistema
A partir da coleta feita e da observação embasadas nos dados referentes às chegadas dos
clientes, ao número de caixas dispostos no local e também ao tempo usado para cada
atendimento, foi definida a modelagem analítica do sistema.
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
3
Segundo Shamblin (1979), grifo dos autores, no modelo M/M/2 tanto a chegada como o
atendimento podem ser marcovianos, seguindo a Distribuição de Poisson (recomendado
para a chegada) ou a Distribuição Exponencial Negativa (recomendado para o processo de
atendimento), entretanto, não foram consideradas apenas as equações referentes às sugestões
dos artífices citados, visto que neste trabalho também será avaliada a utilização da
distribuição de Erlang no processo de atendimento do público.
Para Bolch (1998), “a Notação de Kendall é usada para descrever os sistemas de filas
elementares.”. O sistema observado no estabelecimento é testado no modelo M/M/2 e no
modelo M/Ek/2, onde existe uma única fila para dois atendentes.
Para ambos os processos (chegada e atendimento) se faz necessário o cálculo de:
, seja a taxa média de chegada ao sistema, calculada por:
Onde:
D – Quantidade de Pessoas que chegaram à fila em um intervalo de tempo;
– Ritmo de chegada;
, seja a taxa média de atendimento, calculada por:
Onde:
e – São os tempos de cada intervalo utilizado para o atendimento;
– Frequência Observada para o respectivo intervalo;
Para Hillier (1988) e Bolch (1998), é necessário que o sistema fique em “equilíbrio”, para que
a teoria de filas seja aplicada, considerando sempre ou , fazendo com que o
sistema sempre consiga atender a todos os clientes, pois caso contrário haverá teoricamente a
tendência a formação de uma fila infinita.
A taxa de utilização do sistema pode ser definida como:
4.1 Distribuição de Poisson
A Distribuição de Poisson consiste em uma curva matemática usada na estatística e na
simulação de resultados para representar a probabilidade de que determinado evento ocorraem
um intervalo especificado, sendo antes a probabilidade média conhecida. Na aplicação em
filas é utilizada pelo fato da chegada de pessoas acontecer independe e aleatoriamente no
tempo.
Segundo Prado (2009), a equação de Poisson pode ser escrita na seguinte fórmula:
4.2 Distribuição Exponecial Negativa
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
4
Segundo Moreira (2010) e Prado (2009), a distribuição exponencial negativa é a
correspondente da distribuição de Poisson, porque também se baseia em função de
probabilidade em torno de uma média conhecida. Pode-se dizer que um fenômeno segue as
duas distribuições ao mesmo tempo, Poisson e Exponencial Negativa.
Pelo fato da distribuição estar feita em um intervalo de tempo o resultado calculado através da
integração da equação exponencial negativa, tendo como resultado:
4.3 Distribuição de Erlang
Para Hillier (1998), “a distribuição de tempo de serviço exponencial supõe um grau muito
elevado de variabilidade...”, enquanto a função de Erlang gera uma menor incerteza pelo fato
de encontrar-se em um campo intermediário de variação, onde está localizada a maioria de
serviços reais, ou seja, entre a variação zero e a variação da distribuição exponencial
encontrasse a Distribuição de Erlang.
Ferreira (1998) indica a fórmula da distribuição de Erlang:
Contudo, como a mesma está sendo usada em intervalo de tempo, a equação deve ser utilizada
como a média do intervalo. A utilização da equação desta forma (tempo como média do
intervalo) foi resultado da análise dos autores, sendo esta a que melhor se adequou.
4.4 Teste do Chi-Quadrado
Chi-Quadrado, simbolizado por , é um teste de hipóteses que tem como objetivo, encontrar
um valor da dispersão para duas variáveis nominais, e avaliar a associação existente entre
variáveis quantitativas. O princípio básico deste método é comparar as possíveis divergências
entre as frequências observadas e esperadas para certo evento. Pode-se dizer que dois grupos
se comportam de forma semelhante se as divergências entre as frequências observadas e as
esperadas em cada categoria forem muito próximas à zero. Utiliza-se o teste para verificar se
a frequência com que um determinado acontecimento observado em uma amostra varia
significativamente ou não da frequência com que ele é esperado.
Abaixo segue fórmula do teste Chi-quadrado:
Onde:
o = frequência observada para cada classe
e = frequência esperada para aquela classe
Com o resultado de , calcula-se o grau de liberdade e de posse destes dois valores,
determina-se o máximo erro admissível, que neste trabalho será considerado 5%. Utilizando
os 3 valores ( , G.L. e 5%) encontra-se o valor tabelado, que representa o valor máximo que
pode assumir para o respectivo erro. Caso o resultado seja abaixo do da tabela, a
distribuição é aderente a realidade com erro inferior a 5%.
5. Coleta de dados
5.1 Chegada de clientes
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
5
Os dados de chegada foram coletados entre 6:20h e 7:30h da manhã, desconsiderando 5
observações, de maneira que foram observadas quantos clientes chegava a fila em intervalos
fixos e consecutivos de 1 min, obtendo-se assim a frequência de chegada na fila para cada
minuto, os resultados destas anotações podem ser visualizadas e analisadas na tabela abaixo:
Tabela 1: Frequências Observadas para a chegada de clientes na fila
A partir da Tabela 1 é feito o cálculo da Frequência Calculada, que são os resultados obtidos a
partir da distribuição de Poisson para = 0,92307 .
Tabela 2: Frequências calculadas pela distribuição de Poisson (chegada)
Observado o valor do erro pode-se concluir que o resultado é satisfatório, pois o resultado foi
abaixo do mínimo recomendado para distribuição com grau de liberdade 4 a 5% de
significância, que seria um erro de 9,488. Ou seja, o teste revela que a distribuição de Poisson
para = 0,92307 clientes/min é válido.
Partindo dessa validação o gráfico foi traçado e como resultado obteve-se:
Figura 1: Sobreposição dos dados reais de chegada e a distibuição de Poisson
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
6
O Figura 1 mostra que a distribuição de Poisson se adequa com satisfação a realidade
dapanificadora.
5.2 Atendimento de clientes
Para o atendimento os seus respectivos dados foram obtidos através da cronometragem dos
tempos que cada funcionário gastava para atender um pedido. Como os tempos foram muito
variados, os mesmo foram distribuidos em invervalos de tempo igualmente espaçados.
Para que fossem testadas as distribuições para atendimento, o mesmo foi considerado em
intervalos de tempo de 0min a 7min, dividos em 7classes, onde foram obtidas as frequências
para todos os intervalos, conforme apresenta-se na Tabela 3.
Considerando que OP I (Operador I) e OP II (Operador II).
Tabela 3: Frequência observada e relativa para o atendimento
Com os resultados de e foi feita a comparação entre a taxa de chegada média e a taxa
de atendimento média, usando para isso o resultado multiplicado por 2 (número de
atendentes), bem como o cálculo da taxa de utilização:
Observando que a taxa de utilização do sistema é de 65,3%, ou seja, o sistema encontra-se
estável.
5.2.1 Aplicação da distribuição exponencial negativa
Apartir da tabela 3 foi aplicada à equação da distribuição exponencial negativa para ambos os
funcionarios utilizando suas respectivas taxas de atendimento, gerando como resultado atabela
a seguir.
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
7
Tabela 4: Frequências calculadas pela distribuição Exponencial Negativa (atendimento)
Analisando o resultado do teste do Chi-quadrado e comparando, com os resultados para grau
de liberdade 5 a 5% de significância, obtêm-se o erro máximo admissível de 11,07. Portanto
os erros foram inferiores aos recomendados confirmando a aderência da função à realidade,
para ambos os operadores.
Graficamente, tem-se:
Operador I
Operador II
Figura 2: Sobreposição dos dados reais deatendimento e a distibuição Exponencial Negativa
5.2.2 Aplicação da distribuição de Erlang1
A partir da Tabela 3 foi aplicada à distribuição de Erlang1. Como o modelo proposto neste
tópico adimite vários valores para o fator “k” que modela o formato da curva criada foi
utilizado o valor de k = 1, pois para valores superiores a este, os resultados estavam
desproporcionais a realidade, para ambos os funcionários utilizando suas respectivas taxas de
atendimento, gerando como resultado a seguinte tabela.
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
8
Tabela 5: Sobreposição dos dados reais de atendimento e a distibuição de Erlang1
Analisando o resultado do teste do Chi-quadrado e comparando, com os resultados para grau
de liberdade 3 a 5% de significância, obtêm-se o erro máximo admissível de 11,07. Portanto
como os erros foram inferiores aos recomendados a função adere com satisfação à realidade,
para ambos os operadores, resultados estes iguais ao anterior.
Graficamente, tem-se:
Operador I
Operador II
Figura 3: Sobreposição dos dados reais de atendimento e a distibuição Erlang1
5.3 Simulação com picos de atendimento diferentes
Como parte deste estudo, foi feita a comparação das duas distribuições citadas em situações
onde os picos de atendimento não estão nos intervalos iniciais, mas sim, em faixas
intermediarias de intervalo.
Na tabela abaixo pode ser observada as alterações para a simulação de novos picos de
atendimento, onde o momento de pico de 19 anotações entre 0 e 1 min foi alterado para o
intervalo de 1 a 2 min, mudando os resultados obtidos, para:
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
9
Tabela 6: Frequências
Observadas para a simulação
Desenvolvendo a mesma ideia da tabela 4, distribuição exponencial negativa, porém com os
valores da simulação, tem-se:
Tabela 7: Frequências calculadas pela
Distribuição Exponencial Negativa para a
simulação
Observando o resultado do
teste de aderência conclui-se
que a mesma distribuição
não é mais viável para representa a realidade idealizada, já que o erro é consideravelmente
superior ao valor do erro limite para a significância de 5%. No gráfico é notável a diferença
existente.
Operador I
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
10
Operador II
Figura 4: Sobreposição dos dados reais de atendimento e a Distr. Expo. Negativa para os dados fictícios
Fazendo o cálculo para Erlang k, obteve-se varios tipos de curvas, pelo fato de que cada valor
de k possuir uma curva diferente, e aquela que melhor se adequou aos dados foi Erlang 3, na
qual abaixo encontrasse os resultados:
Tabela 8: Frequências
calculadas pela Distribuição Erlang 3
para a simulação
Comparando os resultados do Chi-
Quadrado da tabela 7 com a tabela 8, é
possível observar a diferença nos
erros, sendo que na distribuição de
Erlang 3 o erro do operrador II é
causado principalmente pelo último intervalo e se caso o mesmo fosse reduzido a 0 o erro
total seria de apenas de 3,3755 e para o operrador I o maior erro gerado foi no penúltimo
intervalo e caso fosse reduzido a zero, semelhante a comparação anterior, o erro total seria de
5,6159. Abaixo segue os respectivos gráficos da distribuição de Erlang 3:
Operador I
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
11
Operador II
Figura 5: Sobreposição dos dados reais de atendimento e a Distr.Erlang 3 para o dado fictícios
6. Análise
Observando os gráficos é possível concluir que a distribuição de Erlang com grau 1 é muito
semelhante à distribuição Exponencial negativa, demonstrando a proximidade dos dois
modelos, para situações com picos de distribuição se concentrando nos primeiros intervalos
de tempo.
Contudo, analisando os gráficos do modelo simulado, conclui-se que a curva de Erlang
melhor se adequou a situação fictícia, mesmo tendo se assemelhado a exponencial negativa
quando . Isso ocorre devido a possibilidade de mutação em sua função pelo fato da
mesma ter o fato de mutação em sua função geratriz, viabilizando a mobilidade (no sentido de
mudança de aspecto) da curva. Sempre tornando-a ajustável a uma situação, dependendo do
grau adotado.
7. Conclusão
Este trabalho foi de grande importância para o aprofundamento sobre as possibilidades que a
Teoria de filas nos proporciona. A modelagem matemática usada neste conteúdo determina a
confiabilidade do estudo, portanto se faz necessário uma boa modelagem para que os
resultados sejam o mais próximo possível da realidade.
A comparação final entre o modelo M/M/c, com chegada pelo modelo de Poisson e
atendimento por Exponecial Negativa, e o modelo M/E1/c, com chegado por Poisson e
atendimento por Erlang de grau 1 (valores superiores para o parâmetro k fazia com que a
curva ficasse altamente discrepante com a realidade, onde quanto maior o valor de k mais a
equação de Erlang tende a uma distribuição normal)mostrou que os modelos são praticamente
idênticos, não havendo diferenças significativas entre os mesmos.
Porém, na simulação, onde mudou-se a taxa de atendimento,fazendo o pico das distribuições
se encontra em intervalos de tempo intermediários, entre 1 e2 min,o modelo M/Ek/c obteve
um melhor desempenho (Erlang a grau 3) quando comparado ao M/M/c. Ou seja, a curva de
Erlang é mais abrangente quando comparada a exponencial negativa, recomendando a sua
utilização para a maioria das situações de atendimento.
Porém, se faz necessário o aprofundamento nas equações de desenpenho de Erlang, visto que
a mesma é dificilmente sitada por autores, na maioria dos casos os mesmo fazem menção a
Erlang para apenas 1 atendente, o que dificulta o aprofundamento no assunto.
Referências
XXXI ENCONTRO NACIONAL DE ENGENHARIA DE PRODUCAO Inovação Tecnológica e Propriedade Intelectual: Desafios da Engenharia de Produção na Consolidação do Brasil no
Cenário Econômico Mundial Belo Horizonte, MG, Brasil, 04 a 07 de outubro de 2011.
12
BOLCH, Gunter / GREINER, Stefan / MEER, Hermann de / TRIVEDI, Kishor S. Queueing Networks and
Markov Chains: Modeling and Performance Evaluation with Computer Science Applications. New York: Wiley-
Interscience, 1998.
FERREIRA, J. O. Simulação de Filas G1/G/m e verificação de aproximações desdás por filas Ph/Ph/m. São
José dos Campos: INPE, 1998.
HILLIER, Frederick S. / LIEBERMAN, Gerald J. Introdução à Pesquisa Operacional. Tradução de Helena L.
Lemos – Rio de Janeiro: Editora Campus / São Paulo: Editora da Universidade de São Paulo, 1988.
MOREIRA, Daniel Augusto. Pesquisa Operacional: Curso Introdutório. 2º Ed. Rev. e Atu. São Paulo: Cengage
Learning. 2010
PRADO, Darci Santos do. Teoria das Filas e da Simulação. Série Pesquisa Operacional – Volume 2. 4º Ed.
Nova Lima: INDG Tecnologia e Serviços Ltda. 2009.
SHAMBLIN, James E. / STEVENS, G.T. JR. Pesquisa Operacional: Uma Abordagem Básica.3º Ed. São
Paulo: Editora Atlas S.A. 1979.
TÔRRES, Oswaldo Fadigas. Elementos da Teoria das Filas. Disponível em:
<http://www16.fgv.br/rae/artigos/2525.pdf>. Acessado em: 10/11/2010.
Top Related