Хотя я не возражаю против практической невозможности (или нецелесообразности) решения этой проблемы методом Монте-Карло с точностью до 6 десятичных знаков, указанных указателем whuber , я думаю, что может быть достигнуто разрешение с шестизначной точностью.
Во-первых, после Glen_b частицы могут обмениваться в стационарном режиме, поэтому достаточно (как при достаточности ) контролировать занятость разных ячеек, поскольку это также составляет марковский процесс. Распределение занятых мест на следующем шаге времени завершено, определяется занятостью в текущий момент времени t . Написание матрицы переходов K определенно нецелесообразно, но моделировать переход просто.т + 1TК
Во-вторых, как отметил Шаббычеф , можно следить за процессом заполнения на 450 нечетных (или четных) квадратах, который остается на нечетных квадратах при рассмотрении только четных времен, т. Е. Квадратовой матрицы Маркова .К2
В- третьих, исходная задача рассматривает только частоту нулевых , после 50 марковских переходов. Принимая во внимание , что начальная точка имеет очень высокое значение для стационарного распределения вероятностей цепи Маркова ( Х ( т ) ) , и при условии , что фокус на одном среднем по всем клеткам, р 0 = 1п^050( Х( т ))можно считать, что реализация цепочки(X(t))в момент времениt=50является реализацией из стационарного распределения вероятностей. Это значительно снижает стоимость вычислений, поскольку мы можем смоделировать непосредственно из этого стационарного распределенияπ, которое является многочленным распределением с вероятностями, пропорциональными 2, 3 и 4 на четном угле, другими ячейками на краю и внутренними ячейками. соответственно.
п^0= 1450Σя = 1450я0( Х( 50 )я)
( Х( т ))t = 50π
Σя = 1450( 1 - πя)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
166,11
Как прокомментировал whuber , оценки должны быть умножены на 2, чтобы правильно ответить на вопрос, следовательно, окончательное значение 332,2137,