Оценки скорости сходимости обобщенных процессов кокса с ненулевым средним и некоторые их применения тема автореферата и диссертации по математике, 01.01.05 ВАК РФ

Артюхов, Сергей Владимирович АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2009 ГОД ЗАЩИТЫ
   
01.01.05 КОД ВАК РФ
Диссертация по математике на тему «Оценки скорости сходимости обобщенных процессов кокса с ненулевым средним и некоторые их применения»
 
Автореферат диссертации на тему "Оценки скорости сходимости обобщенных процессов кокса с ненулевым средним и некоторые их применения"

Московский государственный университет имени М. В. Ломоносова Факультет вычислительной математики и кибернетики

На правах рукописи

АРТЮХОВ Сергей Владимирович

ОЦЕНКИ СКОРОСТИ СХОДИМОСТИ ОБОБЩЕННЫХ ПРОЦЕССОВ КОКСА С НЕНУЛЕВЫМ СРЕДНИМ И НЕКОТОРЫЕ ИХ ПРИМЕНЕНИЯ

01.01.05 - теория вероятностей и математическая статистика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Москва, 2009

003468036

Работа выполнена на кафедре математической статистики факультета вычислительной математики и кибернетики Московского государственного университета имени М. В. Ломоносова

Научный руководитель:

Официальные оппоненты:

доктор физико-математических наук, профессор В. Ю. Королев

доктор физико-математических наук, профессор Ю. С. Хохлов

кандидат физико-математических наук, доцент Е. В. Коссова

Ведущая организация:

Институт проблем информатики РАН

Защита диссертации состоится 29 мая 2009 г. в 11 часов 00 минут на заседании диссертационного совета Д 501.001.44 в Московском государственном университете имени М. В. Ломоносова по адресу: 119991, ГСП-1, Москва, Ленинские горы, МГУ, 2-й учебный корпус, факультет ВМиК, аудитория 685.

С диссертацией можно ознакомиться в библиотеке факультета ВМиК МГУ. С текстом автореферата можно ознакомиться на официальном сайте ВМиК МГУ http://cs.msu.ru в разделе "Наука" - 'Табота диссертационных советов" -"Д 501.001.44".

Автореферат разослан £ I апреля 2009г.

Ученый секретарь диссертационного совета профессор

^ц^, Н. П. Трифонов.

1 Общая характеристика работы Актуальность темы.

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

Пусть N(t) - число скачков случайно блуждающей частицы за период времени [0,i], t > 0. При решении многих практических задач естественно считать, что моменты скачков образуют хаотический точечный случайный процесс на оси времени. Однако этот хаотический случайный процесс может не быть однородным в силу указанных выше представлений. Как известно, наиболее разумными стохастическими моделями неоднородных хаотических точечных процессов являются дважды стохастические пуассоновские процессы, иначе называемые процессалш Кокса (см., например,!2]).

Рассмотрим предложенное Д. Коксом [3] понятие дважды стохастического пуассонов-ского процесса, являющееся естественным обобщением неоднородного пуассоновского процесса.

Пусть N\{t), t > 0, - однородный пуассоновский процесс с единичной интенсивностью, а Л(i), ( > 0, - независимый от Ni(t) случайный процесс, обладающий следующими свойствами: Л(0) = 0, Р(Л (t) < оо) = 1 для любого t > 0, траектории Л (t) не убывают и непрерывны справа. Дважды стохастический пуассоновский процесс N{t), называемый теперь также процессом Кокса, определяется как суперпозиция N\{t) и Л (t):

N(i) = ^(A(i)), t> 0.

В этом случае будем говорить, что процесс Кокса N(t) управляется процессом Л (t). В

1 Королев В. Ю. Вероятностно-статистический анализ хаотических процессов с помощью смешанных гауссовских моделей. - М.: ИПИ РАН, 2007

2 Bening V. Е., Korolev V. Yu. Generalized Poisson Models and Their Applications in Insurance and Finance. - Utrecht: VSP, 2002.

3Cox D. R. Some statistical methods connected with series of events. // J. Roy. Statist. Soc., Ser. B. 1955. Vol. 17. P. 129-164.

f ' L L,

3 О

частности, если процесс A(t) допускает представление

¡

A(t) = J A(r)dr, t > О,

о

в котором A(í) - положительный случайный процесс с интегрируемыми траекториями, то A(í) можно интерпретировать как мгновенную стохастическую интенсивность процесса N(t). Поэтому иногда процесс A(í), управляющий процессом Кокса N(t), будет называться накопленной интенсивностью процесса N(t).

При каждом фиксированном t распределение случайной величины N(t) — Ni (A(f)) является смешанным пуассоновским. Класс смешанных пуассоновских распределений весьма широк и содержит, в частности,

• само пуассоновское распределение (ему соответствует вырожденное смешивающее распределение),

• отрицательное биномиальное распределение (ему соответствует смешивающее гамма-распределение),

• распределение Зихеля, которому соответствует обобщенное обратное нормальное смешивающее распределение,

• бета-пуассоновское распределение при смешивающем бета-распределении,

• распределение Делапорте со сдвинутым гамма-распределением в качестве смешивающего,

• обобщенное распределение Варинга,

• инфекционное распределение Неймана типа А (иначе называемое пуассон-пуассоновским).

Эти и другие примеры смешанных пуассоновских распределений подробно рассмотрены в книгах [4, 5]. В этих книгах, а также в [6] и [2] можно найти детальное описание аналитических и асимптотических свойств процессов Кокса.

В данной диссертации рассматривается частный вариант неоднородного случайного блуждания, в котором число слагаемых в суммах формируется в соответствии с процессом Кокса. Тккие блуждания обычно называют обобщенными дважды стохастическими пуассо-новскими процессами или обобщенными процессами Кокса. Этот случай имеет чрезвычайно важное практическое значение.

Пусть Х-[,Х2,-.. - одинаково распределенные случайные величины, a N{t) - дважды стохастический пуассоновский процесс (процесс Кокса), управляемый процессом A(í), то есть

N(t) = JVi(A(t)), {>0,

где Л (t) - случайный процесс с неубывающими непрерывными справа траекториями, выходящими из нуля, JV](<) - стандартный пуассоновский процесс (однородный пуассоновский

4Grarí(íe/í J. Mixed Poisson Processes. - London: Chapman and Hall, 1997

''Королев О. Ю. Бенинг В. Е., Шоргин С. Я Математические основы теории риска. - М.: ФИЗ-МАТЛИТ, 2007

6 Grandell J. Doubly Stochastic Poisson Process. // Lecture Notes in Mathematics, 1976. Vol. 529. P. 229-234

процесс с единичной интенсивностью), причем процессы JVi(t) и Л(t) стохастически независимы.

Предположим, что при каждом t > 0 случайные величины N(t),Xi,Xi,... независимы. Процесс

называется обобщенным процессом Кокса (при этом для определенности считаем, что ¡Cj=i = 0)- Процессы вида (1) играют очень важную роль во многих прикладных задачах. Достаточно сказать, что при Л(t) = At с А > 0 процесс S(t) превращается в классический обобщенный пуассоновский процесс, широко используемый при моделировании многих явлений в физике, теории надежности, финансовой и актуарной деятельности, биологии и т. д. Большое число разнообразных прикладных задач, приводящих к обобщенным пуассо-новским процессам, описано в книгах [2, 7].

Примеры некоторых задач, сводящихся к обобщенным процессам Кокса, приведены, например, в книгах [2, 7].

Всюду в дальнейшем предполагается, что у случайных величин имеется, по

крайней мере, два первых момента. Обозначим EXi = а О, DXi = и1, 0 < а2 < ос. Как показано, например, в [2, 5], даже в таких предположениях предельные распределения обобщенных процессов Кокса могут иметь произвольно тяжелые хвосты. При этом асимптотическое поведение обобщенных процессов Кокса с нулевым средним принципиально отличается от асимптотического поведения обобщенных процессов Кокса с ненулевым средним. Если в первом случае в качестве предельных законов в оговоренных выше мо-ментных условиях выступают масштабные, то во втором - сдвиговые смеси нормальных законов.

Цель работы

Целью данной диссертации является

• изучение точности асимптотических аппроксимаций распределений обобщенных процессов Кокса, их экстремумов и обобщенных процессов риска с ненулевыми средними сдвиговыми смесями нормальных законов;

• решение задачи оптимизации параметров страховой деятельности, деятельности склада и пункта обмена валют при случайном характере интенсивностей потоков страховых премий и страховых выплат;

• построение гарантированных двусторонних оценок для значения надежности модифицируемых технических или информационных систем при случайной интенсивности потока модификаций.

Научная новизна

Все основные результаты диссертации являются новыми и состоят в следующем;

7Gnedenko В. V., Korolev V. Yu. Random Summation: Limit Theorems and Applications. — Boca Raton, FL: CRC Press, 1996.

m

(1)

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

2. Решена задача оптимизации параметров страховой деятельности при случайном характере интенсивностей потоков страховых премий и страховых выплат.

3. Предложена интерпретация функциональной зависимости спроса и предложения через зависимость интенсивности потока клиентов от параметров спекулятивной деятельности (в частности, от маржи).

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

Методы исследования

В работе используются классические методы теории вероятностей. Базовым теоретическим результатом является теорема 1.1.2, При ее доказательстве используется метод, основанный на представлении конечномерных распределений обобщенных процессов Кокса в виде смесей обобщенных пуассоновских распределений, в которых смешивание производится по соответствующему конечномерному распределению управляющего процесса. Такой подход позволяет использовать известные результаты о точности нормальной аппроксимации для пуассоновских случайных сумм. Используемый метод оказывается достаточно универсальным и может быть успешно применен к построению оценок скорости сходимости распределений некоторых других процессов с ненулевым средним к сдвиговым смесям нормальных законов. В частности, в диссертации с помощью этого же метода указанные оценки построены для распределений экстремумов обобщенных процессов Кокса и обобщенных процессов риска. Затем теорема 1.1.2 сама становится базой для решения оптимизационных задач, рассматриваемых в главах 2-4. Для решения этих задач также используются аналитические и асимптотические свойства пуассоновских случайных сумм.

Теоретическая и практическая значимость

Результаты диссертации имеют теоретический характер. Особенностью представленных в диссертации результатов, отличающей их от предыдущих, является получение универсальной оценки скорости сходимости обобщенных процессов Кокса с ненулевыми средними к сдвиговым смесям нормальных законов. Найденная оценка оказывается применимой и к таким обобщенным процессам Кокса, в которых у управляющего процесса отсутствует дисперсия. Аналогичные универсальные оценки, существенно улучшающие известные ранее, получены для экстремумов обобщенных процессов Кокса и обобщенных процессов риска. Кроме того, впервые получены двусторонние оценки оптимального количества товара в задаче управления запасами при случайном характере интенсивности потока заявок и аналогичные двусторонние оценки начального капитала страховой компании в задаче минимизации издержек при случайном характере интенсивности потоков страховых премий и страховых выплат. Более того, впервые найдены гарантированные доверительные интервалы для нестационарного аналога коэффициента готовности модифицируемых технических и информационных систем (ТИС) с непрерывным временем в рамках экспоненциальной,

логистической и гиперболической моделей изменения надежности ТИС с непрерывным временем.

Апробация работы и публикации.

Основные результаты диссертации опубликованы в 10 печатных работах. Результаты диссертации неоднократно докладывались на научно-исследовательском семинаре «Теория риска и смежные вопросы» на факультете вычислительной математики и кибернетики МГУ, на международном научном семинаре по проблемам устойчивости стохастических моделей в Майори-Салерно (Италия) в сентябре 2005 г., на научных конференциях студентов, аспирантов и молодых ученых «Ломоносов - 2005» и «Ломоносов - 2006», а также на научных семинарах кафедры математической статистики факультета вычислительной математики и кибернетики МГУ и Института проблем информатики РАН.

Структура и объем диссертации

Диссертация состоит из введения, четырех глав и списка литературы, содержащего 78 наименований. Общий объем работы составляет 114 страниц. Работа включает 2 таблицы и 15 графиков.

2 Краткое содержание диссертации

Введение содержит общую характеристику работы, описание объектов исследования и основных результатов.

Первая глава посвящена оценкам точности аппроксимации распределений обобщенных процессов Кокса с ненулевыми средними сдвиговыми смесями нормальных законов. Здесь также решаются задачи построения оценок скорости сходимости как распределений экстремумов обобщенных процессов Кокса, так и распределений обобщенных процессов риска. Эти оценки являются существенным уточнением и обобщением известных оценок (см., например, [2, 8]).

Случай управляющих процессов с конечной дисперсией. Рассмотрим обобщенный процесс Кокса (1). Предположим, что существуют такие числа i. € (0, оо) и s 6 (0, оо), что

EA(í) = Ü, DA(t) = s3í, t > 0.

Стандартную нормальную функцию распределения будем обозначать Ф(х), х € К. Символ => будет обозначать сходимость по распределению.

В книге [2j доказана общая теорема о необходимых и достаточных условиях сходимости одномерных распределений обобщенных процессов Кокса, частным случаем которой является следующий результат.

Теорема 1. Пусть а ф 0. В дополнение к приведенным выше условиям на моменты

Р

управляющего процесса Л(t) предположим, что Л(t) —» оо при f —» оо. Тогда одномерные распределения неслучайно центрированных и нормированных обобщенных процессов Кокса слабо сходятся к распределению некоторой случайной величины Z при t —> оо, то есть

, У*-"*, (i-оо).

8 Королев В. Ю., Соколов И. А. Основы математической теории надежности модифицируемых систем. - М.: ИПИ РАН, 2006.

тогда и только тогда, когда существует случайная величина V такая, что

Л(0 - И

ауД

■V (I оо).

При этом

Р(г <х) = ЕФ(Х,1+ . - У 16 к.

Несложно видеть, что предельная случайная величина ¿Г допускает представление „2,2 1-1/2

г1

1 +

(а2 + о2)(

■Х +

л/ (а2 + а2)( + а2«2

• V,

где ЛТ - случайная величина со стандартным нормальным распределением, независимая от случайной величины V.

В приводимой ниже теореме изучается вопрос о скорости сходимости в теореме 1. В книге [2] приведены некоторые оценки скорости сходимости в теореме 1, но они справедливы для довольно узкого класса распределений предельной случайной величины V, довольно громоздки, содержат трудно вычисляемые характеристики и неудобны для анализа и применения. Здесь будет приведена оценка скорости сходимости в теореме 1 в традиционных терминах.

Прежде чем сформулировать основной результат, введем дополнительные обозначения:

р3 = Е|Хх|3, ¿з-

/Н = вир

(о2 + ст2)3/2 (х) - ЕФ^ц/1 +

, Я(х) = Р

5(«) - аН

V анУ

< х

(а2 + <г2К

<т2)е;

■ вир

Теорема 2. Пусть /З3 < оо, Е|У| < оо. Тогда справедлива оценка

Рь:

VI е£(0,1) \ ^(1 - ()1 I

где С0 < 0.7005,

<3{е) = тах|^,

П + 1

(1 + уГ^-ё) \/2тхе{\ - б)

Следствие 1. Пусть в дополнение к условиям теоремы 5 с&мейстпво случайных вели-

II /(>0

равномерно интегрируемо. Тогда справедлива оценка

Л<д( + 1=. ш£ (-^ + £(1 + <3(0)1

дД «6(0,1) 1 ^/(1 - е)£ ¿V« V/■

Замечание 1. Пусть, к примеру, процесс Л(() является пуассоновским с интенсивностью I = а2. В таком случае распределение случайной величины при каждом { > 0 является так называемым пуассон-пуассоновским или (инфекционным) распределением Неймана типа А. Класс таких распределений введен Ю. Нейманом в работе [9] в связи с некоторыми задачами из области бактериологии и энтомологии. В этом случае величина также имеет порядок поскольку в таком случае справедлива оценка

причем случайная величина V имеет стандартное нормальное распределение, см., например, [5] и, следовательно, распределение случайной величины X также является стандартным нормальным. В таком случае имеет место оценка

Pt<

Замечание 2. Другим примером управляющего процесса Л(i), при котором величина Д( также имеет порядок 0(i-1/2) является случай, когда A(i) является гамма-процессом Леви, приращение которого на единичном интервале имеет гамма-распределение с параметром масштаба Л = i/s2 и параметром формы а = £2/s2. В таком случае распределение случайной величины N{t) при каждом i > 0 является отрицательным биномиальным (см., например, [5]). При этом случайная величина V имеет стандартное нормальное распределение, и, следовательно, распределение случайной величины Z является стандартным нормальным.

Замечание 3. Если случайный процесс N(t) является однородным пуассоновским (с интенсивностью I), то, очевидно, s = 0. При этом логично считать, что распределение случайной величины V является вырожденным в нуле и At = 0. Поэтому в таком случае оценка, приведенная в теореме 2, естественно переходит в оценку скорости сходимости распределений пуассоновских случайных сумм к нормальному закону, приведенную, например, в работе [10] с константой Со, уточненной в работе [и].

Замечание 4. При известных значениях основных параметров можно легко найти численные оценки величины

£6(0,1) I /(1-с)7 f

)l}

Найдем, например, значение данного выражения при условии, что Со = 0.7005, а = 5, а = 10,1 = 4, з = 2, = (что характерно для случая стандартно нормально распределенной случайной величины V), £з = 1,5108. Тогда с помощью численной оптимизации по е получаем

сб"о,1) 1 V(1 - е)е IV

2,3882.

eNeyman J. On a new class of "contagious" distribution, applicable in entomology and bacteriology.

// Annals of Mathematical Statistics, 1939. Vol. 10. P. 35-57.

10Шевцова И. Г. Уточнение структуры оценок скорости сходимости в центральной предельной теореме для сумм случайных величин. // Дис. канд. физ.-матем. наук. - М.: МГУ, 2006.

11 Шевцова И. Г. Уточнение верхней оценки абсолютной постоянной в неравенстве Берри-Эссеена. //Труды молодых ученых факультета вычислительной математики и кибернетики МГУ. - М.: изд-во факультета вычислительной математики и кибернетики МГУ. 2008. С. 101-110.

Следовательно, для данных значений параметров оценка скорости сходимости имеет вид

2,3882

Pi < А( + :

уД

Случай управляющих процессов с бесконечной дисперсией. Предположение о существовании дисперсии управляющего процесса Л(<) в теореме 1 и теореме 2 не является критическим. Можно доказать аналоги этих результатов для того случая, когда предполагается лишь существование математического ожидания процесса Л(С).

По аналогии с обозначениями, введенными выше, положим

Pt = sup

PfS(t)-at

Ч Vi

) VvV + tf2)/

: SUp

< v ) -P(V<v}

Теорема 3. Пусть /З3 < oo, EA(f) < оо при каждом t > О, E|V| < оо. Тогда справедлива оценка

' " ~ Л(t)-t

Vi <e(o,l) I VT^l £

Vi

Замечание 5. В качестве примера ситуации, на которую распространяется действие теоремы 3, но к которой не применимы теоремы 1 и 2, можно привести такую, в которой

A(t) = шах{0, VtV + t} + ^ {—^Vt - l),

где 2 < а < 3, а V - случайная величина с плотностью

а + 1

р[х) =

2(1*1 + 1)" '

iSR.

Несложно проверить, что EA(i) = t при любом i > 0, но второй момент случайной величины A(i) бесконечен вследствие бесконечности второго момента случайной величины V (и, следовательно, бесконечен и второй момент обобщенного процесса Кокса S(t), управляемого таким процессом Л(<)). При этом, однако, как легко видеть,

Л(t) - t V~t

при t —> оо. Этот случай является наглядной иллюстрацией очень интересного и нетривиального факта: в отличие от классической теории суммирования, для сумм со случайным числом слагаемых (в частности, для обобщенных процессов Кокса) с бесконечной дисперсией существование нетривиальных слабых пределов возможно и при нормировке порядка t1/2, являющейся «стандартной» в классической теории лишь для сумм с конечной дисперсией.

Метод, использованный при доказательстве теоремы 2, оказывается достаточно универсальным и позволяет получить оценки, аналогичные приведенным в теореме 2, как для экстремумов обобщенных процессов Кокса (S(t) - maxo<T<tS(t),), так и для обобщенных процессов риска {R[t) = cA(t) — Ху). Этим оценкам посвящены заключительные разделы главы 1.

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

Процесс извлечения спекулятивной прибыли рассматривается на примере пункта обмена валют. Пусть в начальный момент времени пункт обмена валют обладает некоторой суммой денег. Предположим, что у данной организации существует возможность как купить, так и продать любое количество валюты по некоторой цене. Назовем эту цену ценой обмена и пусть она равна c(t) в некоторый момент времени t. Тогда часть денег в начальный момент времени <о = 0 обменивается на валюту по цене c(io) = с* на валютном рынке. Затем обменный пункт выставляет свои цены на покупку и продажу, которые, соответственно, меньше и больше, чем с*. При этом цены покупки и продажи пункта обмена валют должны быть определены таким образом, чтобы к концу отчетного периода получить максимально возможную прибыль от проведенных операций.

Разницу между ценой покупки и ценой продажи будем называть маржей (margin). Отклонение цены покупки от цены обмена будем называть спрэдом покупки (buy spread), а отклонение цены обмена от цены продажи - спрэдом продажи (sell spread).

Опишем формально рассматриваемую модель. Определим при т > 0 случайные процессы

N+(т) N-(t)

М(т) = и + с+ £ Л-+-С- £ XJ

3=1 з=г

и

N+{ г) Л-(г)

G(r) = „-

;=■1 ¿=i

причем ]Cj=i(') = 0. Здесь и > 0 имеет смысл начального капитала пункта обмена валют, величина v > 0 определяет начальный объем имеющейся валюты, положительные числа с+ и с" имеют смысл цены продажи и цены покупки валюты обменным пунктом, соответственно, неотрицательные случайные величины и - это количество валюты, проданной ¿-тому клиенту и купленной у l-того клиента, соответственно, а целочисленные случайные процессы Лг+(т) и N~(t) определяют количество клиентов, пришедших в пункт обмена валют для того, чтобы купить или продать валюту, Везде далее предполагается, что все перечисленные случайные величины и процессы являются независимыми, а последовательности {ЛГ+}„>1 и {Х,Г}п>1 состоят из одинаково распределенных случайных величин. Таким образом, процессы М(т) и G(r) характеризуют капитал компании и объем имеющейся валюты в некоторый момент времени.

Везде далее будем представлять числа с+ и с~ в виде

с+ = с" + <5+, с" = с* - Г,

где <5+ > 0 и 5~ > 0 - это спрэд продажи и спрэд покупки соответственно.

Заметим, что при увеличении цены продажи с+ количество клиентов, желающих купить у обменного пункта валюту, в среднем за единицу времени должно уменьшиться или, другими словами, должна уменьшиться интенсивность потока клиентов-покупателей. С другой стороны, очевидно, при приближении цены продажи к цене с*, интенсивность потока клиентов должна увеличиться. Таким образом, вполне разумным предположением является обратная зависимость интенсивности потока клиентов, покупающих валюту, от разницы

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

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

виде:

ЛГ+(г) лг-(т) /Л'-(т) N+{t)

М(г) + c(r)G(r) = и + c(r)u + 6+ ¿ Х++Г ¿ Хг+Д ¿ jf-_ ¿ X¡

3=1 j=i \ J=1

где Д(г) = с(т) - с*

Объектом дальнейшего изучения является распределение случайной величины

1V+(r) JV~(r) / W-(r)

ñ(r) = 6+ ¿ Х+ + Г £ *7 + Д(г) »+ ¿ X7

J=1 j=l \ J=1

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

Целью обменного пункта является определение 6+ и 5~ таким образом, чтобы прибыль от совершенных операций была в некотором смысле максимальной. Для поиска оптимальных спрэдов покупки/продажи были сформулированы следующие задачи:

Задача 1. Определить и 5~ таким образом, чтобы максимизировать ожидаемую прибыль:

ЕЯ —» шах; (2)

¿+,í-

Задача 2. Найти значения <5+ и 6~, обеспечивающие для некоторого г > 0 и близкой к единице величины 7 соотношение

P(ñ > г) > 7, (3)

то есть требуется найти курсы покупки и продажи валюты, которые дают возможность получить заданную прибыль с достаточно большой вероятностью. Кроме того, соотношение (3) дает возможность оценить величину гарантированной прибыли при заданных значениях <5+ и <5~.

В данной главе рассматриваются несколько частных случаев описанной выше общей постановки задачи.

Случай неоднородных потоков клиентов. Предположим, во-первых, что цена обмена не изменяется во времени, то есть с(т) = c*,Vr > 0. Во-вторых, мы будем предполагать, что Лг+(г) и N~(t) являются дважды стохастическими пуассоновскими процессы (процессы Кокса), управляемые процессами Л+(т) и Л~(т) соответственно. Данное предположение дает возможность построения более реалистичных моделей, так как существуют сезонные колебания спроса как на покупку, так и на продажу валюты (например, в будние дни сделок совершается меньше, нежели в выходные), а также существуют различные случайные всплески интенсивности (например, связанные с выходом новостей, изменяющих стоимость валюты относительно рубля).

Величина прибыли компании в некоторый момент временя t будет выражаться формулой

JV+(t) N-{t)

7=1 3=1

Введем новые обозначения

N•>"(1) tr-(t)

J=1 3=1

Тогда

Везде далее для этого случая мы будем предполагать, что <5+ = б" = 5, то есть ситуация на валютном рынке достаточно стабильна.

Решим поставленные ранее задачи 1 и 2 для данного случая. Дли этого дополнительно предположим, что ЕЛ+(() = и ЕЛ-(4) = А_(. Тогда ожидаемая прибыль компании будет равна

ЕЯ(0 = Л- +

При решении задачи 2 возникают вполне очевидные затруднения, поскольку даже для простейших видов распределений случайных величин и и в предположении, что интенсивности потоков клиентов постоянны, то есть Л'+(т), N~(т) являются пуассоновскими процессами, вычисление функции распределения Л(<) в явном виде представляется крайне затруднительным. Поэтому для оценки функции распределения ЩЬ) используется теорема 2 и свойство регулярности равномерной метрики (см. [12]).

Таким образом, получаем верхние и нижние оценки для гарантированной прибыли:

П = (а+А+ + о-А~)<Й + х/[(Д+)2 + (Д~)2] Ь ■ и - 7 + ) .

= (а+А+ + а-\~)И + у/[(Л+)2 + (1'-)2]< • и - 7 - К+ ,

гг

где

= ^[А+((а+)2 + (<т+)2) + (а+)2(5+)2]52, = х/[А-((а-)2 + (а")2) + (а-)2(*г)2]<52.

'е(о,1) I ^(1- е)А+ А+\ £ .е(0,1) I - б)А- А- V б "Л

для 7 £ (0,1) символом 11(7) обозначается 7-квантиль стандартного нормального закона.

Случай однородных потоков клиентов. Как и для случая неоднородного потока клиентов получены двусторонние оценки для гарантированной прибыли в смысле задачи 2. Пусть N'i' и - проекции процессов Лг+(г) и ЛГ~(г) в момент времени t. Предположим, что случайные величины и имеют пуассоновское распределение с параметрами А+ и А-. Также предположим, что цена обмена не изменяется во времени, то есть ф) = с*,\Л > 0.

Тогда

и(:. 7 _ с + Л>з \ + + л

V (а+^ + Л-^Г / v

+ + А"/^) 5 < г < (А+д/ + А-^Г) <5+

+" (1 ■~ 7 + СоЧ;Л"1з/2) + (4)

V (А+/4 + \-р2) !

12 Золотарев В. М. Современная теория суммирования независимых случайных величин. — М.: Наука, 1986., стр. 73-75

при выполнении ограничений:

АУз + А->з ^ А+/4 + А рз ^ „

{X+ßt + X-ßzf2 (А+д+ + А-^)'/2

Этот результат дает возможность определить пределы будущей прибыли при установленной 5 и известных моментах и (1 < к < 3), которые могут быть оценены статистически. Заметим также, что неравенство (4) можно уточнить, имея дополнительную информацию о моментах и ^ (к > 4), применив, например, теорему 1.5.1 из [13].

При поиске оптимального значения 6 в задачах 1 и 2, необходимы дополнительные предположения о виде зависимости среднего количества клиентов, которые хотят купить/продать валюту, от S. Поскольку в различных условиях функция спроса на товар имеет разный вид, в работе рассматриваются несколько модельных примеров зависимостей интенсивностей потоков клиентов А+ и А~ от величины S: случай линейной, обратно-пропорциональной, экспоненциальной и гауссовской зависимости. Для каждого из случаев при конкретных значениях параметров найдены оптимальные значения 5 в смысле задачи 1 и 2 или показано, что оптимумов несколько.

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

Предположим, что в некоторой организации в начальный момент времени t — 0 имеется и единиц некоторого ресурса, На протяжении интервала времени [О, Т] в данную организацию поступают требования двух видов: во-первых, требования на выдачу данного ресурса, во-вторых, требования о принятии на хранение этого же ресурса. Предположим, что S(t) -величина чистых списаний ресурса из организации к моменту времен t. Тогда, и - S(i) - величина текущего остатка ресурса в организации. Требуется найти оптимальное в некотором смысле значение начального остатка в организации.

Для решения поставленной выше задачи введем следующие функций издержек, связанные с величиной и. Пусть ci(t) - издержки, связанные с хранением избыточного резерва. Избыточность понимается в смысле положительности величины E(u — S(t))+ (здесь мы используем стандартное обозначение х+ = тах{0, х}). Пусть Ci(t) - издержки, выраженные в недостатке ресурса. Недостаток понимается в смысле положительности величины E(S(t)-u)+.

Тогда функция суммарных издержек организации имеет следующий вид:

Г т

D(u) = J ci(t) • E(u - S(i))+ dt + j c2(t) ■ E(S(t) - u)+ dt. о 0

Необходимо найти такое uo, что

D(u0)^mmD{u). (5)

хгВенинг В. Е., Королев В. Ю. Введение в математическую теорию риска. — М.: МАКС-Пресс, 2000.

В работе показано, что поиск минимума в (5) эквивалентен отысканию решения уравнения

Т

У Ы0 + сг(0]Р(5(0 <и)сй = У сг(0Л.

При предположении, что функции издержек первого и второго рода не зависят от времени (данное предположение справедливо только для относительно коротких временных интервалов), то есть сх(() = с\, сг(0 = с.%, где С1 и с% - константы и пусть 6 = (С1С+С2) ■ получаем, что для нахождения оптимального щ нужно решить уравнение

Г

±¡9 (5(0<")Й = 5. (6)

о

Далее, это довольно общая постановка задачи применяется сперва для решения задачи об оптимизации резерва склада, потом для оптимизации параметров процесса риска и, наконец, для оптимизации начального капитала пункта обмена валют.

Оптимизации резерва на складе. Пусть XI - количество единиц товара, которое должно быть поставлено (изъято со склада) в соответствии с г-ой заявкой, а N1}) - количество заявок, поступивших в течение интервала времени [0, (0 < / < Т). Тогда суммарное количество товара, которое должно быть поставлено со склада за время [ОД], очевидно, имеет вид (1),

Предположим, что стандартизованная накопленная интенсивность асимптотически нормальна в том смысле, что

уу _ ос)1 где у ^ лг(0,1). (7)

Далее предположим, что

Дг = Бир

< V - Р (V < у)

Тогда двусторонние оценки для оптимального в смысле (5) начального резерва на складе имеют следующий вид

аЛ (т8 - -Д— - 2Кч/т) <и0<а(1п+^г-+2К\/т]. (9)

где

К

ш { «л +2(§!а+в(«л \

«е(о,1) \ ^(1-е)г е\ с 7/

1 = ^, Г> = = ч/Л2 + (а2 + а2)1

Оптимизации параметров процесса риска. Пусть 0 < Ь < Т, - число страховых выплат до момента 4, а - страховая выплата по г-му страховому случаю такие. Тогда

Щ)

5(0 = *>0, (Ю)

^=1

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

m

R(t) = и - S(t) = u + сЩ - Х5 3=1

имеет смысл резерва страховой компании в момент времени t.

Гарантированные оценки для начального капитала страховой компании. Предположим, что верны (7) и (8). Тогда двусторонние оценки для оптимального в смысле (5) начального капитала страховой компании представимы в виде

a'[T-TS+ , + 2KVT| <щ<а' (т-П--* -2KVT) . (11) \ v2-K\L\ ) \ V27ri2 /

где

в' = (в-с)<, & = у/5Щ, ¿ = 2^.

Стратегию функционирования страховой компании разумно выбирать исходя из требования о том, чтобы в конце рассматриваемого периода резерв, имеющийся в распоряжении страховой компании R(t), был достаточно велик. Это требование можно понимать по-разному. Можно, например, задать некоторое число До и потребовать, чтобы

Ей(0 > До (12)

Можно также задать некоторое число 7 е (0,1) (по возможности значение 7 должно быть близким к единице) и До > 0 и потребовать, чтобы

Р (R{t) > До) > 7 (13)

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

Задача 3. Определить значение величин стартового капитала и и ставки страховой премии с, гарантирующие выполнение требования (12) при условии, что средние издержки минимальны.

Решение этой задачи дается неравенством

-1

/До D' \ ( 2КУ1

(14)

Задача 4. Определить значение величин стартового капитала и и ставки страховой премии с, гарантирующие выполнение требования (13) при условии, что средние издержки минимальны.

Решение этой задачи дается неравенством с>а +

Гарантированные оценки для времени достижения желаемого значения резерва. Неравенство (11) позволяет поставить еще две задачи и получить оценки для их решений.

Задача 5. При заданной ставке страховых премий с найти такое То, чтобы для Т >То было выполнение условие (12) при условии, что средние издержки минимальны.

Решение этой задачи дается неравенством

(_L_

\\/2ж(с-

Т > — " 4 Р

+ 2 К}2 + 7^ + 1 - ..+2*1

Т>

- 4(5(c-a)f)

Задача 6. При заданной ставке страховых премий с найти такое То, чтобы для Т > То было выполнено условие 13) при условии, что средние издержки минимальны.

Решение этой задачи дается неравенством

^(тфаЦ + + "itD'VT)2 ((с - + iSRolc - a)l +

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

Пусть в каждый момент времени t надежность системы можно характеризовать параметром P(t) - вероятностью того, что на сигнал, поданный на вход системы в момент f, система отреагирует правильно. По смыслу такая характеристика надежности ближе всего к традиционно используемому коэффициенту готовности. В случайные моменты времени У0 = 0 < П < У2 < ... в систему вносятся (мгновенные) модификации, в результате чего изменяется параметр P(t). Предположим, что траектории процесса P(t) непрерывны справа и кусочно-постоянны, так что P(t) = P(Yj) при t е [Yj,Yj+1), j > 1.

Задача прогнозирования поведения процесса P(t) чрезвычайно важна. Описанная выше очень общая схема может быть переформулировала в терминах, традиционных для столь разных областей знания, как медицина, программирование или менеджмент. Более подробно об этом см. в книгах [и, 7, 8]

Предположим, что точки {УЬ = 0, УьУг,...}, в которые осуществляются модификации системы, образуют дважды стохастический пуассоиовский точечный процесс , управляемый некоторым процессом A(i). Всюду в дальнейшем будет предполагаться, что EA(t) = t. Это предположение можно интерпретировать и как то, что управляющий процесс в среднем пропорционален времени, и (что существенно важно для построения предельных аппроксимаций) как то, что задача параметризована математическим ожиданием управляющего процесса.

Неоднородные экспоненциальные модели с непрерывным временем. Пусть {(Oj,Tjj)}j>i - последовательность независимых одинаково распределенных двумерных случайных векторов таких, что

0< 01 <1, 0<771<1 почти наверное.

ыКоролев В. Ю. Вероятностные модели. Введение в асимптотическую теорию случайного суммирования. — М.: МАКС Пресс, Москва, 1997.

Отметим, что независимость 6j и тц внутри каждой пары равно как и совпадение распределений и Vj внутри каждой пары не предполагается. Однако мы будем считать, что последовательность пар {(0j,Vj)}j>i стохастически независима от точечного случайного процесса Уь Уг, • ■ •

Задав начальную надежность ро, рассмотрим модель, определяемую рекуррентным соотношением

P(Yj+1) = 6j+1P(Yj) + (1 - P(Yj)), j > 0. (16)

Эта модель называется неоднородной экспоненциальной моделью с непрерывным временем. В такой модели случайные величины 6j описывают возможное уменьшение надежности из-за некачественных модификаций, в ходе которых вместо исправления существующих дефектов в систему могут быть внесены новые, в то время как величины rjj описывают повышение надежности за счет исправления дефектов. Частные случаи модели (16) с двухточечными распределениями случайных величин Bj и 7jj и дискретным временем рассматривались в книгах [15, 1б]. В свою очередь, эти частные случаи представляют собой переформулировку в терминах теории надежности одной модели обучаемости, рассмотренной в книге [п].

Обозначим E(?i = 1 - о, Ещ = Ъ.

Тогда верна следующая теорема.

Теорема 4. Для любого t > 0

Из теоремы 4 вытекает, что, если Ее~(а+,>'Л0 —» 0 при t —» оо (что может иметь место, если, к примеру, A(i) —» оо по вероятности при t -+ оо), то в зависимости от знака величины с = (о + 6)ро-Ь ожидаемая надежность системы либо возрастает (если с < 0), либо убывает (если с > 0). Однако в любом случае справедливо

Следствие 2. Пусть Л (t) —> оо по вероятности при f —> оо. Тогда

lim EP(t) = *

t—»oo w a + b

В силу неотрицательности и ограниченности случайных величин 9j из следствия 2 в свою очередь вытекает

Следствие 3. Пусть Л(i) —> оо по вероятности при i —» оо. Соотношение

lim ЕP{t) = 1

i—»oo

имеет место в том и только в том случае, когда P(0i = 1) = 1.

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

'"Волков Л, И., Шишкевич А. М. Надежность летательных аппаратов. — М.:ВШ, 1975.

'"Волков Л. И. Управление эксплуатацией летательных комплексов. — М.: ВШ, 1981.

17Буш Р., Мостеллер Ф. Стохастические модели обучаемости. — М.: ГИФМЛ, 19G2.

Следствие 4. Пусть Л(0 = I, то есть Я(г) - стандартный пуассоновский процесс. Тогда для любого Ь > О

Рассмотрим ситуацию, когда 9} - 1 почти наверное, более подробно. В этом случае соотношение (16) можно переписать в виде

1-Р(У,+1) = (1-^+1)(1-Р(^)), ,">1. (17)

Обозначим ¡оя( 1 - 7/;) = Тогда из (17) получаем

Щ)

1о8(1-Р(())-1оё(1-р0) = ^а- (18)

В правой часта представления (18) стоит обобщенный дважды стохастический пуассоновский процесс (обобщенный процесс Кокса). Это обстоятельство позволяет воспользоваться предельными теоремами для обобщенных процессов Кокса для уточнения асимптотического поведения надежности системы в рамках неоднородной экспоненциальной модели с непрерывным временем.

Предположим, что 0 < = сг2 < оо, и обозначим = а. Заметим, что, так как О < ^ < 1, то а < 0. В дополнение к условию ЕЛ(<) = введенному выше, предположим, что ОЛ(4) = Л для некоторого э £ [0,оо).

Дополнительно предположим, во-первых, что Е|0|3 < оо, и обозначим

/33 = Е|С1|3, £з=(а2Д2)3/2.

Во-вторых, предположим, что выполнено соотношение

sup

Vi

при некотором к 6 (0,оо). Обозначим

ее(0,1) I у/1 - £ у/тте }

Тогда нижняя гарантированная доверительная граница г7(<) для Р(Ь) с коэффициентом доверия 7 имеет вид:

y(i) = 1 - exp {at + \/t[a2(l + s3) + сг2] • u(7 + K/Vt) + log(l - p0)},

бозначается 7-кванти, 11 > О

P(P(i) > z7(t)) > 7.

I nc ГСЯ

P{zW(t)<P(t)<zW(t))>r

где для 7 £ (0,1) символом и(-/) обозначается 7-квантиль стандартного нормального распределения. При этом при каждом { > О

Найдем гарантированную доверительную полосу {(z!f\t), z^\t)) : t > 0} для Pit) с коэффициентом доверия 7, которая понимается в следующем смысле: при каждом t > О

Она задается выражениями

4°(0 = 1 - exp {at + u(i(7 + 1 + 2K/Vt)) \/t[a2(l + s2) + ст2] + log(l - p0)},

4а)(0 = 1 - exp {at - u(|(7 + 1 + 2K/Vt))\Ji[a2(l + s2) + a1] + log(l - p0)}. Неоднородные логистические модели с непрерывным временем. Обозначим

<3=

1 -Р(ЪУ

Пусть 61,62, ■■■ - независимые одинаково распределенные случайные величины. Будем считать, что последовательность {(0}Ь>1 стохастически независима от точечного случайного процесса Уд, У2. ■ - ■ Предположим, что

= з > о. (19)

Эту модель изменения надежности, называемую логистической, можно интерпретировать следующим образом. Если р^- = Р(У}) - вероятность успеха в последовательности испытаний Бернулли, то величина = р_?/( 1 — р?) характеризует ожидаемый номер испытания, заканчивающегося первой неудачей в рассматриваемой последовательности испытаний Бернулли. Таким образом, величина = характеризует ожидаемое время жизни (безотказной работы) системы после ;-ой модификации, как если бы после момента У] в нее не вносились бы никакие изменения, Следовательно, соотношение (19) можно интерпретировать как формализацию того, что каждая модификация системы изменяет ожидаемое время безотказной работы после модификации в случайное число раз. Положим

<2(г) = <?(:Ут), г > о.

Обозначим Е01 = а. Пусть ро - надежность системы в момент t — 0■ Обозначим

Ро

90 =

1 - Ро'

Теорема 5. В рамках неоднородной логистической модели справедливо соотношение

Е<?(() = воЕе'0"1)^", I > 0.

Следствие 5. В рамках неоднородной логистической модели справедливо соотношение

опЕе(в_1'лИ ЕР(*)< г, 1ШЛ, t>0.

к ' 1 + доЕе^-'М')

Следствие 6. Если Лф) - стандартный пуассоновский процесс, то в рамках неоднородной логистической модели при любом г > 0 справедливы соотношения

= ®Ее<—1>*. Е Р(1)<

1 + дое(а-1)('

Обозначим logfy = Xj■ Тогда из (19) получаем соотношение

i k=1

откуда вытекает, что

N(t)

logQ(t) -logq0 =

k=i

Предположим, что Exj = а, 0 < Dxj = с2 < оо и обозначим /З3 = E|xi|3 < оо.

Тогда гарантированная нижняя доверительная граница для надежности системы P(t) в рамках модели (19), которая для заданного коэффициента доверия 76 1) имеет вид

£,(<) =

exp{i7(t)}

1 + ехр {x7(i)}'

где _

х= at + u(l - 7 - K/Vt) ^¿[а2(1 + *2) + <72] + log q0.

При этом при каждом t > О

P(P(i)>i,(t))>7.

Гарантированная доверительная полоса в рамках логистической модели (19) задается следующими выражениями

г(1) = exp {at - ц(§(7 + 1 + 2К/Щ y/t[a2(l + s2) +7Щ + log go} ^ 1 + exp {at - u(i(7 + 1 + 2K/V~t)) ч/ф2(1 + s2) + a2] + log go}'

(2) exp {at + u(l(7 + 1 + 2K/Vt)) уЛ[а2(1 + a2) + a2] + log go}

1 + ехр {аг + и(^(7 + 1 + 2К/уД))^{аг{1 + 82) + сг2] + к^9о} '

Неоднородные гиперболические модели с непрерывным временем. Рассмотрим модель

3(Уя.1) = С?(У)) + ^+1, ] > 0. (20)

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

<3(0 = 4 > 0.

Обозначим Е0] = а. Пусть, как и ранее, ро - надежность системы в момент Ь - 0, 90 = Ро/(1-Ро)-

Из рекуррентного соотношения (20) следует, что

Щ)

<э(0 = ?о + 4>0-

¿=о

Предположим, что /З3 = Е|01|3 < оо

Запишем гарантированную нижнюю доверительную границу для надежности системы Р(0 в рамках модели (20), которая для заданного коэффициента доверия 7 6 1) имеет вид

где _

2^(0 = аЬ + и(1 - 7 - К/чД) 0^(1 +в2)+ <т2] + 9о.

При этом при каждом < > 0

Р(Р(0 > £,(()) > 7-

Гарантированная доверительная полоса в рамках гиперболической модели (20) имеет вид (1) = оЛ - и{\{ч + 1 + 2К/У1))^1р(1 + а2) + а2] 4- до 4 и 1 + оЛ - и{\{п + 1 + 2К/Щ%/<[а2(1 + в2) + о-2] + 9о' 2<2) = оЛ + и(|(7 + 1 + 2К/Щ^{\ + а') + <т»] + до

2-у 1 + а« + и(|(7 + 1 + 2 К/Щ \/Ф2(1 + «2) + <г2) + 90'

Список публикаций автора по теме диссертации

Научные статьи

1. Артюхов С. В., Вазюкипа О. А., Королев В. Ю., Кудрявцев А. А. Модель оптимального ценообразования, основанная на процессах риска со случайными требованиями // Системы и средства информатики, 2005. Специальный выпуск. С. 207-224.

2. Artyukhov S., Bazyukina О., Korolev V., Kudryavtsev A. On optimization in a demand and supply problem. // Transactions of the XXV International Seminar on Stability Problems for Stochastic Models, Maiori/Salerno, Italy, 20-24 September 2005. University of Salerno, P. 25-31.

3. Артюхов С. В. Применение асимптотических свойств пуассоновских случайных сумм в одной задаче оптимизации резерва // Системы и средства информатики, 2006. Специальный выпуск. С. 238-248.

4. Артюхов С, В., Вазюкипа О. А., Королев В. Ю., Кудрявцев А. А., Шевцова И. Г. Об оптимизации спекулятивной прибыли на примере пункта обмена валют // Актуарий, 2008. No 1. С. 50-56.

5. Артюхов С. В., Жалыбина И. Я., Пузановский А. А. О методе оптимизации прибыли маркет-мейкера // Управление финансовыми рисками, 2008. Выи. 2. С. 35-48.

6. Артюхов С. В., Королев В. Ю. Оценки скорости сходимости распределений обобщенных дважды стохастических пуассоновских процессов с ненулевым средним к сдвиговым смесям нормальных законов // Обозрение прикладной и промышленной математики, 2008. Т. 15, Вып. 6. С. 988-998.

7. Артюхов С. В. Королев В. Ю.. Неоднородные рекуррентные модели изменения надежности модифицируемых систем. Непрерывное время, // Информатика и ее применения, 2008. Т. 2. Вып. 4. С. 57-65.

8. Артюхов С. В. Оценки скорости сходимости распределений экстремумов обобщенных процессов Кокса с ненулевым средним к сдвиговым смесям нормальных законов // Информатика и ее применения, 2009. Т. 3. Вып. 1. С. 69-74.

Тезисы докладов и конференций

9. Артюхов С. В., Вазюкина О. А., Кудрявцев А. А. О методе расчета курсов покупки и продажи валюты обменным пунктом. // Материалы международной конференции студентов, аспирантов по фундаментальным наукам «Ломоносов 2005», 2005. Секция ВМиК. С. 6-8.

10. Артюхов С. В. Об определении оптимального количества рублевых и валютных средств в обменном пункте. // Материалы международной конференции студентов, аспирантов по фундаментальным наукам «Ломоносов 2006», 2006. Секция ВМиК. С. 6-8.

Напечатано с готового оригинал-макета

Издательство ООО "МАКС Пресс" Лицензия ИД N 00510 от 01.12.99 г. Подписано к печати 17.04.2009 г. Формат 60x90 1/16. Усл.печ.л. 1,25. Тираж 75 экз. Заказ 201. Тел. 939-3890. Тел./факс 939-3891 119992, ГСП-2, Москва, Ленинские горы, МГУ им. М.В. Ломоносова, 2-й учебный корпус, 627 к.

 
Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Артюхов, Сергей Владимирович

Введение

1 Оценки точности аппроксимации обобщенных процессов Кокса с ненулевыми средними смесями нормальных законов

1.1 Обобщенные процессы Кокса

1.1.1 Случай управляющих процессов с конечной дисперсией.

1.1.2 Случай управляющих процессов с бесконечной дисперсией.

1.2 Экстремумы обобщенных процессов Кокса

1.3 Обобщенные процессы риска.

2 Оптимизация параметров спекулятивной деятельности

2.1 Общая постановка задачи

2.2 Неоднородные потоки клиентов.

2.3 Однородные потоки клиентов.

3 Применение асимптотических свойств обобщенных процессов Кокса в некоторых задачах оптимизации

3.1 Общая постановка задачи

3.2 Применение асимптотических свойств обобщенных процессов Кокса в одной задаче оптимизации резерва.6S

3.3 Применение асимптотических свойств обобщенных процессов Кокса в задаче оптимизации параметров процесса риска.

3.3.1 Гарантированные оценки для начального капитала страховой компании

3.3.2 Гарантированные оценки для оптимальной ставки страховой премии

3.3.3 Гарантированные оценки для времени достижения желаемого значения резерва.

3.4 Применение асимптотических свойств обобщенных процессов Кокса в задаче оптимизации начального капитала страховой компании. Альтернативный подход.

3.5 Применение асимптотических свойств пуассоновских случайных сумм в одной задаче оптимизации резерва.

4 Применение асимптотических свойств обобщенных процессов Кокса в задаче оценки надежности модифицируемых систем. Случай непрерывного времени

4.1 Неоднородные рекуррентные модели изменения надежности модифицируемых систем. Непрерывное время.

4.1.1 Неоднородные экспоненциальные модели с непрерывным временем

4.1.2 Неоднородные логистические модели с непрерывным временем

4.1.3 Неоднородные гиперболические модели с непрерывным временем

 
Введение диссертация по математике, на тему "Оценки скорости сходимости обобщенных процессов кокса с ненулевым средним и некоторые их применения"

Объект исследования

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

Пусть N(t) - число скачков случайно блуждающей частицы за период времени [0, f], t > 0. При решении многих практических задач естественно считать, что моменты скачков образуют хаотический точечный случайный процесс на оси времени. Однако этот хаотический случайный процесс мол-сет не быть однородным в силу указанных выше представлений. Как известно, наиболее разумными стохастическими моделями неоднородных хаотических точечных процессов являются дважды стохастические пуассоновские процессы, иначе называемые процессами Кокса (см., например, (Bening and Korolev, 2002)). Определение процессов Кокса следует предварить кратким экскурсом в описание некоторых свойств однородных и неоднородных пуассоновских процессов. В полном объеме свойства однородных и неоднородных пуассоновских процессов описаны, например, в книгах (Kingman, 1964) и (Кингман, 2007), где также можно найти многочисленные примеры математических и прикладных задач, в которых возникают такие процессы.

Как известно (см., например, (Bening and Korolev, 2002)), пуассоновский процесс, характеризуемый тем обстоятельством, что интервалы времени между событиями потока стохастически независимы и имеют одинаковое показательное распределение, является наилучшей моделью однородного стохастического хаотического потока событий. Напомним, что привлекательность пуассоновского процесса в качестве модели однородного дискретного хаоса обусловлена как минимум двумя обстоятельствами. Во-первых, показательное распределение обладает максимальной дифференциальной энтропией среди всех абсолютно непрерывных распределений вероятностей, сосредоточенных на всей положительной полуоси и имеющих конечное математическое ожидание, а энтропия, как известно, является очень удобной численной характеристикой неопределенности. Во-вторых, точки (события) пуассоновского потока равномерно распределены на оси времени в том смысле, что для любого конечного интервала времени [г1!, to], t\ < to, условное совместное распределение точек пуассоновского потока, попавших в интервал [ii,^2]э при условии, что в этот интервал попало фиксированное число, скажем, п точек, совпадает с совместным распределением вариационного ряда, построенного по независимой однородной выборке объема п из равномерного на [t\,t2] распределения. Равномерное же распределение обладает максимальной дифференциальной энтропией среди всех абсолютно непрерывных распределений вероятностен, сосредоточенных на конечных интервалах, и очень хорошо соответствует общепринятому представлению об абсолютно непредсказуемой ограниченной случайной величине.

Напомним, что, если N\(t) - однородный пуассоновский процесс с некоторой интенсивностью А > 0, то параметр Л, называемый интенсивностью, имеет смысл среднего числа скачков процесса (количества событий наблюдаемого потока) в единицу времени. При этом

Р(Nx(t) = к) = к = 0,1,. (0.0.1)

ЕNx(t) = DN\(t) = At, t> 0, так что для любых s, t > 0

ENx(t + s)-ENx{t)

-= Л = const. s

Из соотношения (0.0.1) легко видеть, что для любого к — 0,1,.

P(Nx(t) = к) = P(Ni(\t) = к), t> 0, (0.0.2) то есть процессы N\(t) и N\(Xt) стохастически эквивалентны.

Однако, как уже было отмечено, в реальной практике хаотические потоки не бывают однородными. Для наглядности, перед тем как определить дважды стохастический пуассо-новский процесс, рассмотрим неоднородный пуассоновский поток событий N*(t), определяемый следующим образом. Пусть А(/) - некоторая положительная функция. Предположим, что функция X(t.) интегрируема, и обозначим t

A(t) = J Х(т)с1т, t > 0.0. о

Неоднородный пуассоновский процесс N*(t), t > 0, определяется как случайный процесс с независимыми приращениями, траектории которого стартуют из нуля (iV*(0) — 0) и

Р(N*(t) = k) = (Лffl , к = 0,1,. (0.0.3)

А;!

Тогда

EN*(t) = A(t), t> 0, так что

Цю ЕN-{t + .) - ЕЛГМ lim Aft + .) - AM = МЛ (>0 0

40 s sio s dt

Другими словами, при таком определении функция А(/) играет роль мгновенной интенсивности процесса в точке t. При этом по аналогии с (0.0.2) из (0.0.3) вытекает, что

P(N*(t) = k) = P(Ni(A(t)) = k), к — 0,1,., t> 0, то есть процессы N*(t) и Ni(A(t)) стохастически эквивалентны.

Наконец, рассмотрим предложенное Д. Коксом (Сох, 1955) понятие дважды стохастического пуассо7ювского процесса, являющееся естественным обобщением неоднородного пуассоновского процесса.

Пусть N\(t), t > 0, - однородный пуассоновский процесс с единичной интенсивностью, а Л(/.), { > 0, - независимый or N\(t) случайный процесс, обладающий следующими свойствами: Л(0) = О, Р(Л(£) < оо) = 1 для любого t > О, траектории Л(f) не убывают и непрерывны справа. Дважды стохастический пуассоновский процесс N(t), называемый теперь также процессом Кокса, определяется как суперпозиция N\(t) и A(t):

В этом случае будем говорить, что процесс Кокса N(£) управляется процессом A(t). В частности, если процесс Л(£) допускает представление в котором Л(t) - положительный случайный процесс с интегрируемыми траекториями, то А(/.) можно интерпретировать как мгновенную стохастическую интенсивность процесса N(t). Поэтому иногда процесс Л(4), управляющий процессом Кокса N(t), будет называться накопленной интенсивностью процесса N(t).

Несложно убедиться, что для процесса Кокса N(t), управляемого процессом A(t), справедливы соотношения

При каждом фиксированном t распределение случайной величины N(t) = N\ (A(i)) является смешанным пуассоновским. Класс смешанных пуассоновских распределений весьма широк и содержит, в частности,

• само пуассоновское распределение (ему соответствует вырожденное смешивающее распределение),

• отрицательное биномиальное распределение (ему соответств) ет смешивающее гамма-распределение) ,

N(t) = Ni(A(t)), t > 0.

0.0.4) о

Е N(t) = ЕЛ (t), D N(t) = ЕЛ (t) + D A(t).

• распределение Зихеля, котором}- соответствует обобщенное обратное нормальное смешивающее распределение,

• бета-пуассоновское распределение при смешивающем бета-распределении,

• распределение Делапорте со сдвинутым гамма-распределением в качестве смешивающего,

• обобщенное распределение Варинга,

• инфекционное распределение Неймана типа А (иначе называемое пуассон-пуассоновским).

Эти и другие примеры смешанных пуассоновских распределений подробно рассмотрены в книгах (Grandell, 1997), (Королев, Бенинг и Шоргин, 2007). В этих книгах, а также в (Grandell, 1976) и (Bening and Korolev. 2002) также можно найти детальное описание аналитических и асимптотических свойств процессов Кокса.

Актуальность темы исследования

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

Теория случайного суммирования довольно хорошо развита (см., например, монографии (Gut, 1988), (Круглов и Королев, 1990), (Gnedenko and Korolev, 1996), (Kalashnikov, 1997),

Bening and Korolev, 2002), (Silvestrov, 2002) и (Klebanov, 2003)), в которых содержится изложение как основ асимптотической теории случайного суммирования, так и описание специальных ее разделов. Не преуменьшая вклад остальных математиков, посвятивших свои работы исследованию тех или иных свойств случайных сумм, упомянем здесь лишь основополагающие работы Г. Роббинса (Robbins, 1948), который в схеме «нарастающих сумм» нашел достаточные условия сходимости распределений нецентрироваиных случайных сумм к масштабным, а неслучайно центрированных случайных сумм - к сдвиговым смесям нормальных законов; P. JI. Добрушина (Добрушин, 1955), обобщившего результаты Роббинса на произвольные случайно индексированные случайные последовательности при специальном выборе центрирующих и нормирующих констант; Б. В. Гнеденко, который, во-первых, совместно со своим учеником X. Фахимом доказал знаменитую теорему переноса, устанавливающую достаточные условия слабой сходимости случайных сумм независимых одинаково распределенных случайных слагаемых в схеме серий (Гнеденко, Фахим, 1969), и, во-вторых, поставил задачу об отыскании необходимых и достаточных условий упомянутой сходимости, первые шаги в решении которой были сделаны его учениками и прежде всего, А. В. Печинкиныы (Печинкин, 1973) (для случая сходимости к нормальному закону) и Д. Саасом (Саас, 1972), (Szasz, 1972) для общего случая; В. М. Круглова, в частности, нашедшего необходимые и достаточные условия слабой компактности случайных сумм (Круглов, 1998); В. Ю. Королева, который, во-первых, совместно с В. М. Кругловым нашел окончательное решение задачи Гнеденко-Сааса о необходимых и достаточных условиях слабой сходимости случайных сумм независимых одинаково распределенных случайных сумм (Korolev and Kruglov, 1998), и, во-вторых уточнил и обобщил упоминавшиеся выше результаты P. JI. Добрушина, указав необходимые и достаточные условиях слабой сходимости суперпозиции произвольных независимых случайных процессов (Korolev, 1996). Последние результаты, в частности, позволили установить необходимые и достаточные условия слабой сходимости распределений обобщенных процессов Кокса и обобщенных процессов риска как с детерминированными, так и со случайными премиями (Bening, Korolev, 2002), (Королев, Бенипг и Шоргин, 2007).

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

Пусть Xi,X'2,. - одинаково распределенные случайные величины. Предположим, что при каждом t > 0 случайные величины N(t), A'i. Х2,. независимы. Процесс

N(t) sa) = J2 xj> (°-0-5)

3=1 назовем обобщенным процессом Кокса (при этом для определенности считаем, что =

0). Процессы вида (0.0.5) играют очень важную роль во многих прикладных задачах. Достаточно сказать, что при Л(t) = Xt с Л > 0 процесс S(t) превращается в классический обобщенный пуассоновский процесс, широко используемый при моделировании многих явлений в физике, теории надежности, финансовой и актуарной деятельности, биологии и т. д. Большое число разнообразных прикладных задач, приводящих к обобщенным пуассо-повским процессам, описано в книгах (Gnedenko and Korolev, 1990) и (Bening and Ivorolev, 2002).

Приведем примеры некоторых задач, сводящихся к обобщенным процессам Кокса.

Пример 0.0.1. Как уже говорилось, пуассоновский точечный процесс в пространстве описывает абсолютно хаотическое распределение точек. Рассмотрим следующее вполне естественное обобщение этого свойства применительно к задаче о суммарной массе тел. Пусть р(х), х € Е3, - плотность (концентрация) каких-либо материальных частиц в пространстве. Предположим, что частицы - это материальные точки, распределенные в пространстве в соответствии с неоднородным трехмерным нуассоновским точечным процессом, так что если N(B) - число точек в некотором (измеримом) множестве В с R3, то

Р (N(B) = k) = ехр{гг}{-А(Д)}(Л(^, fc = 0,1,2,., где Л (В) = jBp(x)dx. Положим A(t) — J^<t p{x)dx. Пусть Xj - масса, сосредоточенная в j-ой точке. Тогда процесс S(t), определяемый соотношением (0.0.5), описывает суммарную массу, содержащуюся в сфере радиуса £. Подобные объекты представляют традиционный интерес в космологии.

Пример 0.0.2. Рассмотрим следующую модель одномерного броуновского движения -теплового движения некоторой частицы, испытывающей соударения с молекулами вещества, заполняющего среду, в которой движется частица. Пусть N(t) - число соударений частицы с молекулами до момента t. Если бы среда и свойства частицы не изменялись во времени, то было бы естественно предполагать, что N(t) является однородным пуассонов-ским процессом. Однако если допустить, что среда или свойства частицы изменяются во времени, то интенсивность столкновений перестает быть постоянной. Поэтому, обозначив перемещение частицы в результате j-го столкновения через Х:1. можно придти к выводу, что процесс S(t), определяемый соотношением (0.0.5), описывает координату частицы в момент t, если она испытывает броуновское движение в неоднородной (случайной) среде.

Всюду в дальнейшем будем считать, что у случайных величин имеется, по крайней мере, два первых момента. Обозначим EXj = a, DXi = <т2, 0 < а2 < оо. Как показано, например, в (Bening and Korolev, 2002) пли (Королев, Бенинг и Шоргин, 2007), даже в таких предположениях предельные распределения обобщенных процессов Кокса могут иметь произвольно тяжелые хвосты. При этом асимптотическое поведение обобщенных процессов Кокса с нулевым средним принципиально отличается от асимптотического поведения обобщенных процессов Кокса с ненулевым средним. Если в первом случае в качестве предельных законов в оговоренных выше моментных условиях выступают масштабные, то во втором - сдвиговые смеси нормальных законов.

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

Несложно убедиться, что если ЕЛ(£) < оо, то

Е S(t) = aEA(i) (0.0.6) и, если дополнительно Е(Л(£))2 < оо, то

D S(t) = (а2 + <т2)ЕЛ(£) + a2DA(£). (0.0.7)

К моменту начала работы над диссертацией были известны некоторые оценки точности приближения распределений обобщенных процессов Кокса с ненулевым средним сдвиговыми смесями нормальных законов (см., например, (Korolev, 2000), (Bening, Korolev, 2002) и (Королев, Соколов, 2008)). Но данные оценки оказываются справедливыми для довольно узкого класса распределений случайной величины, предельной для стандартизованного процесса накопленной интенсивности, довольно громоздки, содержат трудно вычисляемые характеристики и неудобны для анализа и применения. В диссертации получена оценка скорости сходимости в традиционных терминах, справедливая при минимальных ограничениях на моменты слагаемых и без каких бы то ни было ограничений на распределение случайной величины, предельной для накопленной интенсивности. Более того, данные оценки справедливы как для случая, когда у управляющего процесса существует дисперсия, так и для случая, когда у него дисперсии нет. Рассмотрены примеры, в рамках которых найдены численное значение констант, входящих в оценку, и тем самым показано, что с практической точки зрения полученные оценки оказываются вполне приемлемыми. Метод. использованный для получения указанных оценок, оказался довольно универсальным. Так, в диссертации, он практически без каких-либо изменений был применен для получения оценок скорости сходимости распределений экстремумов обобщенных процессов Кокса и обобщенных процессов риска.

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

• оптимизации резерва в управлении запасами при случайной интенсивности потока заявок,

• оптимизации параметров страховой деятельности при случайном характере интен-сивностей потоков страховых премий и страховых выплат,

• оптимизации параметров спекулятивной деятельности,

• прогнозирования надежности модифицируемых систем.

В задаче оптимизации параметров страховой деятельности используется неклассический стоимостной подход, в рамках которого оптимизируются (минимизируются) суммарные издержки страховой компании за некоторый фиксированный период времени. Стоимостной подход к задачам страхования тесно связан с задачами управления запасами и разрабатывался, в частности, в работах А. Кофмана (Кофман, 1966), Г. В. Ротарь (Ро-тарь, 1972а), Е. В. Булинской (Булинская, 2003), Т. Р. Катаева (Катаев, 2003). Решение рассматриваемой задачи сводится к отысканию корпя некоторого уравнения, связанного с распределением суммарного страхового требования. Поскольку точное решение этого уравнения возможно только при полностью известных распределениях страховых требований и их числа, чего на практике, вообще говоря, быть не может, с помощью оценок точности асимптотических аппроксимаций для распределений обобщенных процессов риска в диссертации приводятся двусторонние оценки для решения упомянутого уравнения. Рассмотрен рад альтернативных способов задания издержек.

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

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

Цель работы

Целью данной диссертации является

• изучение точности асимптотических аппроксимаций распределений обобщенных процессов Кокса с ненулевыми средними сдвиг овыми смесями нормальных законов;

• изучение точности асимптотических аппроксимаций распределений экстремумов обобщенных процессов Кокса с ненулевыми средними сдвиговыми смесями нормальных законов;

• изучение точности асимптотических аппроксимаций распределений обобщенных процессов риска сдвиговыми смесями нормальных законов;

• отыскание двусторонних оценок оптимального значения количества продукта в задаче минимизации издержек при управлении запасами со случайной интенсивностью потока заявок;

• решение задачи оптимизации параметров страховой деятельности при случайном характере интенсивностей потоков страховых премий и страховых выплат;

• решение задачи оптимизации некоторых параметров спекулятивной деятельности при случайной интенсивности потока клиентов;

• построение гарантированных двусторонних оценок для значения надежности модифицируемых технических или информационных систем при случайной интенсивности потока модификаций.

Методы исследования

В работе используются классические методы теории вероятностей. Базовым теоретическим результатом является теорема 1.1.2. При ее доказательстве используется метод, основанный на представлении конечномерных распределений обобщенных процессов Кокса в виде смесей обобщенных иуассоновских распределений, в которых смешивание производится по соответствующему конечномерному распределению управляющего процесса. Такой подход позволяет использовать известные результаты о точности нормальной аппроксимации для иуассоновских случайных сумм. Используемый метод оказывается достаточно универсальным и может быть успешно применен к построению оценок скорости сходимости распределений некоторых других процессов с ненулевым средним к сдвиговым смесям нормальных законов. В частности, в диссертации с помощью этого же метода указанные оценки построены для распределений экстремумов обобщенных процессов Кокса и обобщенных процессов риска. Затем теорема 1.1.2 сама становится базой методов решения оптимизационных задач, рассматриваемых в главах 2-4. Для решения этих задач также используются аналитические и асимптотические свойства пуассоновских случайных сумм.

Теоретическая и практическая значимость

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

Краткое содержание диссертации

 
Список источников диссертации и автореферата по математике, кандидата физико-математических наук, Артюхов, Сергей Владимирович, Москва

1. Артюхов С. В., Базюкина О. А., Королев В. Ю., Кудрявцев А. А. Модель оптимального ценообразования, основанная на процессах риска со случайными требованиями // Системы и средства информатики, 2005. Специальный выпуск. М.: ИПИРАН. С. 207-224.

2. Артюхов С. В. Применение асимптотических свойств иуассоновских случайных сумм в одной задаче оптимизации резерва // Системы и средства информатики, 2006. Специальный выпуск. М.: ИПИРАН. С. 238-248.

3. Артюхов С. В., Базюкина О. А., Королев В. Ю., Кудрявцев А. А. Шевцова И. Г. Об оптимизации спекулятивной прибыли на примере пункта обмена нам ют // Актуарий, 2008. т. С. 50-56.

4. Артюхов С. В., Жалыбина И. Я., Пузановский А. А. О методе оптимизации прибыли маркет-мейкера // Управление финансовыми рисками, 2008. Вып. 2. С. 35-48.

5. Артюхов С. В. Королев В. Ю. Неоднородные рекуррентные модели изменения надежности модифицируемых систем. Непрерывное время. // Информатика и ее применение, 2008. Т. 2. Вып. 4. С. 57-65.

6. Артюхов С. В. Оценки скорости сходимости распределений экстремумов обобщенных процессов Кокса с ненулевым средним к сдвиговым смесям нормальных законов //Информатика и ее применение, 2009. Т. 3. Вып. 1. С. 69-74.

7. Артюхов С. В., Базюкина О. А., Кудрявцев А. А. О методе расчета курсов покупки и продажи валюты обменным пунктом. // Материалы международной конференции студентов, аспирантов по фундаментальным наукам «Ломоносов-2005», 2005. Секция ВМиК. С. 6-8.

8. Артюхов С. В. Об определении оптимального количества рублевых и валютных средств в обменном пункте. // Материалы международной конференции студентов, аспирантов по фундаментальным наукам «Ломоносов-2006», 2006. Секция ВМиК. С. 6-8.

9. Бенине В. E., Королев В. Ю. Асимптотическое поведение обобщенных процессов риска // Обозрение промышленной и прикладной математики, 1998. Сер. «Финансовая и страховая математика». Т. 5, Вып. 1. С. 116-133.

10. Бенин г В. Е., Королев В. Ю. Введение в математическую теорию риска. — М.: MAKC-Пресс, 2000.

11. Бенин? В. Е., Королев В. Ю. Обобщенные процессы риска. — М.: МАКС-Пресс, 2000.

12. Бойков А. В. Стохастические модели капитала страховой компании и оценивание вероятности неразорения. // Дис. канд. физ.-матем. наук. — М.: Матем. ин-т им. В. А. Стеклова РАН, 2003.

13. Булинекая Е. В. Теория риска и перестрахование. Часть I. Упорядочивание рисков. — М.: изд-во мех.-матем. ф-та МГУ, 2001.

14. Булинекая Е. В. О стоимостном и надежностном подходе в теории запасов // Кибернетика и системный анализ, 2005. Т. 1. С. 61-81.

15. Буш Р., Мостеллер Ф. Стохастические модели обучаемости. — М.: ГИФМЛ, 1962.

16. Волков Л. И., Шишкевич А. М. Надежность летательных аппаратов. — М.:ВШ, 1975.

17. Волков Л. И. Управление эксплуатацией летательных комплексов. — М.: В1П, 1981.

18. Гнеденко Б. В., Фахим X. Об одной теореме переноса. // ДАН СССР, 1969. Т. 187. No 1. С. 15-17.

19. ГОСТ '27.002-89. Надежность в технике. Основные понятия. Термины и определения.- М., 1989.

20. Добрушин Р. Л. Лемма о пределе сложной случайной функции. // Успехи математических наук, 1955. Т. 10. No 2. С. 157-159.

21. Золотарев В. М. Современная теория суммирования независимых случайных величин. — М.: Наука, 1986.

22. КацДж., МакКормик Д. Л. Энциклопедия торговых стратегий. — М.: Альпина, 2006.

23. Катаев Т.Р. Асимптотические аппроксимации для распределений случайных сумм и некоторые их применения. // Дис. канд. физ.-матем. наук. — М.:.МГУ им. М. В. Ломоносова, 2003.

24. Кащеев Д. Е. (2001) Моделирование динамики финансовых временных рядов и оценивание производных финансовых инструментов, j j Дис. канд. физ.-ма1ем. наук. — Тверь: Тверской гос. университет, 2001.

25. Кингман Дж. Пуассоновские процессы. — М.: МЦНМО, 2007.

26. Королев В. Ю. Вероятностные модели. Введение в асимптотическую теорию случайного суммирования. — М.: МАКС Пресс, Москва, 1997.

27. Королев В. Ю. Прикладные задачи теории вероятностей: модели роста надежности модифицируемых систем. — М.: Диалог-МГУ, 1997.

28. Королев В. Ю. Смешанные гауссовские вероятностные модели реальных процессов.- М.: МАКС-Пресс, 2004.

29. Королев В. Ю. Вероятностно-статистический анализ хаотических процессов с помощью смешанных гауссовских моделей. — М.: ИПИРАН, 2007.

30. Королев В. Ю. Бенинг В. Е., Шоргин С. Я Математические основы теории риска. — М.: ФИЗМАТЛИТ, 2007.

31. Королев В. Ю., Соколов И. А. Основы математической теории надежности модифицируемых систем. — М.: ИПИРАН, 2006.

32. Королев В. К)., Соколов И. А. Математические модели неоднородных потоков экстремальных событий . — М.: ТОРУС Пресс, 2008.

33. Кофман А. Методы и модели исследования операций. — М.: Мир, 1966.

34. Круглое В. М., Королев В. Ю. Предельные теоремы для случайных сумм. — М.: МГУ, 1990.

35. Круглое В. М. Слабая компактность случайных сумм независимых случайных величин. // Теория вероятностей и ее применения, 1998. Т. 43. No. 2. С. 248-271.

36. Лукач Е. Характеристические функции. — М.: Наука, 1979.

37. Молодцов Д. А. Теория мягких множеств. — М.: Альпина, 2002.

38. Мэрфи Дж. Технический анализ фьючерсных рынков: теория и практика. — М.: Диаграмма, 1998.

39. Мэрфи Дж. Визуальный инвестор. Как определять тренды. — М.: Диаграмма, 2004.

40. Новинский П. В., Зограф И. А. Оценка погрешностей результатов измерений. — Л.: Энергоатомиздат, 1991.

41. Петраков Н. Я., Ротарь В. И. Фактор неопределенности и управление экономическими системами. — М.: Наука, 1985.

42. Печинкин А. В. О сходимости к нормальному закону сумм случайного числа случайных слагаемых. // Теория вероятностей и ее применения, 1973. Т. 19. No. 2. С. 380-382.

43. Пиндайк Р. С., Рабинфелъд Д.Л. Микроэкономика. 5-е международное издание — М.: Питер, 2002.

44. Ротаръ Г. В. Одна задача управления резервом. // Теория вероятностей и ее применения, 1972. Т. 17. Вып. 3. С. 597-599.

45. Ротаръ Г. В. Некоторые задачи планирования резерва . // Дис. канд. физ.-матем. наук. М.: ЦЭМИ, Москва, 1972.

46. Ротаръ Г. В. Об одной задаче управления резервами. // Эконом, мате.м. методы, 1976. Т. 12. Вып. 4. С. 733-739.

47. Саас Д. О классах предельных распределений для сумм случайного числа одинаково распределенных случайных величин. // Теория вероятностей и ее применения, 1972. Т. 17. No. 3. С. 424-439.

48. Темное Г. О. Математические модели риска и случайного притока взносов в страховании. // Дис. канд. физ.-матем. наук. — С.-Петербург: С.-Петербургский гос. архитектурно-строительный ун-т, 2004.

49. Шевцова И. Г. Уточнение структуры оценок скорости сходимости в центральной предельной теореме для сумм случайных величин. // Дис. канд. физ.-матем. наук. — М.: МГУ, 2006.

50. Шевцова И. Г. О точности нормальной аппроксимации для распределений пуассонов-ских случайных сумм. // Обозрение прикладной и промышленной математики, 2007. Т. 14. Вып. 1. С. 3-28.

51. Шевцова И. Г. Об абсолютной постоянной в неравенстве Берри-Эссеена. // Сб. статей молодых ученых факультета ВМиК МГУ. — М. Изд-во факультета ВМиК МГУ, 2008. Вып. 5. С. 101-110.

52. Шевцова И. Г. Уточнение абсолютной константы в классическом неравенстве Берри-Эссеена. // Статистические методы оценивания и проверки гипотез. Пермь: Изд-во Пермского государственного университета, 2008. Вып. 21. С. 159-168.

53. Ширяев А. Н. Вероятность. — М.: Наука, 1980.

54. Эмбрехтс П., Клюппелъберг. К. Некоторые аспекты страховой математики. // Теория вероятностей и ее применения, 1993. Т.38. Вып. 2. С.375-416.

55. Bachelier L. ТЬёопе de la speculation. // Ann. Ecole Norm. Sup., 1900. Vol. 17. P. 21-86. (reprinted in: P. H. Coothner (Ed.). The Random Character of Stock Market Prices. — Cambridge, Ma: MIT Press, 1967. P. 517-531).

56. H. Biilmann. Tendencies of development in risk theory. // Centennial Celebration of the Actuarial Profession in North America, 1989. Vol. 2. P. 499-521.

57. Bening V. E., Korolev V. Yu. Generalized Poisson Models and Their Applications in Insurance and Finance. —- Utrecht: VSP, 2002.

58. Bowers N. L. et. al. Actuarial Mathematics. — Itasca, IL: Society of Actuaries, 1986.

59. Cox D. R. Some statistical methods connected with series of events. // J. Roy. Statist. Soc., Ser. B. 1955. Vol. 17. P. 129-164.

60. Essen C.-G. A moment inequality with an application to the central limit theorem. // Shand. Artuarietidiskr, 1956. Vol. 39. P. 160-170.

61. Gnedenko В. V., Korolev V. Yu. Random Summation: Limit Theorems and Applications. Boca Raton, FL: CRC Press, 1996.

62. Grandell J. Doubly Stochastic Poisson Process. // Lecture Notes in Mathematics, 1976. Vol. 529. P. 229-234.

63. Grandell J. Aspects of Risk Theory. — Berlin-Heidelberg-New York: Springer, 1990.

64. Grandell J. Mixed Poisson Processes. — London: Chapman and Hall, 1997.

65. Gtit A. Stopped random walk. — New York: Springer, 1988.

66. Kalashnikov V. Geometric Sums: Bounds for Rare Events with Applications. — Dordrecht-Boston-London: Kluwer Academic Publishers, 1997.

67. Kingman J. F. C. On doubly stochastic Poisson processes // Cambridge Philos. Soc., 1964. Vol. 60. No. 4. P. 923-930.

68. Klebanov L. B. Heavy-taileB Distributions. — Prague: MATFYZ Press, 2003.

69. Korolev V. Yu. A general theorem on limit behavior of superpositions of independent random processes with applications to Cox processes, j j Journal of Mathematical Sciences, 1996. Vol. 81. No. 5. P. 2951-2956.

70. Korolev V. Yu. The asymptotic behavior of extrema of compound Cox process with non zero mean. // Journal of Mathematical Sciences, 2000. Vol. 99. P. 1273-12850.

71. Koroleu V. Yu., Kruglov V. M., A criterion of convergence of nonrandomly centered random suns of independent identically distributed random variables. // Journal of Mathematical Sciences, 1998. Vol. 89. No. 5. P. 1495-1506.

72. Neyman J. On a new class of "contagious" distribution, applicable in entomology and bacteriology. // Annals of Mathematical Statistics, 1939. Vol. 10. P. 35-57.

73. Robbins H. The asymptotic distribution of the sum of a random number of random variables. // Bull. Amer. Math. Soc., 1948. Vol. 54. No. 12. P. 1151-1161.

74. Romanowski M. Random errors in observations and the influence of modulation of their distribution. — Stuttgart: Verlag Konrad Wittwer, 1979.

75. Szasz D. Limit theorems for the distributions of the sums of a random number of random variables. // Annals of Mathematical Statistics, 1972. Vol. 43. No. 6. P. 1903-1913.

76. Silvestrov D. S. Limit Theorems for Randomly Stopped Stochastic Processes. Research Report 2002 1-4, Department of Mathematics and Physics, Vast,eras: Malardalen University, 2002.