Вычисление вероятности разорения страховой компании в случае выплат, имеющих экспоненциальное распределение со сдвигом Текст научной статьи по специальности «Математика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2012

Управление, вычислительная техника и информатика

№ 4(21)

УДК 519.2

Е.В. Капустин

ВЫЧИСЛЕНИЕ ВЕРОЯТНОСТИ РАЗОРЕНИЯ СТРАХОВОЙ КОМПАНИИ В СЛУЧАЕ ВЫПЛАТ, ИМЕЮЩИХ ЭКСПОНЕНЦИАЛЬНОЕ РАСПРЕДЕЛЕНИЕ СО СДВИГОМ

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

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

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

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

Рассмотрим следующую модель страховой компании: поступление капитала в компанию за счет страховых взносов детерминировано и имеет постоянную скорость с, страховые выплаты образуют пуассоновский поток постоянной интенсивности X, размеры страховых выплат независимы и имеют одинаковое распределение с плотностью п(х) и начальным моментом 1-го порядка т1. В литературе такая модель называется классической моделью страховой компании [6-8].

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

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

1. Постановка задачи

(1)

0, если х < х0,

2. Вычисление вероятности разорения страховой компании

Обозначим Р(^) вероятность разорения компании при уровне капитала £. Известно [6, 7], что Р(£) удовлетворяет уравнению

£ +ад

| Р(£ - х)п(х)ёх + | п(х)ск

сР' (£) -ХР(Б) + Х и граничному условию

= 0 (2)

Иш Р(£) = 0. (3)

£ ——^

Применяя операционный метод [9], можно показать [6, 7], что если выполняется условие нормального функционирования компании

с > %т, (4)

то задача (2), (3) имеет единственное ограниченное решение, причем его изображение (преобразование Лапласа)

Р (р) = | Р( £)е“ PSdS

имеет вид

Р (р)=I-------с-кт—. (5)

р Ср -Х(1 -П (р))

где п(р) - изображение п(х).

Чтобы найти вероятность разорения компании Р(£), нужно найти оригинал функции Р(р). Если страховые выплаты имеют экспоненциальное распределение со сдвигом (1), то

т, = х0 +9 ,

п (р) = -

1 + 9р

(4) и (5) принимают вид

с > Х(х0 + 9), (6)

1 (л Х( х0 +9)

Р(р) = --11 I-------—----------—. (7)

р \ с ) XX е рх°

р------+ --

с с 1 + 9р

Раскладывая дробь, стоящую в правой части (7), в ряд по степеням е имеем

рх0

Р(р) =1 -(1 -х(х° +9)+ £(-1)"(АТ-------1 1 е-прх0 I, (8)

р I с Др-а "=1 ' Iс9) (р-а)п+1 (р + Р)п | ' '

X в 1

где а = —, в= -.

с 9

При достаточно большом х0 ряд в (8) сходится в области Яе р > х0 (точка р = а в эту область не входит), то есть правая часть (8) является изображением [9]. Вычисляя обратное преобразование Лапласа [9], получаем

Х( х0 +9)

г+х (-!)” і ^

) = 1 -|1 -

.. Г С-

X

^ "=1

где 1( х) - единичная функция [9]

если х > 0.

1(С-"о) Г .(9)

■{0,

1(х) 1 0, если х < 0.

Заметим, что при фиксированном £ сумма в правой части (9) содержит лишь конечное число слагаемых.

3. Асимптотическое поведение вероятности разорения страховой компании при больших значениях капитала

Исследуем поведение функции (9) при больших значениях £ Известно [7, 9], что при £ — да

Р(£)~ Р(р)ер£, (10)

р=И)

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

Известно [6, 7], что особенность функции (5) в точке р = 0 устранима. Это используется при выводе формулы (5), но может быть легко проверено непосредственно.

Представим (7) в виде

Р(р) = 1 - (с -х(х0 +9))(1 + 9р) (11)

р (1 + 9р)( ср -X) + Хе~ рх°

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

г = pxo, к = —, к2 = X” • х^0 XXo

Тогда знаменатель второй дроби, стоящей в правой части (11), принимает вид

X/(г), где

/(г) = (1 + V)(к2 г -1) + е~г, (12)

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

к2 > 1 + к .

Из курса математического анализа имеем, что если к2 > 1 + к1, то функция (12)

имеет два вещественных нуля, тривиальный г = 0 и отрицательный. Обозначим

этот нуль х1 (см. рис. 1).

f (*)*

xl^—' 0

X

Рис. 1. График функции f (x)

Покажем, что если к2 > 1 + кг, то функция (12) не имеет нулей в полуплоскости Яе г > 0 и полосе х < Яе г < 0 , имеет единственный нуль 2 = 0 на прямой Яе г = 0 и единственный нуль 2 = х1 на прямой Яе 2 = х1.

Достаточно рассмотреть верхнюю полуплоскость 1т г > 0. Пусть г = х + іу . Тогда уравнение / (г) = 0 равносильно системе уравнений

|"к1к2(х2 - у2) + (к2 - к1)х -1 + е~х соє у = 0,

[2к1к2ху + (к2 - к1)у - е~х єіп у = 0.

Докажем, что в области х > х1, у > 0 система (13) решений не имеет.

Пусть х > 0, у > 0 . Представим второе уравнение в системе (13) в виде

поэтому равенство (14) выполняться не может. Следовательно, точка (х, у) не является решением системы (13).

Пусть х < х < 0, у > 0. Тогда f (х) < 0 (см. рис. 1), поэтому

kjk2(х2 - у2) + (k2 - k1)х -1 + e~х cos у < k1k2 х2 + (k2 - k1)х -1 + e- = f (x) < 0 ,

то есть в системе (13) не выполняется первое равенство, точка (х,у) решением системы (13) не является.

Таким образом, наибольшую вещественную часть среди особых точек функ-

х

ции (7), кроме точки р = 0 , имеет точка р0 = —. Это полюс первого порядка, по-

х0

этому [9]

ex (2kjk2 X + k2 - kj) = Siny.

(14)

y

Из предположений и курса математического анализа имеем

Следовательно, при £ ^-да

(с - Х( х0 + 9))( х0 + 9х) х15

р(8)------V------>_о_—ж_о-----V— ех„ , (15)

Хх0 (9 + х0 е Х1) - с(х0 + 29х1)

где х < 0 - отрицательный нуль функции (12).

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

с = 10, Х = 2, 9=1, х0 = 2.

Рис. 2. Вероятность разорения страховой компании и ее асимптотика Заключение

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

ЛИТЕРАТУРА

1. Глухова Е.В., Капустин Е.В. Расчет вероятности разорения страховой компании с учетом перестраховки // Изв. вузов. Физика. 2000. № 4. С.3-9.

2. Глухова Е.В., Капустин Е.В. Расчет вероятности разорения страховой компании с учетом перестраховки при пуассоновском потоке страховых взносов // Статистическая обработка данных и управление в сложных системах: сб. статей. Томск: Изд-во Том. ун-та, 2000. Вып. 2. С. 34-46.

3. Глухова Е.В., Капустин Е.В., Терпугов А.Ф. Расчет вероятности разорения страховой компании с учетом банковского процента // Методы и алгоритмы прикладной математики в технике, медицине и экономике: Материалы междунар. научно-практич. конф. Юж.-Рос. гос. техн. ун-т. Новочеркасск: НАБЛА, 2001. Ч. 3. С. 36-38.

4. Капустин Е.В. Расчет вероятности разорения для модели страховой компании, учитывающей возможность одновременного наступления нескольких страховых случаев // Обработка данных и управление в сложных системах: сб. статей / под ред. А.Ф. Терпугова. Томск: Изд-во Том. ун-та, 2004. Вып. 6. С. 188-195.

5. Капустин Е.В., Михайленко М. С. Модель страховой компании с пуассоновским потоком взносов и с учетом издержек, равномерных по времени // Информационные технологии и математическое моделирование (ИТММ-2010): Материалы IX Всероссийской научнопрактической конференции с международным участием (19-20 ноября 2010 г.). Томск: Изд-во Том. ун-та, 2010. Ч. 1. С. 13-16.

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

7. Panjer H.H., Willmot G.E. Insurance Risk Models. Society of Actuaries, Schaumburg, Ill. 1992. 442 p.

8. ШтраубЭ. Актуарная математика имущественного страхования. М.: Мир, 1988. 146 с.

9. ЛаврентьевМ.А., Шабат Б.В. Методы теории функций комплексного переменного. М.: Наука, 1987. 688 с.

Капустин Евгений Викторович

Филиал Кемеровского государственного университета в г. Анжеро-Судженске

Е-mail: kapustin@asf.ru Поступила в редакцию 5 мая 2012 г.

Kapustin Eugene V. (Anjero-Sudjensk branch of the Kemerovo Staty University). Calculation of the ruin probability of an insurance company when claim size distribution is exponential with a shift.

Keywords: risk model, the ruin probability, asymptotic behavior.

The paper considers the model of the insurance company when premium rate is constant, flow of claims is generated by Poisson process (the classical risk model) and claim sizes have probability density

The ruin probability of the company and its asymptotic behavior for large values of initial capital are obtained by operational calculus.