Sebbene non mi opponga all'impossibilità pratica (o impraticabilità) di una risoluzione Monte Carlo di questo problema con una precisione di 6 decimali evidenziata da whuber , penso che sia possibile ottenere una risoluzione con sei cifre di accuratezza.
In primo luogo, seguendo Glen_b , le particelle sono scambiabili in un regime stazionario, quindi è sufficiente (come in sufficienza ) monitorare l'occupazione delle diverse cellule, poiché anche questo costituisce un processo di Markov. La distribuzione delle occupazioni al successivo passaggio è completata determinata dalle occupazioni al momento attuale . Scrivere la matrice di transizione è decisamente impraticabile, ma simulare la transizione è semplice.t Kt + 1tK
Secondo, come notato da shabbychef , si può seguire il processo di occupazione sui 450 quadrati dispari (o pari), che rimangono sui quadrati dispari se si considerano solo i tempi pari, cioè la matrice quadrata di Markov .K2
In terzo luogo, il problema originale considera solo la frequenza di zero occupazioni, , dopo transizioni di Markov. Dato che il punto di partenza ha un valore molto elevato per la distribuzione di probabilità stazionaria della catena di Markov e dato che si concentra su una media unica su tutte le celle, possiamo considerare che la realizzazione della catena al momento è una realizzazione dalla distribuzione di probabilità stazionaria. Ciò comporta una notevole riduzione dei costi di elaborazione, poiché possiamo simulare direttamente da questa distribuzione stazionaria50(X(t)) p 0=1p^050( X( t ))(X(t))t=50π
p^0=1450∑i=1450I0(X(50)i)
(X(t))t=50π, che è una distribuzione multinomiale con probabilità proporzionali a 2, 3 e 4 nell'angolo pari, rispettivamente altre celle sul bordo e celle interne.
Ovviamente, la distribuzione stazionaria fornisce direttamente il numero previsto di celle vuote come
pari a ,166.1069
∑i=1450(1−πi)450
166.1069
pot=rep(c(rep(c(0,1),15),rep(c(1,0),15)),15)*c(2,
rep(3,28),2,rep(c(3,rep(4,28),3),28),2,rep(3,28),2)
pot=pot/sum(pot)
sum((1-pot)^450)-450
[1] 166.1069
che è abbastanza vicino a un'approssimazione di Monte Carlo di [basato su simulazioni 10,, che hanno impiegato 14 ore sulla mia macchina]. Ma non abbastanza vicino per 6 decimali.166.11
Come commentato da whuber , le stime devono essere moltiplicate per 2 per rispondere correttamente alla domanda, quindi un valore finale di 332.2137,