Processos Estocásticos
Quinta Lista de Exercı́cios
12 de fevereiro de 2014
1 Suponha que um organismo unicelular pode estar somente em dois estágios distintos A ou B. Um indivı́duo
no estágio A passa para o estágio B com uma taxa exponencial α. Um indivı́duo no estágio B se divide
em dois novos indivı́duos de tipo A com uma taxa exponencial β. Defina uma CTMC apropriada para a
população desses organismos e determine os parâmetros para esse modelo (tempo gasto em cada estado e
probabilidades de transição).
Seja NA (t) o número de organismos no estágio A no instante t e, de forma análoga, seja NB (t) o número de
organismos no estágio B. Assim, podemos tomar os pares ordenados hNA (t), NB (t)i para representar os estados
da CTMC. Isto é, o PE {hNA (t), NB (t)i, t ≥ 0} é uma CTMC com parâmetros
vhn,mi = αn + βm
αn
Phn,mi;hn−1,m+1i =
αn + βm
βn
Phn,mi;hn+2,m−1i =
αn + βm
.
2 Considere duas máquinas que são reparadas por um único técnico. A máquina i funciona por um tempo
exponencial com taxa λi antes de quebrar (com i = 1, 2). Os tempos de reparo (para quaisquer das duas
máquinas) são exponenciais com taxa µ. É possı́vel analisar esse sistema como um PNM? Se sim, quais são
os parâmetros? Se não, seria possı́vel analisá-lo?
Este não é um PNM porque somente a informação de quantas máquinas estão funcionando em um dado instante
não é suficiente para se analisar o sistema. Para tal é necessário saber também quais máquinas estão funcionando
em cada instante. Assim, podemos construir uma CTMC com os seguintes estados:
• 0: ambas as máquinas estão funcionando.
• 1: máquina 1 funcionando, 2 com defeito.
• 2: máquina 2 funcionando, 1 com defeito.
• 3: ambas máquinas com defeito, máquina 1 sendo consertada.
• 4: ambas máquinas com defeito, máquina 2 sendo consertada.
Essa CTMC pode ser representada pelo seguinte diagrama:
µ
µ
1
3
λ1
λ2
0
λ2
λ1
µ
2
µ
1
4
A CTMC acima possui os parâmetros tempo gasto em cada estado
v0 = λ 1 + λ 2
v1 = λ1 + µ
v2 = λ2 + µ
v3 = v4 = µ
e probabilidades de transição
P01 =
P10 =
µ
λ1 + µ
λ2
λ1 + λ2
P14 =
P02 =
λ1
λ1 + λ2
λ1
λ1 + µ
P20 =
P31 = P42 = 1
µ
λ2 + µ
P23 =
λ2
λ2 + µ
.
3 Considere um PNM com taxas de chegadas λi = (i + 1)λ e taxas de saı́da µi = iµ, i ≥ 0.
a. Determine o tempo esperado para se sair do estado 0 e chegar no estado 4.
b. Determine o tempo esperado para se sair do estado 2 e chegar no estado 5.
Começando com E[T0 ] =
i = 1, 2, 3, 4. Assim
1
1
1
µi
= , usamos a identidade E[Ti ] =
+ E[Ti−1 ] para computar E[Ti ] para
λ0
λ
λi
λi
µ
µ1
1
µ 1
1
1
+
E[T0 ] =
+
· =
·
λ1
λ1
2λ 2λ λ
2λ
λ
2
1
µ2
1
2µ 1
µ
µ
1
E[T2 ] =
+
E[T1 ] =
+
·
·
·
=
λ2
λ2
3λ 3λ 2λ
λ
3λ
λ
2
3
3µ 1
µ
µ
1
µ3
1
1
+
·
·
·
E[T3 ] =
+
E[T2 ] =
=
λ3
λ3
4λ 4λ 3λ
λ
4λ
λ
E[T1 ] =
e, generalizando, vemos que
E[Ti ] =
1
·
(i + 1)λ
i
µ
λ
e portanto,
E[T4 ] =
1
·
5λ
4
µ
λ
.
a. Como Ti é o tempo que o processo leva para transitar do estado i para o estado i + 1, o tempo esperado é
E[T0 ] + E[T1 ] + E[T2 ] + E[T3 ].
b. Mesma explicação do item anterior: E[T2 ] + E[T3 ] + E[T4 ].
4 Assume-se que cada indivı́duo de uma população procria com uma taxa exponencial λ e morre com uma taxa
exponencial µ. Além disso, há uma taxa exponencial θ de crescimento da população devido à imigração.
No entanto, a imigração não é permitida se o tamanho da população é maior ou igual a N . Modele essa
situação como um PNM.
Tomando cada estado X(t) como o tamanho da população no instante t, temos um PNM com
λn = nλ + θ,
n<N
λn = nλ,
n≥N
µn = nµ
.
2
5 Uma pequena barbearia, operada por um único barbeiro, pode acomodar no máximo dois clientes ao mesmo
tempo. Clientes em potencial chegam com uma taxa de Poisson de 3 por hora e os tempos de serviço são
VAs exponenciais independentes com média de 1/4 hora.
a. Qual é o número médio de clientes na barbearia?
b. Qual é a proporção de clientes em potencial que entram na loja?
Tomando o número de clientes na barbearia como o estado, temos um PNM com
λ0 = λ1 = 3,
µ1 = µ2 = 4 .
Graficamente temos
λ0
λ1
0
1
µ1
2
µ2
e calculando as equações de fluxo para os estados 0 e 2, vem
3
π0
4
2
3
3
µ2 π2 = λ1 π1 ⇒ π2 = π1 =
π0
4
4
λ0 π0 = µ1 π1 ⇒ π1 =
.
Como π0 + π1 + π2 = 1, podemos resolver o sistema de equações para π0
2
3
π0 = 1
4
16
37
.
2 3
3
30
π1 + 2π2 =
+2
· π0 =
4
4
37
.
3
π0 + π0 +
4
⇒
π0 =
a. O número médio de clientes na barbearia é
b. A proporção de clientes em potencial que entram na loja é
λ(1 − π2 )
9 16
28
= 1 − π2 = 1 −
·
=
λ
16 37
37
.
6 Um centro de atendimento é composto por dois servidores, cada um trabalhando com uma taxa exponencial
de dois serviços por hora. Clientes chegam com uma taxa de Poisson de três por hora. Assuma que a
capacidade do centro é de no máximo três clientes.
a. Que fração dos clientes em potencial entram no sistema?
b. Qual seria o valor do item anterior se houvesse somente um servidor no sistema com uma taxa duas
vezes mais rápida, isto é, µ = 4?
Tomando o número de clientes no centro como o estado, temos um PNM com
λ0 = λ1 = λ2 = 3,
µ1 = 2,
3
µ2 = µ3 = 4 .
Assim, as equações de fluxo se reduzem a
3
π0
2
3
9
π2 = π1 = π0
4
8
27
3
π0
π3 = π2 =
4
32
π1 =
e portanto
−1
3 9 27
32
π0 = 1 + + +
=
2 8 32
143
.
a. A fração de clientes em potencial que entram no sistema é
λ(1 − π3 )
27 32
116
= 1 − π3 = 1 −
·
=
≈ 81.1%
λ
32 143
143
.
b. Com um único servidor trabalhando duas vezes mais rápido temos um PNM com
λ0 = λ1 = λ2 = 3,
µ1 = µ2 = µ3 = 4 .
Agora, as equações de fluxo se reduzem a
3
π0
4
2
3
3
π0
π2 = π1 =
4
4
3
3
3
π3 = π2 =
π0
4
4
π1 =
e portanto
2 3 −1
3
3
3
64
π0 = 1 + +
+
=
4
4
4
175
.
E finalmente, a nova fração de clientes que entram no sistema é
1 − π3 = 1 −
148
27 64
·
=
≈ 84.6%
64 175
175
.
7 Considere um ponto de taxi onde taxis e clientes chegam de acordo com processos de Poisson com respectivas taxas de um e dois por minuto. Um taxi fica sempre em espera, independente do número de taxis já
parados no ponto. Entretanto, um cliente que chega e não encontra um taxi disponı́vel vai embora, isto é,
não existe uma fila de espera de clientes.
a. Calcule o número médio de taxis esperando.
b. A proporção de clientes que chegam e conseguem um taxi.
Sejam os estados denotados pelo número de taxis esperando. Assim, temos um PNM com λn = 1 e µn = 2. Note
que essa CTMC corresponde a um modelo de fila M/M/1, cujas métricas já são conhecidas.
a. A média de taxis esperando é o tamanho médio da fila =
1
.
µ−λ
b. A proporção de clientes que chegam e conseguem um taxi é a proporção de clientes que chegam e encontram
ao menos um taxi esperando. A taxa de chegada desses clientes é 2(1 − π0 ). A proporção dessas chegadas
4
é portanto
2(1 − π0 )
λ
λ
1
= 1 − π0 = 1 − 1 −
= =
2
µ
µ
2
.
8 Para uma fila M/M/1, calcule
a. o número esperado de chegadas durante um perı́odo de serviço; e
b. a probabilidade de que nenhum cliente chegue durante um perı́odo de serviço.
a. Seja S uma VA indicando o tempo de serviço. Como a taxa de serviço é exponencial, sabemos que E[S] =
1/µ. O número esperado de chegadas buscado corresponde então a E[λS] = λE[S] = λ/µ. (O segundo
passo das equações é justificado pelo fato de λ ser uma constante com relação ao tempo de serviço, pois
chegadas e saı́das são eventos independentes.)
b. Novamente tomamos S como uma VA indicando o tempo de serviço. Buscamos a probabilidade condicional
de haver 0 chegadas dado que o perı́odo de serviço é S. Mas essa probabilidade é exatamente igual à
probabilidade o servidor terminar antes de uma nova chegada (caso contrário o sistema muda de estado e o
µ
.
tempo S é “resetado” – tente entender o motivo). Assim, a probabilidade pedida é
λ+µ
9 As máquinas de uma fábrica quebram com uma taxa exponencial de 6 por hora. A fábrica emprega apenas
um técnico que conserta as máquinas com uma taxa exponencial de 8 por hora. O custo causado pela
produção perdida quando há máquinas com defeito é de $10 por hora por máquina. Qual é o custo médio
causado por máquinas defeituosas?
Este problema pode ser modelado por uma fila M/M/1 com λ = 6 e µ = 8. O custo médio é dado por
$10 por hora por máquina × número médio de máquinas quebradas.
Mas o número médio de máquinas quebradas é exatamente L, o tamanho da fila, cujo valor já foi calculado:
L=
λ
6
= =3 .
µ−λ
2
Portanto, o custo médio é de $30 por hora.
10 Considere um sistema M/M/1 onde clientes chegam com taxa λ e são servidos com taxa µ. No entanto,
assuma que em qualquer momento que o servidor estiver ocupado existe uma probabilidade dele quebrar,
levando o sistema a parar. A probabilidade de quebra é descrita por uma taxa exponencial α. Quando o
sistema para, todos os clientes que estavam no sistema partem e não são mais permitidas chegadas até que
o defeito for consertado. O tempo de reparo é exponencialmente distribuı́do com taxa β.
a. Defina os estados apropriadamente.
b. Descreva as equações de balanço de fluxo.
a. O estados são n (n ≥ 0) e b. O estado n indica que há n clientes no sistema e o estado b que uma quebra
ocorreu. A CTMC que modela o sistema descrito é dada pelo diagrama abaixo.
5
b
β
α
λ
0
α
λ
1
µ
α
α
λ
λ
n
2
µ
µ
···
µ
b. As equações de fluxo são
απ0 = µπ1 + βπb
(λ + µ + α)πn = λπn−1 + µπn+1
βπb = α(1 − π0 ) .
6
n+1
···
Download

Resolução_Lista05