close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2010
Управление, вычислительная техника и информатика
№ 4(13)
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
УДК 519.865
К.И. Лившиц, Я.С. Бублик
РАСПРЕДЕЛЕНИЕ УСЛОВНОГО ВРЕМЕНИ
ДО РАЗОРЕНИЯ СТРАХОВОЙ КОМПАНИИ ПРИ ДВАЖДЫ
СТОХАСТИЧЕСКОМ ПОТОКЕ СТРАХОВЫХ ВЫПЛАТ1
Найдена производящая функция условного времени до разорения страховой
компании при дважды стохастическом потоке страховых выплат и малой нагрузке страховой премии.
Ключевые слова: условное время до разорения, вероятность разорения,
дважды стохастический поток, нагрузка страховой премии
В [1] было получено выражение для вероятности разорения страховой компании на бесконечном временном интервале при дважды стохастическом потоке
страховых выплат и малой нагрузке страховой премии. Однако более исчерпывающей характеристикой деятельности компании является распределение условного времени до ее разорения при условии, что разорение происходит [2]. Через
распределение условного времени до разорения может быть выражена, в частности, вероятность разорения на конечном временном интервале. В настоящей работе находится распределение условного времени до разорения страховой компании
при дважды стохастическом потоке страховых выплат и малой нагрузке страховой премии.
1. Математическая модель страховой компании
Как и в [1] , будем считать, что интенсивность потока страховых выплат λ(t )
является однородной цепью Маркова с непрерывным временем и n состояниями
λ(t ) = λ i [3]. Переход из состояния в состояние задаётся матрицей инфинитезимальных характеристик Q = ⎡ q ij ⎤ ранга n − 1 . Таким образом, переход из со⎣ ⎦
стояния i в состояние j за малое время ∆t имеет вероятность
Pij (∆t ) = qij ∆t + ο(∆t ) , i ≠ j ,
(1)
Pii (∆t ) = 1 + qii ∆t + ο(∆t ) , i = 1, n ,
где qij ≥ 0 при i ≠ j и
n
∑ qij = 0 .
(2)
j =1
1
Работа выполнена в рамках аналитической ведомственной целевой программы «Развитие научного
потенциала высшей школы» (2009 – 2010 годы), проект № 4761.
К.И. Лившиц, Я.С. Бублик
16
Обозначим Pi ( t ) = P {λ ( t ) = λi } , i = 1, n . Если управляющая цепь является неразложимой, то существуют финальные вероятности
πi = lim Pi ( t ) ,
t →∞
которые являются решением системы уравнений
n
∑ q ji π j = 0 ;
(3)
π1 + π2 + ... + πn = 1 .
(4)
j =1
Обозначим через λ 0 среднюю интенсивность потока страховых выплат в стационарном режиме
n
λ 0 = ∑ λ i πi .
(5)
i=1
Будем считать далее, что страховые выплаты являются независимыми случайными величинами с плотностью распределения ψ ( x ) , средним значением
M { x} = a и вторым моментом M {x 2 } = a2 .
Наконец, в соответствии с классической моделью страховой компании [4] будем считать, что страховые премии поступают непрерывно во времени с постоянной скоростью c , так что за время ∆t приращение капитала за счет страховых
премий равно c∆t . Пусть S ( t ) – средний капитал компании в момент времени t .
Как показано в [1], при t 1
S ( t ) = S ( 0) + ( c − λ0 a ) t .
(6)
Из выражения (6) следует, что при t 1 капитал компании в среднем монотонно возрастает, если
c = (1 + θ ) λ 0 a ,
(7)
где θ > 0 . При θ < 0 компания разоряется. Параметр θ , как и в классической модели, – нагрузка страховой премии.
2. Производящие функции условного времени
Пусть ( Ω, F , F = ( Ft )t ≥ 0 P ) – вероятностное пространство, на котором определены траектории процесса S ( t ) изменения капитала компании. Пусть в начальный момент времени капитал компании равен s и значение интенсивности
λ ( t ) = λ i . Разобьем все возможные траектории процесса S ( t ) , выходящие из
этой точки на два класса: {Sω ( t ) , ω ∈ Ωi ( s )} – траектории, приводящие к разоре-
{
}
нию, и Sω ( t ) , ω ∈ Ωi ( s ) – траектории, приводящие к выживанию. Пусть ti ( s, ω)
– время до разорения на траектории, приводящей к разорению. Обозначим
Φ i ( s, u ) =
∫
Ωi ( s )
e −uti ( s ,ω) P ( d ω)
(8)
Распределение условного времени до разорения страховой компании
и пусть
Pi ( s ) =
∫
P ( d ω)
17
(9)
Ωi ( s )
– вероятность разорения на бесконечном интервале при условии, что в начальный
момент времени капитал равен s и интенсивность потока равна λ i . Тогда
ϕi ( s, u ) ϕi ( s, u ) =
Φ i ( s, u )
Pi ( s )
(10)
есть производящая функция моментов условного времени до разорения при условии, что в начальный момент времени капитал равен s и интенсивность потока
равна λ i .
Функции Φ i ( s, u ) должны удовлетворять граничным условиям
Φ i ( s, 0 ) = Pi ( s ) ,
lim Φ i ( s, u ) = 0,
s →∞
(11)
второе из которых вытекает из того, что при s → ∞ область интегрирования в (8)
Ωi ( s ) → ∅ .
Пусть Fi ( t , s ) – функция распределения условного времени до разорения при
условии, что в начальный момент времени капитал равен s и интенсивность потока равна λ i . Если Pi ( s, t ) – вероятность разорения страховой компании за время t , то
Pi ( s, t ) = Fi ( s, t ) Pi ( s ) ,
так как для разорения за время, не превосходящее t , компания должна разориться
и время до разорения должно быть не больше чем t . Поэтому определение любой
из вероятностей Pi ( s, t ) или Fi ( t , s ) определяет и вторую вероятность.
Для вывода уравнений, которым должны удовлетворять функции Φ i ( s, u ) ,
рассмотрим два соседних момента времени t и t + ∆t . За время ∆t капитал компании изменится на величину ∆s и
ti ( s, ω) = ∆t + t j ( s, ω) ,
(12)
где номер j соответствует значению интенсивности λ ( t ) в момент времени t + ∆t .
Усредняя соотношение (12), будем иметь
Φ i ( s, u ) = e −u∆t M ∆s , j {Φ j ( s + ∆s, u )} .
(13)
При принятой модели за время ∆t могут произойти следующие события:
1. С вероятностью (1 − λi ∆t )(1 + qii ∆t ) + o ( ∆t ) интенсивность потока не меняется, страховые выплаты не производятся.
2. С вероятностью λ i ∆t ψ ( x ) dx (1 + qii ∆t ) + o ( ∆t ) интенсивность потока не меняется и производится случайная страховая выплата размера x .
3. С вероятностью (1 − λi ∆t ) qij ∆t + o ( ∆t ) происходит изменение интенсивности потока с λ i на λ j , страховая выплата не производится.
Остальные события имеют вероятность o ( ∆t ) .
К.И. Лившиц, Я.С. Бублик
18
Используя формулу полной вероятности, получим из (13)
s
⎡
Φ i ( s, u ) = e−u∆t ⎢(1 − ( λi − qii ) ∆t ) Φ i ( s + c∆t , u ) + λ i ∆t ∫ Φ i ( s + c∆t − x ) ψ ( x ) dx +
⎢⎣
0
∞
⎤
+λi ∆t ∫ ψ ( x ) dx + ∑ qij ∆t Φ j ( s + c∆t , u ) ⎥ + o ( ∆t ) ,
⎥⎦
i≠ j
s
где учтено, что при s < 0 ti ( s, ω) = 0. Считая функции Φ i ( s, u ) дифференцируемыми и переходя к пределу при ∆t → 0 , получим
c
s
n
∂Φ i ( s, u )
= ( λi + u ) Φ i ( s, u ) − ∑ qij Φ j ( s, u ) − λ i ∫ Φ i ( s − x, u )ψ ( x ) dx −
∂s
j =1
0
∞
−λi ∫ ψ ( x ) dx.
(14)
s
Получить точное решение систем уравнений (14) не удается. Поэтому рассмотрим далее асимптотический случай, когда нагрузка страховой премии θ 1 .
Решение системы уравнений (14) будем искать в виде
⎛ θs u ⎞
Φ i ( s, u ) = A ( u , θ ) fi ⎜
, ,θ⎟,
(15)
⎝ 1 + θ θ2 ⎠
где A ( u, θ ) и fi ( z , u , θ ) – некоторые пока не определенные функции. В силу произвольности функции A ( u, θ ) можно считать, что
n
∑ πi fi ( 0, u, θ ) = 1.
(16)
i =1
Также будем считать, что функции fi ( z , u , θ ) по крайней мере дважды дифференцируемы по z и равномерно непрерывны по u и θ , а также, что существует
конечный предел
(
)
lim A θ2u , θ ≠ 0 .
θ→0
Подставляя (15) в уравнения (14) и сделав замену переменной
u
θ2
θ
s = z,
1+ θ
= v получим уравнения относительно функций fi ( z , v, θ )
(
)
n
λ 0 aθfi ( z , v, θ ) = λ i + θ2 v fi ( z , v, θ ) − ∑ qij f j ( z , v, θ ) −
j =1
∞
θ
⎛
⎞
−λi ∫ fi ⎜ z −
x, v, θ ⎟ ψ ( x ) dx + r ( θ ) ,
1
+
θ
⎝
⎠
0
(17)
где
r ( θ ) = λi
∞
∞
λi
⎛ z (1 + θ )
2 ⎞
Φ
−
x
,
θ
v
ψ
x
dx
−
(
)
⎜
⎟
i
∫ ⎝ θ
∫ ψ ( x ) dx .
⎠
A θ2 v, θ z (1+θ )
z (1+θ )
(
θ
)
θ
(18)
Распределение условного времени до разорения страховой компании
19
( )
Можно показать, что при сделанных предположениях r ( θ ) = o θ2 .
Обозначим
fi ( z , v ) = lim fi ( z , v, θ ) .
(19)
θ→ 0
Переходя в (17) к пределу при θ → 0 , получим, что
n
∑ qij f j ( z, v ) = 0 .
(20)
j =1
Так как по условию Rang ⎡⎣ qij ⎤⎦ = n − 1 , то из сравнения систем уравнений (2) и
(20) получаем, что
fi ( z , v ) = f ( z , v ) ∀i ,
(21)
где f ( z , v ) – неизвестная пока функция.
Представим теперь функции fi ( z , v, θ ) в виде
fi ( z , v, θ ) = f ( z , v ) + Bi ( z , v ) θ + o ( θ ) .
(22)
θ
⎛
⎞
Подставляя разложения (22) в уравнения (17), раскладывая f ⎜ z −
x, v ⎟ и
⎝ 1+ θ
⎠
θ
⎛
⎞
Bi ⎜ z −
x, v ⎟ в ряд Тейлора и ограничиваясь членами, имеющими порядок θ ,
⎝ 1+ θ
⎠
получим, переходя к пределу при θ → 0 , что
n
∑ qij B j ( z, v ) = ( λi − λ0 ) af ( z, v ) .
(23)
j =1
n
Так как одновременно
∑ πi qij = 0
n
и
i =1
∑ πi (λi − λ0 ) = 0 , то система уравнений
i=1
(23) совместна и имеет ранг (n − 1) , как и система уравнений (3).
−1
Пусть матрица R = ⎡⎣ Rij ⎤⎦ = ⎡⎣ qij ⎤⎦ , i, j = 1, n − 1 . Тогда решение системы (23)
имеет вид [1]
n −1
Bk ( z , v ) = Bn ( z , v ) + ∑ Rkj ( λ j − λ 0 ) a f ( z , v ) .
(24)
j =1
Представим теперь функции fi ( z , v, θ ) в виде
( )
fi ( z , v, θ ) = f ( z , v ) + Bi ( z , v ) θ + Ci ( z , v ) θ2 + o θ2 .
(25)
θ
⎛
⎞
Подставляя разложения (25) в уравнения (17), раскладывая f ⎜ z −
x, v ⎟ ,
⎝ 1+ θ
⎠
θ
θ
⎛
⎞
⎛
⎞
Bi ⎜ z −
x, v ⎟ и Ci ⎜ z −
x, v ⎟ в ряд Тейлора и ограничиваясь членами,
⎝ 1+ θ
⎠
⎝ 1+ θ
⎠
имеющими порядок θ2 , получим, переходя к пределу при θ → 0 и учитывая (23),
К.И. Лившиц, Я.С. Бублик
20
n
∑ qij C j ( z, v ) − ( λi − λ0 ) aBi ( z, v ) + λi
j =1
a2
f ( z , v ) + λ i a f ( z , v ) − vf ( z , v ) = 0 .
2
(26)
Умножая соотношения (26) на πi , суммируя и учитывая (4), будем иметь
n
−∑ πi ( λ i − λ 0 ) aBi ( z , v ) + λ 0
i =1
a2
f ( z , v ) + λ 0 a f ( z , v ) − vf ( z , v ) = 0 .
2
(27)
Учитывая (24), из (27) получим
A1 f ( z , v ) + A2 f ( z , v ) − vf ( z , v ) = 0 ,
A1 =
где
n −1
n −1
λ 0 a2
− a 2 ∑ ( λ i − λ 0 ) πi ∑ Rij ( λ j − λ 0 ),
2
i =1
j =1
(28)
A2 = λ 0 a .
(29)
Откуда
f ( z , v ) = U1 ( v ) e X1 ( v ) z + U 2 ( v ) e X 2 ( v ) z ,
где
X1 ( v ) =
− A2 − A22 + 4 A1v
,
2 A1
X 2 (v) =
− A2 + A22 + 4 A1v
2 A1
(30)
– корни характеристического уравнения дифференциального уравнения (28).
Будем считать, что A1 > 0 . Можно показать, что это условие выполняется по
крайней мере при n = 1,3 . Так как Φ i ( s, u ) Pi ( s ) – производящая функция,
то Φ i ( s, u ) Pi ( s ) ≤ 1 . Re u и, следовательно, Re v > 0 . Отсюда следует, что
f ( z , v ) → 0 при z → ∞ . Поэтому U 2 ( v ) = 0 . Наконец, из условия (16) получаем,
что U1 ( v ) = 1 . Таким образом,
f ( z, v ) = e
−
A2 + A22 + 4 A1v
z
2 A1
Φ i ( s, u ) = A ( u , θ ) e
и
⎛ u ⎞ θ
X1 ⎜ 2 ⎟
s
⎝ θ ⎠1+θ
(31)
+ O (θ) .
(32)
При выводе соотношения (32) неявно предполагалось, что в уравнениях (14)
s ≠ 0 . Для определения функции A ( u, θ ) рассмотрим теперь уравнения (14) при
s = 0 . Умножая уравнения системы (14) на πi и складывая уравнения, при s = 0
будем иметь
n
c ∑ πi
i =1
∂Φ i ( 0, u ) n
= ∑ πi ( λ i + u ) Φ i ( 0, u ) − λ 0 .
∂s
i =1
(33)
Подставляя в (33) выражения (32), получим, что
−1
⎡
⎛ u ⎞⎤
A ( u, θ ) = λ 0 ⎢λ 0 + u − λ 0 aθx1 ⎜ 2 ⎟ ⎥ .
⎣
⎝ θ ⎠⎦
Таким образом, окончательно получаем, что
−1
(34)
⎛ u ⎞ θ
X1 ⎜ 2 ⎟
s
⎡
⎛ u ⎞⎤
Φ i ( s, u ) = λ 0 ⎢λ 0 + u − λ 0 aθX1 ⎜ 2 ⎟ ⎥ e ⎝ θ ⎠1+θ + O ( θ ) .
⎣
⎝ θ ⎠⎦
(35)
Распределение условного времени до разорения страховой компании
21
Из соотношений (35) и (9) получаем теперь выражение для вероятностей разорения на бесконечном интервале:
⎡ aA
Pi ( s ) = ⎢1 + 2
A1
⎣
θA2
−1
s
⎤ −
θ ⎥ e (1+θ ) A1 + O ( θ ) .
⎦
(36)
s
Пусть n = 1 и ψ ( s ) =
1 −a
e
– показательное распределение. Так как при
a
этом A1 = λ 0 a 2 , A2 = λ 0 a , то выражение (32) принимает вид
θs
1 − (1+θ ) a
e
.
1+ θ
Построенная аппроксимация приводит, таким образом, в этом случае к истинной вероятности разорения [ 2] .
Моменты условного времени до разорения определяются соотношениями
P (s) =
tik ( s ) =
( −1)k ∂ k Φ i ( s, u )
Pi ( s )
∂u k
(37)
u =0 .
Вычисляя производные, получим, что при θ 1 среднее значение условного
времени
s+a
(38)
ti1 ( s ) =
+ O (1) .
λ 0 aθ
а дисперсия условного времени
2A (s + a)
Di ( s ) = 1 3 3 + O (1) .
(39)
θ A2
Отметим, что выражение (38) является точным при n = 1 и показательном
распределении выплат [ 2] .
3. Плотность распределения условного времени
при нулевом начальном капитале
Из соотношений (35) и (36) вытекает, что производящая функция условного
времени до разорения (10) при нулевом начальном капитале
1 + 2α
1
ϕi ( 0, u ) =
+ O (θ) ,
λ 0β χ + p + p
где
p=
u
λ 02β
+
α2
λ 02β
,
2
α=
aA2 θ
,
2 A1
β=
a2
,
A1
χ=
λ 0 (1 + α ) β − α 2
λ 02β2
.
Плотность распределения условного времени до разорения при нулевом начальном капитале gi ( t , 0 ) определится как обратное преобразование Лапласа от
ϕi ( 0, u ) . Или
gi ( t , 0 ) =
λ 02 (1 + 2α ) β
2 π
t
∫
0
x
(t − x )
3
2
e
−
( λ 0β x + 2 α ( t − x ) ) 2
−λ 0 x
4β( t − x )
dx + O ( θ ) .
(40)
К.И. Лившиц, Я.С. Бублик
22
При n = 1 и экспоненциальном распределении страховых выплат выражение
(40) принимает вид
g (t, 0) =
λ λ (1 + θ ) t
2 π
∫
0
x
(t − x )
3
2
e
−λx −
λ ( x +θ( t − x ) )2
4( t − x )
dx + O ( θ ) .
(41)
Истинная плотность распределения p(t , 0) в случае экспоненциального распределения страховых выплат имеет вид [ 2]
(
)
p(t , 0) = 1 + θe −λ ( 2+θ )t I 1 2λ 1 + θt ,
где I1 ( z ) – модифицированная функция Бесселя первого порядка. Зависимость
плотностей распределения gi ( t , 0 ) и p ( t , 0 ) от t приведена на рис. 1. Параметр
θ = 0,5 .
g, p
λ=1
10–1
g
λ=5
p
10–2
10–3
0
1
3
2
t
4
Рис. 1
4. Плотность распределения условного времени
при неограниченно возрастающем начальном капитале
Будем теперь считать, что при θ → 0 начальный капитал компании
s → ∞ согласованно с уменьшением θ . Более точно, будем считать, что
3
lim θs ( θ ) = ∞ , но
θ→ 0
lim θ 2 s ( θ ) = 0 . Обозначим
θ→0
m=
1
,
θA2
σ=
2 A1
3
θ
A23
,
и введем величину
z=
t − ms
σ s
.
γ=
2 A1
A2
Распределение условного времени до разорения страховой компании
23
Производящая функция величины z
m
ϕi , z ( u , s ) = e σ
su
⎛ u
⎞
ϕi ⎜
,s⎟
⎝σ s ⎠
или
⎛ aA ⎞
2 γ 2 θs
⎪⎧
⎛ u ⎞ θs ⎪⎫
λ 0 ⎜1 + 2 θ ⎟ exp ⎨ γu θs +
+ X1 ⎜ 2
⎬
⎟
1+ θ
A1 ⎠
⎝ θ σ s ⎠ 1 + θ ⎪⎭
⎪⎩
⎝
ϕi , z ( u , z ) =
+ O (θ).
u
⎛ u ⎞
λ0 +
− λ 0 aθX1 ⎜ 2
⎟
⎝θ σ s ⎠
σ s
Далее,
u
⎛ u ⎞
2
2
X1 ⎜ 2
=
⎟ = −γ − γ 1 + 2
⎝θ σ s ⎠
γ θs
⎡
u
1 u2
1
u3
⎛ 1 ⎞⎤
= −γ 2 ⎢ 2 +
−
+
+ o⎜
⎟⎥ .
2
3
⎝ θs θs ⎠ ⎦
γ θs 2 γ θs 2 γ θs θs
⎣
(42)
(43)
Подставляя разложение (43) в выражение (42) и переходя к пределу при
θ → 0 , получим, что
lim ϕi , z ( u , s ) =
θ→0
Таким образом при θ 1 и θs
ски нормальное распределение.
u2
e2
.
1 случайная величина z имеет асимптотиче-
Заключение
В работе найдена производящая функция условного времени до разорения
страховой компании при дважды стохастическом потоке страховых выплат и дополнительном предположении о малости нагрузки страховой премии. Предложенная методика может быть использована для анализа других моделей страхования при условии, что нагрузка страховой премии считается малой.
ЛИТЕРАТУРА
1. Лившиц К.И., Бублик Я.С. Вероятность разорения страховой компании при дважды стохастическом потоке страховых выплат // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 1 (10). С. 66 –77.
2. Глухова Е.В., Змеев О.А., Лившиц К.И. Математические модели страхования. Томск:
Изд-во ТГУ, 2004. 180 с.
3. Наумов В.А. Марковские модели потоков требований // Системы массового обслуживания и информатика. М.: УДН, 1978. С. 67−73.
4. Panjer H.Y., Willmont G.E. Insurance Risk Models. Society of Actuaries, 1992. 442 p.
Лившиц Климентий Исаакович
Томский государственный университет
Бублик Яна Сергеевна
Филиал Кемеровского государственного университета в г. Анжеро-Судженске
E-mail: kim47@mail.ru yana@asf.ru
Поступила в редакцию 7 июня 2010 г.
1/--страниц
Пожаловаться на содержимое документа