О способах описания скачкообразных процессов

УДК 530.1

О СПОСОБАХ ОПИСАНИЯ СКАЧКООБРАЗНЫХ

ПРОЦЕССОВ

Как известно, описание случайных процессов диффузионного типа базируется на уравнениях второго порядка (Фоккера - Планка для вероятностей или Шредингера для амплитуд вероятности). В работе [1] были рассмотрены обобщения этих уравнений на случай скачкообразных процессов - псевдодифференциальные уравнения Фоккера -Планка и Шредингера (дифференциальные уравнения бесконечного порядка) и приведены точные решения этих уравнений в некоторых простых частных случаях.

Здесь мы рассмотрим связь псевдодифференциального уравнения Фоккера - Планка со стохастическими дифференциальными уравнениями. Это позволит нам построить обобщение процесса Орнштейна - Уленбека, играющего важную роль в физических приложениях (в частности, в теории броуновского движения), на случай конечных скачков.

Напомним, что формулировке теории диффузионных процессов на основе уравнения Фоккера - Планка

отвечает эквивалентная формулировка на основе стохастического дифференциального уравнения (СДУ) Ланжевена

где А(х^) - коэффициент сноса, связанный с действием неслучайной силы, В(х, ¿) интенсивность случайной силы Ь[1) = - производной винеровского процесса. Т.е.

А. Л. Шелепин

Рассматриваются стохастические дифференциальные уравнения и отвечающие им псевдодифференциальные уравнения Фоккера - Планка для скачкообразных марковских процессов. Приводятся точные решения для аналога процесса Орнштейна - Уленбека.

с1х

(2) позволяет дать описание произвольного диффузионного процесса с помощью простейшего непрерывного процесса W(t). Однако как уравнение Фоккера - Планка, так и уравнение Ланжевена со случайной силой L(t) = W'(t) описывают процессы с непрерывными траекториями и не описывают скачкообразные процессы (процессы с разрывными траекториями, например, процессы Коши и Пуассона).

Для описания скачкообразных марковских процессов необходимо использовать более общие уравнения, чем (1) и (2). Так, в уравнении (2) производная винеровского процесса (гауссов белый шум) может быть заменена на производную некоторых простейших скачкообразных процессов. Вместо же уравнения второго порядка (1) надо рассматривать дифференциальное уравнение бесконечного порядка, записываемое обычно в форме разложения Крамерса - Мойала

ОО /_ 1\п n=l n!

Kn= lim f(x - x')np(x,t + At\x',t)dx, At—»0 Дt J

где коэффициенты I\n(x,t) называют моментами перехода [2, 3].

Отметим, что моменты перехода на самом деле тесно связаны с кумулянтами, а не моментами соответствующего процесса, как об этом можно было бы судить исходя из названия. Кумулянты (семиинварианты) определяются как коэффициенты к в разложении логарифма характеристической функции f(k) по степеням А:,

00 (ИЛп

=

п=1 П"

В [4] показано, что либо все моменты перехода, кроме первых двух, равны нулю, либо отличны от нуля все Кп при четных ?г, причем выполняется неравенство

1<1 < Кп-2Кп+2, п > 4, четное.

Аналогично, либо все кумулянты, кроме первых двух, равны нулю, либо имеется бесконечно много отличных от нуля кумулянтов (теорема Марцинкевича, см. [3, 5]). Особенно проста связь между моментами перехода и кумулянтами в случае постоянных коэффициентов Кп. В этом случае (3) для характеристической функции дает

00 (гк)п

&ДМ) = А

п=1 п!

Сопоставляя (4) и (6), получим

««= Кп

Скачкообразные марковские процессы можно, в свою очередь, разделить на два широких класса:

1. все скачки имеют одинаковую амплитуду (естественным для этого случая является дискретное множество возможных состояний);

2. амплитуда различных скачков различна (непрерывное множество возможных состояний).

Рассмотрим коротко первый случай, к которому относятся, например, пуассоновский процесс и случайное блуждание по прямой. Им отвечают уравнения

д1р(х, г) = 7[р(х + /, г) - р(х, г)], дгр(х, г) = 7[р{х + /, I) + р(х -/,*)- 2р{х,«)],

где / - амплитуда скачков (шаг), 7 - частота скачков (вероятность шага в единицу времени). Представляя перемещение из точки х в х + 1 действием оператора сдвига е'9*, можно переписать их в виде дифференциальных уравнений бесконечного порядка

дгр{х, <) = 7

Эти уравнения с постоянными коэффициентами могут быть легко решены путем перехода к уравнениям на характеристические функции (см. [3]). Разлагая их правые части в ряд, получим для четных моментов перехода соответственно Кп = 7/п или Кп = 2-у1п и точное равенство = /\п_2Л'п+2 как предельный случай неравенства (5). Нечетные моменты отвечают за снос. Кумулянты процесса Пуассона связаны с моментами перехода Кп соотношением (7), кп = 7/"(£ — ¿о)-

Роль малого параметра в разложении (8) по степеням дх играет величина скачка I - в пределе частых малых скачков уравнение для случайного блуждания переходит в уравнение диффузии, а уравнение для пуассоновского процесса - в уравнение Лиувилля.

Для описания более сложных скачкообразных процессов можно применять СДУ, основанные, однако, на использовании в качестве базисного уже не винеровского процесса, а простейшего скачкообразного процесса - процесса Пуассона. Так, при рассмотрении дробового шума используется аналогичное (2) СДУ, где флуктационная сила выражается через производную процесса Пуассона, Ь{1) — А^'(^) [3]. Как уже отмечалось, это описание пригодно лишь в том случае, если все скачки имеют одинаковую

фиксированную амплитуду (случайным является время прихода ¿-импульсов, вызывающих скачки).

Перейдем теперь к случаю скачков произвольной величины. Действительно, для многих процессов амплитуда различных скачков неодинакова, например, если мы рассмотрим изменение скорости за счет соударений при тепловом движении. Здесь можно говорить лишь о верхнем пределе на амплитуду скачка, определяемом скоростью теплового движения молекул.

Как подчеркивается в [2], вывод уравнения Фоккера - Планка основан на математи ческом доказательстве Колмогорова, в котором предполагается, что имеется бесконечно много бесконечно малых скачков; на этом же предположении основано обычное рассмотрение диффузии. Поэтому при обобщении подхода Ланжевена естественным шагом представляется замена винеровского процесса И/(<) на более общий с характер-

ной (максимальной) амплитудой скачков А, в пределе частых малых скачков переходящий в винеровский а при А —> оо - в процесс Коши И^^), характеризующийся бесконечной дисперсией (амплитуда скачков не ограничена).

Указанным требованиям удовлетворяет рассмотренный в [1] процесс, который мы ниже будем называть А-процессом и обозначать как И/д

д1Р

где А имеет смысл величины, а с/А - частоты скачков. Характеристическая функция /(к,£) удовлетворяет уравнению

д/(М)/3* = (с/А) (1 - х/Г+№) /(¿,0,

/

Принимая во внимание начальное условие р(х, ¿|х0, ¿о) = ¿>(х — х0), откуда следует, что /(Аг, ¿о) = ехр(г/гх0), получим

/

Плотность вероятности находится с помощью обратного преобразования Фурье р(х,<) = (1/27г) /ехр(—1кх)}{к^)йк,

Р

ТТЛ

где К\{х) - функция Макдональда, х = х — х0, I = £ — ¿0 и для вычисления интеграла использована формула (2.5.39.2) из [7].

Предельным случаям отвечают хорошо известные процессы. При Л —> оо

1 с{

р

с характеристической функцией ехр(—Ы\к\) (подробнее о процессе Коши см. [6]). Ам плитуда скачков в этом случае не ограничена.

В пределе Л —> 0, с — оо и конечном Лс (9) переходит в уравнение диффузии

дгр(х,Цх о, ¿о) = ^\сд1р(х,Цх0,и),

что соответствует гауссову процессу с коэффициентом диффузии Б = Лс (винеровский процесс получается при Лс = 1). Гауссов процесс, таким образом, отвечает пределу бесконечно частых бесконечно малых скачков.

Используя безразмерные параметры £ = (ж — Жо)/Л и т — с{Ь — ¿о)/Л, распределение (12) (мы будем называть его Л-распределением) можно переписать как

= ^е^ + е)-1/2Кг[

ТТЛ

Из (14) видно, что форма графика функции плотности Л-распределения зависит от одного параметра г = — ¿о)/Л; Л задает естественный масштаб по х. Дисперсия Л-распределения сг2 = \с{Ь — ¿о) [1]. В зависимости от значения т можно выделить три характерных случая: а) т << 1, б) т « 1, в) г >> 1. В случае а) имеем распределение Коши с полушириной (половиной расстояния между точками х\, х2, для которых /(ж) = /шаж/2) с(Ь — ¿о); крылья до значений х яз Л соответствуют убыванию по закону 1/х2. В случае б) распределение Коши описывает центральную часть, крылья убывают по степенному закону типа 1/ж", промежуточному между 1 /ж2 и экспоненциальным. В случае в) справедливо нормальное распределение со среднеквадратичным отклонением сг = ^/Лс(* — ¿о) •

Согласно (7) нечетные кумулянты Л-процесса равны нулю (снос отсутствует), а четные к2п = ^Г(п — 1/2)сЛ2п_1(£ — г0)- В частности, первые два момента < >= их = 0, < И^д^о)^^) >= к2 = сЛ(< — ¿о) совпадают с таковыми для винеровского процесса с Б = Лс.

Рассмотрение СДУ для скачкообразных процессов начнем с процесса, описываемого уравнением

с1х

— = А

где - коэффициент сноса, связанный с действием неслучайной силы, а //(¿) =

- производная Л-процесса. Процесс Ь{1) (негауссов белый шум) имеет те же два первых момента, что и гауссов белый шум, < /,(£) >= 0, < Ь(10)Ь{1) >= с\8{1 — ¿0), но в отличие от гауссова белого шума, который характеризуется нулевыми кумулянтами кп при п > 2, здесь четные кумулянты отличны от нуля. Уравнению (15) отвечает ПдФП уравнение вида

дгр

Если коэффициент сноса не зависит от х и Ь, А(х^) = а, то решение (16) отличается от решения (12) свободного уравнения (9) лишь заменой х —> х — а{Ь — ¿0), что полностью аналогично диффузионнму случаю.

Рассмотрим теперь случай А(х^) = 'ух, которому в диффузионном пределе отвечает процесс Орнштейна - Уленбека,

адм) = дх

Ему отвечает уравнение для характеристической функции /

й/

Это уравнение первого порядка может быть решено методом характеристик, для чего запишем вспомогательные уравнения

<И ¿к с?/

/ -гк (с/ А)(1 — \/1 + X2k2)f

Интегрирование первого из них дает кс 11 — С — const, второго - / = [(1 + Vl + X2k2)/2]~c^Xl exp[-(c/A7)v/l + А2к2]Сг, где Сх = Сх{ке~^). Сх как функция ке~* находится из начального условия 0|х, 0) = ¿(ж — xq), /(к, 0) = exp(ikxo). В результате вычислений имеем

ДМ) = ехр

(Vl + Х2к2е- Vl + А2к2 + In 1 + гхо

А7 V 1 + VI +Х2к2е~2^ J

Решение это достаточно громоздко, тем не менее дисперсия (находимая по формуле М(х2) = —д2/(к^)/дк2\к~о) совпадает с таковой для диффузионного процесса Орн-штейна - Уленбека, а2 = ^(1 — ехр(—2^)).

Стационарное решение (которое легко получить и непосредственно из (18)) имеет

вид

/(&, оо) = ехр

А7 V 2

В пределе Л —» 0 и при конечном сЛ (много малых скачков) получаем характеристическую функцию нормального распределения ехр(—сХк2/^), а в пределе Л —> оо характеристическую функцию распределения Коши ехр( —с|/с|/7) с полушириной с/7.

Как известно, процесс Орнштейна - Уленбека дает описание рэлеевской частицы - броуновской частицы, рассматриваемой в более мелкой временной шкале, когда в качестве стохастической функции берется скорость, а не координата частицы. Уравнение Фоккера - Планка для эволюции плотности функции распределения /(и, 2), где V -скорость броуновской частицы массы М, в этом случае записывается в виде [2, в]

|/<М) - 7 + ££)/

Здесь кТ = ту ц/2, где _ средняя скорость теплового движения молекул, а параметр 7 пропорционален частоте столкновений. Это уравнение имеет стационарное решение

/

В частности, при М — тп стационарная функция распределения

/

представляет собой одномерное распределение Максвелла. Несмотря на то, что при М = т мы получили правильный ответ в пределе t —> оо, уравнение (22) применимо лишь для достаточно тяжелых частиц с массой М/т » 1 (изменение скорости за одно столкновение должно быть незначительным). Для молекул же газа изменение скорости за одно соударение имеет тот же порядок, что и сама скорость, и диффузионное приближение (22) не является удовлетворительным.

Однако и рассмотренное выше псевдодифференциальное уравнение (17) или соответствующее ему СДУ (15) также неадекватны задаче (в частности, не дают правильного стационарного распределения (23)). Дело в том, что снос в этих уравнениях, задаваемый коэффициентами А(у^) = вызван действием некоторой внешней неслучайной

силы (или бесконечно частыми и малыми соударениями), а не случайными соударениями, приводящими к скачкообразному изменению скорости, как в рассматриваемом нами случае. Уравнения вида (17) и (15) описывают случай, когда имеется скачкообразная "диффузия" и плавный снос.

Следуя [1], обратимся к рассмотрению уравнения

dtp

Потребуем, чтобы псевдодифференциальное уравнение Фоккера - Планка вида (25) переходило в (22) в диффузионном пределе т/М —> 0. Это требование позволяет определить коэффициенты Л(г;), с(и), и мы получаем псевдодифференциальное уравнение

dtp

где 7 = l/vо - частота столкновений, a / длина свободного пробега. Функция /(х, оо) (23) является, как нетрудно проверить, стационарным решением псевдодифференциального уравнения (26), описывающего скачкообразные процессы и переходящего в (22) при т/М -> 0.

Если в диффузионном пределе (уравнение Фоккера - Планка) снос задается коэффициентом при первой степени дифференциального оператора дх, то в уравнении ПдФП снос может быть введен двояко: либо по-прежнему, как в диффузионном пределе, что отвечает действию неслучайной силы, либо ненулевыми коэффициентами при бесконечном числе нечетных степеней с)х, что отвечает скачкообразным изменениям скорости при соударениях частиц.

Проведенное выше рассмотрение приводит к следующему выводу. В диффузионном пределе вид уравнения Фоккера - Планка или Ланжевена не зависит от того, обусловлен снос действием некоторой внешней неслучайной силы или теми же соударениями, что и диффузия. Если же соударения не являются частыми и малыми, а приводят к конечным скачкообразным изменениям физических величин, то мы получим различные уравнения в зависимости от того, чем обусловлен снос.

Псевдодифференциальные уравнения, описывающие скачкообразные процессы, могут быть записаны в различных видах в зависимости от порядка расстановки операторов х и дх. Выше мы использовали две формы записи - симметричную (25), (26) и разложение Крамерса - Мойала. Разложение Крамерса - Мойала в диффузионном пределе дает уравнение Фоккера - Планка в форме Ито (производные слева от функций

Кп(х,1)), а уравнения (25), (26) - уравнение Фоккера - Планка в форме Стратоновича (симметричной форме).

Весьма интересным является вопрос о соотношении псевдодифференциальных уравнений

<1х = В{х,1)<1\Ух

Такое представление скачкообразных процессов, возможно, упростило бы их исследование, т.к. при решении СДУ (27) можно было бы опираться на известные свойства процесса Ид^). Интерпретация Стратоновича здесь может оказаться предпочтительнее, т.к. уравнение ПдФП в симметричной форме более удобно для нахождения стационарного решения, чем разложение Крамерса - Мойала.

Автор признателен С. А. Решетняку за полезное обсуждение.

ЛИТЕРАТУРА

[1] Ш е л е п и н А. Л. ЯФ, 60, N 2, 265 (1997).

[2] В а н К а м п е н Н. Г. Стохастические процессы в физике и химии. М., Высшая школа, 1990.

[3] Г а р д и н е р К. В. Стохастические методы в естественных науках. М., Мир, 1986.

[4] Р a w u 1 a R. F. Trans. IEEE, IT-3, 33 (1967).

[5] R a j a g о р а 1 А. К., Sudarshan Е. С. G. Phys. Rev., АЮ, 1S52 (1974).

[6] Ф е л л е р В. Введение в теорию вероятностей и ее приложения, т. 2. М., Наука, 1967.

[7] Прудников А. П., Брычков Ю. А., М а р и ч е в О. И. Интегралы и ряды. Элементарные функции. М., Наука, 1981.

[8] Б а л е с к у Р. Равновесная и неравновесная статистическая механика, т. 2. М., Мир, 1978.

Поступила в редакцию 1 марта 1998 г.