автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Исследование математических моделей процессов страхования при нестационарных потоках страховых рисков

ЗМЕЕВ ОЛЕГ АЛЕКСЕЕВИЧ

На правах доописи

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

05.13.18 - «Математическое моделирование, численные методы и комплексы программ»

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

Томск-2005

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

Научный консультант доктор физико-математических наук,

профессор Терпугов А. Ф.

Официальные оппоненты: доктор физико-математических наук,

профессор Данилов Н. Н. доктор физико-математических наук, профессор Кошкин Г. М. доктор физико-математических наук, профессор Топчий В. А.

Ведущая организация: Институт вычислительного моделирова-

Защита состоится 19 мая 2005 года в 10.30 на заседании диссертационного совета Д 212.267.08 при Томском государственном университете по адресу: 634050, г. Томск, пр. Ленина, 36, корп. 2, ауд. 102.

Отзывы на автореферат (в двух экземплярах, заверенные печатью) просьба направлять по адресу: 634050, г. Томск, пр. Ленина, 36, Томский государственный университет, ученому секретарю ЮО. Буровой

С диссертацией можно ознакомиться в научной библиотеке Томского государственного университета.

Автореферат разослан 12 апреля 2005 г.

ния Сибирского отделения Российской академии наук (г. Красноярск)

Ученый секретарь диссертационного совета, доктор технических наук,доцент

Скворцов А. В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность и состояние проблемы

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

Основы современной актуарной математики были заложены работами Ф.Лундберга и Х.Крамера, в которых была предложена и исследована так называемая классическая модель процесса страхования, основанная на следующих предположениях: процесс поступления страховых премий в компанию считается детерминированным, за время / приращение капитала равно с/, где с -количество средств, поступивших в компанию за единицу времени; страховые выплаты - независимые, одинаково распределенные случайные величины; моменты страховых выплат образуют пуассоновский поток. Описание и различного рода исследования в рамках классической модели можно найти в монографиях [Э. Штрауба, Д. Кокса и В. Смита ,Y. H. Panjer и G. E. Willmont, J.Grandell, H.U. Gerber, H. У.Прабху], обзорах [В. И. Роторя и В.Е.Бенинга, П. Эмбрехтса и К. Клюппен-берта, В. Калашникова и Д. Константинидиса]. Основным достоинством классической модели является ее относительная простота, которая позволяет вычислить в явном виде такие характеристики, как вероятности разорения и выживания страховой компании.

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

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

возможность одновременного наступления страховых случаев. Минимизации вероятности разорения путем выбора инвестиционной стратегии посвящены работы Т. А. Белкиной, А. Г. Фроловой, С. В. Чекалиной, А. В. Бойкова и А. В. Мельникова, в которых предполагается возможность как безрисковых, так и рисковых инвестиций.

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

С другой стороны, кроме классических страховых компаний, в Российской Федерации на рынке страховых услуг существуют объекты, в деятельности которых активную роль играет государство. Построению и исследованию математических моделей таких объектов в последние годы посвящен ряд работ, в которых для исследования работы государственных фондов применяются различные методы теории массового обслуживания или идеи классических моделей страхования. Например, в работе Л Ф. Адашкина строится диффузионная аппроксимация для математической модели деятельности фонда социального страхования РФ. В работах А. А. Назарова, И. Р.Гарайшиной, Я. В. Галайко исследуются математические модели фондов пенсионного страхования.

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

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

Основной трудностью в рамках дискретно-событийного метода является правильное определение переменных состояния, необходимых для реализации моделирования с корректной последовательностью событий и получением интересующей статистики. В связи с этой сложностью необходимо отметить работу L. W. Schruben, который в 1993 году предложил метод представления событий с помощью графов, и работы Т. К. Som и R. G. Sargent, в рамках которых этот метод был значительно усовершенствован.

В последнее время в качестве возможного пути реализации имитационных моделей рассматривается объектно-ориентированное моделирование. Основные идеи применения объектно-ориентированного подхода к вопросу построения имитационных моделей можно найти, например, в работах J. A. Levasseur и D. W. Jones, S. D. Roberts. На сегодняшний день фактическим стандартом в области объектно-ориентированного анализа и проектирования (OOA&D) является использование Унифицированного Языка Моделирования (UML, Unified Modeling Language). UML непосредственно унифицирует известные методы Г. Буча, Д. Рамбо (ОМТ) и А. Джекобсона, при этом он обладает гораздо большими возможностями.

Другим фактором, оказывающим принципиальное влияние на современный процесс реализации программных комплексов, являются вопросы, связанные с так называемым повторным использованием кода. Поэтому, в последнее время возрастает значение каркасов приложений, с помощью которых объектно-ориентированные системы можно использовать повторно в максимальной степени. Обычно в рамках каркаса активно используются типовые решения проектирования (паттерны, шаблоны проектирования). Концепция типовых решений была предложена С. Alexander, S. Ishikawa, M. Silverstein, M. Jacobson, I. Fiksdahl-King, S. Angel, а в области разработки программного обеспечения идея применения паттернов подробно изложена в классической монографии Э. Гаммы, Р. Хемла, Р.Джонсона, Дж.Влиссидеса. В качестве примера работы последнего времени, посвященной тематике применения типовых решений при разработке информационных систем, приведем монографию М. Фаулера.

В представленной работе исследуются модели, учитывающие вышеизложенные факты, что, по мнению автора, и определяет ее актуальность. Исследованию деятельности страховой компании при нестационарных потоках страховых рисков посвящены первые три главы настоящей работы. Наиболее близкими по тематике к вопросам, рассматриваемым в этих главах, являются работы Д. Д. Ахмедовой, А. Ф. Терпугова, В. М. Каца, ,К. И. Лившица и А. А. Назарова, А. Ю. Голубина, С. А. Масяйкина, в которых рассматривается влияние расходов на рекламу в деятельности страховой компании, и исследуется конкурентное взаимодействие страховых компаний на общем страховом поле. В четвертой главе настоящей работы предлагается и исследуется математическая модель Фонда социального страхования Российской Федерации, основанная на адаптации классической модели страхования с учетом особенностей деятельности фонда. И, наконец, последние две главы настоящей работы посвящены разработке каркаса приложений имитационного моделирования в рамках дискретно-событийного метода. Применимость каркаса исследуется на основе построения программного комплекса, предназначенного для моделирования работы страховых компаний, рассмотренных в теоретической части диссертации.

Цель работы

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

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

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

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

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

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

6. Построить и изучить математическую модель фонда социального страхования.

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

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

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

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

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

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

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

Достоверность результатов и методика исследований

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

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

На защиту выносятся

Основные научные результаты, полученные автором и выносимые на защиту, состоят в следующем:

1. Разработана математическая модель страховой компании в виде двумерного

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

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

3. Рассмотрено конкурентное взаимодействие двух страховых компаний на общем рынке страховых услуг и построено переговорное множество (множество Парето) для такого взаимодействия.

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

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

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

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

Апробация работы

Основные положения диссертации и отдельные её результаты докладывались и обсуждались на Межрегиональной научно-практической конференции «Наука и образование: пути интеграции» (Анжеро-Судженск, 1998), Научно-теоретической конференции «Образование и наука на пороге третьего тысячелетия» (Барнаул, 1999), Межрегиональной научно-практической конференции «Качество образования и наука» (Анжеро-Судженск, 1999), Межрегиональной научно-методической конференции «Повышение эффективности научных исследований и совершенствование учебного процесса» (Анжеро-Судженск, 2000), Всероссийской научно-практической конференции «Математическое моделирование экономических систем и процессов» (Чебоксары, 2000), II научно-практической конференции «Наука и образование» (Белово, 2000), Международной научно-практической конференции «Методы и алгоритмы прикладной математики в технике, медицине и экономике» (Новочеркасск, 2001), III международной научно-теоретической конференции «Образование и наука в третьем тысячелетии» (Барнаул, 2001), Всероссийской научно-практической конференции «Новые технологии и комплексные решения: наука, образование производство» (Анжеро-Судженск, 2001), Международной научно-методической конференции «Новые информационные технологии в университетском образовании» (Кемерово, 2002), Межрегиональной VI научно-практической конференции «Научное творчество молодежи» (Анжеро-Судженск, 2002), II международной конференции «Проблемы актуарной и финансовой математики» (Минск, 2002), I, II, III Всероссийских конференциях по финансово-актуарной математике и смежным вопросам (Красноярск, 2002, 2003, 2004), IV

Всероссийской конференции с международным участием «Новые информационные технологии в исследовании сложных структур» (Томск, 2002), Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование» (Анжеро-Судженск, 2002), Всероссийской научно-практической конференции «Наука и практика: диалоги нового века» (Анжеро-Судженск, 2003), VIII Всероссийской научно-практической конференции «Научное творчество молодежи» (Анжеро-Судженск, 2004), V Всероссийском симпозиуме по прикладной и промышленной математике (Кисловодск, 2004), 8th Korea - Russia International Symposium on Science and Technology KORUS 2004 (Tomsk, 2004).

Публикации и личный вклад автора

По теме диссертации опубликовано 68 печатных работ, в том числе 20 - в журналах, рекомендованных ВАК для опубликования основных научных результатов докторских диссертаций, 19 работ опубликовано без соавторов.

Диссертационная работа и все результаты, лежащие в ее основе, выполнена при непосредственном участии автора на всех этапах. Ему в большинстве случаев принадлежит постановка задачи исследования, теоретическое описание и анализ полученных результатов. Исключение составляют работы, выполненные автором под руководством А. Ф. Терпугова, монография в соавторстве с Е. В. Глуховой и К. И. Лившицем, в которой автору принадлежит четвертая глава, объем которой составляет третью часть работы, и статьи, опубликованные совместно с Д.Д.Ахмедовой и А.Н.Моисеевым, в которых, по мнению автора, разделить научный вклад не представляется возможным. Остальные работы опубликованы в соавторстве с аспирантами и студентами, которые в различное время работали под руководством автора.

Структура и объем работы

Диссертация состоит из введения, шести глав и заключения, изложенных на 331 странице, включая 64 рисунка, и списка литературы, который содержит 217 наименований. Условно предлагаемую диссертационную работу можно разбить на три части. В первой части, в которую входят главы I—Ш, предлагаются и исследуются математические модели страховой компании в виде двумерного случайного процесса, компонентами которого являются капитал компании и число застрахованных рисков, при различных предположениях относительно числа застрахованных рисков, капитала компании, влияния рекламы и т.д.

В первой главе рассмотрены математические модели страховых компаний в предположениях, что поток входящих рисков является стационарным марковским потоком. Параграф 1.1 посвящен исследованию модели, в которой интенсивность входящего потока линейно зависит от числа уже имеющихся рисков, а страховое поле считается неограниченным. В 1.1.1 относительно рассматриваемой модели делаются следующие предположения: будем описывать состояние страховой компании в момент времени t двумерным случайным вектором {k(t\ S(t)}, где k(t) - число рисков, застрахованных компанией, a S(t) - ее капитал в момент времени t. Изменения капитала и числа застрахованных рисков происходят в следующих случаях:

1. Компания страхует новый риск. Будем предполагать, что поток приходящих рисков - это примитивный поток с параметром Первое слагаемое отражает поток рисков, которые клиенты страхуют в компании по независящим от нее обстоятельствам, а второе - тот факт, что среди людей, не застраховавших свои риски, распространяется информация о страховой компании, происходит неявная реклама компании. Вероятность того, что за время М компания застрахует новый риск, равна Каждый новый риск приносит компании страховую премию !;> размер которой является случайной величиной с функцией распределения ГДг) имоментами М{^} = а, М^2}=а2.

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

Будем считать, что взносы вносятся независимо друг от друга и поэтому за время в компанию поступит такой взнос с вероятностью

3. Страховое время некоторых рисков заканчивается. Будем считать, что каждый риск покидает компанию независимо от поведения других рисков с интенсивностью ц. Тогда за время Д< компанию покинет риск с вероятностью *цД/ + о(Д/).

4. Наконец, наступают страховые случаи. Будем считать, что с каждым клиентом может наступить страховой случай с интенсивностью цп и эти страховые случаи для различных рисков независимы. Тогда на интервале Д/ наступит страховой случай с вероятностью а компания при этом выплатит страховое возмещение в размере ц, которое является случайной величиной с функцией распределения /^(г) имоментами

Целью исследования является рассмотрение статистических характеристик процессов к($ и 5(0, т. е. поведение числа рисков и капитала компании.

В 1.1.2-1.1.7 проведены исследования модели в предположении, что процесс к{{) находится в стационарном режиме. Полученные результаты сформулированы в виде следующих утверждений и теорем:

Теорема 1.1. В предположениях 1) - 4) ^^J - финальные вероятности того, ЧТО ЧИСЛО рисков компании ™тзип / пттг^тт**ттсппт^ст ^ппггипшриистлш

ГТ^ + 'Р* '

к-1 /«О V +

(1.2)

а условие существования стационарного режима имеет вид

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

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

доказана теорема 1.2.

Теорема 1.2. Если

1. существует конечный предел

2. функция Р(х,е) дважды дифференцируема по х,

то функция Р(х) удовлетворяет уравнению Фоккера-Планка вида

При этом константа К = 1/(1 - ф).

Таким образом, показано, что при р-*оо число рисков является

асимптотически нормальной величиной с и

асимптотически нормальной

_ Р

М (и-РГ

Теорема 1.3. В стационарном режиме функция корреляции процесса к($ определяется выражением

йо*(т)=Л4(т)-М2{ф-Н-ехр[-(ц-р4) |т|}. (1.8)

(1-ФГ

Теорема 1.4. В стационарном для числа рисков режиме математическое ожидание капитала М {5,} определяется выражением

(1.14)

где Sa капитал компании в момент I.

Таким образом, условие возрастания среднего капитала имеет вид Цв + Х^с-Цп6>0, если Цв + А^С-Ц^^О, компания разоряется. Нетрудно

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

1 Теорема 1.5. В стационарном для числа рисков режиме дисперсия капитала

О {5,} определяется выражением

(1.17)

Из равенства (1.17) следует, что при больших I дисперсия капитала, так же как и его среднее, растет пропорционально времени I.

Теорема 1.6. Функция корреляции соу{.У(Д(/+/0)} процессов и в стационарном режиме для числа рисков имеет вид

Хц

Теорема 1.7. Функция корреляции процесса Я, ЯоЛМг) в стационарном режиме для числа рисков определяется выражением

cov{S„*(r+i0)} = K-^T.(1-21)

H-Ps

В 1.1.8 получены характеристики k(t) в нестационарном режиме. Теорема 1.8. В нестационарном режиме среднее число рисков left), дисперсия числа рисков Dk(t) и функция корреляции числа рисков Rqkihh) определяются выражениями

ехр[-(ц-р5)/},

+--^-lexpf- (и - Р, >Xl - expf- (н - р, >}),

(1.25)

(1.26) (1.27)

где /0 число рисков, застрахованных а компании в момент /0.

В 1.1.9 рассмотрено поведение капитала страховой компании в нестационарном режиме для числа рисков:

Теорема 1.9. В нестационарном режиме для числа рисков режиме математическое ожидание капитала М{5() и дисперсия капитала £>{•?,}

определяются выражениями

= + + - +

(1.30)

(X + Ii)

Среднее значение капитала компании в нестационарном режиме при больших t растет пропорционально t, аналогично стационарному случаю. Дисперсия капитала при больших / растет, как t.

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

В параграфе 12 рассмотрен случай, когда потенциальный рынок страховых услуг (страховое поле), на котором действует компания, ограничен, т.е. число рисков, которые может застраховать компания, ограничено. Пусть N -максимально возможное число рисков, за бесконечно малый промежуток времени & каждый из N потенциальных рисков может застраховаться с вероятностью ХД/ + о(Д<). Риск не может быть застрахован повторно, пока не истечет срок текущего договора. Величину £(/), по-прежнему, будем считать числом застрахованных в компании рисков. Тогда суммарный поток поступления рисков в компанию будет примитивным потоком с интенсивностью

В отличие от модели, рассмотренной в предыдущем параграфе, вероятность поступления в компанию нового риска за время Д( в этом случае равна По-прежнему в этой ситуации компания получает страховую премию 5» размер которой является случайной величиной с функцией распределения ^(г), М$ = а и Остальные вероятности, связанные с

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

Целью исследования, как и ранее, является рассмотрение статистических характеристик процессов и 8(1), т. е. поведения числа застраховавшихся рисков и капитала компании в стационарном и нестационарном режимах.

В 1.2.2 доказаны следующие утверждения:

Теорема 1.10. Для модели с ограниченным страховым полем Яу - финальные вероятности того, что число рисков компании равно определяются

соотношениями

где

*о = (1 + р)~". _

я, =СМ1 + р)~*>ям V1 = 1> N ,

(1.37)

Теорема 1.11. При N-¥"0 в стационарном режиме распределение числа рисков является асимптотически нормальным с

Далее в 12.3 - 1.2.7 исследована деятельность компании в предположении, что число застрахованных рисков стационарно. Результаты сформулированы в виде следующих теорем:

Теорема 1.12. В стационарном режиме функция корреляции процесса к^) определяется выпажением

^(тК(т)-М2{лЬ-%ехр(-Мт|) (1.42)

(|+Р)

Теорема 1.13. В стационарном для числа рисков режиме математическое ожидание капитала М {5,} определяется выражением

Лф^о+^^-^+ца),

(1.48)

где капитал компании в момент t¡r

Таким образом, условие возрастания среднего капитала и в этом случае имеет вид компания разоряется.

Теорема 1.14. В стационарном для числа рисков режиме дисперсия капитала /З^,} определяется выражением

(1.51)

Теорема 1.15. Функция корреляции соу{.У(Д(?+ /0)} процессов и в стационарном режиме для числа рисков имеет вид

соу{5,Д(/+/0)}=(сХс-а).-^)^-^[1-ехр{-(Х+и>}]. (1.55)

ТеореМа 1.16. Функция корреляции процесса ^о, (^>'2) в стационарном режиме для числа рисков определяется выражением

ЛГХц

Х+д

+ (вЯ.-4цч ♦^Т^гИ^' (1.56)

(Х+ц)

Аналогично параграфу 1.1.9 в 1.2.9 получены характеристики £(/) в нестационарном режиме. В результате доказана следующая теорема:

Теорема 1.17. В нестационарном режиме среднее число рисков дисперсия числа рисков £>*(/) и функция корреляции числа рисков ЛоД/^г) определяются выражениями

ю.

Х+ц

+ (<!--

да'

Х+ц,

ехр{-(Я+ц>}.

(1.59)

(1.60) (1.61)

где /0 - число рисков, застрахованных в компании в момент /0.

В 1.2.10 рассмотрено поведение капитала страховой компании 5(/) в нестационарном режиме для числа рисков:

Теорема 1.18. В нестационарном режиме для числа рисков режиме математическое ожидание капитала М{5,} и дисперсия капитала определяются выражениями

- - е*р(-

(1.65)

~(Ил62 + \С2)+ (Ь - Цл6 + X N1 (. ЗМЛ А7.(ла + ц)

Среднее значение капитала компании в нестационарном режиме растет пропорционально t, аналогично стационарному случаю. Дисперсия капитала при больших / растет, как I.

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

В параграфе 1.3. в отличие от моделей, рассмотренных в 1.1. и 1.2, делается дополнительное предположение относительно поведения капитала компании. Считается, что за время капитал компании увеличится на где г -

ставка банковского процента.

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

Теорема 1.19. Для моделей с учетом банковского процента в стационарном для числа рисков режиме математическое ожидание капитала компании определяется выражением

где 80 - стартовый капитал компании.

Следствие. Математическое ожидание капитала компании в стационарном режиме для модели с неограниченным страховым полем равно

(186)

Для модели с ограниченным страховым полем

= (1,7)

Заметим, что при г = 0 выражения (1.86) и (1.87) превращаются в соотношения (1.14) и (1.48) соответственно.

Теорема 1.20. Дисперсия капитала страховой компании для модели с работающим капиталом в стационарном для числа рисков режиме определяется выражениями

'H-PJ

-1 —+

+ ^ (Г -fr + -2^ИХ+И>)' ^Ц,.2„ .2

для модели с неограниченным страховым полем и

z>{s(0}=[;

т. Vrt-i

Шь + +un6, - Ял, )—— I--1+

п А. + Ц/ 2г

для модели с ограниченным страховым полем, где S0 - стартовый капитал компании.

Наконец, функция корреляции капитала компании определяется выражением R(tut2Ha2+b2

+ ble <e .D{k}---j—j-,

где bx, öj, b2, M{k}, £){i} и 9, параметры соответствующих моделей.

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

1. С вероятностью А,,Д* + о(Д/) в /-ую компанию поступает новый риск. Клиенты, страхующие риски, вносят страховые премии, являющиеся случайными величинами с математическим ожиданием а,. Риски поступают независимо друг от друга.

2. Время страхования некоторых рисков заканчивается, и клиенты покидают компанию. Считается, что каждый клиент уходит из компании с вероятностью цД/ + о(д/) независимо от другого.

3. С вероятностью Х.^Д/+о(Д/) каждый из клиентов компании уплачивает

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

4. Наконец, с вероятностью щ,Д/+о(Д/) у каждого из клиентов компании

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

dt \ ц В дальнейшем комбинацию (Л^ -¿цп)/ц будем обозначать через 5.

(1.102)

В 1.4.2 описана модель взаимодействия страховых компаний. Считается, что интенсивности прихода клиентов в ту или иную компанию X, зависят от величин страховых взносов at и а2, то есть = \1(а],а2). Далее, довольно естественно считать, что эти вероятности обратно пропорциональны величине страховых взносов, т. е. Х{/к2 = аг/щ. Если

/(а„а2)=Х.1(о„а2)+А.2(а|,а2).

Тогда

=-^_/(0„02), 12 —/(а„а2). (1.103)

Д) + ÖJ Ч\ + а2

Вид зависимости f(aha2) от а, и а2 определить достаточно сложно. В настоящей работе мы будем считать, что f{a{,a2) зависит лишь от некоторого параметра р, который, в свою очередь, зависит от в, и аъ так что р = р(а{, а2). Относительно вида функции /(р) естественно выдвинуть следующие предположения:

1./(0)<-Ко;

2. f{p) монотонно убывает с ростом р\

3. lim f(p)= 0;

4. lim pf(p)=0. р->+®

Что касается самой зависимости р(а{, а2), то к ней можно предъявить следующие достаточно естественные требования. Пусть а2 фиксировано, тогда:

1. Если (jj -wo, то р должно равняться а2, так как при этих условиях вся динамика рынка страхования будет определяться именно этой величиной, клиенты просто будут игнорировать первую компанию. Таким образом, сформулируем следующее утверждение:

lim р(а,,а2)=а2.

2. С уменьшением ах р(а,, а2) также должно монотонно убывать, так как у желающих застраховаться в этом случае появляется возможность выбора. Поэтому должно выполняться следующее условие:

За,

3. При а, -»О p(ait а2) также должна стремиться к нулю, так как в этом случае страховаться будут все потенциальные клиенты, следовательно,

lim р(аиа2)= 0.

в|-<0

4. Наконец, р(а1,а2) должно быть симметричной функцией относительно переменных а{ и а2,т. е. р{ах, я2)= р(а2,а,).

По-видимому, достаточно правдоподобной является следующая зависимость р от а, и а2:

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

Р =

-—Eifa

Таким образом, выражения (1.103) можно записать в следующем виде:

(1.105)

тогда

(1.106)

(1.107)

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

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

В параграфе 2.1. описывается модель страховой компании. Предполагается,

что

1. Компания страхует новый риск. Будем предполагать, что поток приходящих рисков - это пуассоиовский поток с переменной интенсивностью с параметром Вероятность того, что за время компания застрахует новый риск, равна Х(/)Д?+о(Д<). Каждый новый риск приносит компании страховую премию размер которой является случайной величиной с функцией

распределения

и моментами

Л/{);} = *,(<), Л42}=*2(/).

Так как величина

первого страхового взноса - это прерогатива компании и может устанавливаться по её усмотрению, то будем считать, что а2С0 зависят от времени и являются

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

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

случайной величиной с функцией распределения ^(г) и моментами М{^} = С1 и

Будем считать, что взносы вносятся независимо друг от друга и поэтому за время Д/ в компанию поступит такой взнос с вероятностью *Х;Д/ + о(Д<).

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

интенсивностью д. Тогда за время А; компанию покинет риск с вероятностью к цДг + о(Д<).

4. Наконец, наступают страховые случаи. Будем считать, что с каждым клиентом может наступить страховой случай с интенсивностью цч и эти страховые случаи для различных рисков независимы. Тогда на интервале Д/ наступит страховой случай с вероятностью кцпМ + о(Д/), а компания при этом выплатит страховое возмещение в размере я, которое является случайной величиной с функцией распределения /^(г) и моментами Л/{г|} = А1( Л/{г|2]=62.

В параграфе 2.2. предполагается, что параметр потока - непрерывная функция или функция, имеющая разрывы только первого рода. В 2.2.1. определяются выражения для характеристик капитала и числа застрахованных рисков в случае неограниченного страхового поля.

Теорема 2.1. В предположениях 1)-4) математическое ожидание £,(?)=Л/дисперсия Дг(')и функция ковариации числа застрахованных рисков С^определяются выражениями

(г) = *('о>""('~'о) + '^(г^'Лк, (2.1)

(2.2) (2.3)

где к((0) число рисков в начальный момент времени.

Теорема 2.2. В предположениях 1) - 4) математическое ожидание £](/) и дисперсия капитала компании £>х(<) определяются выражениями

.о

где А, капитал компании в момент и.

о ^

- +(1- г'») )к{уУ"А>

(2.21)

о

(Ь,

(2.22)

о

где, для простоты = 0.

Теорема 2.3. В предположениях 1) - 4) функция ковариации капитала и числа застрахованных рисков Ся(<) определяется выражением

о

¡й^-н^Ат** ДОиУж

(2.35)

В 2.2.2. определены выражения для тех же что и в 2.2.1 характеристик капитала и числа рисков в случае ограниченного страхового поля (теоремы 2.4-2.6).

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

ki{t) = k{t0)e "» +wp.(z> » dz,

Iо

DA)- j

-2vif-i)-2 [X(x)&

dz,

(2.39)

(2.40)

чияСш) -МО-11- /»■(*)<&

С*(М2)=А(тт(м2)> , (2.41)

где 4(/0) число рисков в начальный момент времени.

Теорема 2.5. Для модели с ограниченным страховым полем математическое ожидание дисперсия капитала £>$(/) определяются выражениями:

I I I

+ (2.54)

0 0 о

где 50 - капитал в начальный момент времени,

^(<)= Л'}а2(«МмЦ + (б2Ип +с2Х;)}*1(я)Л| -

о о

- + -СЪТ ] \ск{2,и\1и

о очо

dz

0^0

dz +

-2j ¡C^uhiuMuhizjk^u

оЧо

dz,

(2.55)

где, для простоты, Sa = 0.

Теорема 2.6. Для модели с ограниченным страховым полем, функция ковариации капитала и числа застрахованных рисков Csk (t) определяется выражением

С« (') = - )/е4 (/,«V" - /С* «>», («>.(«v«. (2.71) о о

где Ct(t,u) выражается формулой (2.41).

В 2.2.3. предполагается, что л,(/) влияет на интенсивность потока входящих рисков, так как уменьшение a,(i) приводит к увеличению Щ и, наоборот, при увеличении a,(t) число желающих застраховаться уменьшается. Поэтому в общем случае предполагается, что X зависит от а^/) и от

Предположение 5). Для конкретизации будем считать, что

(2.74)

где Г(0) = 1, Г(+ со) = 0 и Г(а,) монотонно убывает с ростом д] (г), ^ о ) имеет смысл максимальной интенсивности потока входящих рисков.

Теорема 2.7. В предположениях 1) - 5) оптимальное управление а{(у), О^уйТ, доставляющее максимум среднему значению капитала в конце рассматриваемого промежутка, определяется из уравнения

р (Д1М) м

В параграфе 2.3 предполагается, что параметр входящего потока Х(г) -непрерывный в средне квадратичном стационарный случайный процесс.

Предположение 6). Пусть характеристики Х(/):

й(2) = Л0(г)+Г, НтЛо(г)=0.

Г-И1

В 2.3.1. и 2.3.2 определяются выражения для характеристик числа застрахованных рисков и капитала для неограниченного страхового поля.

Теорема 2.8. В предположениях 1) - 4), 6) математическое ожидание, дисперсия и функция ковариации числа застрахованных рисков определяются выражениями

(2.81) (2.81)

Г*

М о Ц

~ , ч I 1, ,1 ^'г-'Л1'1'^ , ,

2Ц

I Я^Уск*'— (2.82)

Теорема 2.9. В предположениях 1) - 4), 6) математическое ожидание и дисперсия £>$(/) капитала компании определяются выражениями

Мх {5, (/)} = 50+1 /в, (у)Л + ,

о ^

о I м

II1 _1 М Ц оо V ^ >

^ о ^

" / \ i

JJ/f0(«-v)a)(v)(l-e-il('-,'')t/^v+(l-e-,J')J ¡R^u-v^ydu dv , (2.86) .0 0 -°Ло

где S0 - значение капитала в начальный момент времени.

В 2.3.3 предполагается, что величина öj , обозначающая средние первоначальные взносы является функцией, зависящей от X0(f).

Предположение 7). Пусть а, =ai(X0(<)). Рассмотрим случай, когда интенсивность потока рисков >.(/), страхующихся в компании, связана с Х0(/) соотношением

Выражение (2.95) имеет следующий смысл: Х0(/) есть интенсивность потока рисков, желающих застраховаться, но реально каждый риск страхуется с вероятностью F(ö|) зависящей, естественно, от величины взноса alt таким образом, интенсивность потока рисков, которые будут застрахованы, равна X.(i) = Х,0 (i)F(a,). Желание управлять этой интенсивностью изменяя величину страхового взноса ах, находит свое отражение в том, что а, из константы превращается в функцию от X0(f), т.е. а, = а,(Х0(<)).

Теорема 2.11. В предположениях 1) — 4), 7) оптимальное управление ß[ = £7|(X0(i)), реализующее максимальный темп роста среднего значения капитала определяется как решение уравнения

F'(a,)=0. (2.100)

V ^ /

Теорема 2,12. Пусть Х0(м) = X, и X0(v)sX2 имеют совместную плотность распределения вероятностей

р(Х, ,Х2, w) = р0(Х, Д2, w)+ р(Х,)р(Х2), lim

тогда в предположениях 1) - 4), 7) асимптотическое поведение дисперсии капитала определяется выражением

Km^W-Wobfro WeifroMo +

» Л

1

+— И

jp(X0)X0F(fll(X0)>iX0

"г

00

+ 4

с.Ч-^п 7г

//{>-1^2«. М^! )№, Л2 )Л, А 2. (2.104)

о о

Теорема 2.13. В предположениях 1) - 4), 7) в случае экспоненциально распределенных первоначальных страховых взносов оптимальное управление страховыми взносами, реализующее максимум среднего темпа роста капитала при ограничении на флукхуацию этого капитала, то есть

1ип— ' ' => тах при Ьт—= О0, определяется из уравнения

(X,))+ Га, (X,(X,))+

+ /

И

м ; I ц )

= 0. (2.109)

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

В параграфе 3.1 дается описание модели. Считается, что в момент времени I состояние страховой компании характеризуется ее капиталом 5(/) и числом застрахованных рисков &(?). Далее будем считать, что компания отчисляет часть своего капитала на рекламную компанию, так что на интервале + на ей проведение выделяется а(*)$(*)Д/ денег, где а(г) - доля капитала, выделяемая на рекламу в единицу времени.

1. С вероятностью (Х0 + а(/)Х,.У(г))Дг+о(Д<) в компанию придёт новый риск и, страхуясь, внесет страховую плату которая является случайной величиной с функцией распределения /^(дс). Слагаемое а(<)Х,£(/) описывает приток новых рисков, обусловленный рекламой, параметр Х0 определяет интенсивность прихода рисков, не находящихся под воздействием рекламы; параметр ^ определяет эффективность рекламы.

2. С вероятностью цч£(()Д/ + о(Д/) произойдёт страховой случай и компания выплатит страховое возмещение в размере п, которое является случайной величиной с функцией распределения ^(я). Будем считать, что страховые случаи

независимы.

3. С вероятностью

поступит очередной взнос от

застраховавшихся в компании рисков, который является случайной величинои с функцией распределения /^(х). Предполагается, что взносы осуществляются

независимо друг от друга.

4. Наконец, с вероятностью (й(/)Д/ + о(Д/) прекращается действие страхового договора и один из рисков покинет компанию независимо от других рисков.

В параграфе 3.2 исследуется модель для случая неограниченного страхового поля. В 3.2.1 находятся выражения математического ожидания капитала и числа рисков компании:

(3.5)

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

В 3.2.2 определяются условия эффективности рекламы, и результаты сформулированы в виде следующей теоремы.

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

а+

И

>1, Х.|Я>1 + -

(3.8)

Расшифруем экономический смысл полученных условий эффективности рекламы. Рассматривая первое условие, перепишем его в виде

В последнем неравенстве слева

• первый сомножитель Х^аБ - это интенсивность дополнительного потока рисков, который возникает из-за действия рекламы,

• сомножитель а + \сХ^ -¿>цп)Д1 - это средний доход, который получает компания от каждого застрахованного риска,

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

Смысл второго критерия менее очевиден. Для того, чтобы его получить, представим второе неравенство из (3.8) в следующем виде:

а^оУххЗЧ^.

В этом случае левая часть неравенства аХ^аБ - это дополнительный доход, которые получает компания за счет первоначальных взносов от рисков, привлеченных рекламой, справа а5 - это, как и ранее, расходы на рекламу в единицу времени, а последнее слагаемое запишем в виде - это потери

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

Следствие. При выполнении первого условия (3.8) реклама будет эффективна при любом а > 0. Второе условие (3.8) может быть выполнено лишь при оно даёт ограничение на величину параметра а, определяющего

количество денег, выделяемых на рекламу. Имеем

И

а>-

V-1

(3.11)

То есть вкладывать деньги в рекламу стоит лишь тогда, когда параметр, соответствующий доле денег на рекламу, превышает некоторую величину.

В 3.2.3 решается следующая задача оптимального управления: пусть страховая компания начинает свою работу в момент времени t = 0 и её цель провести рекламную программу таким образом, чтобы к концу рассматриваемого промежутка времени [О, T] темп прироста её среднего капитал стал максимальным. Критерий оптимальности определяется следующим образом:

Максимум определяется по вцду ait), О ¿ИТ с учетом ограничений S0- капитал компании в начальный момент времени. Заметим, что это лишь одна из постановок задач на оптимизацию рекламной деятельности. Для решения задачи рассматриваемый промежуток времени [0,Г] делится на три: причем каждому из промежутков соответствует своё

значение a(f):

Система относительно средних капитала и числа рисков получена в 3.2.1 Система для ковариаций капитала и числа рисков выводится в данном параграфе и имеет следующий вид

dt

+2 (<лс

0"2>

а/

Решение системы находится для каждого из промежутков, где параметр а(/) изменяет свои значения. На первом промежутке решение имеет вид

Ц Ц

(3.35)

(3.36)

-2и<

М ц"1

+ С3,

(3.37)

I ^ 2ц

константы находятся из начальных условий Dy(0,0)= 0, Dt(0,0)=0, Cw(0,0)=0. Далее для каждого из решений, соответствующих второму и третьему промежутку константы находятся из условий сшивания на границе М {5(Г0,0)} = М {S(T0, а)}, Л/'{5(Г0,0)} = Л/'{5(Г0,а)} (для решения на втором промежутке) и

(Г„0) .D^a^D^.O) (на третьем). Поэтому в выражениях для среднего значения и дисперсии капитала в момент времени Т присутствуют константы, зависящие от моментов Т0 (включение рекламы) и Г, (выключение рекламы). При выполнении условий Т0 > О и Т0 <ТХ<Т условие максимума по по Г0 и Тх имеет вид

dU{T0,Tl)r

дТ0 SUM 5Г,

Далее задача сводится к отысканию с помощью численных методов такого момента Г,, при котором функция £/(0,7]) имеет максимум.

В параграфе 3.3 исследуется модель страховой компании с учетом отчислений на рекламу в случае ограниченного страхового поля. В 2.3.1 выводится и решается система относительно средних значений капитала и числа рисков аналогично случаю неограниченного страхового поля. В 2.3.2 рассматривается стационарный режим страховой компании. Находится точка покоя системы

= 0,

(3.69)

► со определяются условия, при которых эта точка является

в асимптотике при N ■ устойчивой.

В 3.3.2 определяются условия эффективности рекламы (теорема 3.2) Теорема 3.2. Для модели с учетом расходов на рекламу в случае ограниченного страхового поля реклама эффективна, если выполнено хотя бы одно из условий

(3.77)

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

В параграфе 4.2 строится диффузионная аппроксимация деятельности фонда. Основной характеристикой состояния фонда является его капитал 8(1) в момент времени /. С этим капиталом происходят следующие изменения:

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

2. Фонд выделяет часть своих средств на социальные программы. Мы будем считать, что эти средства также выделяются непрерывно во времени, однако скорость их выделения с (5) зависит от величины капитала £ в данный момент времени.

Величина с0-с*(5) в дальнейшем обозначена как с(8). Таким образом, с(5) есть скорость изменения капитала за счет детерминированных расходов и она зависит от величины капитала 8. Именно в наличии слагаемою с (5) и зависимости с(8) от 8 и заключается отличие данной модели от классической.

3. Происходят страховые выплаты, Будем считать, что поток страховых выплат является пуассоновским потоком постоянной интенсивности Я,Л,и сами страховые выплаты \ являются независимыми одинаково распределенными случайными величинами с функцией распределения и начальными

моментами

Кроме того, будем считать, что достижение порога S(t) = 0 не приводит к разорению фонда, и даже при S(i) < 0 он продолжает функционировать, только происходят задержки по страховым выплатам.

Заметим, что процесс S(t) - марковский процесс. Обозначим P{S, t)=Pr{SiS(t)sS+dS}/dS.

Теорема 4.1. Если

1. плотность распределения вероятностей P(S, /) дифференцируема по t,

2. произведение c(S)P(S, t) дифференцируемо по S,

I»

3. несобственный интеграл Jp(S + и, t)dF^(u) конечен,

о

то плотность распределения вероятностей P(S, t) удовлетворяет прямому уравнению Колмогорова следующего вида

^ = (4.1)

dt es 0J

Решение уравнения (4.1) в общем виде представляет определенные математические проблемы, но сделав ряд дополнительных предположений, его можно найти, воспользовавшись методом асимптотического анализа марковизируемых систем. Сформулируем эти дополнительные предположения, пусть процесс S(t) принимает достаточно большие значения, тогда выберем бесконечно малый параметр е > 0 таким образом, чтобы при в -»0 процесс s2S(t) не равен тождественно нулю.

C(S)= c,(e2s), te2 = т, Se2 = х(т)+Еу, i?(S, <)= п(у,т,г), (4.3)

Теорема 4.2. Если

1. функция с, (х) дифференцируема,

2. существует конечный предел НтП(у,т,Е) = П(у,т),

е-»0

3. функция П (у, т) дифференцируема по т и дважды по у, то функция П(у, т) удовлетворяет уравнению Фоккера-Планка вида

t) с т)} t ка2 д2П(у, т) (4 5)

dt ду 2 ду2

При этом функция дг(т) определяется как решение дифференциального уравнения

äc,(;c(t))-Xöi. (4.6)

at

Из полученного уравнения (4.5) следует, что процесс

у^ЖМ,

Е

в пределе при е-»0 по распределению сходится к частному случаю процесса Орнштейна-Уленбека у(т), удовлетворяющему стохастическому дифференциальному уравнению вида

^¡Xa2chv(x), (4.8)

где w(t) - стандартный винеровский процесс.

Из уравнений (4.6) и (4.8) вытекает следующее утверждение. Теорема 4,3. Если функция с, (г) дважды дифференцируема, то случайный процесс

г(тМт)+еу(т) (4.9)

удовлетворяет равенству вида

dz(x) = (с, (z) - Ал, + Ефиз2 dw{x) +^-R2dx,

(4.10)

где R2 = -yVMX 0<9^1.

Таким образом, показано, что процесс S(t) можно аппроксимировать решением стохастического дифференциального уравнения

dS(t) = (c(S)-a{k)dt + Ja^dw(t)- (4.13)

Так как фонд не стремится к неограниченному накоплению капитала, то существует стационарное (финальное) распределение вероятностей p(S) капитала S.

p(S) = С • ехр

агк J

(4.14)

где константа С находится из условия нормировки.

Далее в 4.2.3 находится плотность вероятностей величины капитала 5 при фиксированном релейном управлении. В 4.2.4 определяются вероятностные характеристики деятельности фонда при фиксированном управлении: вероятность неплатежеспособности фонда и вероятность выделения денег на социальные расходы - и предложена процедура, позволяющая определить параметры управления, которые обеспечивают заданные вероятностные характеристики работы фонда. В 4.2.5 и 4.2.6 аналогичная задача решается для релейно-гистерезисной процедуры управления капиталом.

В параграфе 4.3. получены точные формулы для плотности вероятностей капитала фонда в стационарном режиме и предложена процедура, позволяющая определить параметры управления, которые обеспечивают заданные вероятностные характеристики работы фонда, в предположении, что страховые выплаты распределены по экспоненциальному закону, когда поток страховых выплат является пуассоновским потоком постоянной интенсивности X, и сами страховые выплаты являются независимыми одинаково распределенными случайными величинами с экспоненциальным распределением р^(х) = ехр{-х/а}/а, х>0.

В параграфе 4.4. решена задача нахождения вероятностных характеристик капитала фонда при произвольном законе управления его капиталом, но в отсутствии гистерезиса при этом управлении.

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

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

В параграфе 5.2. рассмотрены элементы дискретно-событийной имитационной модели. В 5.2.1 приведено формальное описание метода моделирования, а в 5.2.2 выполнена функциональная декомпозиция элементов дискретно-событийной модели, проиллюстрированная соответствующей диаграммой состояний.

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

В параграфе 5.4. рассмотрены вопросы, связанные с проектированием уровня бизнес-логики в рамках архитектуры каркаса. В 5.4.1 приведено решение, позволяющее использовать диаграммы активности и диаграммы состояний иМЬ для документирования метода представления событий с помощью графов, предложенного ЬЖ 8сЬгиЬеп. Далее, в 5.4.2. рассмотрена задача организации моделей, имеющих сложную иерархическую структуру, с применением паттерна компоновщик. Для инициализации компоновщика используется типовое решение строитель. Приведены примеры организации объектов, формирующих модель, в рамках компоновщика. В 5.4.3 рассмотрена проблема взаимодействия с библиотечными пакетами, для решения которой предлагается использовать паттерны стратегия и адаптер. Наконец, в 5.4.4 описан порядок инициализации и разработки имитационной модели с помощью каркаса.

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

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

Публикации

Результаты диссертации опубликованы в 68 работах, основными из которых,

по мнению автора, являются следующие:

1. Глухова Е.В., Змеев ОА, Лившиц К.И. Математические модели страхования. Томск: Изд-во Том. ун-та, 2004.178 с.

2. Змеев ОА Модель функционирования страховой компании при конечном числе возможных клиентов.// Известия вузов. Физика. № 4. 1999, с. 34-39.

3. Змеев ОА Расчет характеристик времени разорения страховой компании для моделей с интенсивностью входного потока, зависящей от числа имеющихся рисков.//Известия вузов. Физика. №4.2000, с. 10-16.

4. Змеев О.А. Математические модели функционирования страховой компании с учетом банковского процента.// Известия вузов. Физика. №1.2001, с. 19-25.

5. Ахмедова Д.Д., Змеев ОА Оптимизация расходов на рекламу при деятельности страховой компании.// Известия вузов. Физика, № 6.2001, с. 3-7.

6. Змеев О.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний.// Известия вузов. Физика, № 2. 2002, с. 24-28.

7. Ахмедова Д.Д., Змеев ОА, Терпугов А.Ф. Оптимизация деятельности страховой компании с учетом расходов на рекламу.// Вестник Том. гос. ун-та. №275,2002, с. 181-184.

8. Змеев ОА, Лезарев А.В. Шаблон объектного проектирования для реализации функциональности процесса моделирования в имитационных моделях систем массового обслуживания.// Вестник Том. гос. ун-та. № 275,2002, с. 108-111.

9. Змеев ОА, Моисеев А.Н. Шаблон диаграммы компонентов информационной системы корпоративного уровня.// Вестник Том. гос. ун-та. № 275, 2002, с. 130-132.

10. Ахмедова Д.Д., Змеев О.А. Нахождение характеристик страховой компании с нестационарным входящим потоком и ограниченными числом рисков.// Вестник Том. гос. ун-та Приложение №1 (I), 2002. С. 3-8.

11. Змеев ОА, Моисеев А.Н., Новиков Д.В. К вопросу проектирования уровня хранения в виде ООРДБ.// Вестник Том. гос. ун-та. Приложение № 1 (II), 2002. С. 363-367.

12. Змеев ОА Математическая модель фонда социального страхования с детерминированными расходами на социальные программы (диффузионное приближение).// Известия вузов. Физика, № 3.2003, с. 83-87.

13. Змеев ОА Математическая модель фонда социального страхования со случайными расходами на социальные программы (диффузионное приближение).// Известия вузов. Физика, № 3.2003, с. 88-93.

14. Змеев О.А. Математическая модель деятельности фонда социального страхования при экспоненциальных страховых выплатах.// Вестник Том. гос. ун-та, 2003, №280. С. 130-135.

15. Змеев ОА, Моисеев А.Н. Сравнительный анализ некоторых методов O-R-преобразования.// Вестник Том. гос. ун-та, 2003, № 280. С. 263-271.

16. Вальц О.В., Змеев О.А. Диффузионная аппроксимация модели фонда социального страхования с релейно-гистерезисным управлением капитала // Известия вузов. Физика, 2004. № 2. - С. 26-31.

17. Змеев ОА Деятельность фонда социального страхования при релейно-гисте-резисном управлении капиталом.// Математическое моделирование, 2004, т.16,№2.-С.43-53.

18. Вальц О.В., Змеев ОА Исследование модели фонда социального страхования.// Обозрение прикладной и промышленной математики. Т 11. Вып. 2.2004.-С. 311-312.

19. Вальц О.В., Змеев ОА Математическая модель деятельности фонда социального страхования при экспоненциальных страховых выплатах и случайными расходами на социальные программы // Вестник Томского государственного университета. № 284,2004. С. 37-41.

20. Войтиков К.Ю., Змеев ОА, Моисеев А.Н., Якушев А.А. Архитектура надстраиваемых приложений клиент/сервер с обобщенным протоколом передачи данных // Вестник Томского государственного университета. № 284, 2004.С. 169-173.

21. Змеев О А, Приступа А.В. Разработка объектно-ориентированного программного комплекса имитационного моделирования систем массового обслуживания // Вестник Томского государственного университета. № 284, 2004. С. 174-176.

22. Змеев О.А. Модель функционирования страховой компании при интенсивности входящего потока, зависящего от числа клиентов // Математическое моделирование. Кибернетика. Информатика. Сб. статей. -Томск: Изд-во Том. ун-та, 1999. - С. 67-73.

23. Змеев О.А. Определение вероятности разорения страховой компании для модели с интенсивностью входящего потока, зависящей от числа клиентов // Статистическая обработка данных и управление в сложных системах. Сб. статей- Томск: Изд-во Том. ун-та, 1999. - Вып. 1. - С. 57-66.

24. Змеев О.А. Определение вероятности разорения страховой компании для модели с конечным числом возможных клиентов // Статистическая обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд~во Том. ун-та, 1999.-Вып. 1.-С. 66-75.

25. Змеев О А, Терпугов А.Ф. Модель страховой компании с ограничением на число клиентов с учетом банковского процента // Математическое моделирование экономических систем и процессов. Всероссийская научно-практиче-ская конференция. Тезисы докладов. Чебоксары. 2000, с 60-63.

26. Змеев ОА, Змеева Е.Е. Определение вероятностных характеристик времени разорения страховой компании при условии, что оно произойдет // Статистическая обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2000. - Вып. 2. - С. 70-78.

27. Змеев ОА, Терпугов А.Ф. Расчет характеристик времени разорения страховой компании для модели с конечным числом возможных рисков // Экономика, технология, предпринимательство. Сб. статей. Томск. 2000, с. 6067.

28. Змеев О.А., Терпугов А.Ф. Математическая модель функционирования страховой компании с учетом банковского процента.// Методы и алгоритмы прикладной математики в технике, медицине и экономике. Международная научно-практическая конференция. Часть 3. Тезисы докладов. Новочеркасск. 2001, с. 33-36.

29. Ахмедова Д.Д., Змеев ОА Модель функционирования страховой компании с учетом расходов на рекламу при ограниченном числе клиентов // Статистическая обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2001. - Вып. 3. - С. 3-14.

30. Змеев ОА Нестационарный режим для математической модели функционирования страховой компании с интенсивностью входного потока, зависящей от числа клиентов.// Статистическая обработка данных и управление в сложных системах. Сб. статей. - Томск: Изд-во Том. ун-та, 2001. - Вып. 3.-С. 45-53.

31. Войтиков К.Ю., Змеев О А, Моисеев А.Н. Реализация классификаторов на сервере приложений в трехзвенной архитектуре «клиент/сервер» // «Новые информационные технологии в университетском образовании». Тезисы международной научно-методической конференции. Кемерово, 2002, с. 120122.

32. Ахмедова Д.Д., Змеев ОА Исследование математической модели страховой компании при нестационарном потоке рисков // Математические методы в финансах и экономике. Материалы второй международной конференции «Проблемы актуарной и финансовой математики». Минск, 2002, с. 5-10.

33. Змеев О.А., Масяйкин С.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний // Математические методы в финансах и экономике. Материалы второй международной конференции «Проблемы актуарной и финансовой математики». Минск, 2002, с. 20-24.

34. Ахмедова Д.Д., Змеев О.А. Нахождение характеристик работы страховой компании с нестационарным входящим потоком клиентов // Труды I Всероссийской ФАМ конференции. Ч. II. Красноярск, 2002, с. 20-27.

35. Ахмедова Д.Д., Змеев О А Нахождение характеристик капитала страховой компании с нестационарным входящим потоком, параметр которого случайная функция.// Материалы Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование», Анжеро-Судженск, 2002, с 23-26.

36. Змеев ОА, Лезарев А.В. Функциональные требования для систем имитационного моделирования систем массового обслуживания.// Материалы Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование», Анжеро-Судженск, 2002, с 128-130.

37. Змеев ОА, Моисеев АН., Новиков Д.В. Решение проблемы переполнения «родительских» таблиц в ООРБД.// Материалы Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование», Анжеро-Судженск, 2002, с 130-131.

38. Ахмедова Д.Д., Змеев ОА Математическая модель функционирования страховой компании с входящими рисками в виде пуассоновского потока событий с переменной интенсивности.// Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2002. - Вып. 4. -С. 3-12.

39. Войтиков К.Ю., Змеев ОА, Моисеев А.Н. Объектный подход к проблеме проектирования подсистемы нормативно-справочной информации.//

Обработка данных и управление в сложных системах. Сб. статей- Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 13-20.

40. Змеев ОА Построение переговорного множества при конкурентном взаимодействии двух страховых компаний.// Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2002. - Вып. 4. -С. 33-47.

41. Змеев ОА, Моисеев А.Н., Новиков Д.В. Применение метаданных при проектировании уровня хранения в системах, использующих ООРБД // Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 48-53.

42. Войтиков К.Ю., Змеев ОА, Моисеев А.Н. К вопросу о реализации классификатора в объектно-ориентированной информационной системе // Вестник Кем. гос. ун-та № 2(10). Кемерово. 2002, с. 167-177.

43. Войтиков К.Ю., Змеев ОА, Моисеев А.Н. Основные функциональные требования к подсистеме «Брокер объектных запросов» в рамках унифицированного процесса разработки программного обеспечения.// Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2003. - Вып. 5. - С. 3-13.

44. Змеев О А Математическая модель деятельности фонда социального страхования с детерминированными расходами на социальные программы при релейно-гистерезисном управлении капиталом.// Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2003. -Вып. 5.-С. 42-56.

45. Змеев ОА, Моисеев АН., Якушев А.А. Каркас брокера распределенной объектной базы данных // Материалы Всероссийской научно-практической конференции «Наука и практика: диалоги нового века» Часть 3. Томск: Изд-во «Твердыня», 2003. С. 90-92.

46. Змеев ОА, Приступа А.В. Классификация коммерческих систем имитационного моделирования // Материалы Всероссийской научно-практической конференции «Наука и практика: диалоги нового века» Часть 3. Томск: Изд-во «Твердыня», 2003. С. 93-95

47. Ахмедова Д.Д., Змеев О.А. Оптимальное управление первоначальным страховым взносом при деятельности страховой компании // Труды II Всероссийской ФАМ конференции. Ч. II. Красноярск, 2002, с. 21-24.

48. Змеев О.А. Диффузионное приближение в математических моделях фонда социального страхования // Труды II Всероссийской ФАМ конференции. Ч. II. Красноярск, 2002, с. 80-85.

49. Змеев ОА, Приступа А.В. Применение паттернов проектирования при построении систем имитационного моделирования // Материалы VIII Всеросс. научн.-практ. конф. «Научное творчество молодежи». Томск: Изд-во ТГУ, 2004. С. 34-36.

50. Змеев ОА, Приступа А.В. Диаграммы состояния UML как способ представления графа событий имитационной модели системы массового обслуживания.// Обработка данных и управление в сложных системах. Вып. 6. - Томск: Изд-во ТГУ, 2004. - С. 76-81.

51. Pristupa А. V., Zmeyev О. A. 8th Korea - Russia International Symposium on Science and Technology KORUS 2004, Tomsk, RUSSIA, June 26 - July 3, 2004, Vol. 1. P. 141-144.

USO

ГГ.

19 МДЙ 2005 \ ïli

К v У

ВВЕДЕНИЕ

Глава I. Модели страховых компаний при марковском стационарном потоке входящих рисков

1.1. Марковская модель с неограниченным страховым полем

1.1.1. Описание модели

1.1.2. Распределение числа рисков в стационарном режиме

1.1.3. Функция корреляции числа рисков в стационарном режиме

1.1.4. Математическое ожидание капитала компании в стационарном режиме для числа рисков

1.1.5. Дисперсия капитала компании в стационарном режиме для числа рисков

1.1.6. Ковариация капитала компании и числа застрахованных рисков в стационарном режиме для числа рисков

1.1.7. Функция корреляции капитала компании в стационарном режиме для числа рисков

1.1.8. Среднее, дисперсия и функция корреляции для числа рисков в нестационарном режиме

1.1.9. Поведение капитала страховой компании в нестационарном режиме для числа рисков

1.2. Марковская модель с ограниченным страховым полем

1.2.1. Описание модели

1.2.2. Распределение числа рисков в стационарном режиме

1.2.3. Функция корреляции числа рисков в стационарном режиме

1.2.4. Математическое ожидание капитала компании в стационарном режиме для числа рисков

1.2.5. Дисперсия капитала компании в стационарном режиме для числа рисков

1.2.6. Ковариация капитала компании и числа застрахованных рисков в стационарном режиме

1.2.7. Функция корреляции капитала компании в стационарном режиме

1.2.8. Среднее, дисперсия и функция корреляции для числа рисков в нестационарном режиме

1.2.9. Поведение капитала страховой компании в нестационарном режиме

1.3. Марковские модели с учетом банковского процента

1.3.1. Дополнительные предположения

1.3.2. Математическое ожидание капитала компании при условии, что число рисков стационарно

1.3.3. Дисперсия капитала компании при условии, что число рисков стационарно

1.3.4. Функция корреляции капитала компании при условии, что число рисков стационарно

1.4. Конкурентное взаимодействие двух страховых компаний в рамках марковских моделей

1.4.1. Модель страховой компании

1.4.2. Модель взаимодействия двух компаний

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

1.4.4. Интервал изменения значений параметра

1.4.5. Построение переговорного множества для отрицательных значений 5 115 Резюме

Глава II. Модели страховых компаний при марковском нестационарном потоке входящих рисков

2.1. Описание модели страховой компании

2.2. Параметр входящего потока - детерминированная функция

2.2.1. Характеристики страховой компании при неограниченном числе страховых рисков

2.2.2. Характеристики страховой компании при ограниченном числе страховых рисков

2.2.3. Оптимальное управление средними страховыми взносами J4Q

2.3. Параметр входящего потока - случайная функция

2.3.1. Характеристики числа рисков страховой компании при неограниченном числе страховых рисков

2.3.2. Характеристики капитала страховой компании при неограниченном числе страховых рисков

2.3.3. Характеристики страховой компании при неограниченном числе страховых рисков в случае зависимости функции средних первоначальных 157 взносов от интенсивности входящего потока рисков

Резюме >

Глава III. Математическая модель и управление деятельностью страховой компании с учетом расходов на рекламу

3.1. Модель страховой компании

3.2. Исследование деятельности страховой компании при неограниченном числе страховых рисков

3.2.1 Динамика капитала и числа застрахованных рисков

3.2.2 Условия эффективности рекламы

3.2.3 Управление денежными средствами, вкладываемыми в рекламную программу страховой компании

3.3. Исследование деятельности страховой компании на ограниченном страховом поле

3.3.1. Поведение капитала и числа рисков

3.3.2. Исследование стационарного режима в страховой компании

3.3.3. Условия эффективности рекламы 218 Резюме

Глава IV. Исследование математической модели фонда социального страхования Российской Федерации

4.1. Описание объекта моделирования/

4.2. Диффузионноя аппроксимация деятельности фонда

4.2.1. Математическая модель деятельности фонда

4.2.2. Построение диффузионного приближения

4.2.3. Релейное управление капиталом

4.2.4. Определение параметров управления

4.2.5. Релейное гистерезисное управление капиталом

4.2.6. Определение параметров управления

4.3. Исследование математической модели деятельности фонда при релейном гистерезисном управлении капиталом и экспоненциально распределенных 235 страховых выплатах

4.3.1. Особенности математической модели деятельности фонда

4.3.2. Релейно-гистерезисное управление капиталом

4.3.3. Стационарная плотность вероятностей величины капитала

4.3.4. Вероятностные характеристики

4.3.5. Релейное управление

4.4. Исследование математической модели деятельности фонда при произвольном законе управления капиталом и экспоненциально 248 распределенных страховых выплатах

4.4.1. Особенности математической модели деятельности фонда

4.4.2. Плотность вероятностей величины капитала

4.4.3. Вероятностные характеристики работы фонда

4.4.4. Временные характеристики деятельности фонда

4.4.5. Релейное управление капиталом 259 Резюме

Глава V. Проектирование каркаса приложений имитационного моделирования смо дискретно-событийным методом

5.1. Языки и среды имитационного моделирования. Достоинства и недостатки

5.2. Дискретно-событийные имитационные модели

5.2.1. Формальное описание метода моделирования

5.2.2. Элементы дискретно-событийной модели и их организация

5.3. Проектирование базовой архитектуры каркаса

5.3.1. Основные идеи для архитектуры приложений, построенных на базе каркаса

5.3.2. Взаимодействие базовых классов на верхнем уровне архитектуры

5.3.3. Порядок инициализации приложения

5.4. Проектирование архитектуры на уровне предметной области

5.4.1. Диаграммы состояний и активности UML как способ представления графа событий дискретно-событийного метода

5.4.2. Реализация моделей, имеющих сложную иерархическую структуру

5.4.3. Организация взаимодействия с библиотечными пакетами

5.4.4. Реализация механизмов синхронизации и обработки событий

5.4.5. Порядок инициализации модели 290 Резюме

Глава VI. Реализация имитационных моделей

6.1 Диаграммы активности для моделей страховых компаний

6.2 Представление моделей в виде компоновщика

6.3 Программное обеспечение

6.3.1 Общая характеристика программы

6.3.2 Системные требования

6.3.3 Инсталляция

6.3.4 Запуск программы

6.3.5 Окончание работы с программой

6.3.6 Работа с моделями

6.4 Результаты моделирования 313 Резюме

АКТУАЛЬНОСТЬ РАБОТЫ

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

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

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

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

Считается, что первыми работами по математической теории страхования являются работы Ф. Лундберга и X. Крамера, в которых была предложена и исследована так называемая классическая модель процесса страхования [138, 169, 207]. Классическая модель страховой компании, благодаря ее относительной простоте, позволяет вычислить в явном виде вероятности разорения и выживания страховой компании, выработать рекомендации по определению необходимого начального капитала и назначению страховых премий. В то же время классическая модель не отражает многие черты деятельности страховых компаний в реальной жизни. Развитию и уточнению классической модели посвящено большое количество работ по математической теории страхования, однако, остается еще много проблем, требующих дополнительного исследования. К числу малоизученных можно отнести, например, проблемы:

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

• учет нестационарности и случайности потоков входящих рисков;

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

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

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

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

ЦЕЛЬ РАБОТЫ

Целью данной работы является:

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

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

3. Рассмотреть конкурентное взаимодействие двух страховых компаний на рынке страховых услуг.

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

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

6. Построить математическую модель фонда социального страхования.

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

СОСТОЯНИЕ ПРОБЛЕМЫ

Как уже было отмечено выше, считается, что первыми работами по математической теории страхования являются работы Ф. Лундберга и X. Крамера, в которых была предложена и исследована так называемая классическая модель процесса страхования. Описание и различного рода исследования в рамках классической модели можно найти, например, в монографиях Э. Штрауба [155], Д. Кокса и В. Смита [117], Y. Н. Panjer и G.E. Willmont [207], J. Grandell [186], H.U. Gerber [185], H. У. Прабху [136], обзорах В.И. Роторя и В.Е. Бенинга [138], П. Эмбрехтса и К. Клюппенберта [160], В.Калашникова и Д. Константинидиса [105]. Классическая модель страховой компании базируется на следующих предположениях [184]: процесс поступления страховых премий в компанию считается детерминированным, за . время t приращение капитала равно ct, где с - количество средств, поступивших в компанию за единицу времени; страховые выплаты - независимые, одинаково распределенные случайные величины; моменты страховых выплат образуют пуас-соновский поток. Таким образом, величина страховых требований, поступивших в компанию за время t, образует сложно-пуассоновский процесс [150]. Основным достоинством классической модели является ее относительная простота, которая позволяет вычислить в явном виде такие характеристики, как вероятности разорения и выживания страховой компании.

В публикациях, посвященных изучению классической модели, в основном исследуются вероятности разорения и выживания страховой компании и принципы выбора нагрузки страховой премии (нагрузки безопасности). Из работ последнего времени отметим, например, работы V.M. Malinovskii [200, 201], в которых рассматривается вероятность разорения на конечном интервале, В.Е. Бенинга и В.Ю. Королева [24], в которой рассматривается вероятность разорения при малой нагрузке страховой премии, J. Grandell [187, 188], в которой исследуются простые аппроксимации вероятности разорения. S. Asmussen [164] исследовал адаптивные процедуры оценки вероятности разорения. В работах О.П. Виноградова [35], Ю.Д. Григорьева и А.В. Куклина [56, 57] рассматриваются возможности построения верхних и нижних границ для вероятности разорения. В работе В.В. Калашникова и Г.Ш. Цициашвили [106] строятся оценки для вероятности разорения при наличии больших выплат, а в работе К.Г. Гунченко [59] - прямые методы оценки вероятности разорения для различных видов распределений суммарного риска. Возможности перестрахования больших рисков рассматриваются в работах Г.А. Медведева [129], Ю.Д. Григорьева [54, 55], Е.В. Глуховой и Е.В. Капустина [51], Е.В. Булин-ской [27], И.В. Черепановой [154] и Л.Д. Шона [156, 157].

В большинстве работ последнего времени рассматриваются более сложные модели, обобщающие классическую модель. В рамках этих работ, процесс поступления страховых премий в компанию также считается случайным процессом. Так, например, в работе К.И. Лившица [120] находятся вероятность разорения и условное время до разорения для случая, когда страховые премии, поступающие в компанию, образуют пуассоновский процесс. В работах М. А. Маталыцкого, Т. В. Романюк [127, 128] страховая компания рассматривается как некоторая система массового обслуживания В работах К.И. Лившица и Л. Ю. Сухотиной [121, 122] рассмотрены характеристики страховой компании при малой нагрузке страховой премии для пуассоновской модели и модель страховой компании с учетом сезонных изменений. В работах В.Е. Бенинга и В.Ю. Королева [22, 23] исследуется случай, когда моменты страховых выплат образуют процесс Кокса (дважды стохастический пуассоновский процесс). Неоднородный поток страховых выплат рассматривается в работе О.П. Виноградова [34].

Большое внимание уделяется также проблемам, связанным с возможностью страховой компании использовать имеющиеся в ее распоряжении свободные средства для получения дополнительной прибыли и уменьшения тем самым вероятности разорения. Е.В. Глуховой и Е.В. Капустиным [50] рассчитывались вероятности выживания страховой компании при размещении части средств на депозитных вкладах, а в [107] учитывается возможность одновременного наступления страховых случаев. Минимизации вероятности разорения путем выбора инвестиционной стратегии посвящены работы Т.А. Белкиной, А.Г. Фроловой, С.В. Чекалиной [20, 21], А.В. Бойкова и А.В. Мельникова [26], в которых предполагается возможность как безрисковых, так и рисковых инвестиций.

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

B.М. Каца, К.И. Лившица и А.А. Назарова [109-116, 195], А.Ю. Голубина [185],

C.А. Масяйкина [125, 126], в которых рассматривается влияние расходов на рекламу на деятельность страховой компании и исследуется конкурентное взаимодействие страховых компаний на общем страховом поле.

В ряде работ рассматриваются более сложные, по сравнению с предыдущими, модели. Например, Н. Schmidli [212] рассматривает возможность одновременного инвестирования и перестрахования. Ф. Еникеева и В. Калашников [60] и D.G. Коп-stantinides, Q.H Tang, G.Sh.Tsitsiashvilii [196, 217], исследовали модели риска с инфляцией. В работе П. Эмбрехтса [159] прослеживается связь между актуарным и финансовым подходами к расчету страховых премий. Применение методов теории чувствительности к задачам страхования и теории финансов рассматривается в работе R. Norberg [208]. Применению франшизы, которая может получить большое распространение в связи введением обязательного автомобильного страхования, посвящена работа Ю.Д. Григорьева и И.Ю. Хекало [58].

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

Построению и исследованию математических моделей таких объектов в последние годы посвящен ряд работ, в которых для исследования работы государственных фондов применяются различные методы теории массового обслуживания или идеи классических моделей страхования применяются с учетом особенности работы таких фондов. Например, в работе Л.Ф. Адашкина [2] строится диффузионная аппроксимация для математической модели деятельности фонда социального страхования РФ. В работах А.А. Назарова, И.Р. Гарайшиной, Я.В. Галайко [41, 43, 44, 46-48, 132] исследуются математические модели фондов пенсионного страхования.

В четвертой главе настоящей работы предлагается и исследуется математическая модель Фонда социального страхования Российской Федерации, основанная на адаптации классической модели страхования с учетом особенностей деятельности фонда. С математической точки зрения вопросы, рассмотренные в этой главе, сводятся к задаче управления так называемым процессом разорения. Описание и обзор основных результатов различного рода исследований по этой тематике можно найти, например, в седьмой главе монографии Л. Такача [145]. К публикациям последнего времени, посвященным исследованию процессов разорения, можно отнести работы А.Т. Семенова [140 - 143].

Последние две главы настоящей работы посвящены разработке и исследованию имитационных моделей процессов страхования, рассмотренных в предыдущих частях работы. По утверждению А. М. Лоу и В. Д. Кельтона [1], имитационное моделирование является одним из самых распространенных методов исследования операций и теории управления. Согласно исследованию Gupta [189] эта технология исследования сложных систем занимает второе место, после математического программирования.

В рамках настоящей работы для реализации имитационного моделирования использован метод дискретно-событийного моделирования [1, 166]. К сожалению, необходимо отметить, что в отечественной литературе этот термин не является устоявшимся. В качестве примера приведем только ряд учебных пособий последних лет, в которых приведено описание этого метода. Например, в учебнике В. И. Вар-фоломеева, С. В. Назарова [33] метод называется дискретно-стохастической моделью (стохастическим автоматом) или Р-схемой, в работе Ю. И. Рыжикова [139] -дискретной моделью, наконец, в учебном пособии Б. Г. Ослина [133] все задачи имитационного моделирования сводятся к указанной схеме.

Основной трудностью в рамках дискретно-событийного метода является правильное определение переменных состояния, необходимых для реализации моделирования с корректной последовательностью событий и получением интересующей статистики. В связи с этой сложностью необходимо отметить работу L.W. Schruben [213], который в 1983 году предложил метод представления событий с помощью графов, и работы Т.К. Som и R.G. Sargent [211, 216], в рамках которых этот метод был значительно усовершенствован. Альтернативный метод формального описания СМО предложен, например, в работе Н.Н. Лябах, М.А. Бутакова [123].

В последние годы значительно вырос интерес к распределенным технологиям имитационного моделирования. Вопросы, связанные с применением распределенных вычислений в рамках дискретно-событийного метода, можно проследить в работах Chandrasekaran U., Sheppard S. [173], J. Misra [203], R. M. Fujimoto [182, 183], D.M. Nicol [204], R.L. Bagrodia [165]. Перечислим основные направления, которые рассматриваются в рамках распределенной реализации метода. Например, идеи, связанные с реализацией многопроцессорной обработки имитационной модели, рассмотрены в работе S. Sheppard et al. [216], в работе J.C. Comfort [177] рассматривается обработка списка событий по принципу «главный-подчиненный». "Другим подходом к распределению имитационной модели между различными процессорами является декомпозиция исходной модели на подмодели, выполнение которых передается различным вычислительным ресурсам. Этот метод распределенного моделирования был изначально разработан К. М. Chandy и J. Misra [174 -176], в более поздних работах J. Misra [203] рассмотрены проблемы технической реализации этого подхода. Альтернативная концепция виртуального времени, связанная с распределением подмоделей между параллельными вычислительными процессами, реализованная через механизм изменения шкалы времени рассмотрена у D.R. Jefferson [192].

Специальные исследования по эффективности распределенного моделирования описаны в работах S. Lavenberg, R. Muntz, В. Samadi [197], J.C. Comfort [177]. В работе P. Heidelberger [190] рассматривается влияние распределенного моделирования на статистическую эффективность со смешанными результатами. Решение вопроса о декомпозиции имитационной модели для распределения между параллельными вычислительными процессами можно найти в монографии М. S. Shanker, R. Padman, W. D. Kelton [214]. Анализ конкретных приложений для распределенного и параллельного моделирования приведен в работах D.M. Nicol, М.М. Johnson, A.S. Yoshimura [205] и C.D. Carathers, В. Topol, R. М. Fujimoto, V. Sunderam [172].

В связи с быстрым развитием технологий World Wide Web необходимо отметить ряд работ последнего времени, связанные с исследованиями самых разнообразных возможностей применения этих технологий для расширения возможностей моделирования. Так, например, в работах P. A. Fishwick [179, 180] исследуется широкий спектр вопросов, связанных с использованием механизмов клиент-сервер для увеличения производительности, распространения имитационных моделей и результатов их выполнения. Общий обзор подходов к моделированию, основанный на использовании веб-узлов и базирующийся на характерные примеры, дан в работе P. Lorenz, Н. Dorwarth, К. С. Ritter, T.J. Schriber [199].

В последнее время в качестве альтернативы или возможного пути реализации для распределенного подхода рассматривается объектно-ориентированное моделирование. Фактически объектно-ориентированное моделирование берет свое начало от объектно-ориентированного языка SIMULA, который появился в начале 60-х годов прошлого века. В качестве источников, в которых рассмотрены вопросы объектно-ориентированного моделирования, отметим J.A. Levasseur [198] и D.W. Jones, S. D. Roberts [193].

В рамках настоящей работы разработка объектно-ориентированных имитационных моделей осуществляется с применением Унифицированного Языка Моделирования (UML, Unified Modeling Language). UML является прямым потомком методов объектно-ориентированнога анализа и проектирования (OOA&D), которые появились в конце 80-х начале 90-х годов прошлого века. UML непосредственно унифицирует известные методы Г. Буча [167], Д. Рамбо (ОМТ) [209, 210] и А. Джекобсона [191], при этом он обладает гораздо большими возможностями. Язык UML стандартизирован консорциумом OMG (Object Management Group) и в настоящее время является стандартом OMG. Подробное описание истории возникновения, целей и задач этого средства визуализации, специфицирования, проектирования и документирования программных систем можно найти в книге Г. Буча, Д. Рамбо, А. Джекобсона [28], краткое изложение основных принципов языка приведено в работах К. Скотта [144] и М. Фаулера [180]. Предложенный в пятой главе работы каркас для разработки систем имитационного моделирования, основан на концепции типовых решений проектирования (паттернов, шаблонов проектирования), которая была предложена С. Alexander, S. Ishikawa, М. Silverstein, М. Jacob-son, I. Fiksdahl-King, S. Angel в [163], в области разработки программного обеспечения идея применения паттернов подробно изложена в классической монографии Э. Гаммы, Р. Хемла, Р. Джонсона Дж. Влиссидеса [42]. В качестве примера работы последнего времени, посвященной тематике применения типовых решений при разработке информационных систем, приведем монографию М. Фаулера [149]. Эти работы образуют методологическую основу решений, которые рассматриваются в последних главах настоящей работы.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе рассмотрены математические модели страховых компаний в предположениях, что поток входящих рисков является стационарным марковским потоком. Первый параграф главы посвящен исследованию модели, в которой интенсивность входящего потока линейно зависит от числа уже имеющихся рисков, а страховое поле считается неограниченным. В 1.1.1 относительно рассматриваемой модели делаются следующие предположения: Будем описывать состояние страховой компании в момент времени t двумерным случайным вектором {£(/), <S(f)}, где k(t) - число рисков, застрахованных компанией, a S(t) - ее капитал в момент времени t. Изменения капитала и числа застрахованных рисков происходят в следующих случаях:

1. Компания страхует новый риск. Будем предполагать, что поток приходящих рисков - это примитивный поток с параметром Х+р^ k{t). Первое слагаемое в последнем выражении отражает поток рисков, которые клиенты страхуют в компании по независящим от нее обстоятельствам, а второе - тот факт, что среди людей, не застраховавших свои риски, распространяется информация о страховой компании, происходит неявная реклама компании. Вероятность того, что за время At компания застрахует новый риск, равна + k)At+o(At). Каждый новый риск приносит компании страховую премию , размер которой является случайной величиной с функцией распределения и моментами М{Е} = а, М^21= а2

2. Будем считать далее, что по каждому из застрахованных рисков регулярно с интенсивностью выплачивается взнос в размере С,, который является случайной величиной с функцией распределения F^(z) и моментами М= с и

М^2 j= с2. Будем считать, что взносы вносятся независимо друг от друга и поэтому за время At в компанию поступит такой взнос с вероятностью kX^At + o(Af).

3. Страховое время некоторых рисков заканчивается. Будем считать, что каждый риск покидает компанию независимо от поведения других рисков с интенсивностью ц. Тогда за время At компанию покинет риск с вероятностью k^At + o(At).

4. Наконец, наступают страховые случаи. Будем считать, что с каждым клиентом может наступить страховой случай с интенсивностью |ип и эти страховые случаи для различных рисков независимы. Тогда на интервале At наступит страховой случай с вероятностью k\i^kt + o(kt), а компания при этом выплатит страховое возмещение в размере г|, которое является случайной величиной с функцией распределения Fn(z) и моментами М{г\} = Ь, м|п2]=Ь2.

Целью исследования является рассмотрение статистических характеристик процессов k(t) и S(t), т. е. поведение числа рисков и капитала компании.

В 1.1.2-1.1.7 проведены исследования модели в предположении, что процесс k(t) находится в стационарном режиме. Полученные результаты сформулированы • в виде следующих утверждений и теорем:

Теорема 1.1. В предположениях 1) - 4) финальные вероятности пj определяются соотношениями 1 да А-1 + ' пМй.

--Д"**-^. 0-2) а условие существования стационарного режима имеет вид р^ /ji < 1.

Так как полученные выражения для финальных вероятностей являются достаточно громоздкими и «неудобными» для дальнейшего использования, то было рассмотрено так называемое диффузионное приближение, когда процесс k(t) аппроксимируется диффузионным процессом. Для нахождения такого приближения использовался асимптотический метод анализа марковизируемых систем. С учетом следующих обозначений доказана теорема 1.2 р = \j\i, ф = Ре/ц. - = s2, B2j = к+ех, -П:=Р(х,е). Р 8

Теорема 1.2. Если

1. существует конечный предел lim Р(х, е) = Р(х),

Е->0

2. функция P{x,z) дважды дифференцируема по х, то функция Р(х) удовлетворяет уравнению Фоккера-Планка вида dx 1— ф dx

При этом константа к = l/(l - ср).

Таким образом, показано, что число рисков является асимптотически нормальной величиной при р-> со с М{к) = —= - ^ и D{k\ = , ^ „ = —. w 1-Ф ц-р w (i-ф)2 f

Теорема 1.3. В стационарном режиме функция корреляции процесса k(t) определяется выражением

Rok {ф-^-ехр{- Ox-pJ |т|}. (1.8)

1-ф)

Теорема 1.4. В стационарном для числа рисков режиме математическое ожидание капитала M{St} определяется выражением

M{St]rSQ+-^-{\xa + Xc>c-^b), (1.14) где SQ капитал компании в момент tQ.

Таким образом, условие возрастания среднего капитала имеет вид + 0, если \ш + \^с-\1цЬ<0, компания разоряется. Нетрудно заметить, что полученное условие имеет достаточно естественное экономическое обоснование, в среднем доходы компании должны превышать ее расходы.

Теорема 1.5. В стационарном для числа рисков режиме дисперсия капитала £>{£,} определяется выражением

-t+

2 Я.ц (ар, +с\ -bp ^-(ехр{-(ц-р4)}-l). (1.17)

Из равенства (1.17) следует, что при больших t дисперсия капитала, так же как и его среднее, растет пропорционально времени t.

Теорема 1.6. Функция корреляции соу{5,Д(/ + /0)} процессов St и k(t) в стационарном режиме для числа рисков имеет вид l^-Pa

Теорема 1.7. Функция корреляции процесса S, R^s{tx,t2) в стационарном режиме для числа рисков определяется выражением 1 op. -ь^ +е-М е-Мь-ч)г]9 (122) и-pj

В 1.1.8 получены характеристики k(t) в нестационарном режиме. В результате доказана следующая теорема:

Теорема 1.8. В нестационарном режиме среднее число рисков дисперсия числа рисков Dk(t) и функция корреляции числа рисков Rokih^) определяются выражениями ш= v-h

1.25) ц + Р^ ц-р4

1л ехр{- (ц-pJKl-ехр{- (ц - р4 >}),

1.26) (1.27)

RQAhh)= At(min(/j,t2))exp{-(|i-p^)\t2 где /0 число рисков, застрахованных в компании в момент /0 •

В 1.1.9 рассмотрено поведение капитала страховой компании S(t) в нестационарном режиме для числа рисков:

Теорема 1.9. В нестационарном режиме для числа рисков режиме математическое ожидание капитала M{St} и дисперсия капитала D{St} определяются выражениями

M{St} = So +-i— (яр4 -vnb + Хр] о ~~

H-Pj аР4 - \хцЬ + \с 1 /0--ехр(- (ц - Р4)/),

Ц-Р

1.30)

5 У

1 1 А, fe } =-а2 + + \С21 iQ (p ka-\LAb+\cf

Г, х 1 2 1

Ф-Pj 1° о 1 "зИЛ 2 M

X + lx) t X ЗА. ^ +

In---— H- (p р-ц^Ь+Х^У'

Ц-Р А.(А. + 2Ц + Р6) V + lo~

Ц-Р

ZJ ■

Ц-Р

2HJ

5 У

4-Pj

2 +ц) +

In (1.31)

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

Среднее значение капитала компании в нестационарном режиме при больших t растет пропорционально t, аналогично стационарному случаю. Дисперсия капитала при больших t растет, как /•.

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

20 начального числа рисков /0 от стационарного значения. Именно это слагаемое и описывает переходный режим.

Во втором параграфе рассмотрен случай, когда потенциальный рынок страховых услуг (страховое поле), на котором действует компания, ограничен, т.е. число рисков, которые может застраховать компания, ограничено. Пусть N - максимально возможное число рисков, за бесконечно малый промежуток времени At каждый из N потенциальных рисков может застраховаться с вероятностью ХД/ + о(Д/). Риск не может быть застрахован повторно, пока не истечет срок текущего договора. Величину k(t) по-прежнему будем считать числом застрахованных в компании рисков. Тогда суммарный поток поступления рисков в компанию будет примитивным потоком с интенсивностью

В отличие от модели, рассмотренной в предыдущем параграфе, вероятность поступления в компанию нового риска за время At в этом случае равна (N-k)XAt+o[At). По-прежнему в этой ситуации компания получает страховую премию , размер которой является случайной величиной с функцией распределения = a и Л/^2 }= д2-Остальные вероятности, связанные с приходом и уходом рисков, аналогичны модели с бесконечным страховым полем, которую мы рассмотрели выше.

Целью исследования, как и ранее, является рассмотрение статистических характеристик процессов k{t) и S(t), т. е. поведения числа застраховавшихся рисков и капитала компании в стационарном и нестационарном режимах.

В 1.2.2 доказаны следующие утверждения:

Теорема 1.10. Для модели с ограниченным страховым полем финальные вероятности Tij определяются соотношениями пj = CjNpJ(1 + p)"N,для V j = IN , (1.37) где р=Х/ц.

Теорема 1.11. При N-><x> в стационарном режиме распределение числа рисков является асимптотически нормальным с

1 j 1+р Х + ц

D{k} =

N р NK]i

О +Р)2 (X + nJ2'

Далее в 1.2.3 - 1.2.7 исследована деятельность компании в предположении, что число застрахованных рисков стационарно. Результаты сформулированы в виде следующих теорем:

Теорема 1.12. В стационарном режиме функция корреляции процесса k(t) определяется выражением W ■- М2 {kh-^%: ехр(- (1.42)

О+Р)

Теорема 1.13. В стационарном для числа рисков режиме математическое ожидание капитала M{St} определяется выражением

M&^So+^^c-^b + iia), (1.48)

1+р где S0 капитал компании в момент /0.

Таким образом, условие возрастания среднего капитала и в этом случае имеет вид \ьа + \^с-\хцЬ>Ъ. Если ца + с-цпb <0, компания разоряется. Нетрудно заметить, что полученное условие имеет достаточно естественное экономическое обоснование, в среднем доходы компании должны превышать ее расходы.

Теорема 1.14. В стационарном для числа рисков режиме дисперсия капитала Z^S",} определяется выражением

Np ——/+ 1+р (ак-Ь^ +с^)2^Н-(ехр{-(^ + ц>}-1). (1.51)

Теорема 1.15. Функция корреляции cov^^^ + Zq)} процессов S, и k(t) в стационарном режиме для числа рисков имеет вид соv{St, k{t +10 )} = -ак-Ьцц [1 - ехр {- (к + ц)}]. (1.55)

Х + ц)

Теорема 1.16. Функция корреляции процесса S, RqXH^i) в стационарном режиме для числа рисков определяется выражением

Rs{hh) = a2\i + b2\xл + с2Х^ +

NX\i {aX-b^ +сХ (1.56)

Аналогично параграфу 1.1.9 в 1.2.9 получены характеристики k{t) в нестационарном режиме. В результате доказана следующая теорема:

Теорема 1.17. В нестационарном режиме среднее число рисков k{t), среднее от числа рисков в квадрате k}{t) и функция корреляции числа рисков Rokih-fi) определяются выражениями

NX 4 w Х + ц

Х + ц ехр{-(?, + ц)}.

1.59)

А. + ц) ла+ц

А. + ц

1о~

NX A. + JI ехр{- (X + ц)}(1 - ехр{- (X + ц)}), (1.60)

Rok ('i Jti) = Dk (minfc ,t2 ))exp{- (X + ц)|* 2 - tx |}. (1.61) где /0 число рисков, застрахованных в компании в момент t0.

В 1.2.10 рассмотрено поведение капитала страховой компании S(t) в нестационарном режиме для числа рисков:

Теорема 1.18. В нестационарном режиме для числа рисков режиме математическое ожидание капитала M{St) и дисперсия капитала £>{£,} определяются выражениями

M{S,} = S0 + ^(v"'а%~is^'o -+ D{st} = (nn62 - + hc2 (<o -:

1.65)

NX fa-\LAb+X^cf

NX-5X + 5\i(. NX

-3 I lo о2

NX

X + tf б(>. + ц)3 i, X + \i) 3(^+ц): + - ■+\c2 )+ fa ~ Vx\b + \cfx

NX + ii) x<

NX lO

3NX

NX(NX + y) 2

NX

1(цnb2 - Xa2 + X&U , {Xa-\iAb + X^cf

NX(NX + 2\i-X) .2 ~lo 1

X-[L

X + »Y "(Х + ц) fa -\1ЦЬ+X^cf(i0

NX X + |i

2X-5NX-2\i e-2^fa-\iAb + Xtcf\iQ

NX X + \i

2(X + \if 3 NX-X+ii

2(X + iif e-^fa-^b + X^cf 3(X + [if

2 NX(NX + \i) 2NX-X+ji

NX

X + \i,

1.66) k + ц)2 ^ + H Как и ранее, в нестационарном.режиме математическое ожидание, дисперсия и функция корреляции числа застрахованных в компании рисков имеют более сложный вид, чем в стационарном режиме.

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

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

В параграфе 1.3. для моделей, рассмотренных в 1.1. и 1.2 , исследуется более сложная модель поведения капитала компании. Считается, что за время At капитал компании увеличится на rS(t)At, где'г - ставка банковского процента.

Как показано выше, в ситуациях, когда входящий поток рисков является стационарным марковским потоком, исследование моделей необходимо начинать с определения вероятностных характеристик процесса k(t). Однако из найденных ранее соотношений следует, что характеристики процесса k{t) могут быть определены независимо от капитала компании S(t). Поэтому дополнительно были исследованы только характеристики капитала компании. (1.3.2-1.3.4).

Теорема 1.19. Для моделей с учетом банковского процента в стационарном для числа рисков режиме математическое ожидание капитала компании определяется выражением

5,(()=S„-е" +Sl±h^Ml.Lr' -l), (1.82) Г где Sq - стартовый капитал компании.

Следствие. Математическое ожидание капитала компании в стационарном режиме для модели с неограниченным страховым полем равно

1-86,

Для модели с ограниченным страховым полем

87)

Заметим, что при г = О выражения (1.86) и (1.87) превращаются в соотношения (1.14) и (1.48) соответственно.

Теорема 1.20. Дисперсия капитала страховой компании для модели с работающим капиталом в стационарном для числа рисков режиме определяется выражениями

2rt

D {£(*)} = Ха2 + С2 + \ьцЬ2 + ^а2 )

2 г

-1

• + NXji (г-(Х + ф2г'+r + x + [l 2 С • --—--t~z—:-ТТ\--г6 oq для модели с неограниченным страховым полем и

V А п \~2rt

NX

D (S(/)} = Ша2 + (х^с2 - Ха2 ) е -1

2 г с •

2 NX\i (г-(Х + ф2г1 -2 геИ^»)> +Г + Х \х + »)2 r{r2-(X^f J для модели с ограниченным страховым полем, где S0 - стартовый капитал компании.

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

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

1. С вероятностью А,Д* + <?(Д/) в компанию поступает новый риск. Клиенты, страхующие риски, вносят страховые премии, являющиеся случайными величинами с математическим ожиданием а. Риски поступают независимо друг от друга.

2. Время страхования некоторых рисков заканчивается, и клиенты покидают компанию. Считается, что каждый клиент уходит из компании с вероятностью цД/ + <з(Д/) независимо от другого.

3. С вероятностью X^At + o(At) каждый из клиентов компании уплачивает дополнительный взнос, который является случайной величиной с математическим ожиданием с.

4. Наконец, с вероятностью цл At + o(At) у каждого из клиентов компании независимо от других наступает страховой случай, выплата по которому является случайной величиной с математическим ожиданием b.

В 1.4.2 описана модель взаимодействия страховых компаний. Считается, что вероятности прихода клиентов в ту или иную компанию Xi зависят от величин страховых взносов и а2, то есть' Xt =Я,;(а1,а2)- Далее довольно естественно г г где 6j, а2, Ь2, М{к}, D{k) и 0, параметры соответствующих моделей. считать, что эти вероятности обратно пропорциональны величине страховых взносов, т. е. Хх/Х2 =аг!а\ • Если f(ax, а2) = Хх(ах, а2)+Х2(ах, а2).

Тогда f(ax,a2),X2= —— f(ax,a2). (1.103) ах + a2 ci\ + a2

Вид зависимости f(ax,й2)отй| и a2 определить достаточно сложно. В настоящей работе мы будем считать, что f(aXi а2) зависит лишь от некоторого параметра р, который, в свою очередь, зависит от ах и а2, так что р = р(ах, а2). Относительно вида функции f(p) естественно выдвинуть следующие предположения:

1. У(0)<+оо;

2. f(p) монотонно убывает с ростом р;

3. lim /(/?) = 0;

00

4. lim pf(p)= 0. р—>+00

Что касается самой зависимости р(ах,а2), то к ней можно предъявить следующие достаточно естественные требования. Пусть а2 фиксировано, тогда:

1. Если ах -» +оо, то р должно равняться а2, так как при этих условиях вся динамика рынка страхования будет определяться именно этой величиной, клиенты просто будут игнорировать первую компанию. Таким образом, сформулируем следующее утверждение: lim р(ах,а2)=а2.

Я,->+00

2. С уменьшением ах р(ах,а2) также должно монотонно убывать, так как у желающих застраховаться в этом случае появляется возможность выбора. Поэтому должно выполняться следующее условие: ) дах

0.

3. При ах -»0 р(ах, а2) также должна стремиться к нулю, так как в этом случае страховаться будут все потенциальные клиенты, следовательно, liin р(ах, а2 ) = 0. 27

4. Наконец, р(а1,а2) должно бьггь симметричной функцией относительно переменных и а2, т. е. р(а{, а2)= р{а2, ах).

По-видимому, достаточно правдоподобной является следующая зависимость р от ах и а2:

Л = (1.104) р ах а2 удовлетворяющая всем вышеперечисленным условиям и напоминающая среднее геометрическое с некоторым параметром v > 0. В данной работе будем считать, что

Р=г ^ . (1.105)

K+oJ

Таким образом, выражения (1.103) можно записать в следующем виде:

Я., = ^ (1.106) tfj + а2 а2 тогда ikt«i]/Wi =£ik±l2)/0,) (1107) a j + д2 а\ + а2

Цель каждой из компаний состоит в том, чтобы, выбирая величину средней страховой премии а/5 максимизировать величину капитала компании. С математической точки зрения получившаяся задача представляет собой кооперативную игру двух лиц с ненулевой суммой. В основе решения получившейся игры лежит построение переговорного множества (множества Парето), на котором происходит согласование стратегий игроков. В работе предложен алгоритм построения переговорного множества для различных значений v и 5, рассмотрено несколько конкретных примеров его построения.

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

В 2.1. описывается модель страховой компании. Предполагается, что

1. Компания страхует новый риск. Будем предполагать, что поток приходящих рисков - это пуассоновский поток'с переменной интенсивностью с параметром Вероятность того, что за время At компания застрахует новый риск, равна А,(г)Д/+о(Д/). Каждый новый риск приносит компании страховую премию размер которой является случайной величиной с функцией распределения и моментами М= ах(/), м{|2}= а2{$). Так как величина первого страхового взноса - это прерогатива компании и может устанавливаться по её усмотрению, поэтому будем считать, что ах (/), a2(t) зависят от времени и являются величинами, которыми можно управлять. Также предположим, что ax(t\ a2(t) могут быть функциями, имеющими разрывы первого рода.

2. Будем считать далее, что по каждому из застрахованных рисков регулярно с интенсивностью Х^ выплачивается'взнос в размере^, который является случайной величиной с функцией распределения F^(z) и моментами М{<^} = с, и с2. Будем считать, что взносы вносятся независимо друг от друга и поэтому за время At в компанию поступит такой взнос с вероятностью kX^At + o{At).

3. Страховое время некоторых рисков заканчивается. Будем считать, что каждый риск покидает компанию независимо от поведения других рисков с интенсивностью [I. Тогда за время At компанию покинет риск с вероятностью к цА/ + o{At).

4. Наконец, наступают страховые случаи. Будем считать, что с каждым клиентом может наступить страховой случай с интенсивностью и эти страховые случаи для различных рисков независимы. Тогда на интервале At наступит страховой случай с вероятностью k\i^At + o{At), а компания при этом выплатит страховое возмещение в размере г|, которое является случайной величиной с функцией распределения F4{z) и моментами = м\г\г\=Ь2.

В параграфе 2.2. предполагается, что поток входящих рисков пуассоновский и его параметр X(t) - непрерывная функция или функция, имеющая разрывы только первого рода. В 2.2.1. определяются выражения характеристик капитала и числа застрахованных рисков в случае неограниченного страхового поля.

Теорема 2.1. В предположениях 1)-4) математическое ожидание kx(t)= M{k{tf$, дисперсия Dk{t)u функция ковариации числа застрахованных рисков Ск (tx, t2) определяются выражениями кх (0=k(tQ У^ + )x(zy^dz, в

2.1)

2.2) (2.3) где k(t0) число рисков в начальный момент времени.

Теорема 2.2. В предположениях 1) - 4) математическое ожидание Sx(t) и дисперсия капитала компании Ds(t) определяются выражениями

Я, (О = S0 + )ах (u)x(u)du - ЬЫ^к х х

Jx(v)(l - + (l - ) )x(vyvdv

2.21) где S0 капитал компании в момент tQ,

Ds{t)=\a2{z)k(z)dz + о dz.

2.22) где, для простоты S0 = 0.

Теорема 2.3. В предположениях 1) -4) функция ковариации капитала и числа застрахованных рисков Csk (t) определяется выражением ft о I cSk(t)=(ci\ ~^iMti) 1 \uX(t-uy^du+te^' \\{uyudu\ (2.35) lo -oo J

В 2.2.2. определены выражения для тех же что и в 2.2.1 характеристик капитала и числа рисков в случае ограниченного страхового поля (теоремы 2.4-2.6).

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

-v{t-tQ\-\\{x)dx tf / ^ l(0 = *('o>

-n\\(z}> dz, i . . \

Dk(t)= J

-2ii(l-z)-2\X(x)dx dz, max (4,12) min((i,'2)

2.39)

2.40)

2.41)

С4(/„/2) = Д(тт(/„/2)> где k(t0) число рисков в начальный момент времени.

Теорема 2.5. Для модели с ограниченным страховым полем математическое ожидание S^t), дисперсия капитала Ds(t) определяются выражениями it t 5, (/) = S0 + jAfa, (u)\(u)du - J&-, (mJs, (u)k(u)du - (ь,цп - сх\(u)du, (2.54) о о о где S0 - капитал в начальный момент времени.

I t

Ds (/) = N ja2 (u)X(u)du + + ) JAj (u)du -n n tfz

- J&! (u]a2 {u)k{u)du + 2(^^-^j2! \ck(z,u)du 0 040 tfz 2(Vn - )J \ck{z,u\a\(u)x(u)+ ax(z)x(z)\lu nV n dz

0V0 r tfz 2 J JQ (z, u)a{ (u)X(u)ax (z)X(z)du oU dz,

2.55) где, для простоты S0 = 0.

Теорема 2.6. Для модели с ограниченным страховым полем, функция ковариации капитала и числа застрахованных рисков CSk (t) определяется выражением t t

-h\)\ck{t,u)du-\ck{t,u)ax{u)x{u)du, (2.71) n n где Ck(t,u) выражается формулой (2.41).

В 2.2.3. предполагается, что ax(t\ a2(t) влияют на интенсивность потока входящих рисков, так как уменьшение a{(t) приводит к увеличению X(t) и, наоборот, при увеличении ax(f) число желающих застраховаться уменьшается. Поэтому в общем случае предполагается, что X зависит от ах (t) и от t. Предположение 5). Для конкретизации будем считать, что

X(t)=F(ai(t))X0(t), (2.74) где F(o) = 1, F(+ со) = 0 и F(a{) монотонно убывает с ростом ax{t), XQ(t) имеет смысл максимальной интенсивности потока входящих рисков.

Теорема 2.7. В предположениях 1) - 5) оптимальное управление ax(v),

О < v < Т, доставляющее максимум среднему значению капитала в конце рассматриваемого промежутка, определяется из уравнения

Следствие. При Т -» оо уравнение (2.78) приобретёт вид

2.79) что и определяет асимптотическое значение tfj(v).

В параграфе 2.3 предполагается, что параметр входящего потока X(t) - непрерывный в средне квадратичном стационарный случайный процесс. Предположение 6). Пусть характеристики X(t):

M{x(t)}=7i,

R(z)=R0(z)+X2, lim R0 (z) = 0.

В 2.3.1. и 2.3.2 определяются выражения для характеристик числа застрахованных рисков и капитала для неограниченного страхового поля.

Теорема 2.8. В предположениях 1) - 4), 6) математическое ожидание, дисперсия и функция ковариации числа застрахованных рисков определяются выражениями мхт=~,

Dk(t) = ^\R0(v)e^dv + K

2.81) (2.81)

2ц 2ц |(!J(i|

Теорема 2.9. В предположениях 1) - 4), 6) математическое ожидание ^.{^(О) и дисперсия Ds(t) капитала компании определяются выражениями мх {5! (0} = +1 [ах (v)dv + xt, t .м J t

Ds{t)=l\a2{z)dz ъ2\1ц+с2\+-(ъх\1ц-сх\} n -f -1) + )]r0{и ~ v)ax(u)ax{v)dudv ^ .00 l^g-HW+Ig^+Jul 2 ), / \(\ 11+ o Л И dx +

0 ^ oo ^ )r0(xj^l -г**- + )'\R0 (u - v)ax (v)(l - e-^-u))iudv + (l - e* ) J | J/?0 (и - v)ax (uyvdu

0 0 -ooVO dx + У dv (2.86) где S0 - значение капитала в начальный момент времени.

В 2.3.3 предполагается, что величина ах, обозначающая средние первоначальные взносы является функцией, зависящей от X0(t).

Предположение 7). Пусть al = cii(X0(t)). Рассмотрим случай, когда интенсивность потока рисков страхующихся в компании, связана с X0(t) соотношением b(') = *o('M*iM')))- (2.95)

Выражение (2.95) имеет следующий смысл: A.0(f) есть интенсивность потока рисков, желающих застраховаться, но реально каждый риск страхуется с вероятностью-F^j) зависящей, естественно, от величины взноса ах, таким образом, интенсивность потока рисков, которые будут застрахованы, равна X(t) = Х0 {t)F(al). Желание управлять этой интенсивностью изменяя величину страхового взноса ах, находит свое отражение в том, что ах из константы превращается в функцию от

М')>т-е- а\ =fli(M0)

Теорема 2.11. В предположениях 1)-4), 7) оптимальное управление а\ ~ а\ (^о (0)' реализующее максимальный темп роста среднего значения капитала определяется как решение уравнения.

F(ax)+{ax + (2.100)

I V

Теорема 2.12. Пусть Л.0(м) = A,j. и A.0(v) = Х2 имеют совместную плотность распределения вероятностей р(Хх ,Х 2, w) = р0(Хх ,Х 2, w)+ р(Хх )р(Х2), lim p0(Xl,X2,w)=0,

W-M тогда в предположениях 1) - 4), 7) асимптотическое поведение дисперсии капитала определяется выражением

D (V) 00 lim = \p(X0)a2{X0)XQF{ax{X0))dX0 + 1 — И

->oo t f

AJ p(xo)X0F(al(X0))dX0

- 2 J Jtfj (A,,)ax (A,2)X{X2F(a{(A,j ))F(aj(A,2l,X2)dXldX

2 +

0 0 2

V J uu uu

J jXlX2F(al (k, ))F(a, (X2 Л2 о 0 (2.104) 00

Теорема 2.13. В предположениях 1) - 4), 7) в случае экспоненциально распределенных первоначальных страховых взносов оптимальное управление страховыми взносами, реализующее максимум среднего темпа роста капитала при ограничении на флуктуацию этого капитала, то есть

Mx{s{f)} lim — w' => max, при

->oo t

HmM.fl,.

->00 t определяется из уравнения li

I ^ /

4 а^Щ'М*

F'biM

J J H 2IXV

F'(ax(Xx))+ 2

AF{ai{Xx))+ax{Xx)Fiax{X\)))

Al k2F(ax (X2 ,X 2 )dX: 0. (2.109)

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

В параграфе 3.1 дается описание модели. Считается, что в момент времени t состояние страховой компании характеризуется её капиталом S{t) и числом застрахованных рисков k(t). Далее будем считать, что компания отчисляет часть своего капитала на рекламную компанию, так что на интервале [/, t + Д*] на её проведение выделяется a(t)s(t)At денег, где а(/) - доля капитала, выделяемая на рекламу в единицу времени.

1. С вероятностью (Я0 + a(t)kiS(t))At + o{At) в компанию придёт новый риск и, страхуясь, внесёт страховую плату , которая является случайной величиной с функцией распределения F^(x). Слагаемое atyfasty) описывает приток новых рисков, обусловленный рекламой, параметр А,0 определяет интенсивность прихода рисков, не находящихся под воздействием рекламы; параметр определяет эффективность рекламы.

2. С вероятностью fin&(/)A/ + o(A/) произойдёт страховой случай и компания выплатит страховое возмещение в размере rj, которое является случайной величиной с функцией распределения Fn(x). Будем считать, что страховые случаи независимы.

3. С вероятностью X^k(t)At + o(At) поступит очередной взнос С, от застраховавшихся в компании рисков, который является случайной величиной с функцией распределения F^(x). Предполагается, что взносы осуществляются независимо друг от друга.

4. Наконец, с вероятностью \xk{t)At + o(At) прекращается действие страхового договора и один из рисков покинет компанию независимо от других рисков.

В параграфе 3.2 исследуется модель для случая неограниченного страхового поля. В 3.2.1 находятся выражения математического ожидания капитала и числа рисков компании:

Щ=С1еУ + С2е™ + аЯ°

Y1Y2

С2(у2 + ха у,у2

3.5) константы определяются из начальных условий S(o)=S0,k(o) = к0.

В 3.2.2 определяются условия эффективности рекламы, и результаты сформулированы в виде следующей теоремы.

Теорема 3.1. Для модели с учетом расходов на рекламу в случае неограниченного страхового поля реклама эффективна, если выполнено хотя бы одно из условий f с\-Ь\1цу а +

1, + (3.8) ч И /

Расшифруем экономический смысл полученных условий эффективности рекламы. Рассматривая первое условие, перепишем его в виде

A.,a S а+—2-L ч aS. И

В последнем неравенстве слева первый сомножитель X{aS - это интенсивность дополнительного потока рисков, который возникает из-за действия рекламы, сомножитель а+\ск^ -Ь[1ЦУ[1 - это средний доход, который получает компания от каждого застрахованного риска,

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

Смысл второго критерия менее очевиден. Для того, чтобы его получить, представим второе неравенство из (3.8) в следующем виде: aX1<zS'>aS' + |xS'.

В этом случае левая часть неравенства ctkxaS - это дополнительный доход, которые получает компания за счет первоначальных взносов от рисков, привлеченных рекламой, справа aS - это, как и ранее, расходы на рекламу в единицу времени, а последнее слагаемое запишем в виде \ik(t)S/k(t) - это потери капитала, связанные с завершением срока договоров страхования в единицу времени. Следовательно, экономический смысл второго критерия эффективности рекламы можно сформулировать в следующем виде: дополнительный доход от первоначальных взносов, полученных за счет рекламы, больше чем расходы на ее проведение и потерь капитала, связанных с завершением срока действия некоторых договоров страхования.

Следствие. При выполнении первого условия (3.8) реклама будет эффективна при любом а>0. Второе условие (3.8) может быть выполнено лишь при A.jd! > 1; оно даёт ограничение на величину параметра а, определяющего количество денег, выделяемых на рекламу. Имеем

7-^-7. (3.11)

То есть вкладывать деньги в рекламу стоит лишь тогда, когда параметр, соответствующий доле денег на рекламу, превышает некоторую величину.

В 3.2.3 решается следующая задача оптимального управления: пусть страховая компания начинает свою работу в момент времени t = О и её цель провести рекламную программу таким образом, чтобы к концу рассматриваемого промежутка времени [0,Г] темп прироста её среднего капитал стал максимальным. Критерий

M{S(T)}-S0 оптимальности определяется следующим образом: и = —, , . . - => max. Макл/^да симум определяется по виду а(/), 0<t<T с учетом ограничений 0<а0; S0- капитал компании в начальный момент времени. Заметим, что это лишь одна из постановок задач на оптимизацию рекламной деятельности. Для решения задачи рассматриваемый промежуток времени [0,г] делится на три: [0,Г0], [7"0,Z*j ] и \ГХ ,Г], причем каждому из промежутков соответствует своё значение а(/): а (0 =

О, 0</<Г0, aQ,T0<t<Tx, О, Тх <t<T.

Система относительно средних капитала и числа рисков получена в 3.2.1 Система для ковариаций капитала и числа рисков выводится в данном параграфе и имеет следующий вид dt 2[с\ - Ьцц )Csk (/)+ (с2\ + Ь2цп )M{k(t)}+Х0а2, iM{k(t)}~ 2\ФЩ}+ {*}+ X0> at

К "f-n. Xxaa(t )M )}+ X 0a + (X x aa(t)- a(t)- (t). (3.30)

Решение системы находится для каждого из промежутков, где параметр а(/) изменяет свои значения. На первом промежутке решение имеет вид

Dk(t,0)=Cxe~2* + Dxe~*

3.35)

Ся(/,0)=(сЯ,с+ Ae^j+Cj^ +{сХ^-Ь»ц + (3.36) И / И

Ds(t,0) =

2Х о(С2хс +д2ц) t—-z^)2 te-»" + jj. M* e

- 2C2 - ) 2A - Z>fin f V- \i2

Сх(сХ^-Ъ\1ц} Cx(cX^-b\i4Y C3,

3.37) / константы находятся из начальных условий Ds(0,0)=0, £^(0,0)=0? C5>t(0,0) = 0.' Далее для каждого из решений, соответствующих второму и третьему промежутку константы находятся из условий сшивания на границе M{s(T0 ,0)} = M{s(T0, a)}, Л/'{5(Г0,0)} = A/'{iS(r0,a)} (для решения на втором промежутке) и Ds(Tx,a)=Ds(Tx,0), CSk(Тх,а)= CSk(Тх,0) ,Dk(Тх,а) = D(TX,0) (на третьем). Поэтому в выражениях для среднего значения и дисперсии капитала в момент времени Т присутствуют константы, зависящие от моментов Т0 (включение рекламы) и Тх (выключение рекламы). При выполнении условий Г0 > 0 и TQ<TX<T условие максимума по по Т0 и Тх имеет вид dU(T0,Tx)l дТ0 dU{T0> 7j)

ЭГ,

Страховая компания включает рекламу в нулевой момент времени. Далее задача сводится к отысканию с помощью численных методов такого момента Тх, при котором функция C/(0,7j) имеет максимум.

В параграфе 3.3 исследуется модель страховой компании с учетом отчислений на рекламу в случае ограниченного страхового поля. В 2.3.1 выводится и решается система относительно средних значений капитала и числа рисков аналогично случаю неограниченного страхового поля. В 2.3.2 рассматривается стационарный режим страховой компании. Находится точка покоя системы

S9 к = yv—----7-Г-1=, (З.оУ) а(А,0 + - + a\i)+ л/D в асимптотике при N -»со определяются условия, при которых эта точка является устойчивой.

В 3.3.2 определяются условия эффективности рекламы (теорема 3.2) Теорема 3.2. Для модели с учетом расходов на рекламу в случае ограниченного страхового поля реклама эффективна, если выполнено хотя бы одно из условий Щ а +

1,

Х,а>1 + -- (3-77) a

Во второй части диссертации (глава IV) предлагается и исследуется математическая модель работы фонда социального страхования при различных предположениях относительно способа страховых выплат. В отличие от обычных страховых компаний в задачу фонда входит не только оплата страховых случаев (временная нетрудоспособность, пособия по беременности и родам и т.д.), но и систематические выплаты по реализации региональных и отраслевых программ по охране здоровья работников, санаторно-курортному лечению, обслуживанию детей и т.д. нисиьпижля ГОСУДАРСТВЕННАЯ БИБЛИОТЕКА

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

В параграфе 4.2 строится диффузионная аппроксимация деятельности фонда. Основной характеристикой состояния фонда является его капитал S(t) в момент времени t. С этим капиталом происходят следующие изменения:

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

2. Фонд выделяет часть своих средств на социальные программы. Мы будем считать, что эти средства также выделяются непрерывно во времени, однако скорость их выделения c*(S) зависит от величины капитала S в данный момент времени.

Величину Cq-c*(S) мы в дальнейшем будем обозначать как c(S). Таким образом, c(S) есть скорость изменения капитала за счет детерминированных расходов и она зависит от величины капитала S. Именно в наличии слагаемого с* (S) и зависимости c(S) от S и заключается отличие данной модели от классической.

3. Происходят страховые выплаты, Будем считать, что поток страховых выплат является пуассоновским потоком постоянной интенсивности X, и сами страховые выплаты Е, являются независимыми одинаково распределенными случайными величинами с функцией распределения F^ (t) и начальными моментами

М© = й,ИМ{Й = А2.

Кроме того, будем считать, что достижение порога S(t) = 0 не приводит к разорению фонда, и даже при S(t) < 0 он продолжает функционировать, только происходят задержки по страховым выплатам.

Заметим, что процесс S(t) - марковский процесс. Обозначим P(S,/) = Pr{S < S(t)<S + dS}/dS.

Теорема 4.1. Если

1. плотность распределения вероятностей P(S, t) дифференцируема по t,

2. произведение c(S)P(S, /) дифференцируемо по S,

00

3. несобственный интеграл J/^S + и, t)dF^ (и) конечен, о ' то плотность распределения вероятностей P(S, t) удовлетворяет прямому уравнению Колмогорова следующего вида

4.1)

Решение уравнения (4.1) в общем виде представляет определенные математические проблемы, сделав ряд дополнительных предположений, его можно найти, воспользовавшись методом асимптотического анализа марковизируемых систем. Сформулируем дополнительные предположения, пусть процесс S(t) принимает достаточно большие значения, если выберать бесконечно малый параметр s > 0 таким образом, чтобы при е О процесс s2S(t) не равен тождественно нулю. c(5) = c1(s25'), т, Ss2=x(x)+ey, -P(S,t) = U(y,x,s), (4.3)

Теорема 4.2. Если

1. функция Cj(x) дифференцируема,

2. существует конечный предел НтП(у, х, е)=П(у, т),

• е-»0

3. функция П(у, т) дифференцируема по т и дважды по у, то функция П(у, т) удовлетворяет уравнению Фоккера-Планка вида gfei = -сЩ8^' Т»т) дх ду 2 ду

При этом функция х(х) определяется как решение дифференциального уравнения

4.6) ах

Из полученного уравнения (4.5) следует, что процесс s в пределе при е—>0 по распределению сходится к процессу у(х), удовлетворяющему стохастическому дифференциальному уравнению вида dy(x) = с[ (x(x))ydx + Jka^dw(x), (4.8) где w(t) - стандартный винеровский процесс.

Из уравнений (4.6) и (4.8) вытекает следующее утверждение.

Теорема 4.3. Если функция сх (х) дважды дифференцируема, то случайный процесс z\

М =х(т)+Еу(т)

4.9) удовлетворяет равенству вида

4.10) где R2 = -у2с"(еу9), 0 £ 0 £ 1.

Таким образом, показано, что процесс S(t) можно аппроксимировать решением стохастического дифференциального уравнения

Так как фонд не стремится к неограниченному накоплению капитала, то существует стационарное (финальное) распределение вероятностей p(S) капитала S. где константа С находится из условия нормировки.

Далее в 4.2.3 находится плотность вероятностей величины капитала S при фиксированном релейном управлении. В 4.2.4 определяются вероятностные характеристики деятельности фонда при фиксированном управлении: вероятность неплатежеспособности фонда и вероятность выделения денег на социальные расходы - и предложена процедура, позволяющая определить параметры управления, которые обеспечивают заданные вероятностные характеристики работы фонда. В 4.2.5 и 4.2.6 аналогичная задача решается для релейно-гистерезисной процедуры управления капиталом.

В параграфе 4.3. получены точные формулы для плотности вероятностей капитала фонда в стационарном режиме и предложена процедура, позволяющая определить параметры управления, которые обеспечивают заданные вероятностные характеристики работы фонда, в предположении, что страховые выплаты распределены по экспоненциальному закону, когда поток страховых выплат является пу-ассоновским потоком постоянной интенсивности X, и сами страховые выплаты ^ являются независимыми одинаково распределенными случайными величинами с экспоненциальным распределением р^(х) -ехр{-х/а}/а, х>0. dS{t) = (c(S) - a{k)dt + Ja^kdw(t).

4.13)

4.14)

В 4.4. решена задача нахождения вероятностных характеристик капитала фонда при произвольном законе управления его капиталом, но в отсутствии гистерезиса при этом управлении.

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

В параграфе 5.2. рассмотрены элементы дискретно-событийной имитационной модели. В 5.2.1 приведено формальное описание метода моделирования, а в 5.2.2 выполнена функциональная декомпозиция элементов дискретно-событийной модели, проиллюстрированная соответствующей диаграммой состояний (рис. 5.2).

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

В параграфе 5.4. рассмотрены вопросы, связанные с проектированием уровня бизнес-логики в рамках стандартной трехслойной архитектуры. В 5.4.1 приведено решение, позволяющее использовать диаграммы активности и диаграммы состояний UML для документирования метода представления событий с помощью графов, предложенного L.W. Schruben. Далее в 5.4.2. рассмотрена задача организации моделей, имеющих сложную иерархическую структуру с применением паттерна компоновщик. Для инициализации компоновщика используется типовое решение строитель. Приведены примеры организации объектов, формирующих модель, в рамках компоновщика. В 5.4.3 рассмотрена проблема взаимодействия с библиотечными пакетами, для решения которой предлагается использовать паттерны стратегия и адаптер. Наконец, в 5.4.4.описан порядок инициализации и разработки имитационной модели с помощью каркаса.

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

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

ОСНОВНЫЕ НАУЧНЫЕ РЕЗУЛЬТАТЫ

Основные научные результаты, полученные автором и выносимые на защиту, состоят в следующем:

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

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

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

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

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

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

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

МЕТОДИКА ИССЛЕДОВАНИЙ

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

ТЕОРЕТИЧЕСКАЯ ЦЕННОСТЬ РАБОТЫ

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

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

ПРАКТИЧЕСКАЯ ЦЕННОСТЬ РАБОТЫ

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

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

ПУБЛИКАЦИИ ПО РАБОТЕ

Содержание работы опубликовано в 68 печатных работах. Основное содержание работы отражено в следующих публикациях в журналах, рекомендованных ВАК для опубликования основных научных результатов докторских диссертаций:

1. Змеев О.А. Модель функционирования страховой компании при конечном числе возможных клиентов.// Известия вузов. Физика. № 4.1999, с. 34-39.

2. Змеев О.А. Расчет характеристик времени разорения страховой компании для моделей с интенсивностью входного потока, зависящей от числа имеющихся рисков.// Известия вузов. Физика. №4.2000, с. 10-16.

3. Змеев О.А. Математические модели функционирования страховой компании с учетом банковского процента.// Известия вузов. Физика. №1.2001, с. 19-25.

4. Ахмедова Д.Д., Змеев О.А. Оптимизация расходов на рекламу при деятельности страховой компании.// Известия вузов. Физика, № 6.2001, с. 3-7.

5. Змеев О.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний.// Известия вузов. Физика, № 2. 2002, с. 24-28.

6. Ахмедова Д.Д., Змеев О.А., Терпугов А.Ф. Оптимизация деятельности страховой компании с учетом расходов на рекламу.// Вестник Том. гос. ун-та. № 275, 2002, с. 181-184.

7. Змеев О.А., Лезарев А.В. Шаблон объектного проектирования для реализации функциональности процесса моделирования в имитационных моделях систем массового обслуживания.// Вестник Том. гос. ун-та. № 275,2002, с. 108-111.

8. Змеев О.А., Моисеев А.Н. Шаблон диаграммы компонентов информационной системы корпоративного уровня.// Вестник Том. гос. ун-та. № 275, 2002, с. 130-132.

9. Ахмедова Д.Д., Змеев О.А. Нахождение характеристик страховой компании с нестационарным входящим потоком и ограниченными числом рисков.// Вест- -ник Том. гос. ун-та. Приложение №1 (I), 2002. С. 3-8.

47

10. Змеев О.А., Моисеев А.Н., Новиков Д.В. К вопросу проектирования уровня хранения в виде OOP ДБ.// Вестник Том. гос. ун-та. Приложение №1 (II), 2002. С. 363-367.

11. Змеев О.А. Математическая модель фонда социального страхования с детерминированными расходами на социальные программы (диффузионное приближение).// Известия вузов. Физика, № 3. 2003, с. 83-87.

12. Змеев О.А. Математическая модель фонда социального страхования со случайными расходами на социальные программы (диффузионное приближение).// Известия вузов. Физика, № 3. 2003, с. 88-93.

13. Змеев О.А. Математическая модель деятельности фонда социального страхования при экспоненциальных страховых выплатах.// Вестник Том. гос. ун-та,

2003, №280. С. 130-135.

14. Змеев О.А., Моисеев А.Н. Сравнительный анализ некоторых методов O-R-преобразования.//Вестник Том. гос. ун-та, 2003, № 280. С. 263-271.

15. Вальц О.В., Змеев О.А. Диффузионная аппроксимация модели фонда социального страхования с релейно-гистерезисным управлением капитала // Известия вузов. Физика, 2004. № 2. - С. 26-31.

16. Змеев О.А. Деятельность фонда социального страхования при релейно-гисте-резисном управлении капиталом.// Математическое моделирование, 2004, т. 16, № 2. -С.43-53.

17. Вальц О.В., Змеев О.А. Исследование модели фонда социального страхования.// Обозрение прикладной и промышленной математики. Т 11. Вып. 2.

2004.-С. 311-312.

18. Вальц О.В., Змеев О.А. Математическая модель деятельности фонда социального страхования при экспоненциальных страховых выплатах и случайными расходами на социальные программы.// Вестник Том. гос. ун-та. № 284, 2004. С. 37-41.

19. Войтиков К.Ю., Змеев О.А., Моисеев А.Н., Якушев АЛ. Архитектура надстраиваемых приложений клиент/сервер с обобщенным протоколом передачи данных.//Вестник Том. гос. ун-та. № 284,2004. С. 169-173.

20. Змеев О.А., Приступа А.В. Разработка объектно-ориентированного программного комплекса имитационного моделирования систем массового обслуживания.// Вестник Том. гос. ун-та. № 284,2004. С. 174-176.

Частично результаты, отражающие основное содержание работы, опубликованы в монографии

Глухова Е.В., Змеев О.А., Лившиц К.И. Математические модели страхования. Томск: Изд-во Том. ун-та, 2004. 178 с. В ней автору принадлежит глава 4.

Кроме этого, по работе имеются следующие публикации:

1. Змеев О.А., Змеева Е.Е. Модель функционирования страховой компании при конечном числе возможных клиентов // Наука и образование: пути интеграции. Часть 2. Тезисы докладов. Анжеро-Судженск. 1998, с. 32-33.

2. Змеев О.А., Терпугов А.Ф. Модель функционирования страховой компании при интенсивности потока, зависящего от числа клиентов // Наука и образование: пути интеграции. Часть 2. Тезисы докладов. Анжеро-Судженск. 1998, с. 33-34.

3. Змеев О.А. Модель функционирования страховой компании при интенсивности входящего потока, зависящего от числа клиентов // Математическое моделирование. Кибернетика. Информатика. Сб. статей. - Томск: Изд-во Том. унта, 1999.-С. 67-73.

4. Змеев О.А., Терпугов А.Ф. О некоторых подходах к вопросам моделирования деятельности страховых компаний // Образование и наука на пороге третьего тысячелетия: научно-теоретическая конференция. Тезисы докладов. Барнаул. 1999, с. 88-90.

5. Змеев О.А. Определение вероятности разорения страховой компании для модели с интенсивностью входящего потока, зависящей от числа клиентов // Статистическая обработка данных и управление в сложных системах. Сб. статей- Томск: Изд-во Том. ун-та, 1999. - Вып. 1. - С. 57-66.

6. Змеев О.А. Определение вероятности разорения страховой компании для модели с конечным числом возможных клиентов // Статистическая обработка данных и управление в сложных системах. Сб. статей - Томск: Изд-во Том. ун-та, 1999. - Вып. 1. - С. 66-75.

7. Змеев О.А., Змеева Е.Е. Расчет характеристик времени разорения страховой компании для моделей с интенсивностью входного потока, зависящей от числа имеющихся рисков // Качество образования и наука. Тезисы докладов. Анжеро-Судженск. 1999, с. 32-33.

8. Змеев О.А., Терпугов А.Ф. Расчет характеристик времени разорения страховой компании для моделей с конечным числом возможных рисков // Качество образования и наука. Тезисы докладов. Анжеро-Судженск. 1999, с. 33-34.

9. Змеев О.А. Математическая модель страховой компании с «работающим» капиталом.// Повышение эффективности научных исследований и совершенствование учебного процесса. Ч. I. Межрегиональная научно-методическая конференция. Тезисы докладов. Анжеро-Судженск. 2000.

10. Змеев О.А., Терпугов А.Ф. Модель страховой компании с ограничением на число клиентов с учетом банковского процента // Математическое моделирование экономических систем и процессов. Всероссийская научно-практическая конференция. Тезисы докладов. Чебоксары. 2000, с 60-63.

11. Змеев О.А., Змеева Е.Е. Определение вероятностных характеристик времени разорения страховой компании при условии, что оно произойдет // Статистическая обработка данных и управление в сложных системах. Сб. статей-Томск: Изд-во Том. ун-та, 2000. - Вып. 2. - С. 70-78.

12. Змеев О.А., Терпугов А.Ф. Расчет характеристик времени разорения страховой компании для модели с конечным числом возможных рисков // Экономика, технология, предпринимательство. Сб. статей. Томск. 2000, с. 60-67.

13. Змеев О.А., Моисеев А.Н., Погудин А.А. Использование объектно-ориентированного проектирования для модификации модели базы данных в среднем звене многоуровневой информационной системы.// II научно-практическая конференция «Наука и образование». Тезисы докладов. Белово. 2001, с. 44-46.

14. Змеев О.А., Терпугов А.Ф. Математическая модель функционирования страховой компании с учетом банковского процента.// Методы и алгоритмы прикладной математики в технике, медицине и экономике. Международная научно-практическая конференция. Часть 3. Тезисы докладов. Новочеркасск. 2001, с. 33-36.

15. Змеев О.А. Математические модели страхования в виде дважды стохастического случайного процесса.// Образование и наука в третьем тысячелетии. Третья международная научно- теоретическая конференция. Часть 1. Тезисы докладов. Барнаул. 2001, с. 46.

16. Ахмедова Д.Д., Змеев О.А. Модель функционирования страховой компании с учетом расходов на рекламу при ограниченном числе клиентов // Статистическая обработка данных и управление в сложных системах. Сб. статей - Томск: Изд-во Том. ун-та, 2001. - Вып. 3. - С. 3-14.

17. Змеев О.А. Нестационарный режим в математической модели функционирования страховой компании с учетом банковского процента/// «Новые технологии и комплексные решения: наука, образование производство». Материалы всероссийской научно-практической конференции. Часть II. (Математика). Тезисы докладов. Анжеро-Судженск. 2001, с. 44-45.

18. Змеев О.А. Нестационарный режим для математической модели функционирования страховой компании с интенсивностью входного потока, зависящей от числа клиентов.// Статистическая обработка данных и управление в сложных системах. Сб. статей - Томск: Изд-во Том. ун-та, 2001. - Вып. 3. - С. 4553.

19. Змеев О.А., Масяйкин С.А. Расчет характеристик страховой компании с учетом расходов на рекламу.// «Новые технологии и комплексные решения: наука, образование производство». Материалы Всероссийской научно-практической конференции. Часть II. Тезисы докладов. Анжеро-Судженск. 2001, с. 36-37.

20. Змеев О.А., Моисеев А.Н., Новиков Д.В. Технология хранения объектов в реляционной базе данных.// «Новые технологии и комплексные решения: наука, образование производство». Материалы всероссийской научно-практической конференции. Часть VI. Тезисы докладов, г. Анжеро-Судженск. 2001, с. 19-20.

21. Войтиков К.Ю., Змеев О.А., Моисеев А.Н. Реализация классификаторов на сервере приложений в трехзвенной архитектуре «клиент/сервер» // «Новые информационные технологии в университетском образовании». Тезисы международной научно-методической конференции. Кемерово, 2002, с. 120-122.

22. Гарайшина И.Р., Змеев О.А. Расчет математического ожидания капитала страховой компании с учетом банковского процента в нестационарном режиме // Сборник трудов межрегиональной VI научно-практической конференции «Научное творчество молодежи». Тезисы докладов. Анжеро-Судженск. 2002, с. 41-43.

23. Ахмедова Д.Д., Змеев О.А. Исследование математической модели страховой компании при нестационарном потоке рисков // Математические методы в финансах и экономике. Материалы второй международной конференции «Проблемы актуарной и финансовой математики». Минск, 2002, с. 5-10.

24. Змеев О.А., Масяйкин С.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний // Математические методы в финансах и экономике. Материалы второй международной конференции «Проблемы актуарной и финансовой математики». Минск, 2002, с. 20-24.

25. Ахмедова Д.Д., Змеев О.А. Нахождение характеристик работы страховой компании с нестационарным входящим потоком клиентов // Труды I Всероссийской ФАМ конференции. Ч. II. Красноярск, 2002, с. 20-27.

26. Ахмедова Д.Д., Змеев О.А. Нахождение характеристик капитала страховой компании с нестационарным входящим потоком, параметр которого случайная функция.// Материалы Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование», Анжеро-Судженск, 2002, с 23-26.

27. Змеев О.А., Лезарев А.В. Функциональные требования для систем имитационного моделирования систем массового обслуживания.// Материалы Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование», Анжеро-Судженск, 2002, с 128-130.

28. Змеев О.А., Моисеев А.Н., Новиков Д.В. Решение проблемы переполнения «родительских» таблиц в ООРБД.// Материалы Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование», Анжеро-Судженск, 2002, с 130-131.

29. Ахмедова Д.Д., Змеев О.А. Математическая модель функционирования страховой компании с входящими рисками в виде пуассоновского потока событий с переменной интенсивности.// Обработка данных и управление в сложных системах. Сб. статей - Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 3-12.

30. Войтиков К.Ю., Змеев О.А., Моисеев А.Н. Объектный подход к проблеме проектирования подсистемы нормативно-справочной информации.// Обработка данных и управление в сложных системах. Сб. статей - Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 13-20.

31. Змеев О.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний.// Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2002. - Вып. 4. -С. 33-47.

32. Змеев О.А., Моисеев А.Н., Новиков Д.В. Применение метаданных при проектировании уровня хранения в системах, использующих ООРБД // Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 48-53.

33. Войтиков К.Ю., Змеев О.А., Моисеев А.Н. К вопросу о реализации классификатора в объектно-ориентированной информационной системе // Вестник Кем. гос. ун-та № 2 (10). Кемерово. 2002, с. 167-177.

34. Ахмедова Д. Д., Змеев О. А. Оптимальное управление первоначальным страховым взносом при деятельности страховой компании.// II Всероссийская ФАМ конференция: Тезисы докладов. Красноярск, 2003, с. 57.

35. Змеев О.А. Диффузионное приближение в математической модели фонда социального страхования со случайными расходами на социальные программы.// II Всероссийская ФАМ конференция: Тезисы докладов. Красноярск, 2003, с. 62.

36. Войтиков К.Ю., Змеев О.А., Моисеев А.Н. Основные функциональные требования к подсистеме «Брокер объектных запросов» в рамках унифицированного процесса разработки программного обеспечения.// Обработка данных и управление в сложных системах. Сб. статей - Томск: Изд-во Том. ун-та, 2003. -Вып. 5.-С. 3-13.

37. Змеев О.А. Математическая модель деятельности фонда социального страхования с детерминированными расходами на социальные программы при ре-лейно-гистерезисном управлении капиталом.// Обработка данных и управление в сложных системах. Сб. статей.- Томск: Изд-во Том. ун-та, 2003. - Вып. 5.-С. 42-56.

38. Ахмедова Д.Д., Змеев О.А. Характеристики числа рисков и капитала страховой компании при зависимости среднего первоначального взноса и средних страховых выплат от времени.// Материалы Всероссийской научно-практической конференции «Наука и практика: диалоги нового века» Часть 3. Томск: Изд-во «Твердыня», 2003. С. 23-24.

39. Змеев О.А., Моисеев А.Н., Якушев А.А. Каркас брокера распределенной объектной базы данных // Материалы Всероссийской научно-практической конференции «Наука и практика: диалоги нового века» Часть 3. Томск: Изд-во «Твердыня», 2003. С. 93-95.

40. Змеев О.А., Приступа А.В. Классификация коммерческих систем имитационного моделирования // Материалы Всероссийской научно-практической конференции «Наука и практика: диалоги нового века» Часть 3. Томск: Изд-во «Твердыня», 2003. С. 93-95

41. Ахмедова Д.Д., Змеев О.А. Оптимальное управление первоначальным страховым взносом при деятельности страховой компании // Труды II Всероссийской ФАМ конференции. Ч. II. Красноярск, 2002, с. 21-24.

42. Змеев О.А. Диффузионное приближение в математических моделях фонда социального страхования // Труды II Всероссийской ФАМ конференции. Ч. И. Красноярск, 2002, с. 80-85.

43. Змеев О.А., Приступа А.В. Применение паттернов проектирования при построении систем имитационного моделирования // Материалы VIII Всеросс. научн.-практ. конф. «Научное творчество молодежи». Томск: Изд-во ТГУ, 2004. С. 34-36.

44. Змеев О.А., Приступа А.В. Диаграммы состояния UML как способ представления графа событий имитационной модели системы массового обслуживания.// Обработка данных и управление в сложных системах. Вып. 6. - Томск: Изд-во ТГУ, 2004. - С. 76-81.

45. Ахмедова Д.Д., Змеев О.А. Модель страховой компании с нестационарным потоком входящих рисков и переменной интенсивностью наступления страховых случаев.// Ill Всероссийская ФАМ конференция: Тезисы докладов. Красноярск, 2004, с. 14.

46. Вальц О.В., Змеев О.А. Диффузионная аппроксимация модели фонда социального страхования с релейно-гистерезисным управлением капиталом.// III Всероссийская ФАМ конференция: Тезисы докладов. Красноярск, 2004, с. 17.

47. Pristupa А. V., Zmeyev О. A. 8th Korea - Russia International Symposium on Science and Technology KORUS 2004, Tomsk, RUSSIA, June 26 - July 3, 2004, Vol. l.P. 141-144.

АПРОБАЦИЯ РАБОТЫ

Основные положения диссертации и отдельные её результаты докладывались и обсуждались на:

1. Межрегиональной научно-практической конференции «Наука и образование: пути интеграции» Анжеро-Судженск. 1998.

2. Научно-теоретической конференции «Образование и наука на пороге третьего тысячелетия». Барнаул. 1999.

3. Межрегиональной научно-практической конференции «Качество образования и наука». Анжеро-Судженск. 1999.

4. Межрегиональной научно-методической конференции «Повышение эффективности научных исследований и совершенствование учебного процесса». Анжеро-Судженск. 2000.

5. Всероссийской научно-практической конференции «Математическое моделирование экономических систем и процессов». Чебоксары. 2000.

6. II научно-практической конференции «Наука и образование». Белово. 2000.

7. Международной научно-практической конференции «Методы и алгоритмы прикладной математики в технике, медицине и экономике». Новочеркасск. 2001.

8. III международной научно-теоретической конференции «Образование и наука в третьем тысячелетии». Барнаул, 2001.

9. Всероссийской научно-практической конференции «Новые технологии и комплексные решения: наука, образование производство». Анжеро-Судженск. 2001.

10. Международной научно-методической конференции «Новые информационные технологии в университетском образовании». Кемерово, 2002.

11. Межрегиональной VI научно-практической конференции «Научное творчество молодежи». Анжеро-Судженск. 2002.

12. II международной конференции «Проблемы актуарной и финансовой математики». Минск, 2002.

13. I, II, III Всероссийской ФАМ конференции. Красноярск, 2002,2003,2004.

14. IV Всероссийской конференции с международным участием «Новые информационные технологии в исследовании сложных структур». Томск, 2002.

15. Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование». Анжеро-Судженск, 2002.

16. Всероссийской научно-практической конференции «Наука и практика: диалоги нового века». Анжеро-Судженск, 2003.

17. VIII Всероссийской научно-практической конференции «Научное творчество молодежи». Анжеро-Судженск, 2004.

18. V Всероссийском симпозиуме по прикладной и промышленной математике. Кисловодск,2004.

19. 8th Korea - Russia International Symposium on Science and Technology KORUS 2004, Tomsk, RUSSIA, 2004.

БЛАГОДАРНОСТИ

Автор хотел бы выразить огромную благодарность доктору физ.-мат. наук, профессору А.Ф. Терпугову, доктору технических наук, профессору А.А. Назарову, доктору технических наук, профессору К.И. Лившицу за неоценимую помощь в работе над диссертацией; кандидату технических наук А. Н. Моисееву, кандидатам физико-математических наук Д. Д. Ахмедовой (Даммер) и Е.Е. Змеевой, совместная работа с которыми привела к ряду результатов, изложенных в настоящей диссертации; аспирантам и студентам, которые в разные годы работали под моим руководством и выполняли различного рода исследования, связанные с тематикой настоящей диссертации: О.В. Вальц, К.Ю. Войтикову, И.Р. Гарайшиной, А.В. Ле-зареву, С.А. Масяйкину, Д.В. Новикову, А.А. Погудину, А.В Приступе, А.А. Якушеву.

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

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

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

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

1. Аверилл М. Jloy, В. Кельтон Д! Имитационное моделирование. Классика CS 3-е изд. - СПб.: Питер; Киев: Издательская группа BHV, 2004. - 847 с.

2. Ахмедова Д.Д. Модель функционирования страховой компании с учетом расходов на рекламу // Студент и научно технический прогресс. Раздел экономика. Тезисы докладов. Часть 1 г. Новосибирск - 2000 г. С. 6-8.

3. Ахмедова Д.Д. Оптимизация расходов на рекламу при деятельности страховой компании. // Новые технологии и комплексные решения: наука, образование, производство. Тезисы докладов. 4 2.- Анжеро-Судженск, 2001. С 3-4.

4. Ахмедова Д.Д. Функционирования страховой компании с учетом расходов на рекламу // Повышение эффективности научных исследований и совершенствование учебного процесса. Тезисы докладов. 4 1.- Анжеро-Судженск, 2000. -С 4-6.

5. Ахмедова Д.Д., Змеев О.А. Модель страховой компании с нестационарным потоком входящих рисков и переменной интенсивностью наступления страховых случаев.// III Всероссийская ФАМ конференция: Тезисы докладов. -Красноярск, 2004 г. С. 14.

6. Ахмедова Д.Д., Змеев О.А. Нахождение характеристик работы страховой компании с нестационарным входящим потоком клиентов// Труды I Всероссийской ФАМ конференции. Ч. II. Красноярск, 2002. С. 20-27.

7. Ахмедова Д.Д., Змеев О.А. Оптимальное управление первоначальным страховым взносом при деятельности страховой компании.// II Всероссийская ФАМ конференция: Тезисы докладов. Красноярск, 2003. С. 57.

8. Ахмедова Д.Д., Змеев О.А. Оптимальное управление первоначальным страховым взносом при деятельности страховой компании// Труды II Всероссийской ФАМ конференции. Ч. II. Красноярск, 2002. С. 21-24.

9. Ахмедова Д.Д., Змеев О.А. Оптимизация расходов на рекламу при деятельности страховой компании// Известия вузов. Физика, № 6. 2001, С. 3-7.

10. Ахмедова Д.Д., Змеев О.А., Терпугов А.Ф. Оптимизация деятельности страховой компании с учетом расходов на рекламу.// Вестник Том. гос. ун-та. № 275, 2002, С. 181-184.

11. Ахмедова Д.Д., Терпугов А.Ф. Математическая модель функционирования страховой компании с учетом расходов на рекламу // Известия вузов. Физика. -2001. -№ 1.- С. 25-28.

12. Баруча-Рид А.Т. Элементы теории Марковских процессов и их приложения. -М.: Изд-во «Наука», 1969. 512 с.

13. Белкина Т.А., Фролова А.Г., Чекалина С.В. Исследование динамической модели страхования: различные инвестиционные стратегии и вероятность разорения // Обозрение прикладной и промышленной математики. 2001. - Т. 8. -Вып. 2.-С. 534-535.

14. Белкина Т.А., Фролова А.Г., Чекалина С.В. О минимизации вероятности разорения при выборе инвестиционных стратегий, не использующих заимствований денежных средств // Обозрение прикладной и промышленной математики. 2002. - Т. 9. - Вып. 2. - С. 333.

15. Бенинг В.Е., Королев В.Ю. Асимптотические разложения для квантилей обобщенных процессов Кокса и некоторые их приложения к задачам финансовой и актуарной математики // Обозрение прикладной и промышленной математики. 1998. - Т. 5. - Вып. Г. - С. 22-43.

16. Бенинг В.Е., Королев В.Ю. Асимптотическое поведение обобщенных процессов риска // Обозрение прикладной и промышленной математики. 1998. - Т. 5. - Вып. 1.- С. 117-133.

17. Бенинг В.Е., Королев В.Ю. Асимптотическое разложение для вероятности разорения в классическом процессе риска при малой нагрузке безопасности // Обозрение прикладной и промышленной математики. 2000. - Т. 7. - Вып. 1. -С. 177-179.

18. Беркович В. Чудин А. Практика применения паттернов проектирования// RSDN Magazine. 2003. - №3.

19. Бойков А.В., Мельников А.В. О вероятности неразорения при воздействии экономической среды и инвестирования на рынке Башелье // Обозрение прикладной и промышленной математики. 2002. - Т. 9. - Вып. 2. - С. 340-341.

20. Булинская Е.В. Теория риска и перестрахование ч.1. Упорядочивание рисков. ~М.: ЦПИ, 2001.

21. Буч Г., Рамбо Д., Джекобсон А. Язык UML. Руководство пользователя. М.: ДМК, 2000.-432 с.

22. Вальц О.В., Змеев О.А. Диффузионная аппроксимация модели фонда социального страхования с релейно-гистерезисным управлением капиталза.// Известия вузов. Физика, 2004. № 2. С. 26-31.

23. Вальц О.В., Змеев О.А. Диффузионная аппроксимация модели фонда социального страхования с релейно-Гистерезисным управлением капиталом.// III Всероссийская ФАМ конференция. Тезисы докладов. Красноярск, 2004. С. 17.

24. Вальц О.В., Змеев О.А. Исследование модели фонда социального страхования.// Обозрение прикладной и промышленной математики-2004.- Т 11. Вып. 2.-С. 311-312.

25. Вальц О.В., Змеев О.А. Математическая модель деятельности фонда социального страхования при экспоненциальных страховых выплатах и случайными расходами на социальные программы// Вестник Том. гос. ун-та. № 284, 2004, С. 37-41.

26. Варфоломеев В.И. Назаров С.В. Алгоритмическое моделирование элементов экономических систем. М.: «Финансы и статистика», 2004. - 264 с.

27. Виноградов О.П. Вероятность разорения страховой компании, когда интервалы между выплатами имеют неодинаковые показательные распределения // Теория вероятностей и ее применения. 1998. - Т. 43. - Вып. 2.

28. Виноградов О.П. Об одном элементарном методе получения оценок вероятности разорения // Обозрение прикладной и промышленной математики. 1998. -Т. 5.-Вып. 1.-С. 134-139.

29. Войтиков К.Ю., Змеев О.А., Моисеев А.Н. К вопросу о реализации классификатора в объектно-ориентированной информационной системе.// Вестник Кем. гос. ун-та №2 (10). Кемерово. 2002, С. 167-177.

30. Войтиков К.Ю., Змеев О.А., Моисеев А.Н. Объектный подход к проблеме проектирования подсистемы нормативно-справочной информации.// Обработка данных и управление в сложных системах. Сб. статей Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 13-20.

31. Войтиков К.Ю., Змеев О.А., Моисеев А.Н., Якушев А.А. Архитектура надстраиваемых приложений клиент/сервер с обобщенным протоколом передачи данных.// Вестник Том. гос. ун-та. № 284, 2004, С. 169-173.

32. Галайко Я.В., Назаров А.А. Исследование числа лиц, застрахованных в Пенсионном фонде Российской Федерации при нестационарном входящем потоке // Вестник Том. гос. ун-та. 2003. - № 280. - С. 103-109.

33. Гамма Э., Хемл Р., Джонсон Р. Влиссидес Дж. Приемы объектно-ориентированного проектирования. Паттерны проектирования. СПб: Питер, 2001. -386 с.

34. Гарайшина И.Р. Исследование процесса изменения накопленного капитала пенсионного фонда // Материалы VII Всероссийской научно-практической конференции «Научно творчество молодежи». Томск: Изд-во Том. ун-та, 2004.-С. 22-24.

35. Гарайшина И.Р. Математические модели процесса изменения накопленного капитала пенсионного фонда // Обработка данных и управление в сложных системах: Сб. статей. Томск: Изд-во Том. ун-та, 2004. - Вып. 6. - С. 39-48.

36. Гарайшина И.Р., Назаров А.А. Исследование математической модели изменения страхового капитала Пенсионного фонда // Вестник Том. гос. ун-та. -2003.-№280.-С. 109-112.

37. Гарайшина И.Р., Назаров А.А. Исследование математической модели процесса изменения страхового капитала Пенсионного фонда при стационарном пу-ассоновском входящем потоке страховых взносов. // Известия вузов. Физика, 2004.№2.-С.44-54.

38. Гарайшина И.Р., Назаров А.А. Исследование математической модели процесса изменения накопленного капитала пенсионного фонда при нестационарном пуассоновском входящем потоке// Вестник Том. гос. ун-та. 2004. - № 284. -С. 35-37.

39. Глухова Е.В., Змеев О.А., Лившиц К.И. Математические модели страхования. Томск: Томск: Изд-во Том. ун-та, 2004, 180 с.

40. Глухова Е.В., Капустин Е.В. Расчет вероятности выживания страховой компании с учетом банковского процента при пуассоновском потоке взносов // Известия вузов. Физика. -2001. -№ 6. С. 7-12.

41. Глухова Е.В., Капустин Е.В. Расчет вероятности разорения страховой компании с учетом перестраховки // Известия вузов. Физика. 2000. - № 4. - С. 3-9.

42. Гнеденко Б.В. , Коваленко И.Н. Введение в теорию массового обслуживания. -М.: Наука, 1987.-336 с.

43. Гнеденко Б.В. Курс теории вероятностей. М.: Наука, 1969. - 400 с.

44. Григорьев Ю.Д. Две теоремы об оптимальном уровне удержания для моделей перестрахования Борча и Штрауба. //Труды III Всероссийской ФАМ конференции. Ч. 2. Красноярск, 2004. - С. 94-100.

45. Григорьев Ю.Д. Дележи и оптимизация страхового риска. // IV Всероссийская ФАМ конференция. Тезисы докладов. Красноярск, 2005. - С. 24.

46. Григорьев Ю.Д., Куклин А.В. Вычисление нижних оценок вероятностей разорения в случае логнормального распределения размеров выплат// Математические модели природы и общества: Труды конференции ММПО. -Красноярск, 2002. С. 40-45.

47. Григорьев Ю.Д., Куклин А.В. К вопросу о вычислении вероятности разорения страховой компании в классической модели риска //Труды I Всероссийской ФАМ конференции. Ч. 2. Красноярск, 2002. - С. 99-105.

48. Григорьев Ю.Д., Хекало И.Ю. Что такое оптимальная франшиза? //Математические модели природы и общества: Труды конференции ММПО. -Красноярск, 2002.-С. 51-55.

49. Гунченко К.Г. Прямые методы оценки вероятности разорения страховой компании. // Обозрение прикладной и промышленной математики.-2004.- Т 11. Вып. 2.-С. 322-323. ;

50. Еникеева Ф., Калашников В. Модель риска с инфляцией //Обозрение прикладной и промышленной математики. 1998. - Т. 5. - Вып. 1. - С. 140-147.

51. Змеев О.А. Деятельность фонда социального страхования при релейно-гисте-резисном управлении капиталом.// Математическое моделирование, 2004, т. 16, № 2. -С.43-53.

52. Змеев О.А. Диффузионное приближение в математических моделях фонда социального страхования.// Труды II Всероссийской ФАМ конференции. Ч. И. Красноярск, 2002, С. 80-85.

53. Змеев О.А. Диффузионное приближение в математической модели фонда социального страхования со случайными расходами на социальные программы.// II Всероссийская ФАМ конференция: Тезисы докладов. Красноярск, 2003.-С. 62.

54. Змеев О.А. Математическая модель деятельности фонда социального страхования при экспоненциальных страховых выплатах.// Вестник Том. гос. ун-та, 2003,№280.-С. 130-135.

55. Змеев О.А. Математическая модель фонда социального страхования с детерминированными расходами на социальные программы (диффузионное приближение).// Известия вузов. Физика, № 3. 2003. С. 83-87.

56. Змеев О.А. Математическая модель фонда социального страхования со случайными расходами на социальные программы (диффузионное приближение).// Известия вузов. Физика, № з. 2003, С. 88-93.

57. Змеев О.А. Математические модели страхования в виде дважды стохастического случайного процесса// Образование и наука в третьем тысячелетии. Третья международная научно- теоретическая конференция. Часть 1. Тезисы докладов. Барнаул. 2001. С. 46.

58. Змеев О.А. Математические модели функционирования страховой компании с учетом банковского процента// Известия вузов. Физика. №1.2001. С. 19-25.

59. Змеев О.А. Модель функционирования страховой компании при интенсивности входящего потока, зависящего от числа клиентов// Математическое моделирование. Кибернетика. Информатика. Сб. статей Томск: Изд-во Том. унта, 1999.-С. 67-73.

60. Змеев О.А. Модель функционирования страховой компании при конечном числе возможных клиентов Известия вузов. Физика. № 4. 1999, с. 34-39.

61. Змеев О.А. Определение вероятности разорения страховой компании для модели с конечным числом возможных клиентов// Статистическая обработкагданных и управление в сложных системах. Сб. статей Томск: Изд-во Том. ун-та, 1999.-С. 66-75.

62. Змеев О.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний// Известия вузов. Физика, № 2. 2002, С. 24-28.

63. Змеев О.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний// Обработка данных и управление в сложных системах. Сб. статей Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 3347.

64. Змеев О.А. Расчет характеристик времени разорения страховой компании для моделей с интенсивностью входного потока, зависящей от числа имеющихся рисков. Известия вузов. Физика. № 4.2000, с. 10-16.

65. Змеев О.А., Змеева Е.Е. Модель, функционирования страховой компании при конечном числе возможных клиентов// Наука и образование: пути интеграции. Часть 2. Тезисы докладов. Анжеро-Судженск. 1998, с. 32-33.

66. Змеев О.А., Змеева Е.Е. Расчет-характеристик времени разорения страховой компании для моделей с интенсивностью входного потока, зависящей от числа имеющихся рисков// Качество образования и наука. Тезисы докладов. Анжеро-Судженск. 1999, с. 32-33.

67. Змеев О.А., Лезарев А.В. Шаблон объектного проектирования для реализации функциональности процесса моделирования в имитационных моделях систем массового обслуживания.// Вестник Том. гос. ун-та. № 275, 2002, с. 108-111.

68. Змеев О.А., Моисеев А.Н. Сравнительный анализ некоторых методов 0-R-преобразования.// Вестник Том. г;ос. ун-та, 2003, № 280. С. 263-271.

69. Змеев О.А., Моисеев А.Н. Шаблон диаграммы компонентов информационной системы корпоративного уровня.// Вестник Том. гос. ун-та. № 275, 2002, с. 130-132.

70. Змеев О.А,, Моисеев А.Н., Новиков Д.В. Применение метаданных при проектировании уровня хранения в системах, использующих ООРБД.// Обработка данных и управление в сложных системах. Сб. статей Томск: Изд-во Том. ун-та, 2002. - Вып. 4. - С. 48-53.

71. Змеев О.А., Приступа А.В. Применение паттернов проектирования при по-строении+ систем имитационного моделирования.// Материалы VIII Всеросс. научн.-практ. конф. «Научное творчество молодежи». Томск: Изд-во ТГУ, 2004.-4.1. С. 34-36.

72. Змеев О.А., Приступа А.В. Разработка объектно-ориентированного программного комплекса имитационного моделирования систем массового обслуживания.// Вестник Том. гос. ун-та. № 284,2004. С. 174-176.

73. Змеев О.А., Терпугов А.Ф. Модель функционирования страховой компании при интенсивности потока, зависящего от числа клиентов// Наука и образование: пути интеграции. Часть 2. Тезисы докладов. Анжеро-Судженск. 1998, с. 33-34.

74. Змеев О.А., Терпугов А.Ф. О некоторых подходах к вопросам моделирования деятельности страховых компаний// Образование и наука на пороге третьего тысячелетия: научно-теоретическая конференция. Тезисы докладов. Барнаул. 1999, с. 88-90.

75. Змеев О.А., Терпугов А.Ф. Расчет характеристик времени разорения страховой компании для моделей с конечным числом возможных рисков// Качество образования и наука. Тезисы докладов. Анжеро-Судженск. 1999, с. 33-34.

76. Змеев О.А., Терпугов А.Ф. Расчет характеристик времени разорения страховой компании для модели с конечным числом возможных рисков// Экономика, технология, предпринимательство. Сб. статей Томск. 2000, с. 60-67.

77. Калашников В., Константинидис Д. Вероятность разорения // Фундаментальная и прикладная математика. 1996. - Т. 2. -№ 4. - С. 1005-1100.

78. Калашников В.В., Цициашвили Г.Ш. Асимптотически точные двухсторонние оценки вероятности разорения при наличии больших выплат //Обозрение прикладной и промышленной математики. 1998. - № 5. - Вып. 1. - С. 66-82.

79. Кац В.М. О некоторых подходах к вопросам моделирования деятельности страховых компаний // Рыночная экономика России в XXI веке: Сб. трудов. -Томск: Изд-во «SPRINT», 2001. ^ С. 238-240.

80. Кац В.М. Основные принципы моделирования деятельности страховых компаний в условиях рынка // Вторая областная научно-практическая конференция студентов, аспирантов и молодых ученых: Сб. статей. Томск: Изд-во «Графика», 2001. - С. 86-87.

81. Кац В.М., Лившиц К.И. Влияние расходов на рекламу на характеристики страховой компании // Известия вузов. Физика. 2001. - Т. 44, № 1. - С. 2833.

82. Кац В.М., Лившиц К.И. Конкурентное взаимодействие двух страховых компаний на ограниченном страховом рынке // Вестник Томск, гос. ун-та. Приложение. -2002. № 1 (I). - С. 163-166.

83. Кац В.М., Лившиц К.И. О конкурентном взаимодействии двух страховых компаний на ограниченном страховом рынке // Обозрение прикладной и промышленной математики. 2002. - Т. 9. - Вып. 2. - С. 389.

84. Кац В.М., Лившиц К.И. Характеристики страховой компании при малой нагрузке страховой премци // Вестник Томск, гос. ун-та. Приложение. 2002. -№1(1).-С. 159-163.

85. Кац В.М., Лившиц К.И., Назаров А.А. Исследование нестационарных беско-нечнолинейных систем массового обслуживания и их применение к анализу экономико-математических моделей // Вестник Томск, гос. ун-та. 2002. - № 275.-С. 189-192.

86. Кокс Д., Смит В. Вероятность разорения. М: Советское радио, 1967. -350 с.

87. Ларман К. Применение UML и шаблонов проектирования, введение в объектно-ориентированный анализ и проектирование.: Пер. с англ.: Уч. пос. — М.: Издательский дом «Вильяме», 2001. 496 с.

88. Лезарев А.В. Имитационная модель деятельности автотранспортного предприятия.// Материалы Всероссийской научно-практической конференции «Информационные технологии и математическое моделирование», Анжеро-Судженск, 15 ноября 2002, с 204-206.

89. Лившиц К.И. Вероятность разорения страховой компании для пуассоновской модели//Известия вузов. Физика. -1999. Т. 42.-№ 4. - С. 28-33.

90. Лившиц К.И., Сухотина Л.Ю. Математическая модель страховой компании с учетом сезонных изменений // Обработка данных и управление в сложныхгсистемах: Сб. статей. Томск: Изд-во Том. ун-та, 2004. - Вып. 6. - С. 118— 126.

91. Лябах Н.Н., Бутакова М.А. Алгоритм формального описания структуры системы массового обслуживания. // Обозрение прикладной и промышленной математики. 2004. - Т. 11. - Вып. 2. - С. 366-367.

92. Льюс Р.Д., РайфаХ. Игры и решения. М.: Изд-во ин. лит., 1961. - 642 с.

93. Масяйкин С.А. Построение переговорного множества при конкурентном взаимодействии двух страховых компаний //Математические модели природы и общества: Труды конференции ММПО -2002. -Красноярск, 2002. С. 133138.

94. Маталыцкий A.M., Романюк Т.В. Исследование стохастической модели обработки разнотипных исков страховой компании // Материалы второй международной конференции «Проблемы актуарной и финансовой математики» Минск 2002.-С. 70-75.

95. Медведев Г.А. Математические модели финансовых рисков. -Минск, БГУ, 2001.-291 с.

96. Назаров А.А. Асимптотический анализ марковизируемых систем. Томск: Изд-во Том. ун-та, 1991. - 158 с.

97. Ослин Б.Г. Имитационное моделирование систем массового обслуживания -Томск.: Изд-во ТПУ, 2003. 106 с.

98. Оуэн Г. Теория игр. -М.: Мир, 1971. 230 с.

99. Первозванский А.А. Курс автоматического управления. М.: Наука, 1986. -615с.

100. Прабху Н.У. Стохастические процессы теории запасов. М: Мир, 1984. - 184 с.

101. Радюк Л.Е., Терпугов А.Ф. Теория вероятностей и случайных процессов. -Томск: Изд-во Том. ун-та, 1988. 174 с.

102. Ротарь В.И., Бенинг В.Е. Введение в математическую теорию страхования// Обозрение прикладной и промышленной математики. -1994. Т. 1, вып. 5. -С.699-779.

103. Рыжиков Ю.И. Имитационное моделирование. Теория и технологии. СПб.: «Корона», 2004.-384 с.

104. Семенов А.Т. Граничные задачи и процессы разорения // Труды II Всероссийской ФАМ конференции. Часть И. Красноярск, 2003. - С. 154-160.

105. Семенов А.Т. Фундаментальное творчество Вальда и процессы разорения // Труды II Всероссийской ФАМ конференции. Часть II. Красноярск, 2003. - С. 161-166.

106. Семенов А.Т. Асимптотика вероятности разорения в модели Спарре Андерсона. // Труды III Всероссийской ФАМ конференции. Часть II. Красноярск, 2004.-С. 202-206.

107. Семенов А.Т. Асимптотические/разложения вероятности разорения для процесса Леви. // IV Всероссийская ФАМ конференция. Тезисы докладов. Красноярск, 2005. - С. 29-30.

108. Скотт К. Основы UML М. Мир, 1999.

109. Такач JI. Комбинаторные методы в теории случайных процессов. М.: Изд-во «Мир», 1971.-263 с.

110. Тенцер А. База данных хранилище объектов// Компьтер-пресс. 2001. - № 8.

111. Терпугов А.Ф. Математика рынка ценных бумаг. Томск: Изд-во. HTJI, 2004. - 164 с.

112. Терпугов А.Ф. Экономико-математические модели. Томск: Изд-во Том. унта, 1999. - 185 с.

113. Фаулер М. Архитектура корпоративных программных приложений.: Пер. с англ. М.: Издательский дом «Вильяме», 2004. - 544с.

114. Феллер В. Введение в теорию вероятностей и ее приложения / Пер. с анг.: В 2 т.-М.: Мир, 1967.-Т. 1.-498 с.

115. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисленияМ.: Наука, 1966. Т. I. - 608 с. .

116. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисленияМ.: Наука, 1966. Т. II. -800 с.

117. Харкевич А.Д., Лившиц Б.С., Пшеничников А.П. Теория телетрафика. М.: Связь, 1979.-224 с.

118. Черепанова И.В. Математические методы программы перестрахования. // IV Всероссийская ФАМ конференция. Тезисы докладов. Красноярск, 2005. - С. 35-36.

119. Штрауб Э. Актуарная математика имущественного страхования. -Цюрих, 1988,-147 с.

120. Шон Л.Д. Оптимальный уровень удержания при пропорциональном перестраховании, минимизирующий верхнюю границу вероятности разорения. // IV Всероссийская ФАМ конференция. Тезисы докладов. Красноярск, 2005. - С. 27-28. '

121. Шон Л.Д. Выбор оптимального уровня удержания в схеме эксцедентного перестрахования. // IV Всероссийская ФАМ конференция. Тезисы докладов. -Красноярск, 2005. С. 28-29. >

122. Эльсгольц Л.Э. Дифференциальгые уравнения и вариационное исчисление -М: Изд-во. Наука, 1969. 158 с.

123. Эмбрехтс П. Актуарный и финансовый подходы к расчету стоимости в страховании //Обозрение прикладной и промышленной математики. -1996. — Т. 5. -Вып. 1.-С. 6-22.

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

125. Язык бизнеса. Термины / Под ред. В. А. Чжена. М.: ИПК «Шарк», 1995. -736 с.

126. A.Buss. Component Based Simulation Modeling. Proceedings of the 2000 Winter Simulation Conference, ed., J.A. Joinses, R.R. Barton, K. Kang and P.A. Fishwick. Institute of Electrical and Electronics Engineers, Piscataway, New Jersey, U.S.A.

127. Alexander C., Ishikawa S., Silverstein M., Jacobson M., Fiksdahl-King I., Angel S. A Pattern language. Oxford University Press, New York, 1977.

128. Asmussen S. On the Ruin Problem for Some Adapted Premium Rules. Probabilistic Analysis of Rare Events: Theory .and Problems of Safety, Insurance and Ruin. -Riga, 1999,-P. 3-15.

129. Bagrodia R.L. Perils and Pitfalls of Parallel Discrete-Event Simulation, 1996 Winter Simulation Conference, San Diego, p. 136-143, 1996.

130. Banks J, Carson J.S, Nelson B.L. Discrete-Event System Simulation, 2d ed., Prentice-Hall, Upper Saddle River, New Jersey, 1996. p. 128-138.

131. Booch G. Object-Oriented Analysis and Design with Applications. Benjamin/Cum-mings, Redwood City, CA, 1994. Second Edition.

132. Bowden R.O. The Spectrum of Simulation Software, HE Solution, 30, May 1998: p. 44-46.

133. Bower N., Gerber H., Hickman J. et al. Actuarial Mathematics. Shaumburg: Society of Actuaries, 1997.

134. Buschmann F., Meunier R., Rohnert H., Sommeral P., Stal M. System of patterns. West Sussex, England: Wiley, 1996.

135. CACI Products Company: MODSIM III User Manual, La Jolla, California, 1997.

136. Carathers C.D., Topol В., Fujimoto R.M., Sunderam V. Visualizing Parallel Simulations in Network Computing Environments: A Case Study, Proc. 1997 Winter Simulation Conference, Atlanta, p. 110-117, 1997.

137. Chandrasekaran U., Sheppard S. Discrete Event Distributed Simulation A Survey, Proc. Conference on Methodology and Validation, Orlando, p. 32-37, 1987.

138. Chandy K.M., Misra J. Asynchronous Distributed Simulation via a Sequence of Parallel Computations, Commun. Assoc. Comput. Mach., 24: 198-206, 1981.

139. Chandy K.M., Misra J. Distributed Deadlock Detection, Assoc. Comput. Mach. Trans. Computer Systems, 1: 144-156, 1983.

140. Chandy K.M., Misra J. Distributed Simulation: A Case Study in Design and Verification of Distributed Programs, IEEE trans. Software Eng., SE-5:440-452, 1979.

141. Comfort J.C. The Simulation of a Master-Slave Event-Set Processor, Simulation, 42: 117-124, 1984.

142. Deutsch L. Peter Design reuse and frameworks in the Smalltalk-80 system. In Ted J. Biggerstaff and Alan J. Perlis, editors, Software Reusability, Vol. II: Applications and Experience, Addison-Wesley Reading, MA, 1989, h. 57-71.

143. Fishwick P. A. Web-Based Simulation, Proc. 1997 Winter Simulation Conference, Atlanta, p. 100-102, 1997.

144. Fishwick P.A. Web-Based Simulation: Some Personal Observations, Proc. 1996 Winter Simulation Conference, San Diego, p. 772-779,1996.

145. Fowler M. UML Distilled, Reading, MA: Addison-Wesley, 1997.

146. Fujimoto R.M. Parallel and Distributed Simulation, in Handbook of Simulation, J. Banks, ed., John Wiley, New York, 1998.

147. Fujimoto R.M. parallel Discrete Event Simulation: Will the Field Survive? INFORMS J. Cornput., 5: 213-230, 1993.

148. Gerber H.U. An Introduction to Mathematical Risk Theory. Wharton School. University of Pennsylvania, 1979. 164 p.

149. Golubin A.Y. Design of a Pareto-optimal insurance policy in the individual riskrmodel. // VI International Petrozavodsk Conference « Probabilistic methods in Discrete Mathematics». Abstracts of Communications. Part 1.2004. P. 271 -272.

150. Grandell J. Aspect of Risk Theory. New York: Springer-Verlag, 1991.

151. Grandell J. Risk and geometric sums // Applied Stochastic Models and Inform. Processes. Petrazavodsk, 2002. P. 52-53.

152. Grandell J. Simple Approximation of Ruin Probabilities. Probabilistic Analysis of Rare Events: Theory and Problems of Safety, Insurance and Ruin. Riga, 1999, -P. 47-51.

153. Gupta U.G. Using Citation Analysis to Explore the intellectual base, Knowledge dissemination and Research Impact of Interfaces 91970-1992), Interfaces, 27: 85101, 1997.

154. Heidelberger P. Discrete Event Simulations and Parallel processing: Statistical Properties, SIAM J. Statist. Comput., 9: 1114-1132, 1988.

155. Jacobson I., Christerson M., Jonsson P., Overgaard G. Object-Oriented Software Engineering A Use Case Driven Approach. Addison-Wesley, Wokingham, England, 1992.

156. Jefferson D.R. Virtual Time, Assoc. Comput. Mach. Trans. Programming Languages and System, 7:404-425, 1985.

157. Jones D.W., Roberts S.D. Object-Oriented Simulation. Handbook of Simulation, J. Banks ed., John Wiley, New York <1998).

158. Joohnson R.E., Foote brain Designing reusable classes. Journal of Object-Oriented Programming, 1(2): June/July 1988. p. 22-35.

159. Katz V.M., Livshits K.I. Optimization of Advertising Expenses in the Functioning of an Insurance Company // Applied Stochastic Models and Inform. Processes. Petrazavodsk, 2002. P. 82-84.

160. Konstantinides D.G., Tang Q.H., Tshitsiashvili G.Sh. Two-sides Bounds for Ruin Probability under Constant Interest Force // Applied Stochastic Models and Inform. Processes. Petrazavodsk, 2002. P. 98-101.

161. Lavenberg S., Muntz R., Samadi B. Performance Analysis of Rollback Method for Distributed Simulation, Performance '83, A.K. Agrawala and S. K. Tripathi, eds., 1983.

162. Levasseur J.A. The Case for Object-Oriented Solution, OR/MS Today August 1996, 23: p. 65-67.

163. Lorenz P., Dorwarth H., Ritter K.C., Schriber T.J. Towards a Web Based Simulation Environment, Proc. 1997 Winter Simulation Conference, Atlanta, p. 13381344, 1997.

164. Malinovskii V.K. Approximations and upper bounds an probabilities of large deviations in the problem of ruin within finite time // Scand. Actuarial J., 1996. P. 124147.

165. Malinovskii V.K. Probabilities of Ruin when the Safety Loading Tends to Zero. Laboratory of Actuarial Mathematics University of Copenhagen. Working Paper № 153, 1998.-36 p.

166. Marti J. Object-Oriented Modeling and Simulation with MODSIM III, CACI Products Company, La Jolla, California, 1999.

167. Misra J. Distributed Discrete-Event Simulation, Computing Surveys, 18: 39-65, 1986.

168. Nicol D.M. Principles of Conservative Parallel Simulation, Proc. 1996 Winter Simulation Conference, San Diego, p. 128-135,1996.

169. Nicol D.M., Johnson M.M., Yoshimura A.S. The IDES Framework: A Case Study in development of a Parallel Discrete-Event Simulation System, Proc. 1997 Winter Simulation Conference, Atlanta, p. 93-99, 1997.

170. Norberg R. Sensitivity Analysis in Insurance and Finance // Applied Stochastic Models and Inform. Processes. Petrazavodsk, 2002. P. 121-124.

171. Panjer H.Y., Willmont G.E. Insurance Risk Models. Society of Actuaries, 1992. -442 p.

172. Pristupa A.V., Zmeyev O.A. 8th Korea Russia International Symposium on Science and Technology KORUS 2004, Tomsk, RUSSIA, June 26 - July 3,2004, Vol. 1. P. 141-144.

173. Rumbaugh J. The life of an object model: How the object model changes during development. Journal of Object-Oriented Programming, 7 (1): 24-32, March/April, 1994.

174. Rumbaugh J., Blaha M., Premerlani W., Eddy F., Lorenson W. Object-Oriented Modeling and Design. Prentice Hall, Englewood Cliffs, 1991.

175. Sargent R.G. Event Graph Modeling for Simulation with an Application to Flexible Manufacturing System. Management Sci., 34, 1988: p. 1231-1251.

176. Schmidli H. Asymptotics of ruin probabilities for controlled risk processes in the small claim case // Applied Stochastic Models and Inform. Processes. Petrazavodsk, 2002.-P. 141-142.

177. Schruben L.W. Simulation Modeling with event Graphs. Commun. Assoc. Comput. Mach., 26, 1983.: p. 957-963.

178. Shanker M.S., Padman R., Kelton W.D. Efficient Distributed Simulation through Dynamic Load balancing, HE Trans., 2000.

179. Sheppard S., Yuong R. E., Chandrasekaran U., Krishnamurthi M. Three Mechanisms for Distributing Simulation, Рос. 12th Conference ofNSF Production research and technology program, Madison, Wisconsin, p. 67-70, 1985.

180. Som Т.К., Sargent R.G. A Formal Development of Event Graphs as an Aid to Structured and Efficient Simulation Programs. ORSA J. Comput., 1, 1989, p. 107125.

181. Tshitsiashvili G.Sh. Quality Properties of Risk Models Under Stochastic Interest Force // Applied Stochastic Models and Inform. Processes. Petrazavodsk, 2002. P. 149-152.1. АКТоб использовании результатов диссертационной работы О.А. Змеева

182. Исследование математических моделей страхования при нестационарныхпотоках страховых рисков» (

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

184. Оценка потенциального рынка автострахования в рамках страхового поля, ограниченного территорией г. Анжеро-Судженска;

185. Исследование возможности оказания дополнительных страховых услуг, связанных с медицинским обслуживанием населения города.

186. Предложенные в работе решения позволили менеджменту компании профессионально оценить текущую ситуацию и принять соответствующие решения о предоставлении населению соответствующих услуг.2112.04

187. ООО «Губернская страховая компания Кузфа^са», v директор представительства по г. Анжеро*Сужденску,-Г

188. ГОСУДАРСТВЕННОЕ УЧРЕЖДЕНИЕ -КУЗБАССКОЕ1. РЕГИОНАЛЬНОЕ ОТДЕЛЕНИЕ

189. ФОНДА СОЦИАЛЬНОГО СТРАХОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ1. ФИЛИАЛ №1

190. Об использовании результатов диссертационной работы О.А.Змеева «Исследование математических моделей страхования при нестационарных потоках страховых рисков»

191. Использование этих моделей позволило оценить

192. Качественные и количественные параметры деятельности фонда за 20022003 годы

193. Провести оптимизацию процесса планирования расходов на финансирование социальных программ на 2004 год

194. Предложенные в работе решения благотворно повлияли на качество обслуживания населения города и заслуживают дальнейшего внимания.652470, Кемеровская область, г. Анжеро-Судженск, ул. Ленина, 32 тел/факс (253) 6-11-66на №от1. АКТ1. В.Е.Герасимович