Научная статья на тему 'РАСПРЕДЕЛЕНИЕ ВРЕМЕНИ ДО НАЧАЛА ФИНАЛЬНОЙ ОСТАНОВКИ ДИФФУЗИОННОГО ПОЛУМАРКОВСКОГО ПРОЦЕССА НА ИНТЕРВАЛЕ С НЕДОСТИЖИМЫМИ ГРАНИЦАМИ'

РАСПРЕДЕЛЕНИЕ ВРЕМЕНИ ДО НАЧАЛА ФИНАЛЬНОЙ ОСТАНОВКИ ДИФФУЗИОННОГО ПОЛУМАРКОВСКОГО ПРОЦЕССА НА ИНТЕРВАЛЕ С НЕДОСТИЖИМЫМИ ГРАНИЦАМИ Текст научной статьи по специальности «Математика»

CC BY
21
4
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕПРЕРЫВНЫЙ ПОЛУМАРКОВСКИЙ ПРОЦЕСС / ПОЛУМАРКОВСКАЯ ЦЕПЬ / ФИНАЛЬНЫЙ ИНТЕРВАЛ ПОСТОЯНСТВА / ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ / ПРЕОБРАЗОВАНИЯ ЛАПЛАСА / НЕДОСТИЖИМАЯ ГРАНИЦА / ИНТЕГРАЛЬНОЕ ПРЕДСТАВЛЕНИЕ

Аннотация научной статьи по математике, автор научной работы — Харламов Борис Павлович, Расова София Станиславовна

Рассматривается одномерный процесс с непрерывными траекториями, заданный на неотрицательной полуоси, обладающий однородным марковским свойством относительно момента первого выхода из любого открытого интервала (полумарковский процесс). За исключением слова «непрерывные» это определение совпадает с определением полумарковского процесса с кусочно-постоянными траекториями. Непрерывные полумарковские процессы служат математической моделью многих физических, биологических и социальных явлений. Диффузионность процесса состоит в том, что вероятность первого выхода на любую из двух границ симметричной окрестности начальной точки процесса стремится к 1/2 при стремлении диаметра этой окрестности к нулю. Исследуется распределение момента начала финального интервала постоянства выборочной траектории процесса. Так называется бесконечный интервал постоянства, определение которого опирается на вид полумарковских переходных производящих функций процесса. Получено представление производящей функции этого распределения в интегральном виде. Подынтегральное выражение этого представления объясняет смысл квадратичного члена разложения полумарковской производящей функции процесса по степеням диаметра симметричной окрестности начальной точки процесса при стремлении этого диаметра к нулю. А именно, траектории процесса не имеют финального интервала постоянства тогда и только тогда, когда коэффициент этого квадратичного члена равен нулю.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Харламов Борис Павлович, Расова София Станиславовна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

TIME DISTRIBUTION FROM ZERO UP TO BEGINNING OF THE NAL STOP OF SEMI-MARKOV DI USION PROCESS ON INTERVAL WITH UNATTAINABLE BOUNDARIES

A one-dimensional process with continuous trajectories on non-negative semi-axis is considered. The process has the Markov property with respect to the first exit time from any open interval (semi-Markov process). This process is called to be diffusion if probability for its first exit point from any symmetric neighborhood of its initial point across any boundary tends to 1/2 as length of this neighborhood tends to zero. Time distribution from zero up to beginning of the final interval of constancy is investigated. This distribution depends on semi-Markov transition generating functions of the process. Representation for Laplace transform of this distribution is obtained in an integral form. The integrand of this representation explains sense of quadratic members of Tailor decomposition of a semi-Markov transition generating function by powers of diameter of symmetric neighborhood of the process initial point. Namely trajectory of the process has no any final interval of constancy if and only if coefficient of such a quadratic member is equal to zero.

Текст научной работы на тему «РАСПРЕДЕЛЕНИЕ ВРЕМЕНИ ДО НАЧАЛА ФИНАЛЬНОЙ ОСТАНОВКИ ДИФФУЗИОННОГО ПОЛУМАРКОВСКОГО ПРОЦЕССА НА ИНТЕРВАЛЕ С НЕДОСТИЖИМЫМИ ГРАНИЦАМИ»

УДК 519.2 Вестник СПбГУ. Математика. Механика. Астрономия. 2022. Т. 9 (67). Вып. 3

МБС 60040

Распределение времени до начала финальной остановки диффузионного полумарковского процесса на интервале с недостижимыми границами

Б. П. Харламов, С. С. Расова

Институт проблем машиноведения РАН,

Российская Федерация, 199178, Санкт-Петербург, Большой пр. В. О., 61

Для цитирования: Харламов Б. П., Расова С. С. Распределение времени до начала финальной остановки диффузионного полумарковского процесса на интервале с недостижимыми границами // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2022. Т. 9(67). Вып. 3. С. 517-526. https://doi.org/10.21638/spbu01.2022.312

Рассматривается одномерный процесс с непрерывными траекториями, заданный на неотрицательной полуоси, обладающий однородным марковским свойством относительно момента первого выхода из любого открытого интервала (полумарковский процесс). За исключением слова «непрерывные» это определение совпадает с определением полумарковского процесса с кусочно-постоянными траекториями. Непрерывные полумарковские процессы служат математической моделью многих физических, биологических и социальных явлений. Диффузионность процесса состоит в том, что вероятность первого выхода на любую из двух границ симметричной окрестности начальной точки процесса стремится к 1/2 при стремлении диаметра этой окрестности к нулю. Исследуется распределение момента начала финального интервала постоянства выборочной траектории процесса. Так называется бесконечный интервал постоянства, определение которого опирается на вид полумарковских переходных производящих функций процесса. Получено представление производящей функции этого распределения в интегральном виде. Подынтегральное выражение этого представления объясняет смысл квадратичного члена разложения полумарковской производящей функции процесса по степеням диаметра симметричной окрестности начальной точки процесса при стремлении этого диаметра к нулю. А именно, траектории процесса не имеют финального интервала постоянства тогда и только тогда, когда коэффициент этого квадратичного члена равен нулю.

Ключевые слова: непрерывный полумарковский процесс, полумарковская цепь, финальный интервал постоянства, дифференциальное уравнение, преобразования Лапласа, недостижимая граница, интегральное представление.

1. Полумарковский метод. В настоящей работе рассматриваются так называемые полумарковские процессы диффузионного типа, переходные производящие функции которого удовлетворяют дифференциальному уравнению

±Г + А(х)Г-В( А,*)/ = 0, (1)

© Санкт-Петербургский государственный университет, 2022 https://doi.org/10.21638/spbu01.2022.312

где А(х) — непрерывно дифференцируемая функция, Б(Л, х) — неотрицательная функция, непрерывная по второму аргументу и имеющая вполне монотонную производную по первому аргументу Л > 0 (см. [1, с. 167]).

Полумарковским процессом мы называем случайный процесс, обладающий однородным марковским свойством относительно момента первого выхода из любого открытого множества (см., например, [1] и [2]). В этом кратком определении подразумевается, что процесс принимает значения в топологическом пространстве и что его область определения допускает односторонние сдвиги. Только в этом случае можно корректно применять термин «момент марковской регенерации» — марковский момент, относительно которого выполняется однородное марковское свойство. Класс полумарковских процессов диффузионного типа привлекает к себе внимание как благодаря прикладной стороне, т. е. многочисленным приложениям, так и благодаря математической стороне, т. е. специальным методам исследования, генерирующим новые задачи.

Область применения таких процессов — моделирование физических, биологических и социальных явлений — та же, что и у марковских диффузионных процессов, которые составляют более узкий класс. Но полумарковское моделирование возможно и тогда, когда марковское свойство представляется слишком обременительным, например когда нужно принимать в расчет замену времени, не сохраняющую марковское свойство процесса. Такая ситуация возникает, например, в хроматографии (см. [3]). В этом популярном методе исследования предполагается, что при движении частиц анализируемого вещества в длинном и тонком сосуде исходный броуновский характер движения превращается в движение с остановками — конечными интервалами постоянства, которые возникают в результате пребывания частиц на стенках сосуда. Движение происходит в результате чередования интервалов адсорбции и моментов десорбции.

Специальная методика исследования возникла благодаря невозможности использовать для немарковского процесса аппарат теории марковских процессов. Поэтому для вычисления даже стандартных функционалов приходится применять нестандартный набор полумарковских переходных функций. Однако существуют задачи, когда такой полумарковский подход оказывается полезным и для марковских процессов. Особенно наглядно польза полумарковского подхода проявляется в задачах о замедленном отражении и вообще — в задачах о поведении процесса на границах между различными средами.

В настоящей работе рассматривается случайный процесс, выборочные траектории которого имеют предел на бесконечности. Частным случаем существования такого предела является наличие так называемой финальной остановки процесса — интервала постоянства бесконечной длины. Особый интерес представляет ситуация, когда не существует никакого сигнала к началу этого бесконечного интервала — он начинается спонтанно, как результат некоторого случайного выбора. Такой момент начала финальной остановки не является марковским моментом относительно натуральной фильтрации.

Нас будет интересовать также значение процесса после начала финальной остановки (см. [4]). Эта точка фазового пространства имеет большое прикладное значение. Достаточно заметить, что в результате таких остановленных движений накапливаются горы осадочной породы на Земле. Недостаток концептуальной модели такого физического процесса состоит в том, что она не является марковской.

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

2. Случайный процесс и его атрибуты. Каждый случайный процесс рассматривается далее как случайный элемент измеримого пространства V функций £ : [0, то) ^ К, непрерывных справа и имеющих предел слева в каждой точке Ь > 0. Предполагается, что на V определены метрика Скорохода (точнее — Стоуна— Скорохода) и борелевская сигма-алгебра Т подмножеств множества V, совпадающая, как известно, с сигма-алгеброй, порожденной натуральной фильтрацией (Тг). Количественное описание случайного процесса задается его вероятностной мерой или согласованным семейством вероятностных мер.

Популярным представлением случайного процесса является функция от двух аргументов X(Ь,£) = Х^£). Такое обозначение обычно сокращают до X(Ь) = Хг. Это означает, что случайный процесс может быть представлен семейством функций (Хг) одного аргумента, где Хг(£) = £(Ь). Семейство функций (Хг) согласовано с семейством операторов сдвига вг, где вг отображает V в себя так, что при любых неотрицательных в и Ь справедливо (Х3вг)(£) = Хг+3(£). Отсюда следует, что

= е^(£).

Обозначим через Т множество всех измеримых неубывающих функций т. На множестве т < то определим функцию Хт и оператор вт, где

Хт (£) = Хт ю(£), вт (£) = вт

а также функцию

т + Т1 = т + Т\вт,

где т"1 € Т. Заметим, что операция + ассоциативна, но не перестановочна. При т + т1 < то из определения следует

вт+т1 = вт1 вт, Хт+т1 = Хт1 вт . (2)

Докажем последнее равенство. Имеем

Хт+т1 (£) = Хт ю+т^т «)(0 = Хт1вт ш(вт (£)) = Хт1 (вт (£)),

что и требовалось доказать.

Рассмотрим одномерный случайный процесс Х(Ь) (Ь > 0) со значениями в некотором интервале До = (а0,Ь0). Пусть (Рх) (х € Д0) — семейство распределений (вероятностных мер) процесса, индексированное начальными точками процесса: Рх(Х (0)= х) = 1.

Обозначим через а а оператор первого выхода из множества Д С До. Таким образом, рассматривается отображение аА : V ^ [К+] = [0, то], где

ал(0 = М{I > 0: ф) 0 Д};

и, по определению, ад(£) = то, если множество в фигурных скобках пусто. Известно, что а а является марковским моментом относительно натуральной фильтрации (см. [6, с. 194]).

3. Полумарковский процесс. Для марковского момента т стандартным образом определяем сигма-алгебру Тт всех событий, предшествующих этому моменту.

Марковское свойство относительно марковского момента т проще всего определить в терминах математических ожиданий относительно семейства вероятностных мер Рх. Момент т называется моментом марковской регенерации, если

Ех(М2вт); т < то) = Ех{1± Ех(т)Ш; т < то),

где /1 — Тт-измеримая и /2 — ^"-измеримая функции.

Семейство вероятностных мер Рх называется согласованным, если у него существует хотя бы один нетривиальный момент марковской регенерации.

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

Известно, что для такого процесса любая конечная итерация таких моментов первого выхода также является моментом марковской регенерации. Пусть Д = (а, Ь), где ао < а < Ь < Ьо. Обозначим

9а(А,х) = Ех(ехр(-\а а ); а а < то, X (а а) = а),

На(А,х) = Ех(ехр(-\аА); а а < то, X (аА) = Ь),

где А > 0. Это так называемые переходные производящие функции процесса. Из полумарковского свойства процесса следует, что

9а2(а,х) = ЙА1'(а,х)§а2(А, с) + ьа1'(а,х)§а^(А,й), (3)

на2 (А, х) = дА1 (а,х)на2 (а,с) + на^ (а,х)на1' (А,Л), (4)

где Д1 = (с,3), Д2 = (а, Ь) и Д1 С Д2.

Рассмотрим ((£) = т^о > 0 : (Ш > 1о) £(1) = £{Ьо)). Такой случайный момент называется моментом начала финального интервала постоянства. Он не является марковским моментом, но является терминальным моментом. Это означает, что для любого т <еТ на множестве т < то справедливо равенство

С = т + С-

3.1. Диффузионный полумарковский процесс. Пусть С — множество всех непрерывных функций £ <EV с областью определения [0, то).

Рассмотрим конечный интервал До = (ао, Ьо) с недостижимыми границами для полумарковского процесса с начальной точкой х внутри этого интервала. Недостижимость означает, что для любой начальной точки х процесса справедливо равенство Рх(а а0 = то) = I.

Определение. Непрерывный полумарковский процесс называется диффузионным полумарковским процессом, если при r ^ 0

9{x-r,x+r){\x) = t¿(1 ~A(x)r - В(Х, х)г2) + о(г2),

h(x-r,x+r){\x) = + A{x)r - B(\x)r2) +o(r2)

равномерно по Л > 0 и по x в каждом интервале (a, b), где ao < a <b <bo.

Из работы [7] следует, что процесс является диффузионным полумарковским процессом тогда и только тогда, когда его переходные производящие функции удовлетворяют дифференциальному уравнению (1).

Условие 1. При любом y G (ao,bo) выполняются равенства

h(a0,y)(\ У) = 9(y,b0)(\ У) = 1

Условие 2. При любом y G До существуют производные

h(ao,y)(Л, У ) > 0 9(уМ(Л,У+) < 0

Из определения функций Н(а0,ь)(Л,у) и 9(а,ь0)(Л,У) следует, что первая из них не убывает на интервале (a, bo), а вторая не возрастает на интервале (a, bo). Таким образом, условие 2 предписывает существование производных на соответствующих концах интервалов и их знаки.

Очевидно, что если выполняются условия 1 и 2, то при r ^ 0 справедливы равенства

\а,у)(Л,У — r) = 1 — h(a,y)(\y-)r + o(r),

9(y,e)(\ У + r) = 1 + 9(y,e)(\ У+)г + o(r)■

Теорема. Пусть переходные производящие функции диффузионного полумарковского процесса удовлетворяют дифференциальному уравнению (1) и выполняются условия 1 и 2. Тогда производящая функция начала финального интервала постоянства процесса имеет интегральное представление

Ех(ехр(-Л(); ( < ж) =

= Г 2B(0,y)h(atX)(\,y) dy + rb° 2B(0,y)g(Xtb)(X,y) L0h[a y]{\,y-)-g[y b]{\,y+) У Jx h'{a y]{X,y-) - g'{y b]{X,y+) У'

3.2. Доказательство теоремы. Пусть {c1,...,cm} (m > 1) — конечная последовательность точек таких, что ao < ci < ■ ■■ < cm < bo. При этом ck — ck-i = L/(m + 1), где L = bo — ao, 1 < k < m, co = ao. Такая последовательность точек порождает равномерное разбиение интервала Д0 на m +1 частей.

Это равномерное разбиение определяет так называемое правильное покрытие интервала Д0 интервалами Дk (1 < k < m), где Дk = (ck-1, ck+1), co = ao, cm+1 = bo.

Пусть x G Дк при некотором k. Предыдущее правильное покрытие порождает последовательность пар (pn,X(pn)) (n > 0), где po = 0, X(po) = x, а другие члены (если есть) определяются правилом: 1) если pk < ж, то pk+1 = pk + , 2) если

Рк < то и рк+1 < то, то X(рк+1) = X(аАк+1 )0рк, 3) если рк < то и рк+1 = то, то X(рк+\.) не определено (последовательность пар определена до номера к).

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

Обозначим через (т начало финального интервала постоянства этой полумарковской цепи. Очевидно, что Рх(Ст < С, 11тСт = 0 = 1.

Пусть в таково, что 1 < в < т, св = х. Рассмотрим {Дая, Дя, Дья} — «неправильное» покрытие интервала Д, состоящее из трех областей: Дая = (ао,ся), Дя = (ся-1,ся+1), Дья = (ся,Ьо). Точки св-1, ся, ся+1 назовем начальными точками этих интервалов.

Пусть то — момент первого достижения процессом точки св. Предыдущее неправильное покрытие порождает последовательность пар (тптп)), п > 0, где то = 0, X(то) = ся, а другие члены (если есть) определяются правилом:

1) если тк < то и X(тк) = Ся-1, то тк+1 = тк + аАаз,

2) если тк < то и X(тк) = Ся+1, то тк+1 = тк + аАЪв,

3) если тк < то и X(тк) = с8, то тк+1 = тк + аАв,

4) если тк = то, то тк+1 = то и X(тк+1) не определено (последовательность пар определена до последнего номера к, при котором тк < то).

Благодаря полумарковскому свойству исходного процесса эти пары определяют полумарковскую цепь с финальной остановкой, вложенную в исходный полумарковский процесс. Эту цепь можно рассматривать как случайный процесс с кусочно-постоянными (разрывными) траекториями. Однако начало финальной остановки этой полумарковской цепи в общем случае не совпадает с Ст, а число возможных состояний равно трем: ся-1, ся, ся+1.

Распределение этой полумарковской цепи (заданной с точностью до начального состояния) определяется матрицей математических ожиданий Q(А) порядка 3 х 3 переходных производящих функций процесса:

( 0 НА(ая)(А,Ся-1) 0

Q(А) = I дА(в)(А,се) 0 На(я)(А, Ся)

\ 0 дА(Ья)(А,Ся+1) 0

Заметим, что ся-1 и ся+1 — начальные точки интервалов Дая и Дья соответственно. Поэтому моменты первого выхода из этих интервалов возможны только через одну и ту же точку св.

Очевидно, что при то < то справедливо равенство РСз (то = 0, т1 = аА3) = 1. Относительно меры РСд справедливо утверждение: если т1 < то, то X(т1) = ся-1 или X(т1) = ся+1, а также определено т2 = т1 + аАая или т2 = т1 + аАЬя соответственно. В обоих этих случаях справедливо равенство X(т2) = ся, когда т2 < то, что соответствует матрице Q2(А), а именно, распределение случайной величины т2 определяется интегралом

Q222(А) = Ес3 (ехр(-Ат2); т2 < то, Xт) = С2),

где ^П2(А) — элемент матрицы Qn(А) (п > 1), стоящий на пересечении второй строки и второго столбца. В данном случае Q22 = 0 и Q22(А) = Вя(А), где

DS(X) = hAas (X,cs-1) gAs (\,cs) + hAs (X,cs) gAb(X,cs+i). Очевидно, что DS(X) < 1 при X > 0.

Заметим, что после попадания второго компонента вложенной полумарковской цепи в состояние cs применяется правило построения вложенной полумарковской цепи с этим начальным состоянием. Учет этого свойства позволяет записать равенство

3.3. Регулярная последовательность равномерных разбиений исходного интервала. Назовем последовательность равномерных разбиений регулярной, если каждое следующее разбиение в этой последовательности получено в результате деления на две равные части элементов предыдущего разбиения.

Для нас важно то, что число точек разбиения в регулярной последовательности равномерных разбиений стремится к бесконечности, поэтому расстояния между соседними точками разбиений стремятся к нулю. Заметим, что каждая из внутренних точек разбиения cs = y (с меняющимся номером, но с фиксированным положением внутри исходного интервала) находится посередине каждого интервала из некоторой последовательности интервалов с убывающими длинами, а также что в правильном покрытии каждая точка из несчетного числа точек покрывается дважды. Исключение составляют только счетное множество границ интервалов покрытия и две краевые области разбиения интервала Д, длины которых стремятся к нулю.

3.3.1. Определение членов, зависящих от начальной точки. При т0 < ж разделим траекторию процесса на две части: одну — на интервале от нуля до т0 и вторую — от то до бесконечности. После этого мы используем терминальное свойство момента начала финальной остановки:

Ех(ехр(-Х(); ( < то) =

= Ех(ехр(-Хто)(ехр(-Х()вТо); то < то, (0То < то, X(()вТо = ев). Пользуясь полумарковским свойством процесса, получаем

Ех(ехр(-Х(); ( < то) = Ех(ехр(-Хто); то < то) ЕСв (ехр(-ХС); С < то). (5) 3.3.2. Определение члена, зависящего от точки е3. Имеем равенства Ес3 (ехр(-Х(т); < то, X((т) = ев) =

X (т2п) = Q22W = d:(x).

Z = то + Z,

откуда следует равенство

= Ecs (ехр(—Хтп); Тп < ж, X(тп) = ci, Тп+i = ж) =

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

п=0

Ecs (exp(-XT2n); Т2п < ж, X(Т2п) = Cs, T2n+i = ж)

п=о

Pcs (Tl = ж) I 1 + Ecs (ехр(-Хт2п); Т2п < ж, X (Т2п) = Cs)

п=1

Получаем формулу для полумарковской цепи:

P ÍT1 = оо)

ECs (exp(-ACm); Cm < с», X{Qm) = cs) = D {X) ■ (6)

3.4. Интегральное представление. Рассмотрим отношение

PCs(t 1 = 00) 1 — Ds(\) '

(7)

Пусть y = cs и r = L/(m +1). Из определения диффузионного полумарковского процесса следует, что

9AS (0, у) = giy-r,y+r)( О, У) = \( 1 " " В( 0, у)г2 + о(г2)),

(О, у) ее fe(y_ri!/+r)(0, у) = i (1 + А(у)г - В(0, Уу + о(г2))

равномерно по y € А = (а, 6), где А(ж), В (А, ж) — коэффициенты исходного дифференциального уравнения. Отсюда числитель определяется как

py (ti = ж) = 1 - h{y_r,y+r)(0,y) - g(y-r,y+r)(0,y) = В(0, уУ + оУ~). В знаменателе имеем

Ds(A) = Ь(ао,Св)(А cs-l ) gAs (А cs) + g(cs,b0)(A cs + l) hAs (А cs).

Пусть cs — cs-i = cs+i — cs = r. Полагая

h(ao,cs)(A cs-1) = h(a0,cs)(A cs) — rh(ao (\ cs) + o(r),

g(.cs,bo)(A cs+i) = g(cs,bo)(A cs) + ^(csM^ cs) + o(r), где первые члены этих формул Тейлора равны единице, получаем

^-Ds{X) = \{h[a^v){\y-)-g[vM{\y+))r + o{r).

Следовательно,

2B(0,y)r2 + o(r2) 2B(0,y)r + o(r)

(Ка0,у)(Х У-) - д'ум^ У+))Г + 0(г) (к[а0,у)(Х,У-) - 9(у,,Ь0)(Х,У+))г + 0(Г) При этом дробь в последнем выражении при г ^ 0 имеет предел, равный

_2В(0,у)_

к[а0,у)^У-) ~ 9[уМ{\У+У

Учитывая множители, зависящие от х, суммируя последнее выражение по в, где 1 < I < т, и переходя к пределу при т ^ то, получаем

Ex(exp(-\(); ( < œ) = lim Ex(exp(-A(2m); (2™ < œ) =

_ fx 2B{0,y)hM{X,y) [b° 2B(0,y)g(xM)(X,y)

Ja0 h(a,y)(X'V-) - g'(y,b0)(X,y+) Jx h(a,y)(X,V-) - 9(yM)(X,V+)' Теорема доказана. □

Следствие 1. Если выполняются неравенства

_h(a0,x)(X,y)_ >о _9{х,Ьо){\у)_ >0

к(а у)(Х,У-) - 9(у,Ь0)(.Х,У+) ' к(а, у)(Х,У-) - 9(у,Ь0)(Х,У+)

почти для всех у € До, то исходный диффузионный полумарковский процесс не имеет бесконечного интервала постоянства тогда и только тогда, когда В(0, х) = 0 почти для всех х € Д0.

Литература

1. Харламов Б. П. Непрерывные полумарковские процессы.. Санкт-Петербург, Наука (2001).

2. Harlamov B. P. Continuous semi-Markov processes. London, ISTE, Wiley (2008).

3. Harlamov B. P. Stochastic model of gas capillary chromatography. Communications in ¡Statistics — ¡Simulation and Computation, no. 41, 1023—1031 (2011).

4. Харламов Б. П. Финальное распределение диффузионного процесса с остановкой. Записки научных семинаров ПОМИ, № 431, 209-241 (2014).

5. Дынкин Е. Б. Марковские процессы.. Москва, ФМ (1963).

6. Гихман И. И., Скороход А. В. Теория случайных процессов. Т. 2. Москва, Наука (1973).

7. Расова С. С., Харламов Б. П. Об одном локальном свойстве одномерного линейного дифференциального уравнения второго порядка. Дифференциальные уравнения и процессы управления, №2, 65-79 (2021).

Статья поступила в редакцию 24 января 2022 г.;

доработана 1 марта 2022 г.; рекомендована к печати 3 марта 2022 г.

Контактная информация:

Харламов Борис Павлович — д-р физ.-мат. наук; b.p.harlamov@gmail.com Расова София Станиславовна — канд. физ.-мат. наук; s.rasova@gmail.com

Time distribution from zero up to beginning of the final stop of semi-Markov diffusion process on interval with unattainable boundaries

B. P. Harlamov, S. S. Rasova

Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences, 61, Bolshoy pr. V. O., St Petersburg, 199178, Russian Federation

For citation: Harlamov B. P., Rasova S. S. Time distribution from zero up to beginning of the final stop of semi-Markov diffusion process on interval with unattainable boundaries. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2022, vol. 9(67), issue 3, pp. 517-526. https://doi.org/10.21638/spbu01.2022.312 (In Russian)

A one-dimensional process with continuous trajectories on non-negative semi-axis is considered. The process has the Markov property with respect to the first exit time from any open interval (semi-Markov process). This process is called to be diffusion if probability for

its first exit point from any symmetric neighborhood of its initial point across any boundary tends to 1/2 as length of this neighborhood tends to zero. Time distribution from zero up to beginning of the final interval of constancy is investigated. This distribution depends on semi-Markov transition generating functions of the process. Representation for Laplace transform of this distribution is obtained in an integral form. The integrand of this representation explains sense of quadratic members of Tailor decomposition of a semi-Markov transition generating function by powers of diameter of symmetric neighborhood of the process initial point. Namely trajectory of the process has no any final interval of constancy if and only if coefficient of such a quadratic member is equal to zero.

Keywords: continuous semi-Markov process, semi-Markov chain, final interval of constancy, differential equation, Laplace transformation, unattainable boundary, integral representation.

References

1. Harlamov B.P. Continuous semi-Markov processes. St Petersburg, Nauka Publ. (2001). (In Russian)

2. Harlamov B.P. Continuous semi-Markov processes. London, ISTE, Wiley (2008).

3. Harlamov B. P. Stochastic model of gas capillary chromatography. Communications in Statistics — Simulation and Computation, no. 41, 1023—1031 (2011).

4. Harlamov B. P. Final distribution of a diffusion process with stop. Zapiski nauchnykh seminarov POMI, no. 431, 209-241 (2014). (In Russian)

5. Dynkin E. B. Markov processes. Moscow, FM Publ. (1963). (In Russian)

6. Gikhman 1.1., Skorokhod A. V. Theory of random processes. Vol. 2. Moscow, Nauka Publ. (1973). (In Russian)

7. Rasova S. S., Harlamov B. P. On one local property of one-dimension linear differential equation of the second order. Differential equations and control processes, no. 2, 65-79 (2021). (In Russian)

Received: January 24, 2022 Revised: March 1, 2022 Accepted: March 3, 2022

Authors' information:

Boris P. Harlamov — b.p.harlamov@gmail.com Sofia S. Rasova — s.rasova@gmail.com

i Надоели баннеры? Вы всегда можете отключить рекламу.