discrete-time Markov chain, esempio

[[Stato della pagina: Bozza]]

Questo pagina presenta un esercizio ispirato da [1], esercizio che si risolve attraverso una semplicissima DTMC (discrete-time Markov chain), con risoluzione transiente.

Previsioni meteo

Supponete di essere in una cittadina il cui sistema meteo può essere così modellizzato:

  • Se oggi è una giornata soleggiante, allora c'è il 60% delle possibilità che sia bel tempo anche domani, 30% che sia nuvoloso e 10% che piova.
  • Se oggi è una giornata nuvolosa, allora c'è il 40% delle possibilità che domani sia soleggiato, il 45% che rimanga coperto e 15% che inizi a piovere.
  • Se oggi è una giornata piovosa, allora c'è il 15% delle possibilità che domani sia una giornata soleggiata, il 60% che sia coperta e 25% che rimanga piovosa.

Dunque, se venerdì piove, quanto è la possibilità (in percentuale) che lunedì sia nuvoloso?

Soluzione

La catena di Markov utilizzabile per questo esercizio può essere rappresentata come segue,
Visualizzazione grafica della dtMc
La matrice delle transizioni,P, si ricava banalmente:

(1)
P = \left[ \begin{array}{ccc} .6 & .3 & .1 \ .4 & .45 & .15 \ .15 & .6 & .25 \end{array} \right]

Passo 0: venerdì

Anche se la catena di Markov è un processo casuale, i dati iniziali sono - ovviamente - deterministici (oggi è venerdì). Dunque, al passo 0, il vettore di probabilità di occupazione degli stati risulta:

(2)
\pi (0) = \left( 0 , 0 , 1 \right)

Passo 1: da venerdì a sabato

Semplicemente,

(3)
\pi (1) = \pi (0) P = \left( 0, 0, 1 \right) \left[ \begin{array}{ccc} .6 & .3 & .1 \ .4 & .45 & .15 \ .15 & .6 & .25 \end{array} \right] = \left( .15,.6,.25 \right)

Passo 2:da sabato a domenica

Continuiamo con il procedimento in modo iterativo:

(4)
\pi (2) = \pi (1) P = \left( .15,.6,.25 \right) \left[ \begin{array}{ccc} .6 & .3 & .1 \ .4 & .45 & .15 \ .15 & .6 & .25 \end{array} \right] = ( .3675,.465,.1675 )

Passo 3: da domenica a lunedì

E finalmente l'ultimo passaggio:

(5)
\pi (3) = \pi (2) P = (.4316, .42, .1484)

Dunque per lunedì è prevista una giornata nuvolosa al 42%1.

blinkbitsBlinkListblogmarksco.mmentsconnoteadel.icio.usDe.lirio.usdiggFarkfeedmelinksFurlLinkaGoGoMa.gnoliaNewsVineNetvouzRawSugarRedditscuttleShadowsSimpySmarkingSpurlTailRankWistsYahooMyWebFacebook

page_revision: 15, last_edited: 1193865857|%e %b %Y, %H:%M %Z (%O ago)
Unless stated otherwise Content of this page is licensed under Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License