Аннотация:
Пусть {Sn,n⩾0} – случайное блуждание, приращения которого принадлежат без центрирования области притяжения устойчивого распределения индекса α, т.е. существует такой процесс {Yt,t⩾0}, что Snt/an⇒Yt, t⩾0, при n→∞ для некоторых нормирующих констант an. Предполагая, что S0=o(an) и Sn⩽φ(n)=o(an), мы доказываем ряд условных предельных теорем для распределения случайной величины Sn−m, предполагая, что m=o(n) и min0⩽k⩽nSk⩾0. Эти теоремы дополняют утверждения, установленные Ф. Каравенной и Л. Шамоном в 2013 г. Полученные результаты используются при исследовании размера популяции критического ветвящегося процесса, эволюционирующего в неблагоприятной случайной среде.
Библиография: 28 названий.
Ministry of Science and Technology (MOST) of China
G20221740071
Исследование В. А. Ватутина выполнено в МЦМУ МИАН при финансовой поддержке Минобрнауки России (соглашение № 075-15-2022-265), а также Ministry of Science and Technology of the People's Republic of China (грант № G20221740071). Исследование К. Донга выполнено при поддержке Ministry of Science and Technology of the People's Republic of China (грант № G20221740071). Исследование Е. Е. Дьяконовой выполнено в МЦМУ МИАН при финансовой поддержке Минобрнауки России (соглашение № 075-15-2022-265).
Исследование условных распределений случайных блужданий при условии их положительности или неотрицательности имеет долгую историю (см., например, статьи [4], [6]–[12], [14], [17], [18], [20], [25]). Интерес к доказательству различных принципов инвариантности для условных процессов, чему и посвящено большинство указанных работ, обусловлен не только естественным развитием теории случайных блужданий, но и широким использованием таких результатов в теории ветвящихся процессов, эволюционирующих как в постоянной, так и в случайной средах (см. [2], [19], [26]), в статистической физике, в частности, для изучения случайных полимеров (см. [13]), а также в других областях.
В настоящей работе мы также анализируем свойства случайных блужданий, рассматриваемых при условии их неотрицательности. Наше исследование мотивировано работой Ф. Каравенны и Л. Шамона [11]. Авторы статьи [11] показали, что если распределение шага случайного блуждания принадлежит (без центрирования) области притяжения устойчивого закона, то при условии неотрицательности такого блуждания на отрезке длины n распределение соответствующим образом нормированной экскурсии длины n, порожденной этим случайным блужданием, сходится к распределению экскурсии устойчивого процесса Леви, рассматриваемого при условии его неотрицательности на отрезке [0,1]. Мы дополняем этот результат исследованием поведения подходящим образом нормированной траектории случайного блуждания в левой окрестности n−m экскурсии, где m=o(n). Оказалось, что в зависимости от скорости стремления отношения m/n к нулю в пределе возникает три различных распределения. Эти три режима приводят к трем различным скоростям роста числа частиц в критических ветвящихся процессах, эволюционирующих в неблагоприятной случайной среде и рассматриваемых при условии их невырождения к далекому моменту времени n. Наши утверждения для критических ветвящихся процессов в случайной среде дополняют соответствующие теоремы, установленные в работах [26], [27] и [28], в которых при условии невырождения таких процессов к моменту n были исследованы распределения числа частиц в них в моменты времени m=o(n) и m=[nt], 0<t⩽1.
Статья организована следующим образом.
В § 2 мы вводим основные понятия, описываем наши основные условия, накладываемые на случайные блуждания, напоминаем некоторые известные локальные предельные теоремы и доказываем вспомогательные утверждения для случайных блужданий при условии их неотрицательности.
В § 3 содержится несколько условных предельных теорем для экскурсий решетчатого случайного блуждания при условии его неотрицательности.
В § 4 мы доказываем соответствующие условные предельные теоремы для экскурсий случайных блужданий, шаг которых имеет абсолютно непрерывное распределение.
В § 5 мы устанавливаем условную предельную теорему для последовательности случайных величин, сходящихся почти наверное.
В § 6 мы применяем полученные результаты для случайных блужданий при условии их неотрицательности для описания распределения размера популяции критического ветвящегося процесса, эволюционирующего в неблагоприятной среде при условии его выживания к далекому моменту времени.
§ 2. Обозначения и предположения
Для сохранения единства изложения и удобства ссылок мы будем в основном придерживаться обозначений и предположений, введенных в работе [11].
В последующем изложении через C1,C2,… мы будем обозначать некоторые абсолютные константы, не обязательно совпадающие в различных формулах.
Введем множества N:={1,2,…} и N0:=N∪{0}. Для двух положительных последовательностей {cn,n∈N} и {dn,n∈N} будем, как обычно, писать:
cn∼dn, если limn→∞(cn/dn)=1;
cn=o(dn), если limn→∞(cn/dn)=0;
cn=O(dn), если limsupn→∞(cn/dn)<∞;
cn≍dn, если 0<lim infn→∞(cn/dn)⩽lim supn→∞(cn/dn)<∞.
Напомним, что положительная последовательность {cn,n∈N} (или действительная функция c(x)) называется правильно меняющейся на бесконечности с индексомγ∈R, что обозначается как cn∈Rγ или соответственно c(x)∈Rγ, если cn∼nγl(n) (соответственно c(x)∼xγl(x)), где l(x) – медленно меняющаяся функция, т.е. такая положительная функция, что l(cx)/l(x)→1 при x→∞ для любого фиксированного c>0.
– подмножество в пространстве R2. Для пары (α,β)∈A и случайной величины X мы будем использовать запись X∈D(α,β), если распределение величины X принадлежит (без центрирования) области притяжения устойчивого закона с плотностью g(x)=gα,β(x), x∈(−∞,+∞), и характеристической функцией
с независимыми и одинаково распределенными приращениями. Всюду далее в этом параграфе мы предполагаем, что случайное блуждание (S={Sn,n∈N0},P) удовлетворяет следующим условиям.
Условие A1. Приращения Xn, n=1,2,…, случайного блуждания S принадлежат D(α,β) с |β|<1.
Это означает, в частности, что существует такая возрастающая последовательность положительных чисел
an:=n1/αℓ(n),n=1,2,…,
где последовательность ℓ(1),ℓ(2),… медленно меняется на бесконечности, что при n→∞
L{S[nt]an,t⩾0}⟹L{Yt,t⩾0};
символ ⟹ обозначает сходимость по распределению в пространстве D[0,+∞) с топологией Скорохода, а процесс (Y={Yt,t⩾0},P) является строго устойчивым, т.е. имеет стационарные независимые приращения и маргинальные распределения, описываемые характеристическими функциями
EeiwYt=Gα,β(wt1/α),t⩾0.
Известно (см., например, [5; с. 380]), что если Xnd=X∈D(α,β) при всех n∈N, то существует предел
limn→∞P(Sn>0)=ρ=P(Y1>0),
где
ρ=12+1παarctg(βtgπα2).
Пусть ΩRW:=RN0 – пространство непрерывных справа функций со скачками в целочисленных точках неотрицательной полуоси, и пусть Ω:=D([0,∞),R) – пространство непрерывных справа и имеющих конечные пределы слева действительных функций на множестве [0,∞). Снабженное топологией Скорохода пространство Ω становится польским пространством с соответствующей σ-алгеброй борелевских множеств. Положим также ΩRWN:=R{0,1,…,N} и Ωt:=D([0,t],R) для t∈(0,∞).
Обозначим через Px закон распределения случайного блуждания, начинающегося в точке x∈R, т.е. вероятностный закон на ΩRW, порождаемый случайным блужданием S+x, где P – закон распределения случайного блуждания S, начинающегося в нуле; а для процесса Y={Yt,t⩾0} с распределением P на Ω будем обозначать через Pa распределение смещенного процесса Y+a для a∈R.
Условие A2. Будем предполагать, что выполнено одно из следующих ограничений.
∙ ((h;c)-решетчатый случай.) Мера P приращения X1 случайного блуждания сосредоточена на решетке c+hZ, где шаг h>0 выбран максимальным (т.е. распределение X1 не сосредоточено на какой-либо решетке c0+h0Z при h0>h и c0∈R). Отметим, что мы можем взять c∈[0,h).
∙ (абсолютно непрерывный случай.) Мера P приращения X1 случайного блуждания абсолютно непрерывна относительно меры Лебега на множестве R, причем существует n∈N такое, что плотность распределения fn(x):=P(Sn∈dx)/dx случайной величины Sn ограничена (и, следовательно, fn(x)∈L∞).
Положим
LN:=min1⩽k⩽NSk.
Законом распределения моста длиныN∈N положительного случайного блуждания, начинающего в точке x∈[0,∞) и заканчивающегося в точке y∈[0,∞), мы будем называть вероятностный закон на множестве ΩRWN, задаваемый равенством
P↑,Nx,y(⋅):=Px(⋅∣LN⩾0,SN=y).
Ясно, что в решетчатом случае для корректности определения закона распределения P↑,Nx,y(⋅) необходимо предположить положительность вероятностей
q+N(x,y):=Px(LN⩾0,SN=y)>0.
Аналогично, в абсолютно непрерывном случае нам нужно обеспечить положительность функции
где K(N−1):={s1>0,…,sN−1>0} – область интегрирования, а f(⋅)=f1(⋅) – плотность приращения X1 случайного блуждания.
Для t∈[0,∞) и a,b∈[0,∞) обозначим через P↑,ta,b закон распределения на множестве Ωt, соответствующий мосту длины t неотрицательного процесса Леви, начинающемуся в точке a и заканчивающемуся в точке b, который неформально может быть определен как
P↑,ta,b(⋅):=Pa(⋅∣Ys⩾0∀s∈[0,t],Yt=b)
(см. п. 6.1 статьи [11], где содержится детальное обоснование корректности такого определения).
В случае α=2, ρ=1/2, т.е. когда процесс Y является стандартным броуновским движением, закон распределения P↑,1a,b(⋅) задает распределение броуновской экскурсии.
Для N∈N определим нормированное отображениеϕN:ΩRWN→Ω1 соотношением
(ϕN(S))(t):=S[Nt]aN,
где {aN,N∈N} – нормирующая последовательность из условия A1.
Используя это определение, обозначим через P↑,Nx,y∘ϕ−1N закон распределения, индуцируемый на множестве Ω1:=D([0,1],R) при помощи распределения P↑,Nx,y и отображения ϕN.
Теорема 1. Пусть a,b∈[0,∞), и пусть {xN,N∈N} и {yN,N∈N} – две неотрицательные последовательности такие, что xN/aN→a и yN/aN→b (в (h;c)-решетчатом случае будем предполагать, что (yN−xN)∈Nc+hZ для всех N∈N). Если выполнены условия A1 и A2, то при N→∞
P↑,Nx,y∘ϕ−1N⟹P↑,1a,b.
Из этой теоремы вытекает, что если a=b=0, распределение X1 решетчато, а m=m(N)=o(N) при N→∞, то для любого z>0
limN→∞PxN(SN−maN⩾z|LN⩾0,SN=yN)=0.
Следовательно, SN−m/aN→0 по вероятности при N→∞, и, таким образом, теорема 1 дает мало информации о распределении случайной величины SN−m в этой ситуации.
В этом параграфе мы укажем такие центрирование и нормирование случайной величины SN−m, которые обеспечивают сходимость преобразованной случайной величины к собственному невырожденному распределению. Оказывается, что вид предельного распределения существенным образом зависит от скорости изменения параметра yN. Более того, чтобы включить в рассмотрение и абсолютно непрерывны случай, мы покажем, что
limN→∞PxN(SN−m−SNam⩽z|LN⩾0,SN⩽yN)=A(z),
где A(z) – собственное невырожденное распределение, вид которого различается в случаях yN/am→0, yN/am→T∈(0,∞), yN/am→∞.
Для доказательства заявленных результатов нам понадобятся новые обозначения, в которых инфимум по пустому множеству есть бесконечность.
Укажем, что для проверки третьего равенства необходимо использовать соотношение
{Sn−Sn−k,k=0,1,…,n}d={Sk,k=0,1,…,n}.
Для x⩾0 ведем функции восстановления
V±(x):=∞∑k=0P(H±k⩽x)=∞∑k=0∞∑n=0P(τ±k=n,±Sn⩽x).
Заметим, что функции V±(x) непрерывны справа и не убывают, причем
V±(0)=∞∑k=0P(H±k=0)=11−ζ.
Далее, если выполнено условие A1, то limn→∞P(Sn>0)=ρ∈(0,1) и, следовательно, случайные величины τ+1 и τ−1 являются собственными, причем (см., например, [21] или [23])
P(τ+1>n)∈R−ρ,V+(x)∈Rαρ,
P(τ−1>n)∈R−(1−ρ),V−(x)∈Rα(1−ρ).
Отметим, что в силу [25; соотношения (15) и (31)] и [11; соотношение (3.18)] существуют такие положительные константы ˆC, C+ и C−, что
Для x⩾0 введем непрерывные слева функции восстановления
V_±(x):=∞∑k=0P(H±k<x)=∞∑k=0∞∑n=0P(τ±k=n,±Sn<x).
Если распределение случайной величины X1 абсолютно непрерывно, то
V_±(x)=V±(x).
Мы неоднократно будем использовать этот факт, ссылаясь на результаты работ [14] и [25].
Далее нам потребуются различные функции вида T(xN, yN), зависящие от действительных параметров xN и yN.
Для фиксированной положительной последовательности an, n=1,2,…, будем для краткости писать
“T(xN,yN)=o(1) или T(xN,yN)∼c>0 равномерно по xN=o(aN), yN=o(am) при min(m,N)→∞”,
имея в виду, что
“T(xN,yN)=o(1) или T(xN,yN)∼c>0 равномерно по xN∈(0,δNaN], yN∈(0,δNam] для любой положительной последовательности δN→0 при min(m,N)→∞”.
Для упрощения последующего изложения напомним несколько локальных предельных теорем, справедливых в решетчатом случае и установленных в [11; § 4].
Пусть g(⋅) – плотность распределения случайной величины Y1 и g+(⋅) – плотность распределения извилины процесса Леви Y в момент t=1 (см. [8]):
∫y0g+(x)dx=P0(Y1⩽y|inf0⩽s⩽1Ys⩾0),y⩾0
(см. лемму 4 в [12], где дано детальное описание свойств этой плотности). Аналогично, пусть g−(⋅) – плотность распределения извилины процесса Леви −Y в момент t=1. Наконец, для a,b∈[0,∞) введем функцию
C(a,b):=Pa(inf0⩽s⩽1Ys⩾0|Y1=b).
Лемма 1. Пусть выполнены условие A1 и условие A2 ((h;c)-решетчатый случай). Тогда при n→∞ для x,y⩾0, где (y−x)∈nc+hZ, верны следующие соотношения:
1) равномерно по x=o(an) и y⩾0
g+n(x,y)=hP(τ−1>n)anV−(x)(g+(yan)+o(1));
2) равномерно по y=o(an) и x⩾0
g+n(x,y)=hP(τ+1>n)anV+(y)(g−(xan)+o(1));
3) равномерно по y=o(an) и x=o(an)
g+n(x,y)=h(1−ζ)g(0)nanV−(x)V+(y)(1+o(1));
4) для любого T>1 равномерно по x,y∈(T−1an,Tan)
g+n(x,y)=hang(y−xanr)C(xan,yan)(1+o(1)).
Отметим, что соотношение (2.12) является следствием принципа инвариантности Лиггетта для мостов (см. [20]) и локальной теоремы Гнеденко (см. [16]).
Введенные выше функции восстановления были определены в терминах характеристик слабых лестничных моментов случайных блужданий. Следует заметить, что и строгие лестничные моменты {ˆτ±k,k⩾0}, и величины {ˆH±k,k⩾0}, определяемые соотношениями ˆτ±0:=0, ˆH±0:=0 и
ˆτ±k:=inf{n>ˆτ±k−1:±Sn>±Sˆτ±k−1},ˆH±k:=±Sτ±k
при k⩾1, также часто используются при исследовании свойств случайных блужданий и ветвящихся процессов в случайной среде (см., например, [1], [3], [14], [25], [27]). Последовательности {ˆH±k,k⩾0} порождают функции восстановления
Для абсолютно непрерывного случая результаты, аналогичные лемме 1, были также получены в [11].
Лемма 2. Пусть выполнены условие A1 и условие A2 (абсолютно непрерывный случай). При n→∞ для x,y⩾0 верны следующие соотношения:
1) равномерно по x=o(an) и y⩾0
f+n(x,y)=P(τ−1>n)anV−(x)(g+(yan)+o(1));
2) равномерно по y=o(an) и x⩾0
f+n(x,y)=P(τ+1>n)anV+(y)(g−(xan)+o(1));
3) равномерно по y=o(an) и x=o(an)
f+n(x,y)=g(0)nanV−(x)V+(y)(1+o(1));
4) для любого T>1 равномерно по x,y∈(T−1an,Tan)
f+n(x,y)=1ang(y−xan)C(xan,yan)(1+o(1)).
Соотношение (2.16) не было явно упомянуто в статье [11]. Однако его несложно вывести из (2.15) и (2.2) в силу симметрии, а именно, рассматривая случайное блуждание −S вместо S, меняя местами x и y и заменяя каждый встречающийся знак “+” на знак “−”.
Соотношение (2.18) является следствием принципа инвариантности Лиггетта для мостов (см. [20]) и локальной теоремы, доказанной Стоуном в [24].
Замечание 1. Отметим, что значение константы C∗ не является универсальным в следующем смысле. Если мы изменим нормирующую последовательность, положив ¯an=can для некоторого c>0, то
Пусть, далее, U(w):=C∗∗wα(1−ρ),w⩾0. Из доказательства теоремы 1.1 в [10] (см. рассуждения между формулами (3.11) и (3.12) и определение (3.1)) следует, что
E[U(Y1)|inf0⩽s⩽1Ys⩾0]=C∗∗∫∞0wα(1−ρ)g+(w)dw=1.
Отсюда, рассматривая распределение извилины процесса −Y в момент t=1 и замечая, что параметр положительности процесса −Y равен 1−ρ, нетрудно вывести в силу соображений симметрии, что
C∗=limn→∞V+(an)P(τ+1>n)=limn→∞ˆV+(an)P(ˆτ+1>n),
и, вводя функцию ¯U(w):=C∗wαρ, получить соотношение
E[¯U(−Y1)|inf0⩽s⩽1(−Ys)⩾0]=C∗∫∞0wαρg−(w)dw=1.
Лемма доказана.
Положим
bn:=1nan=1n1+1/αℓ(n).
Из [3; предложение 2.3] и равенств (2.13) несложно вывести следующее утверждение, дополняющее леммы 1 и 2.
Лемма 4. Пусть выполнено условие A1. Тогда существует такое число C>0, что равномерно по всем x,y⩾0 и n∈N выполняется соотношение
Px(y−1⩽Sn<y,Ln⩾0)⩽CbnV−(x)V+(y),
дающее в (h;0)-решетчатом случае оценку
Px(Sn<y,Ln⩾0)⩽CbnV−(x)y∑z=0V+(z),
а в абсолютно непрерывном случае – неравенство
Px(Sn<y,Ln⩾0)⩽CbnV−(x)∫y0V+(z)dz.
Доказательство. В [3; утверждение 2.3] было показано, что если выполнено условие A1, то существует число C>0 такое, что равномерно по всем x,y⩾0 и всем n
Px(y−1⩽Sn<y,Ln⩾0)⩽CbnˆV−(x)ˆV_+(y).
Теперь соотношение (2.20) следует из (2.13) и неравенств
ˆV−(x)⩽V−(x),ˆV_+(y)⩽V+(y).
Второе и третье утверждения леммы выводятся с помощь суммирования и интегрирования соответственно.
Лемма доказана.
§ 3. Предельные теоремы для решетчатого случая
Следуя [11], в дальнейшем мы будем рассматривать лишь (1;0)-решетчатый случай условия A2, т.е. будем предполагать, что распределение величины X1 сосредоточено на множестве Z и является апериодическим. Общий (h;c)-решетчатый случай требует лишь более громоздких обозначений, в то время как проводимые при этом доказательства не нуждаются в дополнительных соображениях.
Нашей целью является исследование асимптотического поведения вероятностей
PxN(SN−m⩽zam,LN⩾0,SN=yN)
при условии, что max(xN,yN)/aN→0 и m=o(N) при N→∞.
Рассмотрим отдельно три случая: yN/am→0, yN/am→T∈(0,∞) и yN/am→∞.
3.1. Случай yN/am→0
Лемма 5. Пусть выполнены условия A1 и A2 и распределение случайной величины X1 является (1;0)-решетчатым. Если min(m,N)→∞ так, что m=o(N), то для любого z∈(0,∞)
PxN(SN−m⩽zam∣LN⩾0,SN=yN)=A1(z)(1+o(1))
равномерно по xN=o(aN) и yN=o(am), где вследствие (2.6) и леммы 3 функция
Следствие 1. Пусть выполнены условия A1 и A2 и распределение случайной величины X1 является (1;0)-решетчатым. Если min(m,N)→∞ так, что m=o(N), то для любого z∈(0,∞)
PxN(SN−m⩽zam,LN⩾0,SN⩽yN)∼A1(z)PxN(LN⩾0,SN⩽yN)
равномерно по xN=o(aN) и yN=o(am).
Используя суммирование, из (3.9), (3.2) и (3.3) получаем требуемое утверждение.
Лемма 6. Пусть выполнены условия A1 и A2, а распределение случайной величины X1 является (1;0)-решетчатым. Если min(m,N)→∞ так, что m=o(N), то для любого z∈(0,∞)
PxN(SN−m⩽zam∣LN⩾0,SN=j)∼A2(z,jam)
равномерно по xN=o(aN) и j∈[t0am,t1am], где
A2(z,t):=t−αρ∫z0wαρg(t−w)C(w,t)dw.
Доказательство. Воспользуемся снова разложением (3.5). Оценки (3.6) и (3.7) остаются неизменными. Далее, используя (2.12) вместо (3.8), получаем, что при m→∞
Pk(Sm=j,Lm⩾0)∼1amg(j−kam)C(kam,jamr)
равномерно по εam⩽k⩽zam и j∈[t0am,t1am]. Отсюда, учитывая (3.7) и (3.10), заключаем, что
В отличие от п. 3.1 и п. 3.2, где изучалось распределение случайной величины S_{N-m}, здесь мы исследуем распределение разности S_{N-m}-S_{N}.
Лемма 7. Пусть выполнены условия A1 и A2, а распределение случайной величины X_{1} является (1;0)-решетчатым. Если \min (m,N)\to \infty и a_{m}=o(y_{N}), то для любого z\in (-\infty,\infty)
Так как y_{N}/a_{m}\to \infty , то согласно принципу инвариантности для случайных блужданий, распределения приращений которых принадлежат без центрирования области притяжения устойчивого закона,
равномерно по x_{N}=o(a_{N}) и y_{N}=o(a_{N}). Ввиду (2.1) это влечет утверждение леммы.
Следствие 3. Пусть выполнены условия A1 и A2, а распределение случайной величины X_{1} является (1;0)-решетчатым. Если m\to \infty и a_{m}=o(y_{N}), то для любого z\in (-\infty,\infty)
Доказательство. Доказательства утверждений 1)–3) идейно повторяют аргументы соответствующих утверждений следствий 1–3. Поэтому мы проверим лишь справедливость соотношения (4.1).
§ 5. Предельная теорема для последовательностей, сходящихся почти наверное
В этом параграфе нам понадобится условная мера \mathbb{P}_x^{\uparrow}(\cdot), которая порождает случайное блуждание, начинающееся в точке x и рассматриваемое при условии неотрицательности этого блуждания на всей временно́й оси (см., например, [4], [10] и [25]). Эта новая мера задается следующим образом: для x\geqslant 0, всех N\in \mathbb{N} и любого измеримого множества B, принадлежащего \sigma -алгебре, порожденной случайными величинами S_{1},\dots,S_{N}, положим
Теорема 2. Пусть выполнены гипотезы A1 и A2, а H_{1},H_{2},\dots – равномерно ограниченная последовательность случайных величин, согласованная с фильтрацией \widetilde{\mathcal{F}}\,{=}\,\{\widetilde{\mathcal{F}}_{k},\,k\,{\in}\,\mathbb{N}\} и сходящаяся \mathbb{P}_{0}^{\uparrow}-п.н. к случайной величине H_{\infty} при n\to \infty . Предположим, что параметр m=m(n) стремится к бесконечности при n\to \infty так, что m=o(n). Тогда:
где \mathcal{S}'=\{S_{n}',\,n=0,1,2,\dots\} – вероятностная копия случайного блуждания \mathcal{S}, не зависящая от множества \{S_{j},\,j=0,1,\dots,k\} . Согласно следствиям 1, 2 и соотношению (2.21) существуют константы C, C_{1}, C_{2} и C_{3} такие, что для любого фиксированного k и всех n\geqslant k и z>0
Далее, опираясь на следствие 1 и соотношение a_{\lambda m}\sim \lambda^{1/\alpha}a_{m}, справедливое при m\to \infty , заключаем, что при m,n\to \infty и m=o(n)
Устремляя сначала n, а потом и k к бесконечности и применяя теорему о мажорируемой сходимости, мы видим, что для каждого \lambda >1 правая часть предыдущего соотношения стремится к нулю.
Опираясь на этот результат, мы приходим к цепочке равенств
Объединение полученных оценок, доказывает утверждение (5.1) для (1;0)-решетчатого случая.
Для доказательства соотношения (5.2) для (1;0)-решетчатого случая нужно повторить почти дословно аргументы, приведенные при проверке справедливости (5.1), заменяя ссылку на следствие 1 ссылкой на следствие 2.
Для проверки справедливости соотношения (5.3) для (1;0)-решетчатого случая нужно повторить почти дословно аргументы, приведенные при доказательстве (5.1), заменяя ссылку на следствие 1 ссылкой на следствие 3.
Для доказательства утверждения (5.1) в абсолютно непрерывном случае необходимо заменить в приведенных выше рассуждения выражение \sum_{j=0}^{\varphi (n)}V^{+}(j) на \displaystyle\int_{0}^{\varphi (n)}V^{+}(w)\,dw и использовать лемму 8. Аналогичные аргументы позволяют проверить и соотношения (5.2) и (5.3) для случая абсолютно непрерывного распределения случайной величины X_{1}.
В этом параграфе мы применим результаты, полученные для случайных блужданий, к исследованию размера популяции критических ветвящихся процессов, эволюционирующих в неблагоприятной случайной среде. Для описания рассматриваемой модели и детальной постановки задач, которые мы намерены решать, введем пространство \mathcal{M}=\{\mathfrak{f}\} всех вероятностных мер на множестве \mathbb{N}_{0}. Для упрощения обозначений будем отождествлять меру \mathfrak{f}=\{\mathfrak{f}(\{0\}),\mathfrak{f}(\{ 1\}),\dots\} \in \mathcal{M} с соответствующей вероятностной производящей функцией
и не будем делать различия между \mathfrak{f} и f. Оснащенное метрикой, индуцируемой расстоянием по вариации между распределениями, пространство \mathcal{M}=\{\mathfrak{f}\} =\{f\} превращается в польское пространство. Пусть
– последовательность независимых вероятностных копий случайной величины F. Бесконечная последовательность \mathcal{E}=\{F_{n},\,n\in \mathbb{N}\} называется случайной средой.
Последовательность неотрицательных целочисленных случайных величин \mathcal{Z}=\{Z_{n},\,n\in \mathbb{N}_{0}\} , заданных на некотором вероятностном пространстве (\Omega,\mathcal{F},\mathbb{P}), называется ветвящимся процессом в случайной среде (ВПСС), если Z_{0} не зависит от \mathcal{E} и при фиксации среды \mathcal{E} процесс \mathcal{Z} является марковской цепью
для всех n\in \mathbb{N}, z_{n-1}\in \mathbb{N}_{0} и f_{1},f_{2},\ldots\in \mathcal{M}, где \xi_{n1},\xi_{n2},\dots – последовательность независимых одинаково распределенных случайных величин с распределением f_{n}. Таким образом, Z_{n-1} – это размер (n-1)-го поколения ветвящегося процесса, а f_{n} – распределение числа потомков каждой частицы (n-1) -го поколения.
Известно (см. [1; теорема 1.1 и следствие 1.2]), что если выполнены условия B1, B2, то существуют такие константа \theta \in (0,\infty) и последовательность l(1),l(2),\dots, медленно меняющаяся на бесконечности, что при n\to \infty
где \mathcal{Y}^{+}=\{Y_{t}^{+},\,0\leqslant t\leqslant1\} обозначает извилину строго устойчивого процесса \mathcal{Y} индекса \alpha .
Таким образом, если ВПСС является критическим, то при условии, что Z_{n}>0, случайные величины \log Z_{[nt]}, t\in (0,1], и S_{n} (значение сопровождающего случайного блуждания, обеспечивающего выживание популяции к моменту времени n) растут как a_{n}, умноженное на случайные сомножители.
Эти результаты были дополнены в [28] анализом распределений соответствующим образом нормированных случайных величин \log Z_{[nt]}, t\in (0,1], рассматриваемых при условии, что Z_{n}>0 и S_{n}\leqslant \varphi (n), где \varphi (n)\to \infty при n\to \infty так, что \varphi (n)=o(a_{n}). В этом случае ввиду (6.1) событие \{S_{n}\leqslant \varphi(n)\} можно трактовать как неблагоприятное для развития критического ветвящегося процесса в случайной среде, рассматриваемого при условии его невырождения к далекому моменту времени.
Введем обозначение
\begin{equation*}
A_{\mathrm{u.s}}:=\{Z_{n}>0\text{ при всех }n>0\}
\end{equation*}
\notag
для события, обозначающего невырождение процесса за все время его существования, и положим
Заметим, что в силу [27; теорема 1], если выполнены условия B1 и B2, а распределение приращений сопровождающего случайного блуждания абсолютно непрерывно и \varphi(n)=o(a_{n}), то при n\to \infty
Аналогичная асимптотика имеет место и для (1;0)-решетчатого случая при замене \displaystyle\int_{0}^{\varphi (n)}V^{+}(w)\,dw на \sum_{j=0}^{\varphi (n)}V^{+}(j).
Наряду с асимптотическим поведением вероятности невырождения критического ВПСС в следующей теореме описана скорость роста размера популяции в логарифмической шкале.
Теорема 3 (см. [28]). Пусть выполнены условия B1 и B2. Если \varphi (n)\to \infty при n\to \infty так, что \varphi (n)=o(a_{n}), то для любого y\in (0,1]
где \mathcal{Y}^{++}=\{Y_{t}^{++},\,0\leqslant t\leqslant1\} – экскурсия строго устойчивого процесса \mathcal{Y} индекса \alpha .
Имеется существенная разница между порядками нормировок в соотношениях (6.4) и (6.5), показывающая, что должен существовать фазовый переход скорости роста размера популяции в случае, когда мы рассматриваем процесс внутри интервала [n-m,n), m=o(n), m\to \infty . Такой переход действительно имеет место, и целью этого параграфа является доказательство того, что при нормировке случайных величин \log Z_{n-m} величинами a_{m} возникают три различных предельных распределения. Вид этих распределений зависит от того, какое из следующих трех условий выполнено при \min (n-m,m)\to \infty : функция \phi(n) имеет порядок o(m); функция \phi (n) пропорциональна m; m=o(\phi (n)).
Наш основной результат выглядит следующим образом. Напомним, что в решетчатом случае мы условились для простоты изложения ограничиться анализом (1;0)-решетчатого случая.
Теорема 4. Пусть выполнены условия B1 и B2, \min (m,n)\to \infty и m=o(n). Тогда:
1) если \varphi (n)\to \infty при n\to \infty так, что \varphi (n)=o(a_{m}), то для любого z\in (0,\infty)
В [28; теорема 1] было показано, что если r_{1},r_{2},\dots – такая последовательность натуральных чисел, что r_{n}\leqslant n, r_{n}\to \infty и \varphi (n)\to \infty при n\to \infty так, что \varphi (n)=o(a_{n}), то при n\to \infty
для любой точки непрерывности x\in (0,\infty) распределения случайной величины W. Используя (6.7), получаем, что для любого \varepsilon >0 существует M=M(m,n) такое, что
и исследуем отдельно асимптотическое поведение слагаемых в правой части этого равенства при \min(m,n)\to \infty , m=o(n). Ввиду (6.9) нам достаточно проанализировать величину
Начиная с этого момента и вплоть до конца доказательства мы предполагаем, что распределение случайной величины X_{1} абсолютно непрерывно. Для доказательства требуемых утверждений в случае решетчатого распределения случайной величины X_{1} необходимо всюду ниже заменить знак \displaystyle\int на знак \sum .
\mathbb{P}^{\uparrow}-п.н. при n-j\to \infty . Более того, H_{\infty}(k) >0\mathbb{P}^{\uparrow}-п.н. согласно [1; утверждение 3.1]. Далее, для q\in (-\sqrt{\varphi (n)},0] имеем
Теперь мы отдельно рассмотрим случаи \varphi (n)=o(a_{m}), \varphi (n)\sim Ta_{m} и a_{m}=o(\varphi (n)).
1) Допустим, что \varphi (n)=o(a_{m}). В этом случае, полагая H_{n}:=\mathbb{P}(Z_{n-j}>0\mid \mathcal{E}, Z_{0}=k) в (5.1) и вспоминая (6.12), получаем, что для каждого q\in [-\sqrt{\varphi (n)},0]
3) Рассмотрим, наконец, случай, когда m=m(n)\to \infty и a_{m}=o(\varphi (n)) при n\to \infty . Введем обозначение S_{n-m,n}:=S_{n-m}-S_{n} и для z\in (-\infty,\infty) запишем равенство
Несложно проверить, что если во всех соотношениях, находящихся между формулами (6.10) и (6.11), мы заменим S_{n-m} на S_{n-m,n} и S_{n-m-j} на S_{n-m-j,n-j}, то все оценки, стоящие между (6.10) и (6.11), останутся справедливыми. Таким образом,
Авторы выражают благодарность рецензенту, замечания которого позволили улучшить представление изложенных в статье результатов.
Список литературы
1.
V. I. Afanasyev, J. Geiger, G. Kersting, V. A. Vatutin, “Criticality for branching processes in random environment”, Ann. Probab., 33:2 (2005), 645–673
2.
В. И. Афанасьев, “Принцип инвариантности для критического процесса Гальтона–Ватсона, достигающего высокого уровня”, Теория вероятн. и ее примен., 55:4 (2010), 625–643; англ. пер.: V. I. Afanasyev, “Invariance principle for a critical Galton–Watson process attaining a high level”, Theory Probab. Appl., 55:4 (2011), 559–574
3.
V. I. Afanasyev, C. Böinghoff, G. Kersting, V. A. Vatutin, “Limit theorems for weakly subcritical branching processes in random environment”, J. Theoret. Probab., 25:3 (2012), 703–732
4.
J. Bertoin, R. A. Doney, “On conditioning a random walk to stay nonnegative”, Ann. Probab., 22:4 (1994), 2152–2167
5.
N. H. Bingham, C. M. Goldie, J. L. Teugels, Regular variation, Encyclopedia Math. Appl., 27, Cambridge Univ. Press, Cambridge, 1987, xx+491 pp.
6.
E. Bolthausen, “On a functional central limit theorem for random walks conditioned to stay positive”, Ann. Probab., 4:3 (1976), 480–485
7.
A. Bryn-Jones, R. A. Doney, “A functional limit theorem for random walk conditioned to stay non-negative”, J. London Math. Soc. (2), 74:1 (2006), 244–258
8.
L. Chaumont, “Excursion normalisée, méandre at pont pour les processus de Lévy stables”, Bull. Sci. Math., 121:5 (1997), 377–403
9.
F. Caravenna, “A local limit theorem for random walks conditioned to stay positive”, Probab. Theory Related Fields, 133:4 (2005), 508–530
10.
F. Caravenna, L. Chaumont, “Invariance principles for random walks conditioned to stay positive”, Ann. Inst. Henri Poincaré Probab. Stat., 44:1 (2008), 170–190
11.
F. Caravenna, L. Chaumont, “An invariance principle for random walk bridges conditioned to stay positive”, Electron. J. Probab., 18 (2013), 60, 32 pp.
12.
L. Chaumont, R. A. Doney, “Invariance principles for local times at the maximum of random walks and Lévy processes”, Ann. Probab., 38:4 (2010), 1368–1389
13.
F. den Hollander, Random polymers, École d'Été de Probabilités de Saint-Flour XXXVII – 2007, Lecture Notes in Math., 1974, Springer-Verlag, Berlin, 2009, xiv+258 pp.
14.
R. A. Doney, “Local behaviour of first passage probabilities”, Probab. Theory Related Fields, 152:3-4 (2012), 559–588
15.
В. Феллер, Введение в теорию вероятностей и ее приложения, т. 2, Мир, М., 1967, 752 с. ; пер. с англ.: W. Feller, An introduction to probability theory and its applications, т. 2, John Wiley & Sons, Inc., New York–London–Sydney, 1966, xviii+626 с.
16.
Б. В. Гнеденко, А. Н. Колмогоров, Предельные распределения для сумм независимых случайных величин, Гостехиздат, М.–Л., 1949, 264 с. ; англ. пер.: B. V. Gnedenko, A. N. Kolmogorov, Limit distributions for sums of independent random variables, Addison-Wesley Publishing Company, Inc., Cambridge, MA, 1954, ix+264 с.
17.
D. L. Iglehart, “Functional central limit theorems for random walks conditioned to stay positive”, Ann. Probab., 2:2 (1974), 608–619
18.
W. D. Kaigh, “An invariance principle for random walk conditioned by a late return to zero”, Ann. Probab., 4:1 (1976), 115–121
19.
G. Kersting, V. Vatutin, Discrete time branching processes in random environment, Math. Stat. Ser., John Wiley & Sons, London; ISTE, Hoboken, NJ, 2017, xiv+286 pp.
20.
T. M. Liggett, “An invariance principle for conditioned sums of independent random variables”, J. Math. Mech., 18:6 (1968), 559–570
21.
Б. А. Рогозин, “Распределение первого лестничного момента и высоты и флуктуации случайного блуждания”, Теория вероятн. и ее примен., 16:4 (1971), 593—613; англ. пер.: B. A. Rogozin, “The distrbution of the first ladder moment and height and fluctuation of a random walk”, Theory Probab. Appl., 16:4 (1971), 575–595
22.
Е. Сенета, Правильно меняющиеся функции, Наука, М., 1985, 142 с. ; пер. с англ.: E. Seneta, Regularly varying functions, Lecture Notes in Math., 508, Springer-Verlag, Berlin–New York, 1976, v+112 с.
23.
Я. Г. Синай, “О распределении первой положительной суммы для последовательности независимых случайных величин”, Теория вероятн. и ее примен., 2:1 (1957), 126–135; англ. пер.: Ya. G. Sinai, “On the distribution of the first positive sum for a sequence of independent random variables”, Theory Probab. Appl., 2:1 (1957), 122–129
24.
C. Stone, “A local limit theorem for nonlattice multi-dimensional distribution functions”, Ann. Math. Statist., 36:2 (1965), 546–551
25.
V. A. Vatutin, V. Wachtel, “Local probabilities for random walks conditioned to stay positive”, Probab. Theory Related Fields, 143:1-2 (2009), 177–217
26.
V. Vatutin, E. Dyakonova, “Path to survival for the critical branching processes in a random environment”, J. Appl. Probab., 54:2 (2017), 588–602
27.
В. А. Ватутин, Е. Е. Дьяконова, “Критические ветвящиеся процессы, эволюционирующие в неблагоприятной случайной среде”, Дискрет. матем., 34:3 (2022), 20–33; англ. пер.: V. A. Vatutin, E. E. Dyakonova, Critical branching processes evolving in an unfavorable random environment, 2022, 15 с., arXiv: 2209.13611
28.
В. А. Ватутин, Е. Е. Дьяконова, “Размер популяции критического ветвящегося процесса, эволюционирующего в неблагоприятной среде”, Теория вероятн. и ее примен., 68:3 (2023), 509–531; англ. пер.: V. A. Vatutin, E. E. Dyakonova, “Population size of a critical branching process evolving in an unfavorable environment”, Theory Probab. Appl., 68:3 (2023), 411–430
Образец цитирования:
В. А. Ватутин, К. Донг, Е. Е. Дьяконова, “Случайные блуждания, остающиеся неотрицательными, и ветвящиеся процессы в неблагоприятной среде”, Матем. сб., 214:11 (2023), 3–36; V. A. Vatutin, C. Dong, E. E. Dyakonova, “Random walks conditioned to stay nonnegative and branching processes in an unfavourable environment”, Sb. Math., 214:11 (2023), 1501–1533
Vladimir A. Vatutin, Elena E. Dyakonova, “Branching processes under nonstandard conditions”, Stoch. Qual. Control, 39:1 (2024), 1–1
В. А. Ватутин, Е. Е. Дьяконова, “О минимуме случайного блуждания, сосредоточенного на неотрицательной полуоси”, Дискрет. матем., 36:3 (2024), 50–79; V. A. Vatutin, E. E. Dyakonova, “On the prospective minimum of the random walk conditioned to stay nonnegative”, Discrete Math. Appl., 34:6 (2024), 337–362
В. А. Ватутин, К. Донг, Е. Е. Дьяконова, “Некоторые функционалы для случайных блужданий и критические ветвящиеся процессы в экстремально неблагоприятной среде”, Матем. сб., 215:10 (2024), 58–88; V. A. Vatutin, C. Dong, E. E. Dyakonova, “Some functionals for random walks and critical branching processes in an extremely unfavourable random environment”, Sb. Math., 215:10 (2024), 1321–1350