МОДЕЛИРОВАНИЕ В ЭКОНОМИКЕ
УДК 368.1:519
С. Д. Голубев, Л. А. Черная
ПОСТРОЕНИЕ ФУНКЦИИ ПОЛЕЗНОСТИ В ЗАДАЧАХ ИМУЩЕСТВЕННОГО СТРАХОВАНИЯ
Рассматривается метод построения функции полезности индивидуума, принимающего решение застраховать свое имущество на условиях рынка имущественного страхования. Выведена система функциональных уравнений, определяющих функцию полезности, и предложен метод их решения.
Во многих задачах принятия решений в условиях риска, в частности в задачах имущественного страхования, возникает необходимость количественного описания вкусов и предпочтений индивидуума, принимающего решение (ИПР). Такое описание успешно достигается в рамках так называемой линейной теории полезности. В силу условий задач анализа рисковых ситуаций ИПР вынужден сравнивать не отдельные значения доходов (ущербов или других значимых для ИПР показателей управляемой им экономической системы), а распределения этих показателей. В линейной теории полезности доказывается [1; 2, гл. 7, п. 7.3], что сравнение функций распределения Г(х), заданных на множестве Х0 возможных значений анализируемого показателя, с точки зрения системы предпочтений ИПР эквивалентно сравнению ожидаемой полезности, т. е. распределение (х) считается более предпочтительным, чем ^2(х), если
M {U(x) |Fi(x)} = J U(x)dFi(x) >M {U(x) |F(x)} = J U(x)dF2(x),
Xo Xo
(1)
где и(х) — функция полезности, М{-} — математическое ожидание. Представление ожидаемой полезности в форме функционала вида
т. е. в форме линейного функционала, обусловило появление термина "линейная" применительно к теории полезности.
(2)
Благодаря простоте сравнения распределений в терминах ожидаемой полезности методы теории полезности стали использовать в анализе различного рода рисков, в том числе и страховых, однако трудности построения функций полезности для различных ИПР в определенной мере препятствуют их широкому распространению. Тем не менее, уже в настоящее время в ряде задач страхования и перестрахования весьма затруднительно найти альтернативу показателю ожидаемой полезности при сравнении различных схем страховой (перестраховочной) защиты. Можно также ожидать, что по мере совершенствования и упрощения методов построения функций полезности методы теории полезности будут становиться более популярными и использование показателей полезности в задачах управления рисками станет таким же привычным, как использование ожидаемых значений дохода, прибыли, ущерба и т. д. Основные результаты теории полезности кратко резюмируются следующим образом:
— в рамках аксиоматической теории полезности [2, гл. 7, п. 7.7] сформулирована система аксиом, следствием которых является существование скалярной функции и(х) такой, что для любой пары распределений Fl(x) и ^2(х) выполнено отношение предпочтения — ^\(х) предпочтительнее F2(x) — тогда и только тогда, когда соблюдается неравенство в формуле (1), при этом функция полезности определяется однозначно с точностью до положительного линейного преобразования, т. е. если и(х) — функция полезности, то любая другая функция полезности имеет вид
V (х) = хи (х) + С, (3)
где х > 0, а С — произвольная константа;
— разработана диалоговая процедура построения функции полезности, в которой операционист формулирует последовательность "вопросов" к эксперту, устанавливающему в ответах детерминированный эквивалент некоторой лотереи [2, гл. 7, п. 7.5; 3, разд. 2, п. 2.1]; по результатам этой процедуры строится функция полезности в решетчатом (табулированном) виде, которая затем аппроксимируется подходящей аналитической функцией;
— исследована совокупность аналитических функций, которые могут использоваться в качестве функций полезности, в частности, при аппроксимации решетчатых функций, построенных в процессе диалога с экспертом [3, разд. 2, п. 2.4].
Типовая ситуация принятия решения страхователем, т. е. лицом, стремящимся сохранить свое имущество от воздействия разрушительных факторов, в терминах функции полезности формализуется следующим образом.
Пусть страхователь на основании прошлого опыта (данных статистического учета) знает, что на некотором интервале времени (обычно годовом) его имуществу может быть нанесен ущерб Y, который является случайной величиной с функцией распределения FY(y), т. е.
P {Y < y} = Fy(y).
Случайная величина ущерба Y может принимать только неотрицательные значения, что влечет за собой тождество FY (y) = 0 для отрицательных значений y.
Относительно функции распределения страхового ущерба FY (y) будем предполагать, что существует конечное значение y = Ymin, для которого выполняется равенство
Fy (y) = 1, (4)
т. е. Ymin — минимальный корень уравнения (4). Принятие данного предположения (которое практически всегда соблюдается) дает возможность ограничиться рассмотрением интегралов с конечными, точнее, от 0 до Ymin, пределами при исчислении ожидаемой полезности.
Предположим, что на рынке страховых услуг страховые компании предлагают полную страховую защиту (т. е. полное возмещение ущерба страхователя в случае возникновения на этом интервале такого ущерба) за страховую плату G0; тогда страхователь, уже знающий функцию полезности U(x), должен сравнить значения ожидаемой полезности в двух вариантах:
1) при отказе от услуг страховой компании;
2) при принятии предложения страховой компании со страховой платой Go.
В первом варианте поведения ожидаемая полезность составит величину
Y ■
1 min
M {U(x - Y)} = у U(x - y)dFY(y), (5)
o
где x — начальный капитал страхователя.
Во втором случае значение полезности составит величину
U(x - Go). (6)
Для выбора более выгодного варианта поведения страхователю следует сравнить величины (5) и (6). При этом, если окажется, что
Y ■
^ min
U(x - Go) < У U(x - y)dFY(y), (7)
o
от предложения страховщика следует отказаться (это означает, что затребованная страховщиком плата G0 слишком высока). Если
Y ■
1 min
U(x - Go) > J U(x - y)dFY(y), (8)
0
то предложение о страховании имущества целесообразно принять. Если же
Y
^ min
U(x - Go) = I U(x - y)dFY(y), (9)
0
то оба варианта эквивалентны друг другу, т. е. необходим дополнительный анализ для окончательного принятия решения.
Объединяя соотношения (8) и (9), определим необходимое условие
страхования имущества в виде нестрогого неравенства
Y
^ min
U(x - Go) > У U(x - y)dFY(y). (10)
0
Итак, если страхователь руководствуется принципом максимизации ожидаемой полезности, то ему следует принять предложение страховщика, если соблюдается неравенство (10). Однако страхователь не всегда занимается построением своей функции полезности U(x), тем не менее, в силу определенных соображений страхователь может (считает для себя целесообразным) принять предложение о страховании своего имущества на описанных выше условиях. В этом случае, по-видимому, можно утверждать, что гипотетическая функция полезности такого страхователя удовлетворяет условию (10). Более того, условие (10) может быть использовано для построения функции полезности.
Для решения нестрогого неравенства (10) запишем его в эквивалентной форме:
Y
^ min
U(x - Go) = J U(x - y)dFY(y) + C^(x), (11)
0
где ф(x) — неотрицательная функция, ф(x) > 0; C — положительная константа.
Естественно, что функция ф(x) и константа C могут быть выбраны многими способами, причем указанная неоднозначность вытекает из множественности решений исходного соотношения (10).
Отметим, что существует один частный случай, когда множественность решений соотношения (10) устраняется, а именно если G0 — предельно допустимое значение страховой платы, при которой страхователь еще соглашается на предложение страховой компании о защите имущества. В этом частном случае нестрогое неравенство (10) обращается в равенство, т. е. имеет место уравнение (9). С учетом данного замечания необходимую "добавку" C0(x) в правой части неравенств (8) и (10) можно интерпретировать как меру "запаса по полезности", которую предоставляет страховая компания клиентам прежде, чем у страхователей возникнут сомнения в целесообразности страхования. Обычно страховщики, работающие в том или ином сегменте страхового рынка, могут достаточно точно оценить предельное (максимально допустимое) значение Go, выше которого начнется отток клиентов. Поэтому с точки зрения страховщика представляет интерес оценка функции полезности страхователя в некоторой ограниченной окрестности максимально допустимой страховой платы.
Полученное уравнение (11) относится к классу неоднородных интегро-разностных уравнений. Его решение, естественно, можно найти в виде суммы:
U (x) = Uo (x) + AU (x), (12)
где Uo (x) — решение соответствующего однородного интегро-разност-ного уравнения вида
Y ■
^ min
Uo(x - Go) = У Uo(x - y)dFY(y), (13)
o
а AU(x) — поправка к базовому решению уравнения (13).
Чтобы составить уравнение относительно AU(x), подставим выражение (12) в уравнение (11). Проведя необходимые сокращения с учетом соотношения (13), получим
Y ■
^ min
AU(x - Go) = У AU(x - y)dFY(y) + Cty(x). (14)
o
Итак, функция полезности (12) определяется как сумма решений однородного интегро-разностного уравнения (13) и неоднородного интегро-разностного уравнения (14).
Однородное уравнение (13) будем интерпретировать как условие вычисления максимальной страховой платы Go при условии, что функция полезности страхователя есть Uo (x). Поскольку значение страховой платы Go в уравнении (13) не зависит от начального капитала x, то
данное обстоятельство позволяет предположить, что U0 (x) описывается экспоненциальной функцией (характеристическим свойством функции полезности экспоненциального вида является отсутствие влияния начального капитала x на предпочтения лица, принимающего решение, в частности на величину предельной страховой платы G0 [3, п. 2.4]). Сказанное позволяет искать решение уравнения (13) в виде
Uo (x) w -esx, (15)
где s — константа, подлежащая определению.
Чтобы функция полезности (15) отвечала традиционным требованиям, предъявляемым к функциям полезности лица, не склонного к риску, нужно, чтобы параметр s имел вещественное отрицательное значение [3, п. 2.4]. Подставляя выражение (15) в уравнение (13), получаем характеристическое уравнение относительно s:
Y ■
1 min
-es(x-Go) = -esx j e-syö,Fy(y), o
откуда после сокращения получаем
Y
L min
e-sGo = J e-sydFY(y). (16)
o
Чтобы доказать существование отрицательного числа s, обращающего в тождество равенство (16), а также определить условия существования такого числа, исследуем поведение левой и правой частей равенства в области s < 0.
Отметим сначала, что точка s = 0 всегда является простым корнем уравнения (16). Действительно, левая часть уравнения (16) равна 1 при s = 0. Вычислим
YY
^ min ^ min
lim У e-sydFY(y)= I dFY(y) = Fy(Ymin) - Fy(0). (17) oo
Поскольку возможный ущерб Y является неотрицательной величиной, то с учетом непрерывности функции FY(y) слева естественно принять Fy(0) = 0. С учетом свойства (4) функций распределения
имеем Fy (Ymin) = 1 и формула (17) принимает вид
Y
^ min
lim У e-sydFY(y) =1. (18)
o
Итак, левая и правая части уравнения (16) в точке 8 = 0 одновременно становятся равными 1.
Далее, исследуем поведение левой и правой частей уравнения (16) при малых по модулю отрицательных значениях 8. Учитывая аналитичность левой и правой частей уравнения (16), вычислим их производные по 8 в точке 8 = 0. Имеем
д
lim— e~sG° = -G0,
s^o ds
lim j e~sydFy(y) =lim^— j , ">//-\ (/,) oo
Ym
Ym
(18a)
Ym
= — ydFy (y) = - my,
где ту — математическое ожидание случайной величины У. Поскольку страховая плата С0 складывается из средних страховых выплат ту и рисковой надбавки, то всегда имеет место соотношение
G0 > my,
(19)
или
—G0 < —my
(20)
Следовательно, во-первых, при малых по модулю отрицательных значениях 8 имеем неравенство (см. рисунок)
e-sGo > Z(s),
Ym
(21)
где г(8) = J е~зуdFу(у), во-о
вторых, производные левой и правой частей уравнения (16) в точке 8 = 0 отличны между собой, а отсюда следует, что корень 8 = 0 — простой.
Далее предположим, что на страховую плату наложено ограничение вида
Go + а < Ymin,
(22)
Решение характеристического урав- где а — некоторая положительная нения (16)
константа.
Покажем, что при выполнении условия (22), где Ymjn — минимальный корень уравнения (4), существует отрицательный корень уравнения (16). Для этого необходимо доказать, что при достаточно больших по модулю отрицательных значениях s имеет место неравенство
Z(s) > e-sGo. (23)
Следующая цепочка соотношений приводит к требуемому неравенству:
Y ■
1 min
Z(s)= I e-sdFY(y) =
0
Go+a 1min 1min
= J e-sydFY(y)+ У e-sydFY (y) > j e-sydFY (y) >
0 Go+a Go+a
Y
L min
> e-s(Go+a) j dFY(y) = e-s(Go+a)[1 - Fy(Go + a)] > e-sGo. (24)
Go+a
Так как в силу неравенства (22) и определения Ymin выражение в квадратной скобке формулы (24) есть неотрицательная величина, т. е.
1 - Fy (Go + а) > 0, (25)
следовательно, при достаточно больших по модулю отрицательных значениях s имеем
[1 - Fy (Go + a)]e-sa > 1. (26)
Отметим, что ограничительные условия (19) и (22) на выбор G0 естественно интерпретировать как условия согласованности назначения страховой платы с вероятностной характеристикой страхового ущерба. При нарушении указанных условий, например при попытке назначить страховую плату G0 больше, чем Ymln, сформулированная задача определения U0 (x) в форме (15) не будет иметь решения, в частности кривая Z(s) в отрицательной области "не настигнет" e-sGo.
Поскольку при малых по модулю отрицательных значениях s имеет место неравенство (21), а при больших по модулю отрицательных значениях s знак неравенства (23) меняется на противоположный, то в силу непрерывности левой и правой частей характеристического
уравнения (16) существует такое отрицательное значение s = s, в котором
Y ■
1 min
e-sGo = J e-sydFY(y).
0
Отметим, что получаемое в форме (15) решение однородного интегрального уравнения (13) справедливо на всей вещественной оси, т. е. на интервале (-то, +то).
Переходя к решению неоднородного интегро-разностного уравнения (14), применим преобразование Лапласа. Выбор преобразования Лапласа как метода решения неоднородного интегро-разностного уравнения (14) влечет за собой ограничение на область определения "поправки" AU(x). А именно, при использовании одностороннего преобразования Лапласа "поправка" AU(x) определяется только в области положительных значений x. С практической точки зрения этот недостаток оказывается не слишком ограничительным, так как величина поправки, как показывают расчеты для реальных функций распределения FY(y), весьма невелика по сравнению с U0 (x), тем более, она будет несущественной в области x < 0, где U0 (x) быстро возрастает по модулю при x ^ -то. Воспользовавшись тем свойством преобразования Лапласа, что свертке оригиналов соответствует произведение изображений, получаем
AU (p)e-pGo = AU (v)vFy (р) + Ctp(p), (27)
где
AU(p) = J e-pxAU(x)dx,
0
С
Fy(p) = j e-pxFy(x)dx, (28)
0
С
ф(р) = e-px^(x)dx.
0
Разрешим уравнение (27) относительно ДЦ/(р):
АЩр) = ГСФ(Р1 , (29)
- рРу (р) У J
Прежде чем приступить к обращению изображения (29), установим следующий факт: полюсы изображения (29), обусловленные ну-
лями знаменателя, совпадают с корнями характеристического уравнения (16). Чтобы доказать это, преобразуем правую часть уравнения (16), заменив в ней s на р. Имеем
сю сю
I e-pdFy (y) = e-pFY(y)|0 + pj e-pFY(y)dy. (30) 0 0
При Re p > 0 первое слагаемое в формуле (30) обращается в нуль, второе слагаемое в силу формул преобразования Лапласа (28) равно pFY (р). Таким образом, уравнение
e-pG° — pFy(р) = 0 (31)
с точностью до обозначения переменной совпадает с характеристическим уравнением (16). Выше было установлено, что корнями уравнения (31) являются, во-первых, р = 0, во-вторых, некоторая вещественная отрицательная величина р = s.
Сложности обращения (29) состоят, во-первых, в том, что данное выражение не является дробно-рациональной функцией, а относится к классу так называемых мероморфных функций, т. е. таких аналитических функций, которые в качестве особенностей имеют только полюсы, причем их число конечно в любой ограниченной области [6], во-вторых, операционист обычно располагает только графическим изображением функции распределения FY(y), но не имеет аналитического выражения его изображения по Лапласу FY (р). Поэтому при оперировании выражением FY(р) предполагается, что оно определено интегральным преобразованием (28), реализуемым средствами компьютерного моделирования. Чтобы воспользоваться теоремой разложения оригинала при обращении такого рода изображений, следует убедиться в его правильности в правой полуплоскости комплексного переменного р (Re р > 0): данное изображение убывает при |р| не медленнее, чем -. В рассматриваемом случае правильность изображения р
(29) может быть достигнута только за счет соответствующего выбора его оригинала ф(р), поскольку сомножитель (e-pG° — pFY(р))-1 при р стремится к бесконечности. Действительно, lim e-pG° = 0, и
по теореме о предельных значениях оригинала
lim pFY (р) = lim Fy(x) = Fy (0) = 0,
откуда и следует высказанное утверждение.
Из предельного соотношения
lim pFY(p) = lim Fy(x) = 0
следует, что FY (р) при больших значениях р ведет себя как функ-1
ция —, где а > 1.
ра
В анализе поведения FY (р), а следовательно, всего выражения (29) можно продвинуться дальше, если принять дополнительное предположение, что производная оригинала FY(x) существует в точке x = +0, причем, естественно, эта производная положительна. Применяя предельное соотношение к производной функции FY(x), можно записать
lim р2Fy(р) = lim FY(x) = A, A> 0. (32)
Из формулы (32) следует, что FY (р) при больших р ведет себя как
функция —:
р2
A
МР) « - (33)
р2
Теперь с учетом формулы (33) изображение (29) при больших р ведет себя как
Д£7(р) «---т ■ ф(р). (34)
e-pG0 - р- — р2
Покажем, как можно добиться правильности выражения (34) при больших р в частном случае, когда в качестве ф(x) выбрана плотность гамма-распределения вида
0, если х < 0,
ф(х) = р-( (х; а, в ра (35)
——— х°~1е~13х, если х > 0, г(а)
где а и в — параметры распределения, а > 0, в >0.
Изображение по Лапласу оригинала (35) определяется с привлечением формул (28):
сю сю
ф(р) = [ е-рхф(х)д,х = —Ц- [ ха-1е~(р+Р)х(1х. У Г(а) У
о о
Сделав замену переменных у = (р + [3)х, ¿х = —, получим
Р + в
сю
ф(р) = 1 1 [ уа-1е~Чу = = . (36)
^' Г(а) (р + ¡3)а ] у Т(а)(р + р)а (р + р)а V 7
о
Если в формулах (35), (36) параметр а выбран из условия а > 2, то формулу (29) можно записать в следующем виде:
1 С
Аи= [е-^о _ Рру(р)] (р + /З)2 (р + ¡3)а~2' (37)
В полученном выражении первый сомножитель является правильной
мероморфной функцией, т. е. при больших р ведет себя как -. В этом
р
можно убедиться, если воспользоваться формулой (34):
1 1 1
(е-^о _ р . ^ (р + 0)2 (р + _ АР(Р + ^
Р (38)
Второй сомножитель в формуле (37) в качестве оригинала имеет масштабированную плотность гамма-распределения (37) с параметрами а — 2 и в. Отметим здесь, что "отщепление" от ф(р) сомножителя
--— с присоединением его к изображению---- . . обес-
(р + в)2 * * в-Р°0 — р]ру (р)
печивает его правильность, что позволяет непосредственно перейти к
его обращению. Вместе с тем, "отщепление" сомножителя --—,
(р + в )7
где 7 > 2, нецелесообразно, поскольку усложняет вычисление составляющей оригинала, соответствующего кратному корню р = —в.
Рассмотрим процедуру вычисления оригинала, соответствующего изображению
= [е-^о -рЁу(р)] (р + Р)2' (39)
В силу формулы (38) \¥(р) удовлетворяет требованиям, предъявляемым к изображениям в форме мероморфных функций, что позволяет воспользоваться теоремой разложения оригинала с использованием вычетов данной функции, а именно
w(x) = res W(p)epx. (40)
Pk k
Здесь полюсами являются следующие точки:
1) рг = 0;
2) р2 = §;
3) р3 = —в — двукратный полюс.
Вычет в точке pi = 0 вычисляется по известной [6, гл. 6, п. 82] формуле
res W(p)epx = lim-pi p—^0
pe
px
im
(41)
-pGo
e-pxdFY (x)
где выражение р]у (р) заменено в соответствии с формулой (30) его интегральным эквивалентом. Раскрывая неопределенность в точке р = 0 по правилу Лопиталя, получаем
res W(p)e = lim
pi p— 0
1
Y
i min
-G0e-pGo + / e-XäFr(.)
0
1
epx =
ß 2
ß2 (-G0 + mY)
= Ci. (42)
Аналогично вычисляется вычет в полюсе p2 = s:
res W(p)epx = lim
P2 P—s
p - s
Y
L min
e-pGo - f e-pxdFYV
0
epx = 02 esx,
(p + ß )2
(43)
где
C2 =
1
Ym
-G0e-sGo - J e-sxdFY(x) 0
(S + ß )2
Вычет в полюсе р3 = —в вычисляется с учетом его кратности, равной двум [6, гл. 6, п. 82]:
_ д
res W(p)epx = lim — ■ рз p—-e dp
(p + ß )2epx
im
e-PGo - e-pxdFY(x)
(p + ß )2
= Ö3ixe-ßx + C30e-ex, (44)
ßx
где
Сз1 =
Ym
eeGo — / eß*d,FY (x)
C30 =
Ym
—G0eßGo + J eßxxdFY(x) 0
Ym
е^с - / евх(Еу (х) о
Подставляя выражения (42)-(44) в формулу разложения (40), получаем оригинал У)(х) в виде
w(x) = Ci + C2 esx + Cbixe-ßx + C30e-ßx.
-ßx
(45)
Оригинал изображения
ß
а-2
в силу формул (35) и (36) есть,
(р + ß )а-2
очевидно, также плотность гамма-распределения с параметрами а — 2 и ß, т. е.
L
1
в
•а-2
(р + ß)
а2
= pY(x; а — 2, ß) = ф^).
(46)
Поскольку удалось найти оригиналы каждого из сомножителей выражения (37), то можно вернуться к вычислению оригинала изображения Аи(р), воспользовавшись тем известным фактом, что произведению изображений соответствует свертка оригиналов:
AU(x) = C J w(x — y^^y)dy, 0
(47)
где У)(х) определяется формулой (45), а ф\(х) — формулой (46).
Подстановка формул (45), (46) и (35) в соотношение (47) дает следующий результат:
AU (x) = C
Ci фх(у^у + C2 es(x-y^i (y)dy+
+ C31 / (x — y)e-ß(x-y^i(y)dy + C30 / e-ß(x-y)фl(y)dy
. (48)
2
x
x
x
x
x
Вычислим входящие в формулу (48) интегралы:
X X
! фг(у)д,у = I р., (у; а — 2, в № = (х; а — 2, в), (49) о о
где (•) - интегральный закон гамма-распределения с параметрами
а — 2 и в,
x
s x
J e^My)dy = г(а _ 2) J ya-3e-^dy =
00
x
= e*Jp,(y;o - 2,ß + s)dy = e*fy(X; а - 2, ß + s), (5°)
0
где Fy(x; а - 2, ß + S) — функция гамма-распределения с параметрами
а - 2 и ß + S,
(x - y)e-e(x-yVi(y)dy =
0
x x
p ßx f p ßx f
e 1 ' ^ -*-3„-ßyj„._ e „,\„.a-3,
(x - y)ya-3e-ßydy = --- / (x - y)ya~3dy =
Г(а - 2) J K * Г(а - 2)
00
x О
-вт f / \ /X a~3
e' y\ 3/y\ „j/y
0
Г(а - 2) J \ x/ \x/ \x
x 1-- xa_d - xd
Сделав замену переменных — = получим
x
1
^ (ж _ = ц^ц I ха~\1 - =
оо
1
„-вхха-1 р
= Г1^2) (51)
о
Последнее выражение можно преобразовать, воспользовавшись интегралом Эйлера первого рода [7, гл. 6, п. 140]:
B(p, q)=/ tq-i (1 - t)p-idt, (52)
x
x
i
где В(р,д) — бета-функция, выражаемая через гамма-функцию формулой [7, гл. 6, п. 140]
1
С учетом ф°рмУл (52) и (») „нтеграл / (1 - ОГ-"« пре°бразуется
о
к виду
/ ,1 - ОС-Ч = В,2, „ - 2) = ю = , (М)
J Г(а) Г(а)
о
так как Г(2) = 1! = 1 в силу свойств гамма-функции. С учетом приведенных формул имеем
(x - y)e—ß(x—y)^i(y)dy =
_ e~ßxxa~l Г(а - 2) _ xa~le~ßx ~ Г(а - 2) Г(а) ~ Г(а)
e-ß(x-y)^ (y)dy =
= pY(x; а, ß), (55)
x x
e-ßx с e ßx с
1 eßvya~3e~ßvdy = -- / ya~3dy =
Г(а - 2) J * Т(а - 2)
0 0
-ßx 1 xa—2e—ßx
1 xa-2 = x e
Г(а - 2) (а - 2) Г(а - 1)
= pY(x; а - 1, ß). (56)
Подставляя выражения сверток в формулу (48), получаем окончательное выражение "поправки":
Ли(х) = С [01(х; а - 2, в) + С2(х; а - 2, в + §) +
+ Сэ1р7(х; а, в) + Сзо'Рч(х; а - 1, в)] • (57)
Полученные выражения позволяют сформулировать ограничения на выбор параметров плотности гамма-распределения (35), а именно
1
x
e
в дополнение к указанному ранее ограничению а > 2 параметр в следует выбирать из условия
в + § > 0, (58)
так как в противном случае, если окажется, что в + § < 0, выражения (50) и (56) потеряют смысл (в функциях плотности и гамма-распределения параметры могут принимать только положительные значения).
Важная особенность полученного результата (57) состоит в том, что входящие в эту формулу функции плотности и гамма-распределения вычисляются с помощью встроенных функций в программно-математической среде пакетов типа Mathcad, что существенно облегчает компьютерную реализацию описанных алгоритмов.
Другой вывод, вытекающий из полученных результатов, состоит в том, что эвристически принятое описание "запаса по полезности" в форме плотности гамма-распределения (35) оказалось весьма удобным с точки зрения возможности получения расчетных формул в форме стандартных функций (в замкнутой аналитической форме вычисляются свертки (49), (50), (55) и (56)), причем на возможность получения такого рода формул не влияет конкретный характер исходных данных (принятые ограничительные условия на исходные данные рассмотренной задачи не являются обременительными с точки зрения использования предлагаемых вычислительных алгоритмов, легко проверяются и обычно выполняются в практических задачах).
СПИСОК ЛИТЕРАТУРЫ
1. Фон Нейман Дж., Моргенштерн О. Теория игр и экономическое поведение / Под ред. Н.Н. Воробьева. - М.: Наука, 1970. - 707 с.
2. Д е Грот М. Оптимальные статистические решения / Под ред. Ю.В. Линника и А.М. Кагана. - М.: Мир, 1974. - 491 с.
3. Голубин Ю. А. Математические модели в теории страхования: построение и оптимизация. - М.: Анкил, 2001. - 160 с.
4. Н а т а н с о н И. П. Теория функций вещественной переменной. - М.: ГИТТЛ, 1957.- 552 с.
5. Б а у э р с Н., Г е р б е р Х., Джонс Д., Н е с б и т т С., Х и км а н Д ж. Актуарная математика / Под ред. В.К. Малиновского. - М.: Янус-К, 2001.-409 с.
6. Лаврентьев М. А., Ш а б а т Б. В. Методы теории функций комплексного переменного. - М.: Наука, 1973. - 736 с.
7. С м и р н о в В. И. Курс высшей математики для техников и физиков. Т. 3. М.-Л.: ГТТИ, 1933. - 754 с.
Статья поступила в редакцию 07.02.2005