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

36

ВЕСТН. моек. УН-ТА. СЕР.1, МАТЕМАТИКА. МЕХАНИКА. 2023. №6

11. Хаусдорф Ф. Теория множеств. М.: ОНТИ. Гл. ред. техн.-теор. лит-ры, 1937.

12. Ветохин А.Н. Класс Бэра полунепрерывных снизу минорант показателей Ляпунова // Дифференц. уравнения. 2019. 34, № 10. 1313-1317.

Поступила в редакцию 24.05.2023

УДК 519.21

ОБ АСИМПТОТИЧЕСКОМ ПОВЕДЕНИИ ТОЧЕЧНЫХ ПРОЦЕССОВ ВЫХОДОВ ЗА ВЫСОКИЕ УРОВНИ ГАУССОВСКОЙ СТАЦИОНАРНОЙ ПОСЛЕДОВАТЕЛЬНОСТИ

В. И. Питербарг1

Изучается асимптотическое поведение точечных процессов выходов гауссовской стационарной последовательности за уровень, стремящийся к бесконечности медленнее, чем в пуассоновской предельной теореме для числа выходов. Доказана сходимость по вариации таких точечных процессов к маркированному пуассоновскому процессу. Применяются результаты Ю. В. Прохорова о наилучшем приближении распределения Бернулли смесью гауссовского и пуассоновского распределений. Эта задача поставлена А. Н. Колмогоровым в начале 50-х годов прошлого века.

Ключевые слова: гауссовская последовательность, большие выбросы, пуассоновская предельная теорема, сходимость по вариации.

We study the asymptotic behavior of point processes of exits of a Gaussian stationary sequence beyond a level tending to infinity more slowly than in the Poisson limit theorem for the number of exits. Convergence in variation of such point processes to a marked Poisson process is proved. The results of Yu. V. Prokhorov on the best approximation of the Bernoulli distribution by a mixture of Gaussian and Poisson distributions are applied. A. N. Kolmogorov proposed this problem in the early 1950s.

Key words: Gaussian sequence, large excursions, Poisson limit theorem, convergence in variation.

DOI: 10.55959/MSU0579-9368-1-64-6-5

1. Введение. Рассмотрим гауссовскую стационарную последовательность X(k), k € Z, с нулевым средним EX(k) =0 и единичной дисперсией VarX(k) = 1. Обозначим Cov(X(0),X(k)) = EX(k)X(0) = r(k). Пусть B — ст-алгебра ограниченных борелевских множеств на прямой R. Введем на B точечный процесс

Vu(B) := £ I {X(k) >u} , B € B, (1)

kenB

где последовательность натуральных чисел n = n(u) стремится к бесконечности с ростом уровня и. Мы рассматриваем предельное поведение процесса Пи(■) при различных скоростях роста этого уровня к бесконечности. Известно, что при достаточно быстром убывании корреляции r(k) к нулю на бесконечности процесс пи(■) слабо сходится в естественной нормировке к пуассоновскому. Естественная нормировка — это

П = П{и) = Р{Х{1)>УУ Л>0'

1 Питербарг Владимир Ильич — доктор физ.-мат. наук, гл. науч. сотр. каф. теории вероятностей мех.-мат. ф-та МГУ; e-mail: piter@mech.math.msu.su.

Piterbarg Vladimir Ilich — Doctor of Physical and Mathematical Sciences, Chief Researcher, Lomonosov Moscow State University, Faculty of Mechanics and Mathematics, Chair of Probability Theory.

ВЕСТН. МОСК. УН-ТА. СЕР.1, МАТЕМАТИКА. МЕХАНИКА. 2023. №6

37

т.е. среднее число точек в множестве nB асимптотически постоянно. В случае независимых X(k) это, по сути, классическая теорема Пуассона для бернуллиевской последовательности Пи([0,1]). Поскольку для гауссовской стандартной функции распределения Ф выполнено

1 2

Щи) ^ р(и) := 1 - Ф(и) ^ (1 - гГ2)Ф(и), и> 0, где Ф(и) = /2, (2)

V2nu

удобнее использовать эквивалентную нормировку n(u) = А/Ф(и), что в свою очередь эквивалентно соотношению

г-- \ log log п + 1с^(Ад/7г72)

и = ип = V21ogn - ----b 0(l/logn), п ->• оо (3)

V2 log n

(см., например, [1]). Нетрудно проверить, что это соотношение остается в силе и при np(un) — ж вместе с n — ж и un — ж, т.е. при А — ж. В работе [2], по сути, утверждается следующее.

Теорема 1. Пусть для ковариационной функции r(k) последовательности X (k) выполнено соотношение

r(k) log k — 0 при k — ж.

Если для последовательности точечных процессов пи(В), B € B, n = 1, 2,..., последовательность уровней u = un удовлетворяет соотношению (3), то при u — ж имеет место слабая сходимость

Пи(В) ^ ПА(В), (4)

где Па (В), В € B, — пуассоновский точечный процесс с интенсивностью А.

На самом деле Й. Миттал и Д. Илвисакер доказали предельную теорему для максимума гауссовской стационарной последовательности, однако теорема 1 очевидно отсюда следует в силу теоремы Калленберга о сходимости точечных процессов (детали см. в [1] или [3]). В этих монографиях также приводятся необходимые и достаточные условия на r(k) для выполнения соотношения (4). В настоящей работе рассматривается задача об описании предельного поведения точечного процесса Пи (В), когда уровень u по-прежнему стремится к бесконечности вместе с n, но медленнее, чем в соотношении (4). Основываясь на теоремах Ю.В. Прохорова [4] о приближении распределения Бернулли в зависимости от поведения его параметров, мы рассмотрим сначала предельное поведение последовательности nu([0,1]) при различных соотношениях между n и u = u(n) при n,u — ж. В случае независимых X(k) нами дана фактически переформулировка соответствующей теоремы Ю. В. Прохорова в терминах сходимости маркированных точечных процессов. Это возможно, поскольку теоремы Ю.В. Прохорова дают оценки скорости сходимости в метрике Ii (определение сходимости дано ниже соотношением (5)) к пуассоновскому, нормальному или вырожденному распределению для любых соотношений между u и n при их стремлении к бесконечности. Как и в [4], в модели (1) мы рассматриваем различные соотношения стремления u и n к бесконечности, при этом, если u растет медленнее стандартной (3), его интенсивность стремится к бесконечности. В связи с этим при помощи операции прореживания (см. ниже) мы представляем процесс nu(•) и предельный пуассо-новский в виде маркированных точечных процессов, в которых марки (кластеры) также являются точечными процессами. Основной результат работы — теорема 4, доказательство которой основано на теореме Прохорова с последующим применением техники сравнения гауссовских распределений и которая представляет собой теорему о сходимости к маркированному точечному процессу. Марки предельного процесса являются независимыми пуассоновскими процессами. В заключительных замечаниях к работе рассматриваются дальнейшие возможные шаги в этом направлении и возможные применения.

2. Последовательность с некоррелированными значениями [4]. В 1953 г. Ю.В. Прохоров, развивая результаты П. А. Козуляева 1939 г. [5], опубликовал статью о соотношении качества аппроксимации биномиального распределения В(n,p) пуассоновским и нормальным для различных поведений вероятности успеха при росте числа испытаний [4]. Как указано в этой статье, задача о выборе оптимального предельного распределения, т.е. пуассоновского или нормального, была предложена А. Н. Колмогоровым в руководимом им семинаре. В то время как в [5] доказано, что

пределом Б(п,р) при пр ^ то может быть только пуассоновское, нормальное или же вырожденное распределение, Ю.В. Прохоров установил, что если пр ^ то, то и пуассоновское, и нормальное распределения аппроксимируют Б(п,р), но с разной степенью точности.

В этом пункте теорема Ю.В. Прохорова о сходимости к пуассоновскому распределению практически повторяется в терминах изучаемой модели. Рассматриваются большие значения уровня и, т.е. вероятность р(и) мала. В [4] рассмотрен также случай и < 0, т.е. аппроксимация числа событий {X(к) > и}, и < 0, к = 0,..., п. В силу симметрии гауссовского распределения этот случай очевиден и здесь не рассматривается. Для р = р(и) (см. (2)) обозначим

В (к) = ркдп~к, д = 1-р, Щк) = к = 0,..., п,

к Г ' к!

и введем метрику 11 соотношением

п

Р!(и,п) = ^] |Б(к) - П(к)|. (5)

к=0

Теорема 2 [4, теорема 2]. Пусть г (к) = 0 для всех к > 0. Тогда

Р1(и, п) = А^(и) + р(и)0 ^шт(1, (пр(и))-1/2^ , пр ^ то,

где

2

Л1 = -== = 0.483 .... у2пе

Из теоремы 2 следует, что для любого множества Б € В и любой нормировки п распределение величины п«(Б) может быть аппроксимировано пуассоновским при п, и ^ то, даже если пр(и) ^ то.

Пример. Степенная шкала. Пусть для некоторых неотрицательного а и положительного с выполнено соотношение р(и)па ^ с. Тогда при а = 1 имеет место пуассоновская предельная теорема даже для зависимых X(к) — это теорема 1. Беря в соотношении (3) па вместо п, а > 0, получаем соотношение

■ log log п + \og(c^an/2)

/ \ /-л- 9 о ub 'v 1 1V76 V V / 7 х-л/ /л \

и = u(a) = v 2a logn — --.--h 0(1/ logn), n —> то.

л/2а log n

Ниже увидим, что теперь c — не обязательно интенсивность.

Заметим, что из других теорем статьи Ю.В. Прохорова [4] следует, что пуассоновская аппроксимация предпочтительней для уровней с а > 1/3, в то время как для а < 1/3 — нормальная. При этом обе аппроксимации возможны. Далее будут описаны предельные точечные процессы, соответствующие уровням и(а) и основанные на пуассоновской аппроксимации. Заметим, что, если а € (0,1), интенсивность предельного процесса равна бесконечности.

2.1. Прореживание и кластеры. Рассмотрим пуассоновскую аппроксимацию точечного процесса Пи(B), B € B, при np — то, при этом n — то и p — 0. Продолжаем считать, что r(k) = 0 для всех k > 0. Очевидным, но важным для нас следствием теоремы 2 является равномерная оценка этой аппроксимации.

Следствие 1. Пусть r(k) = 0 для всех k > 0. Существует константа C, такая, что для произвольного множества B € B имеет место оценка

Pi(u,n,B) := Y^ |B(k) - n(k)| < Cp(u).

kenB

Действительно, величина O (min(1, (np(u))-1/2) в теореме 2 не превосходит константы. Поскольку интенсивность точечного процесса пи может стремиться к бесконечности, представим его и соответствующий пуассоновский точечный процесс в удобном виде. Выделим из пи точечный процесс п1и(B), B € B, l-точек:

k € Пи(■) ^ Пи(') : max X(k - i) ^ u, X(k) > u.

То есть перед выходом за уровень и в течение времени I выходов не было. В непрерывном времени для процессов с недифференцируемыми траекториями — это естественный путь выделения пересечений [1, 3].

Положим I зависимым от п, I = 1(п), так, чтобы для некоторого Л € (0, то) выполнялось соотношение

n(1 - p)l(n)p = Л + o(1), n ^ то,р ^ 0. (6)

В силу пуассоновской предельной теоремы для любого непустого множества В € B распределение числа ¿(^-точек в пВ, где

п = п{и) =

стремится к пуассоновскому с параметром А|В| (|В| — мера (длина) В). Логарифмируя (6) и используя разложение Тейлора при p — 0, получаем

1{п) = l°gnP ~ lQg(A + °(1)) = р-1 log(np/A)(l + о(1)), „->оо,р->0. (7)

p + o(p)

2.2. Сходимость к маркированному пуассоновскому процессу. Рассмотрим предельное поведение введенного выше l(n)-прореженного точечного процесса

nU(n) (В) = У I ( max X(k - i) < u,X(k) > u I, В € B. «.t^n li=1.-.i(n) J

v ' ^ I i=1,...,l(n) k&nß K ' ' w

Аналогично доказательству теоремы 1 нетрудно показать, что он слабо сходится при u ^ то к пуассоновскому точечному процессу с интенсивностью А. Однако из следствия 1 теоремы 2 следует

сходимость по вариации [6]. Имея в виду (7), обозначим ni := ¿(n)/log np, так что pni ^ 1 при

7 l(n) f \

n ^ то, и сопоставим каждой точке k точечного процесса nU (■) точечный процесс

ZU(в) = у I{X(k + j) >u}, в es.

j €niB

Образно говоря, мы рассматриваем под лупой точки процесса Пи(B), расположенные между соседними ¿(п)-точками, скажем l(n)-пачки.

Таким образом, на алгебре B х Bœ введен маркированный точечный процесс

-П(£)= у Ы max X (k - i) < u,X(k) > u} у I {X(k + j) > u} | , (8)

fcenB [ ^ (n) J jeniBfc J

где обозначено В := (В, Вк) еВх

(Определение маркированного точечного процесса см. в [7, гл. 3].) Поскольку величины X(к) независимы, следующее утверждение является, по сути, переформулировкой теоремы 2: изменится только нумерация счетного числа событий в сумме (5). Заметим при этом, что алгебра В х является представлением алгебры В, адаптированным к операции кластеризации.

Теорема 3. Маркированный точечный процесс 2П(В), В € В х , сходится в метрике ^ (см. (5)) к маркированному пуассоновскому процессу

ПА(В)= ^ пк(Вк), В еВхВ». (9)

к€пА(Б)

Пуассоновский процесс пд(-) на В с интенсивностью Л и пуассоновские процессы ("марки") пк(■) на Вк, к = 1, 2,..., с интенсивностями 1, к = 1, 2,... , независимы в совокупности. Рассмотрим степенную шкалу. Из соотношения (7) следует, что

Кп) (а-Ц 1о§(ст/(Л + °(1))) - 2-У 1ой£2 Л + 0^ 1оё(1-КО) " с А \ \lognyy

В этой шкале l(n) = 0 при a = 1. Кроме того, интенсивность Л растет до бесконечности, если I убывает к нулю, что вполне естественно.

3. Общая гауссовская последовательность. Для перехода к общему случаю последовательности X(k) с ненулевой корреляцией применим технику сравнения гауссовских распределений [1, 3, 8]. Для последовательности действительных чисел Xk, k € Z, обозначим через Au алгебру множеств, порожденную множествами {xk > u}, k € Z. Обозначим для краткости X = {X(k),k = 1,... , n} и Xo = {Xo(k),k = 1,..., n}, где Xo(k) — последовательность гауссовских стандартных независимых величин. Основным инструментом перехода к общему случаю является следующая оценка близости по вариации вероятностных мер на Au, соответствующих последовательностям X и Xo [1, следствие 2.3.2].

Утверждение 1. Пусть |r(k)| < 1, k = 1,... ,n. Для любого и имеет место неравенство

sup |Р(Х G А) - Р(Хо € А)\ < - У (n~fc)|r(fc)l ехр (. (10)

Заметим, что для каждого элементарного события ш € Q множества {nu(B) € K}, nB С (0, n + 1), K С Zo, содержатся в Au. Поэтому стремление к нулю правой части неравенства (10) означает сближение по вариации точечных процессов построенных по последовательности X(k) и по

рассмотренной в предыдущем пункте последовательности с некоррелированными значениями.

Отметим также, что сходимость в норме Ii из этого неравенства, вообще говоря, не следует. Проблема такой сходимости требует отдельного рассмотрения.

Выясним, при каких соотношениях u и n правая часть (10) стремится к нулю, но при этом np — то. Так как случай np — 0 здесь не рассматривается, то, если для некоторого е > 0 выполнено неравенство infn,p np ^ е, из соотношения (3) следует, что

u

lim sup . ^ 1.

n^œ V2l°g n

Обозначим p(k) := sup^k |r(1)|.

Лемма. Пусть выполнены соотношения

r(k)k1-p(1) — 0, k — то, (11)

и _

liminf > v^iî). (12)

n^œ ^/2 log n

Тогда правая часть в (10) стремится к нулю при n — то.

Доказательство. Из условия (11) следует неравенство р(1) < 1. Поэтому неравенство (10) применимо. Обозначим для краткости 7 := 1 — р(1). Возьмем произвольное а € (0,7) и разобьем сумму в (10) на две: до [na] и после [na] + 1, обозначим их Si и £2 соответственно. Очевидно, что

a 1 /и2

Si < пп р( 1) =ехр —-

^1-р2(1) V 1 + Р(1).

Логарифмируя, убеждаемся, что правая часть стремится к нулю при n — то, если

u2 - (1 + р(1))(1 + a) log n —то при n — то. Поскольку 1 + a < 1 + 7 = р(1) < 1, то Е1 — 0 при n — то, если

• , и /1 + р(1)

lim mf —¡= > \ -

n^œ y 2 log n V 2

Далее,

n / 2

n u2

£2 < , = У |r(k)| ехр (--f-— и

Пользуясь условием (11) и суммируя величины к7, получаем, что для любого произвольно малого е > 0 и всех достаточно больших п сумма в правой части этой цепочки не превосходит еп2-7. Отсюда заключаем, что

en2 7 / u2

s2 < . exp -

Vl"P2(l) V 1 + p(KD,

для тех же e и n. Логарифмируя, получаем, что, поскольку e можно брать произвольно малым, Х2 стремится к нулю тогда и только тогда, когда

U2 , ч ,

--(2 — 7)iogn—>00 при п —> то.

1+n-aY к " 6 Н

А это эквивалентно неравенству

liminf — - > 1.

V (2 - 7) log n

Осталось только вспомнить, что 2 — 7 = 1 + р(1). Лемма доказана.

Из этой леммы, утверждения 1, теоремы 3 и следствия 1 вытекает основное утверждение настоящей работы.

Теорема 4. Пусть для ковариационной функции гауссовской стационарной последовательности X(k) выполнено соотношение (11). Пусть натуральное n и уровень u оба стремятся к бесконечности так, что nP(X(1) > u) ^ то, при этом выполнено неравенство (12). Тогда маркированный точечный процесс (8) сходится по вариации к пуассоновскому маркированному процессу (9). Сходимость не медленнее чем max{p(u),7(u,n)}, где через 7(u,n) обозначена правая часть неравенства (10).

4. Замечания. 4-1- Броуновские марки. Если для каких-либо целей потребуется нормировка ¿(^-пачек (марок) Zu(•), которые "сгущаются", т.е. их интенсивность стремится к бесконечности: pni ^ то, то на основе другой теоремы Ю.В. Прохорова [4, теорема 3] можно по этой же схеме доказать сходимость к пуассоновскому маркированному точечному процессу, марки которого — независимые броуновские движения, а также выбрать преимущественное приближение — броуновское или пуассоновское [4, теорема 1].

4.2. Модель Дерриды случайной энергии [9, 10]. Классическая модель случайной энергии Дер-риды — это сумма

Sn(I3) /5 > 0,

k=1

где X(k) — независимые случайные величины. Стандартные задачи для этой модели — исследование предельного поведения средней энергии Sn (e)/N при N ^ то в зависимости от в и предельного поведения распределения нормированной Sn (в)• Часто при этом основную роль играют большие значения величин X(k), как в асимптотическом методе Лапласа, поэтому результаты настоящей работы могут найти здесь применение. Можно также рассмотреть зависимые X(k). Кроме того, теоремы Прохорова позволяют исследовать качество получаемых аппроксимаций, включая асимптотические разложения.

Работа выполнена при частичной поддержке Института системных исследований РАН и Международной лаборатории стохастического анализа и его приложений НИУ "Высшая школа экономики".

Автор приносит благодарность А. В. Булинскому, М.Я. Кельберту и С. А. Молчанову за плодотворные обсуждения.

СПИСОК ЛИТЕРАТУРЫ

п

1. Питербарг В.И. Двадцать лекций о гауссовских процессах. М.: МЦНМО, 2015.

2. Mittal Y., Ylvisaker D. Limit distribution for the maxima of stationary Gaussian processes // Stochast. Process. and their Appl. 1975. 3. 1-18.

3. Piterbarg V.I. Asymptotic methods in theory of Gaussian random processes and fields. Providence: Amer. Math. Soc., Ser. Translations of Mathematical Monographs, 2012. Vol. 148.

4. Прохоров Ю.В. Асимптотическое поведение биномиального распределения // Успехи матем. наук. 1953. 8, № 3. 135-142.

5. Козуляев П.А. Асимптотический анализ одной основной формулы теории вероятностей // Уч. зап. Моск. ун-та. 1939. 15. 179-182.

6. Богачев В.И. Основы теории меры. Т. 1. М.; Ижевск: НИЦ "Регулярная и хаотическая динамика", 2003.

7. Resnick S.I. Extreme values, regular variation, and point processes. N.Y.; Berlin; Heidelberg: Springer-Vеrlag, 1987.

8. Berman S. Sojourns and extremes of stochastic processes. London: CRC Press, Taylor&Francis, 1992.

9. Ben Arous G., Bogachev L., Molchanov S. Limit theorems for sums of random exponentials // Probab. Theory and Related Fields. 2005. 132, N 4. 579-612.

10. Molchanov S., Panov V. Limit theorems for the alloy-type random energy model // Stochastics. 2019. 91, N 5. 754-772 (arXiv:1802.05071v1 [math.PR]).

Поступила в редакцию 16.05.2023

УДК 519

МОДЕЛИ СТРАХОВАНИЯ С ДИСКРЕТНЫМ ВРЕМЕНЕМ

Е. В. Булинская1

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

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

Two discrete-time insurance models are considered. The first model studies non-proportional reinsurance and bank loans. For this model, we establish the optimal control and stability to small fluctuation of parameters and perturbation of random variables distributions describing the model. The second model is dual and the ruin probabilities are compared under assumption that the gains distributions satisfy one of four partial orders.

Key words: discrete-time insurance models, optimal control, stability with respect to small fluctuation of parameters and perturbation of random variables distributions describing the model, ruin probability, stochastic orders.

DOI: 10.55959/MSU0579-9368-1-64-6-6

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

1 Булинская Екатерина Вадимовна — доктор физ.-мат. наук, проф. каф. теории вероятностей мех.-мат. ф-та МГУ, e-mail: ebulinsk@mech.math.msu.su, ebulinsk@yandex.ru.

Bulinskaya Ekaterina Vadimovna — Doctor of Physical and Mathematical Sciences, Professor, Lomonosov Moscow State University, Faculty of Mechanics and Mathematics, Chair of Probability Theory.