В математике случайное блуждание , иногда называемое блужданием пьяницы , — это случайный процесс , описывающий путь, состоящий из последовательности случайных шагов в некотором математическом пространстве .
Элементарным примером случайного блуждания является случайное блуждание по прямой с целыми числами , которая начинается с 0 и на каждом шаге перемещается на +1 или -1 с равной вероятностью . Другие примеры включают путь, прослеживаемый молекулой при ее движении в жидкости или газе (см. Броуновское движение ), путь поиска животного - фуражника или колеблющуюся цену акций и финансовое положение игрока . Случайные блуждания применяются в технике и во многих научных областях, включая экологию , психологию , информатику , физику , химию , биологию , экономику и социологию . Термин «случайное блуждание» был впервые введен Карлом Пирсоном в 1905 году. [1]
Реализации случайных блужданий могут быть получены с помощью моделирования Монте-Карло . [2]
Популярной моделью случайного блуждания является модель случайного блуждания по регулярной решетке, где на каждом шаге местоположение переходит в другой узел в соответствии с некоторым распределением вероятностей. При простом случайном блуждании местоположение может переходить только на соседние узлы решетки, образуя путь решетки . В простом симметричном случайном блуждании по локально конечной решетке вероятности перехода местоположения к каждому из его непосредственных соседей одинаковы. Наиболее изученным примером является случайное блуждание по d -мерной целочисленной решетке (иногда называемой гиперкубической решеткой) . [3]
Если пространство состояний ограничено конечными размерами, модель случайного блуждания называется простым симметричным случайным блужданием с границей , а вероятности перехода зависят от местоположения состояния, поскольку в краевых и угловых состояниях движение ограничено. [4]
Элементарным примером случайного блуждания является случайное блуждание по строке целых чисел, которое начинается с 0 и на каждом шаге перемещается на +1 или -1 с равной вероятностью.
Эту прогулку можно проиллюстрировать следующим образом. Маркер ставится на нулевую позицию на числовой прямой и подбрасывается честная монета. Если выпадает решка, маркер перемещается на одну единицу вправо. Если выпадает решка, маркер перемещается на одну единицу влево. После пяти бросков маркер теперь может находиться на -5, -3, -1, 1, 3, 5. При пяти бросках, трех орлах и двух решках в любом порядке он приземлится на 1. Существует 10 способов приземление на 1 (подбрасывая три орла и две решки), 10 способов приземления на −1 (подбрасывая три решки и две решки), 5 способов приземления на 3 (подбрасывая четыре орла и одну решку), 5 способов приземления на -3 (подбрасывая четыре решки и одну решку), 1 способ приземления на 5 (подбрасывая пять решок) и 1 способ приземления на -5 (подбрасывая пять решок). На рисунке ниже показаны возможные результаты 5 бросков.
Чтобы формально определить это блуждание, возьмите независимые случайные переменные , где каждая переменная равна 1 или −1, с вероятностью 50% для любого значения, и установите и Ряд называется простым случайным блужданием по . Этот ряд (сумма последовательности -1 и единиц) дает чистое пройденное расстояние, если длина каждой части пути равна единице. Ожидание равно нулю . То есть среднее значение всех бросков монеты приближается к нулю по мере увеличения количества бросков. Это следует из свойства конечной аддитивности ожидания:
Аналогичный расчет, использующий независимость случайных величин и тот факт, что , показывает, что:
Это намекает на то , что ожидаемое расстояние перевода после n шагов должно быть порядка . Фактически [5]
Чтобы ответить на вопрос, сколько раз случайное блуждание пересечет границу, если ему разрешено продолжать движение вечно, простое случайное блуждание пересечет каждую точку бесконечное число раз. У этого результата много названий: феномен пересечения уровней , повторение или разорение игрока . Причина такого названия следующая: игрок с конечной суммой денег в конечном итоге проиграет, играя в честную игру против банка с бесконечной суммой денег. Деньги игрока совершят случайное блуждание и в какой-то момент достигнут нуля, и игра будет окончена.
Если a и b — положительные целые числа, то ожидаемое количество шагов до тех пор, пока одномерное простое случайное блуждание, начинающееся с 0, впервые не достигнет b или − a , равно ab . Вероятность того, что это блуждание достигнет b раньше −a , равна , что можно вывести из того факта, что простое случайное блуждание является мартингейлом . И эти ожидания и вероятности попадания могут быть вычислены в общей одномерной цепи Маркова случайного блуждания.
Некоторые из упомянутых выше результатов можно вывести из свойств треугольника Паскаля . Количество различных обходов из n шагов, где каждый шаг равен +1 или -1, равно 2 n . Для простого случайного блуждания каждое из этих блужданий одинаково вероятно. Для того чтобы Sn был равен числу k, необходимо и достаточно, чтобы количество +1 в обходе превышало число −1 на k . Отсюда следует, что +1 должно появиться ( n + k )/2 раза среди n шагов обхода, следовательно, количество проходов, которые удовлетворяют требованиям, равно количеству способов выбора ( n + k )/2 элементов из набора n элементов, [ 6] обозначено . Чтобы это имело смысл, необходимо, чтобы n + k было четным числом, а это означает, что n и k либо оба четные, либо оба нечетные. Следовательно, вероятность равна . Представляя элементы треугольника Паскаля в виде факториалов и используя формулу Стирлинга , можно получить хорошие оценки этих вероятностей для больших значений .
Если для краткости пространство ограничено знаком +, количество способов, которыми случайное блуждание приземлится на любое заданное число с пятью переворотами, может быть показано как {0,5,0,4,0,1}.
Эта связь с треугольником Паскаля продемонстрирована для малых значений n . При нулевых оборотах единственной возможностью будет оставаться на нуле. Однако за один ход есть один шанс приземлиться на -1 или один шанс приземлиться на 1. За два хода маркер с 1 может переместиться на 2 или вернуться к нулю. Маркер в -1 может переместиться в -2 или вернуться к нулю. Следовательно, есть один шанс попасть на −2, два шанса попасть на ноль и один шанс попасть на 2.
Центральная предельная теорема и закон повторного логарифма описывают важные аспекты поведения простых случайных блужданий на . В частности, первое означает, что по мере увеличения n вероятности (пропорциональные числам в каждой строке) приближаются к нормальному распределению .
В качестве прямого обобщения можно рассмотреть случайные блуждания по кристаллическим решеткам (бесконечнократно абелевы накрывающие графы над конечными графами). На самом деле в этой ситуации можно установить центральную предельную теорему и теорему о больших уклонениях. [7] [8]
Одномерное случайное блуждание также можно рассматривать как цепь Маркова , пространство состояний которой задается целыми числами. Для некоторого числа p, удовлетворяющего , вероятности перехода (вероятность P i,j перехода из состояния i в состояние j ) заданы к
Гетерогенное случайное блуждание на каждом временном шаге рисует случайное число, определяющее локальные вероятности прыжка, а затем случайное число, определяющее фактическое направление прыжка. Главный вопрос — вероятность пребывания в каждой из различных площадок после прыжков, причем в пределе этой вероятности, когда она очень велика.
В более высоких измерениях набор случайно пройденных точек обладает интересными геометрическими свойствами. Фактически получается дискретный фрактал , то есть набор, который демонстрирует стохастическое самоподобие в больших масштабах. В небольших масштабах можно наблюдать «зубчатость», возникающую из-за сетки, по которой совершается прогулка. Траектория случайного блуждания — это совокупность посещенных точек, рассматриваемая как набор без учета того, когда блуждание достигло этой точки. В одном измерении траектория представляет собой просто все точки между минимальной высотой и максимальной высотой, достигнутой при прохождении (обе в среднем порядка ).
Чтобы визуализировать двумерный случай, можно представить человека, случайно идущего по городу. Город фактически бесконечен и представляет собой квадратную сетку тротуаров. На каждом перекрестке человек случайным образом выбирает один из четырех возможных маршрутов (включая тот, по которому он изначально шел). Формально это случайное блуждание по множеству всех точек плоскости с целочисленными координатами .
Чтобы ответить на вопрос о том, сможет ли человек когда-нибудь вернуться в исходную отправную точку прогулки, можно использовать двумерный эквивалент проблемы переезда, обсуждавшейся выше. В 1921 году Джордж Полиа доказал, что человек почти наверняка будет совершать двумерное случайное блуждание, но для трех измерений и выше вероятность возвращения к началу координат уменьшается по мере увеличения числа измерений. В трех измерениях вероятность снижается примерно до 34%. [9] Математик Шизуо Какутани, как известно, ссылался на этот результат следующей цитатой: «Пьяный человек найдет дорогу домой, но пьяная птица может заблудиться навсегда». [10]
Другой вариант этого вопроса, который также задал Полья: «Если два человека покинут одну и ту же отправную точку, встретятся ли они когда-нибудь снова?» [11] Можно показать, что разница между их местоположениями (два независимых случайных блуждания) также является простым случайным блужданием, поэтому они почти наверняка встретятся снова в 2-мерном блуждании, но для 3-х измерений и выше вероятность уменьшается с увеличением количество размеров. Пол Эрдеш и Сэмюэл Джеймс Тейлор также показали в 1960 году, что для размерностей, меньших или равных 4, два независимых случайных блуждания, начинающихся из любых двух заданных точек, почти наверняка имеют бесконечное количество пересечений, но для размерностей выше 5 они почти наверняка пересекаются только конечно часто. . [12]
Асимптотическая функция двумерного случайного блуждания при увеличении числа шагов определяется распределением Рэлея . Распределение вероятностей является функцией радиуса от начала координат, а длина шага постоянна для каждого шага. Здесь длина шага предполагается равной 1, N — общее количество шагов, а r — радиус от начала координат. [13]
Винеровский процесс — это стохастический процесс, поведение которого похоже на броуновское движение — физическое явление диффузии мельчайших частиц в жидкости. (Иногда винеровский процесс называют «броуновским движением», хотя, строго говоря, это смешение модели с моделируемым явлением.)
Винеровский процесс — это масштабный предел случайного блуждания в размерности 1. Это означает, что если существует случайное блуждание с очень маленькими шагами, то существует приближение к винеровскому процессу (и, менее точно, к броуновскому движению). Точнее, если размер шага равен ε, нужно пройти прогулку длиной L /ε 2 , чтобы аппроксимировать длину Винера L . Когда размер шага стремится к 0 (и количество шагов пропорционально увеличивается), случайное блуждание сходится к винеровскому процессу в соответствующем смысле. Формально, если B — пространство всех путей длины L с максимальной топологией, и если M — пространство меры над B с нормальной топологией, то сходимость происходит в пространстве M . Точно так же винеровский процесс в нескольких измерениях является пределом масштабирования случайного блуждания в том же количестве измерений.
Случайное блуждание — это дискретный фрактал (функция с целочисленными размерностями; 1, 2, ...), но траектория винеровского процесса — это настоящий фрактал, и между ними существует связь. Например, совершите случайное блуждание, пока не достигнет круга радиуса, умноженного на длину шага. Среднее количество шагов, которые он выполняет, равно r 2 . [ нужна цитата ] Этот факт является дискретной версией того факта, что блуждание винеровского процесса является фракталом хаусдорфовой размерности 2. [ нужна цитата ]
В двух измерениях среднее количество точек, которые одно и то же случайное блуждание имеет на границе своей траектории, равно r 4/3 . Это соответствует тому факту, что граница траектории винеровского процесса представляет собой фрактал размерности 4/3, факт, предсказанный Мандельбротом с помощью моделирования, но доказанный только в 2000 году Лоулером , Шраммом и Вернером . [14]
Винеровский процесс обладает многими симметриями, которых нет у случайного блуждания. Например, блуждание винеровского процесса инвариантно к поворотам, а случайное блуждание - нет, поскольку лежащая в основе сетка не инвариантна (случайное блуждание инвариантно к поворотам на 90 градусов, но винеровские процессы инвариантны к поворотам, например, на 17 градусов). слишком). Это означает, что во многих случаях задачи случайного блуждания легче решить, переведя их в винеровский процесс, решив задачу там, а затем переведя обратно. С другой стороны, некоторые проблемы легче решить с помощью случайных блужданий из-за их дискретного характера.
Случайное блуждание и винеровский процесс могут быть связаны , а именно проявляться в одном и том же вероятностном пространстве зависимым образом, что заставляет их быть достаточно близкими. Простейшей такой связью является вложение Скорохода , но существуют и более точные связи, такие как аппроксимационная теорема Комлоша – Майора – Туснади .
Сходимость случайного блуждания к винеровскому процессу контролируется центральной предельной теоремой и теоремой Донскера . Для частицы, находящейся в известном фиксированном положении при t = 0, центральная предельная теорема говорит нам, что после большого количества независимых шагов в случайном блуждании положение ходока распределяется в соответствии с нормальным распределением полной дисперсии :
где t — время, прошедшее с начала случайного блуждания, — размер шага случайного блуждания, — время, прошедшее между двумя последовательными шагами.
Это соответствует функции Грина уравнения диффузии , которое управляет винеровским процессом, что предполагает, что после большого количества шагов случайное блуждание сходится к винеровскому процессу.
В 3D дисперсия, соответствующая функции Грина уравнения диффузии, равна:
Приравнивая эту величину к дисперсии, связанной с положением случайного блуждающего человека, можно получить эквивалентный коэффициент диффузии, который следует учитывать для асимптотического винеровского процесса, к которому случайное блуждание сходится после большого числа шагов:
Два приведенных выше выражения дисперсии соответствуют распределению, связанному с вектором , который связывает два конца случайного блуждания, в 3D. Отклонение, связанное с каждым компонентом , составляет лишь одну треть от этого значения (все еще в 3D).
Для 2D: [15]
Для 1D: [16]
Случайное блуждание, размер шага которого варьируется в соответствии с нормальным распределением , используется в качестве модели для данных реальных временных рядов, таких как финансовые рынки.
Здесь размер шага представляет собой обратное кумулятивное нормальное распределение , где 0 ≤ z ≤ 1 — равномерно распределенное случайное число, а μ и σ — среднее и стандартное отклонения нормального распределения соответственно.
Если µ не равно нулю, случайное блуждание будет изменяться по линейному тренду. Если v s — начальное значение случайного блуждания, ожидаемое значение после n шагов будет v s + n µ.
В частном случае, когда µ равно нулю, после n шагов распределение вероятностей расстояния перевода определяется как N (0, n σ 2 ), где N () — обозначение нормального распределения, n — количество шагов , а σ — это обратное кумулятивное нормальное распределение, как указано выше.
Доказательство: Гауссово случайное блуждание можно рассматривать как сумму последовательности независимых и одинаково распределенных случайных величин, X i из обратного кумулятивного нормального распределения со средним, равным нулю, и σ исходного обратного кумулятивного нормального распределения:
но у нас есть распределение суммы двух независимых нормально распределенных случайных величин Z = X + Y , которое определяется выражением
В нашем случае µ X = µ Y = 0 и σ 2 X = σ 2 Y = σ 2 дают
Но для гауссовского случайного блуждания это всего лишь стандартное отклонение распределения расстояния перевода после n шагов. Следовательно, если µ равно нулю, и поскольку среднеквадратичное (RMS) расстояние перевода составляет одно стандартное отклонение, существует вероятность 68,27%, что среднеквадратичное расстояние перевода после n шагов упадет между . Аналогично, существует 50% вероятность того, что расстояние перевода после n шагов упадет между .
Число различных сайтов, посещаемых одним случайным блуждающим, тщательно изучалось для квадратных и кубических решеток, а также для фракталов. [17] [18] Эта величина полезна для анализа проблем захвата и кинетических реакций. Это также связано с колебательной плотностью состояний, [19] [20] процессами диффузионных реакций [21] и распространением популяций в экологии. [22] [23]
Скорость передачи информации гауссовского случайного блуждания относительно квадрата расстояния ошибки, т.е. его квадратичная функция искажения скорости , задается параметрически [24]
Как уже упоминалось, диапазон природных явлений, которые были предметом попыток описания с помощью той или иной разновидности случайных блужданий, значителен, особенно в физике [25] [26] и химии, [27] материаловедении , [28] [29] и биологии. . [30] [31] [32] Ниже приведены некоторые конкретные применения случайных блужданий:
Был рассмотрен ряд типов случайных процессов , похожих на чистые случайные блуждания, но в которых простую структуру можно сделать более обобщенной. Чистая структура может быть охарактеризована шагами, определяемыми независимыми и одинаково распределенными случайными величинами . Случайные блуждания могут происходить в различных пространствах, таких как графы , целые числа, вещественная линия, плоскость или векторные пространства более высокой размерности, на искривленных поверхностях или римановых многообразиях более высокой размерности и на группах . Также возможно определить случайные блуждания, которые совершают свои шаги в случайное время, и в этом случае позиция X
тдолжно быть определено для всех моментов времени t ∈ [0, +∞) . Конкретные случаи или пределы случайных блужданий включают модели полета Леви и модели диффузии , такие как броуновское движение .
Случайное блуждание длины k по возможно бесконечному графу G с корнем 0 представляет собой случайный процесс со случайными величинами, такими что и является вершиной, выбранной равномерно случайным образом из соседей . Тогда число — это вероятность того, что случайное блуждание длины k, начинающееся в v , заканчивается в w . В частности, если G — граф с корнем 0 , это вероятность того, что случайное блуждание с шагом -шага вернется к 0 .
Опираясь на аналогию из предыдущего раздела о высших измерениях, предположим теперь, что наш город больше не представляет собой идеальную квадратную сетку. Когда наш человек достигает определенного перекрестка, он с равной вероятностью выбирает между различными доступными дорогами. Таким образом, если на перекрестке семь выходов, человек пройдет к каждому с вероятностью одна седьмая. Это случайное блуждание по графу. Доберется ли наш человек до своего дома? Оказывается, что и в довольно мягких условиях ответ по-прежнему положительный, [42] но в зависимости от графика ответ на вариантный вопрос «Встретятся ли два человека снова?» не может быть, чтобы они встречались бесконечно часто и почти наверняка. [43]
Примером случая, когда человек почти наверняка доберется до своего дома, является случай, когда длины всех блоков находятся между a и b (где a и b — любые два конечных положительных числа). Обратите внимание: мы не предполагаем, что граф плоский , т. е. в городе могут быть туннели и мосты. Одним из способов доказательства этого результата является использование подключения к электрическим сетям . Возьмите карту города и разместите в каждом квартале резистор номиналом 1 Ом . Теперь измерьте «сопротивление между точкой и бесконечностью». Другими словами, выберите некоторое число R , возьмите все точки электрической сети, находящиеся на расстоянии больше R от нашей точки, и соедините их вместе. Теперь это конечная электрическая сеть, и мы можем измерить сопротивление от нашей точки до проводных точек. Возьмем R до бесконечности. Пределом называется сопротивление между точкой и бесконечностью . Оказывается, верно следующее (элементарное доказательство можно найти в книге Дойла и Снелла):
Теорема : граф является переходным тогда и только тогда, когда сопротивление между точкой и бесконечностью конечно. Неважно, какая точка выбрана, если граф связный.
Другими словами, в переходной системе достаточно преодолеть конечное сопротивление, чтобы добраться до бесконечности из любой точки. В рекуррентной системе сопротивление от любой точки до бесконечности бесконечно.
Эта характеристика быстротечности и повторяемости очень полезна и, в частности, позволяет нам проанализировать случай города, нарисованного на плоскости с ограниченными расстояниями.
Случайное блуждание по графу — это особый случай цепи Маркова . В отличие от общей цепи Маркова, случайное блуждание по графу обладает свойством, называемым временной симметрией или обратимостью . Грубо говоря, это свойство, также называемое принципом детального баланса , означает, что вероятности прохождения заданного пути в том или ином направлении имеют между собой очень простую связь (если граф регулярный , они просто равны). Это свойство имеет важные последствия.
Начиная с 1980-х годов, было проведено много исследований по связи свойств графа со случайными блужданиями. Помимо описанного выше подключения к электрической сети, существуют важные связи с изопериметрическими неравенствами , подробнее см. здесь , функциональными неравенствами, такими как неравенства Соболева и Пуанкаре , и свойствами решений уравнения Лапласа . Значительная часть этих исследований была сосредоточена на графах Кэли конечно порожденных групп . Во многих случаях эти дискретные результаты переносятся на многообразия и группы Ли или выводятся из них .
В контексте случайных графов , особенно модели Эрдеша-Реньи , были получены аналитические результаты для некоторых свойств случайных блуждающих. К ним относятся распределение времени первого [44] и последнего попадания [45] пешехода, где время первого попадания определяется моментом, когда пешеход впервые заходит на ранее посещенный участок графа, а время последнего попадания соответствует в первый раз ходок не может выполнить дополнительное движение, не посетив ранее посещенный участок.
Хорошим справочником по случайному блужданию по графам является онлайн-книга Олдоса и Филла. Информацию о группах см. в книге Вёсса. Если ядро перехода само по себе является случайным (в зависимости от среды ), то случайное блуждание называется «случайным блужданием в случайной среде». Когда закон случайного блуждания включает в себя случайность , закон называется отожженным законом; с другой стороны, если закон рассматривается как фиксированный, его называют угасшим законом. См. книгу Хьюза, книгу Ревеса или конспекты лекций Зейтуни.
Мы можем думать о выборе всех возможных ребер с той же вероятностью, что и о локальном максимизации неопределенности (энтропии). Мы также могли бы сделать это глобально — в случайном блуждании с максимальной энтропией (MERW) мы хотим, чтобы все пути были равновероятными, или, другими словами: для каждых двух вершин каждый путь заданной длины одинаково вероятен. [46] Это случайное блуждание имеет гораздо более сильные свойства локализации.
Существует ряд интересных моделей случайных путей, в которых каждый шаг сложным образом зависит от прошлого. Все они более сложны для аналитического решения, чем обычное случайное блуждание; тем не менее, поведение любой модели случайного блуждающего человека можно получить с помощью компьютеров. Примеры включают в себя:
Самоизбегающий обход длины n — это случайный путь из n шагов, который начинается в начале координат, совершает переходы только между соседними узлами в , никогда не посещает узел повторно и выбирается единообразно среди всех таких путей. В двух измерениях из-за самозахвата типичный путь самоизбегания очень короткий, [48] тогда как в более высоких измерениях он выходит за все пределы. Эта модель часто использовалась в физике полимеров (с 1960-х годов).
Случайное блуждание, выбранное для максимизации уровня энтропии , имеет гораздо более сильные свойства локализации.
Случайные блуждания, при которых направление движения в один момент времени коррелирует с направлением движения в следующий раз. Он используется для моделирования движений животных. [53] [54]