CC BY
79
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА 2017 Управление, вычислительная техника и информатика № 39
УДК 519.2
DOI: 10.17223/19988605/39/4
К.И. Лившиц, А.А. Назаров
ПРОСТАЯ АППРОКСИМАЦИЯ ВЕРОЯТНОСТИ РАЗОРЕНИЯ СТРАХОВОЙ КОМПАНИИ ДЛЯ МОДЕЛИ КРАМЕРА-ЛУНДБЕРГА СО СТОХАСТИЧЕСКИМИ ПРЕМИЯМИ
Предлагается и исследуется оценка вероятности разорения страховой компании для модели Крамера-Лундберга со стохастическими премиями. Проводится сравнение предлагаемой оценки как с точными выражениями, так и с известной ранее асимптотической оценкой вероятности разорения.
Ключевые слова: вероятность разорения; модель Крамера-Лундберга со стохастическими премиями; аппроксимация вероятности разорения.
Стандартной задачей актуарной математики является задача вычисления вероятности разорения страховой компании, т.е. вероятности ситуации, когда страховая компания не может исполнять свои финансовые обязательства ввиду отсутствия денежных средств при различных предположениях о потоках, поступающих в компанию страховых премий и страховых выплат, производимых страховой компанией [1-3]. Сложность состоит, как правило, в нахождении явных решений соответствующих инте-гро-дифференциальных уравнений, определяющих вероятность разорения. Поэтому представляет интерес нахождение оценок для вероятности разорения, для построения которых не нужно решать соответствующие уравнения. В настоящей работе предлагается и исследуется оценка вероятности разорения страховой компании для модели Крамера-Лундберга со стохастическими премиями [4, 5].
1. Модель Крамера-Лундберга со стохастическими премиями
Модель Крамера-Лундберга со стохастическими премиями строится при следующих предположениях: предполагается, что поток страховых премий, поступающих в компанию, является пуассонов-ским с интенсивностью X, премии - независимые, одинаково распределенные случайные величины с
плотностью распределения Ф(x) и моментами ak = M{xk} [Там же]. Страховые выплаты также образуют пуассоновский поток с интенсивностью ц, выплаты - независимые, одинаково распределенные случайные величины с плотностью распределения x) и моментами bk = M {xk }. Предполагается
также, что с начала деятельности страховой компании прошло достаточно длительное время, так что поток страховых выплат не зависит от потока страховых премий.
Пусть в момент времени t капитал компании равен S(t). При сделанных выше предположениях
m(t) n(t)
S(t) = S(0) +X x -X У, , (1)
i=1 j=1
где S(0) - начальный капитал, m(t) - число поступивших страховых премий до момента времени t, n(t) - число страховых выплат, xi - величина i -й страховой премии, у, - величина j -й страховой
выплаты. Определим T = minjt :S(t) < 0} и T = да, если S(t) > 0 Vt. Случайная величина T - момент
разорения. Тогда вероятность предельного разорения страховой компании, при условии, что ее капитал в начальный момент равен S,
P(S) = Pr {T < да}. (2)
Можно показать [4, 5], что вероятность разорения Р(5) (2) определяется уравнением
(X + ц) Р(5) = Х|™ Р(5 + х )Ф( х)с1х + " Р( X - х)¥( х )с1х + "[¥( х )^х (3)
0 X
с граничным условием Р(+м) = 0 . Для существования отличного от Р(Х) = 1 решения уравнения необходимо выполнение условия
Ха, = (1+9) "А, (4)
где 9> 0 - нагрузка страховой премии. Точное решение уравнения (3) удается найти лишь в нескольких частных случаях, например при экспоненциальных или гиперэкспоненциальных распределениях страховых премий и выплат [4, 6].
Известны следующие оценки вероятности разорения Р(Х). Оценка вероятности разорения сверху [4, 5]. Обозначим
со со
ЬФ (ю) = | е-ю5 Ф( X , (ю) = | е-ю5 ¥(5 . о о
Тогда если уравнение
имеет положительный корень к, то Во-вторых, при 9 много меньше 1 [7]:
ХЬФ (к) + цЬр (к) = Х + ц
Р(5) < е-к5. (5)
ц ехр(--2"+Ь— 95)
р (5) --Ха2"+цЬ-. (6)
ОО
2"Ь1
-9х
Х+ц-Х|е х°2+"Ь2 Ф(х)ёх о
Точность оценки (6) тем выше, чем меньше значение параметра 9.
2. Простая оценка вероятности разорения
Следуя [8], назовем «простой» оценкой вероятности разорения такую ее аппроксимацию, для построения которой используются только моменты ак и Ьк распределений страховых премий и страховых выплат и не требуется строить оценку решения уравнения (3). Для классической модели Крамера-Лундберга известна оценка вероятности разорения Эе УуЫег [8, 9], которая основана на идее замены истинного процесса риска с произвольным распределением страховых выплат на процесс риска, в котором страховые выплаты имеют экспоненциальное распределение, а параметры аппроксимирующего процесса подобраны так, чтобы три первых момента у истинного и аппроксимирующего процессов риска совпадали. За оценку вероятности разорения при этом принимается вероятность разорения, соответствующая аппроксимирующему процессу. Воспользуемся этой идеей.
Пусть страховые премии и страховые выплаты имеют экспоненциальные распределения
1 - 1 -Ф(5) = -е а, ¥(5) = -еЬ. (7)
а Ь
Тогда вероятность разорения [4]
Ха-цЬ 5
Р(5) = "£±Ь) е"(Х+Ц)аЬ . (8)
(X + ц)а
Или, учитывая (4),
" -95
Р( 5) = "(а + Ь) ~ (Х+") а (X + ц)а
Обозначим
тк(0 = М}, к = 1,4 . (9)
Считая для простоты, что 5"(0) = 0 , несложно показать, что для процесса 5"(^) (2) имеют место соотношения
т^Г) = (Аа1 -уЬ^, (10)
т2(Г) = (Аа2 + уЬ2У + (Аа1 -уЬх)2 Г2, (11)
т3(^) = (Аа3 -уЬ3У + 3(Аа2 + уЬ2)(Аах -уЬ1)^2 + (Аах - уЬ1)313, (12)
3
т4^) = (Аа4 + уЬ4У + (2(Аах - уЬх)(Аа3 - уЬ3) + 3(Аа2 + уЬ2)2 +—(Аа3 - уЬ3)(Аах - уЬх))^2 +
+ 3(Аа2 +уЬ2)(Аа1 -уЬх)213 + (Аах -уЬх)414. (13)
При экспоненциальных распределениях страховых премий и страховых выплат (7) имеем, очевидно,
ак = к!ак, Ьк = к!Ьк.
Если интенсивности потоков страховых премий и страховых выплат в аппроксимирующем процессе риска равны А0 и у0 соответственно, то для совпадения моментов тк (^) двух процессов достаточно выполнения условий
= Аа1 -уЬ = А0а -у0Ь, (14)
d2 = -2-(Аа2 + уЬ2) = А 0 а 2 +у 0Ь 2, (15)
й3 =—(Аа3 -уЬ3) = А0а3 -у0Ь3, (16)
6
d4 = -^(Аа,^ +уЬ4) = А0 а4 + у0Ь 4. (17)
В результате получаем систему уравнений (14)-(17) на определение неизвестных параметров аппроксимирующего процесса а,Ь,А0,у0 . Из уравнений (14), (15)
= , у0 = . (18)
0 а (а + Ь) Ь(а + Ь)
Подставляя (18) в уравнения (16) и (17), получим
41аЬ + 42(а -Ь) = 43, (19)
[41(а - Ь) + 42]аЬ + 42(а - Ь)2 = . (20)
Обозначим
V = а - Ь, и = аЬ . (21)
Из уравнений (19) и (20) получим
4 - 4- 4244 /'-»'-»ч
V = -, и =-(22)
4 22 4143 4 2
Наконец, соотношения (21) дают
V + л/ V2 + 4и , -V + -\/V2 + 4и .„„.
а =-, Ь =-. (23)
2 2
Оценка вероятности разорения будет иметь вид
А0а-у0Ь ^ р(£) = у0(а + Ь) е (А0аЬ
(А0 +у0)а '
где параметры а, Ь, А 0, у0 вычисляются по формулам (18), (23).
Так как параметры а и Ь должны быть положительны, то очевидно, что для применимости построенной оценки должны выполняться условия
d 4 d 2 d 2 d 4
<0, -Ч4-> 0. (25)
d3 d 22 d3 d 22
3. Примеры
Пример 1. Распределения Эрланга страховых премий и страховых выплат. Пусть страховые премии и страховые выплаты имеют распределения Эрланга
Ф(5) = Ае"а, ¥(5) = ^е Р . (26)
а2 Р2
Решение уравнения (3) будем искать в виде
Р(5) = А1е-^ + А2е-^ . (27)
Подставляя выражение (27) в уравнение (3), получим, что для этого должны выполняться условия
A1 , A2
= -1, (28)
PYi-1 PY2-1
=1, (29)
A1 , A2
(PY1-1)2 (РУ2-1)2 X ■ Ц = Х + ц, i = 1,2. (30)
(аУг + 1)2 (PYí -1)2
Таким образом, величины yi должны быть положительными корнями уравнения
f(z) = -^-7 + 7P-E1)2 - ^ = 0. (31)
(az +1)2 (Pz -1)2
Функция f (z) обладает следующими свойствами: f (0) = 0, f '(0) = - (Xa - ц—) < 0, f (—= да,
a
f (р) = да, f (±да) = -(X + ц), откуда следует, что уравнение f (z) = 0 имеет два положительных корня:
у1 е (0, р) и у 2 е (р, да). Решая систему уравнений (28), (29), получаем, что вероятность разорения в нашем случае имеет вид
P(S) = У2(1-Ру1)2 e~y1s - у1(1-Ру2)2 e~y>s . (32)
У 2 -У1 У 2 -У1
На рис. 1 приведены зависимость вероятности разорения P( S) от начального капитала, вычисленная по формуле (32), и ее оценки, построенные по формулам (6) и (24) соответственно. Параметры
a = 0,5, р = 5, ц = 1, X = (1 + 9)Р . Параметр 9 = 0,05; 0,1; 0,5; 1.
a
Пример 2. Постоянные страховые премии и страховые выплаты. Пусть страховые премии и страховые выплаты постоянны:
Ф(S) = 5(S-a), ¥(S) = 8(S-b). Точное решение уравнения (3) получить в этом случае не удается. В таблице приведены оценка вероятности разорения, полученная в результате имитационного моделирования, и оценки вероятности
разорения, построенные по формулам (6) и (24). Параметры a = 1, b = 10, ц = 1, X = (1 + 9) — . Параметр
a
9 = 0,05; 0,1; 1. Объем выборки при моделировании - 50 000. Отметим, что оценка вероятности разорения, полученная путем имитационного моделирования, всегда является заниженной из-за конечности
времени моделирования процесса S(V) (1). Как видно из приведенных данных, при всех значениях 9 оценка (24) оказывается лучше, чем оценка (6). Выигрыш растет с увеличением нагрузки страховой премии.
- 1 0 = 0.05
5 \ N X. \ 9=0.1
5 Хл ч NN 9 = 0.5
ч ч X \ 9=1
1
°о 5 ю 15 5
Рис. 1. Зависимость вероятности разорения от начального капитала (сплошные линии) и ее оценки, полученные по формуле (6) (штриховые линии) и формуле (24) (штрихпунктирные линии) при распределениях Эрланга страховых премий и выплат
Оценки зависимости вероятности разорения от начального капитала при постоянных страховых премиях и страховых выплатах
5 0 1 2 3 4 5 6 7 8 9 10
Результат модел. е = 0,05 0,952 0,948 0,942 0,934 0,93 0,925 0,917 0,907 0,9 0,893 0,884
Оценка (24) 0,939 0,931 0,922 0,914 0,906 0,898 0,891 0,883 0,875 0,867 0,86
Оценка (6) 0,914 0,905 0,897 0,889 0,881 0,873 0,865 0,858 0,85 0,842 0,835
Результат модел. е = 0,1 0,91 0,899 0,891 0,88 0,871 0,858 0,846 0,834 0,817 0,799 0,781
Оценка (24) 0,886 0,871 0,856 0,842 0,828 0,814 0,8 0,786 0,773 0,76 0,747
Оценка (6) 0,836 0,821 0,806 0,792 0,778 0,764 0,75 0,737 0,724 0,711 0,698
Результат модел. е = 1 0,699 0,648 0,619 0,598 0,569 0,543 0,512 0,478 0,439 0,404 0,367
Оценка (24) 0,631 0,589 0,549 0,512 0,477 0,445 0,415 0,387 0,361 0,336 0,314
Оценка (6) 0,445 0,408 0,374 0,343 0,314 0,288 0,264 0,242 0,222 0,203 0,186
Пример 3. Гиперэкспоненциальное распределение страховых премий и выплат. Пусть страховые премии и страховые выплаты имеют гиперэкспоненциальные распределения
т п
Ф(5) = £ Лк акв, 5) = £ Бк ркв-Рк8. (33)
к=1 к=1
Тогда вероятность разорения Р(5) имеет вид [6]:
Р(5) = £Ре ^, (34)
]=1
где - положительные корни уравнения
f (z) = xjr ^Мь--х-ц = 0!
k=i ak + z k=i Pk - z
a Pj - решение системы уравнений
n 1 1
P/=-
м Pk-У/ j Pk
(36)
OS
\ \ 9 = 0.1
\\ \ \\ \
\ ^ \\ \ V \ к
\ \ V" ч ^ X -
0.8
0.6
<u
o.:
9 = 1
i\ 1 \ 1 \ \; \
\i. \ \ \ \ \
\ 4 \
100
150
Рис. 2. Зависимость вероятности разорения от начального капитала (сплошные линии) и ее оценки, полученные по формуле (6) (штриховые линии) и формуле (24) (штрихпунктирные линии) при гиперэкспоненциальных распределениях страховых премий и выплат
На рис. 2 приведены вероятность разорения, вычисленная по формулам (34)-(35) при т = п = 2, и ее оценки, полученные по формулам (6) и (24) при 0 = 0,1 и 0 = 1. Параметры
а = 1, Ь = 10, ц = 1, X = (1 + 0) — . Как видно из рис. 2, в этом случае обе оценки являются плохими, но и в
а
этом примере оценка (24) ведет себя «чуть лучше», чем оценка (6). Формально использовать оценку (24) в рассматриваемой ситуации нельзя, так как параметры аппроксимирующего процесса, вычисленные по формулам (23) и (18), отрицательны. Например, при 0 = 0,1 а = -0,405, Х0 = -16,762. Однако комбинации величин, непосредственно входящие в оценку (24), положительны.
Заключение
В работе получены расчетные формулы, позволяющие вычислить оценку вероятности разорения страховой компании при произвольных распределениях страховых выплат и страховых премий для модели Крамера-Лундберга со стохастическими премиями. Как следует из рассмотренных примеров, во всех рассмотренных случаях точность предлагаемой оценки лучше, чем точность известной ранее оценки вероятности разорения. В принципе тот же подход может быть использован для вычисления вероятностей разорения и для более сложных моделей, например, для случая, когда и моменты поступления страховых премий образуют ММР поток [7].
ЛИТЕРАТУРА
1. Panjer H.Y., Wilmot G.E. Insurance Risk Models. Society of Actuaries, 1992. 442 p.
2. Бауэрс Н., Гербер Х., Джонс Д., Несбит С., Хикман Д. Актуарная Математика. М. : Янус-К, 2001. 656 с.
3. Глухова Е.В., Змеев О.А., Лившиц К.И. Математические модели страхования. Томск : Изд-во Том. ун-та, 2004. 180 с.
4. Livshits K.I. Probability of Ruin of an Insurance Company for the Poisson Model // Russian Physics Journal. 1999. V. 42, No. 4.
P. 394-399.
5. Бойков А.В. Модель Крамера-Лундберга со стохастическими премиями // Теория вероятностей и ее применение. 2002. Т. 47,
вып. 3. С. 549-553.
6. Лившиц К.И., Сухотина Л.Ю. Вероятность разорения страховой компании при гиперэкспоненциальных распределениях
страховых премий и страховых выплат для различных моделей страхования // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2016. № 2 (35). С. 37-45.
7. Лившиц К.И., Бублик Я.С. Вероятность разорения страховой компании при дважды стохастических потоках страховых пре-
мий и страховых выплат // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. № 4 (17). С. 64-73.
8. Grandell J. Simple approximation of ruin probabilities // Insurance: Mathematics and Economics. 2000. V. 26. P. 157-173.
9. De Vylder F.E. A practical solution to the problem of ultimate ruin probability // Scandinavian Actuarial Journal. 1978. P. 114-119.
Лившиц Климентий Исаакович, д-р техн. наук, профессор. E-mail: kim47@mail.ru Назаров Анатолий Андреевич, д-р техн. наук, профессор. E-mail: nazarov.tsu@gmail.com Национальный исследовательский Томский государственный университет
Поступила в редакцию 10 декабря 2016 г.
Livshits Klimenty I., Nazarov Anatoly A. (Tomsk State University, Russian Federation).
Simple approximation of an insurance company ruin probability for Cramer-Lundberg model with stochastic premiums. Keywords: probability of ruin; Cramer-Lundberg model with stochastic premiums; ruin probabilities approximation.
DOI: 10.17223/19988605/39/4
In this paper the calculated formulas which allow to construct the simple approximation for the ruin probability of an insurance company for Cramer-Lundberg model with stochastic premiums are obtained. The approximation is based on the replacement of the true risk process with arbitrary distributions of the insurance premiums and the insurance payments
m(t) n (t)
s (t) = S (0) + X x "X yj,
i=1 j=1
where S(0) is the initial capital, m(t) is the number of received insurance premiums to the time t, n(t) is the number of insurance payments, xt is the value of i -th insurance premiums, y. is the value of j-th insurance payment, ak = M ^xik } , bt = M |yik } , the flow of insurance premiums is Poisson with intensity X, the flow of insurance payments is Poisson with intensity on the risk process with the insurance premiums and the insurance payments having the exponential distributions. At that the premiums flow intensives X0, the payments flow intensives |0, the average values of the premiums a and the payments b in the approximating process are chosen so that the first four moments of the true and approximating risk processes coincided. This requires the fulfillment of the conditions
d =Xaj "|bj = X0a "|0b, d2 = 2(Xa2 +|b2) = X0a2 + |0b2,
d3 = — (Xa3 - |b3) = X0a3 - |0b3, d4 = —(Xa4 + |b4) = X0a4 + |0b4. 6 24
The parameters a, b, X0, |0 are determined by the ratios
d1b + d2 d 2 - d1a
Xq =-, I0 =-,
a(a + b) b(a + b)
v Wv2 + 4u -v Wv2 + 4u
where
di d 4 d 2 dl di d 2d 4
v = -, u =—-=-2- ■
dj d3 d 2 dj d3 d 2
For the rating of the ruin probability is taken the ruin probability of the approximating process
p(S) = + b) e (^0+^0)ab .
(Xo +^o)a REFERENCES
1. Panjer, H.Y. & Wilmot, G.E. (1992) Insurance Risk Models. Society of Actuaries. Schaumburg, IL.
2. Bowers, N.L., Gerber, H.U., Hekman, J.C., Jones, D.A. & Nesbitt, C.J. (1998) Actuarial Mathematics. Society of Actuaries, Michigan.
3. Gluhova, E.V., Zmeev, O.A. & Livshits, K.I. (2004) Matematicheskie modeli strakhovaniya [Mathematical Models of Insurance].
Tomsk: Tomsk State University.
4. Livshits, K.I. (1999) Probability of Ruin of an Insurance Company for the Poisson Model. Russian Physics Journal. 42(4). pp. 394-
399. DOI: 10.1007/BF02509675
5. Boykov, A.V. (2002) Cramer-Lundberg Model with Stochastic Premiums. Theory of Probability and Its Applications. 47(3). pp. 549-
553. DOI: 10.1137/S0040585X9797987
6. Livshits, K. I. & Suhotina, L. Yu. (2016) Ruin Probability of an Insurance Company with Hyperexponential Distribution of Insurance
Premiums and Insurance Payments for Different Insurance Models. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(35). pp. 37-45. (In Russian). DOI: 10.17223/19988605/35/4
7. Livshits, K.I. & Bublic, Ya.S. (2011) Ruin probability of an insurance company under double stochastic flows of insurance premium
and insurance payments. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika -Tomsk State University Journal of Control and Computer Science. 4(17). pp. 64-73. (In Russian).
8. Grandell, J. (2000) Simple approximation of ruin probabilities. Insurance: Mathematics and Economics. 26. pp. 157-173. DOI:
10.1016/S0167-6687(99)00050-5
9. De Vylder, F.E. (1978) A practical solution to the problem of ultimate ruin probability. Scandinavian Actuarial Journal. pp. 114-119.
DOI: 10.1080/03461238.1978.10419484
