Математическая модель процедуры пожизненной ренты Текст научной статьи по специальности «Математика»

УДК 519.872

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕДУРЫ ПОЖИЗНЕННОЙ РЕНТЫ

А.А. Назаров, И.А. Ананина

Томский государственный университет E-mail: ananinaia@sibmail.com

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

Ключевые слова:

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

Queuing system with unlimited number of lines and phases, streams of references, a method of limiting decomposition, a life annuity.

Введение

Одной из основных задач экономико-математического моделирования является исследование финансовых потоков, обладающих существенной нерегулярностью по величине и во времени. Наиболее адекватными математическими моделями этих потоков являются случайные потоки однородных событий [1], интенсивно исследуемые в теории массового обслуживания [2] - науке, возникшей из потребностей анализа сетей связи [3]. Такие потоки требуют существенного обобщения и развития методов их исследования применительно к экономико-математическим моделям.

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

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

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

В данной работе предложена математическая модель процедуры формирования указанных потоков в виде системы массового обслуживания (СМО) с неограниченным числом линий и бесконечным числом фаз обслуживания в каждой линии [4].

Математическая модель процедуры

пожизненной ренты

В качестве математической модели процедуры пожизненной ренты рассмотрим систему массового обслуживания (рис. 1) с неограниченным числом линий и бесконечным числом фаз обслуживания в каждой линии.

На вход системы поступает пуассоновский поток заявок (владельцев недвижимости, желающих заключить договоры пожизненной ренты) с интенсивностью Я. Заявка входящего потока, поступившая в систему, занимает одну из свободных линий (клиент заключает договор пожизненной ренты). Линия считается свободной, если свободны все ее фазы обслуживания. Обслуживание заявки выполняется последовательно на каждой фазе - временных периодах, определяемых договором для регулярной выплаты ренты (ежемесячно). Будем полагать, что продолжительность обслуживания на каждой фазе случайная (в частности детерминированная) с функцией распределения В (х).

Завершив обслуживание на к-й фазе, заявка с вероятностью 1-гк покидает систему, а с вероятностью гк переходит для обслуживания на следующую (к+1)-ю фазу. Вероятности гк определяются функцией дожития [5] в виде:

гк=8(к)/8 (к-1), (1)

где £(к) - функция дожития, то есть вероятность, что человек доживет до возраста, определяемого номером к рассматриваемой фазы обслуживания, а 5(0)=1.

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

Найдем совместную производящую функцию трехмерного случайного процесса {у('), п('), ш(')}

нху^^Щх^'У^}.

Метод предельной декомпозиции

Пуассоновский входящий поток рассматриваемой СМО разделим по полиномиальной схеме с равными вероятностями на N независимых пуас-соновских потоков с интенсивностями Я/N. Для заявок каждого из потоков определим единствен-

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

Реализуем сказанное для нахождения производящей функции Н(х,у ,£,/).

Исследование однолинейной системы с отказами

Рассмотрим однолинейную систему обслуживания с отказами (рис. 2).

Для этой системы рассмотрим аналогично трехмерный случайный процесс {у(/,Л), п(/,Л), т(/,Л)}. Для марковизации этого процесса определим еще две компоненты: к(0 - номер фазы, занятой в линии в момент времени ¡. Если к(0=0, то линия свободна; и(/) - длина интервала от момента / до момента окончания обслуживания на текущей фазе.

Очевидно, что пятимерный случайный процесс {к(0,и(0,у(/,Л),п(/,Л),т(/,Л}} является марковским, поэтому для его распределения вероятностей Р0(у,п,т,1,Л)=Р{к(1)=0,у(Щ=у, п(1,Л)=п,т(1,Л)=т}, Рк(у,п,т,и,1,Л)=

=Р{к(Г)=к,у(1,Л)=у,п(1,Л)=п,т(1,Л)=т,и(Г)<и}

получим систему дифференциальных уравнений Колмогорова

дР0(у, п, т, г, И) Х„. лгч

0 ' =-—Ро (у, п, т, г, И) +

дг N

^дР^п^т—лгИ! а - )

+X - (1 Гк),

-_1 ди

дР1(у, п, т,и, г, И) _ дР1(у, п, т, и, г, И)

дг

ди

дР(у,п,т,0,г,И) ЛГ. X .

--^ , , , , , > + Ро (у -1, п, т, г, И)—В (и),

ди И

дРк(у, п, т,и, г, И) _ дРк(у, п, т, и, г, И)

дг ди

дРк (у, п, т, 0, г, И) +

ди

ди

Обозначив

XXX хуУ"2тРо(у,п,т,г,И) _ Но(х, у,г,г,И),

у_0 п_о т_0

XXX хуУ"2тР-(у,п, т,и,г, И) _ Н- (х, у, ^и, г,И)

у_0 п_0 т_0

из системы (2) получим систему уравнений для функций Н0(х,у,1^,Л) и Нк(х,у,1,и,1,Л)

дНо (х, у, г, г, И ) X

К д(-- = ~^Но (х, у, г, г, И )+

+гX'дНк (x, У, г, 0 г, И) (1 - ),

- _1 ди - '

дН1 (х, у, г, и, г, И) _ дН1 (х, у, г, и, г, И)

дг

ди

дН,(х, у, г,0, г, И) ТТ , „„ X .

--* ' " > + хН о(х, у, г, г, И)— В (и), (3)

ди И

dHk (х, y, z, u, t, N) = öHk (x, y, z, u, t, N) _ Fk (x, y, z, u, t) =

dt du t

дИк (x, y, z,0, t, N) dHk _,(x, y, z,0, t, N) = \{_k(x, У, z, s) + Угк _A _i(x, У, z, s) B(u +t _ s)}ds,

я ^---ä-Гк _1B(u )

du du

0

откуда получим

Так как в начальный момент времени ¿=0 ис- др ( 0 )

ходная СМО с неограниченным числом линий ^ (х, у, г, :) = —'(хуг ' ' -.

свободна, то свободна каждая линия, поэтому на- ди

чальные условия для решения Н0(х,у,1^,Щ, >

Нк(х,у,1,ы,1,Ы) системы (3) определим в виде = ^х\в (: ~ ^^ =^хВ(:),

Но(х,у,г,0,Ш)=\, Нк (хуАиДЛИ. (4) ° д

Найдя решение задачи Коши (3,4), определим ^ (х у г :) = к(x, у, г ,0 ,:)

маргинальную производящую функцию трехмер- к ди ного процесса {у(/Д), пЦ,Щ, ш(1,Щ

Н (х , у, г ,:, И) =М {хУ('И) уи ) г™)} =

I

: yrk_i \h_i(x, У, z, t _ z)dB(z).

= H0 (x, y, z, t, N) + lim У Hk (x, y, z, u, t, N). (5) Из этой РекУРРентной пр°цедуры полУчаем, что

к=1 функции hk(x,y,z,t) имеют вид

Для различных однолинейных СМО процессы hl(x,y,z,i)=hcB(i),

{v(t,N),n(t,N),m(t,N)} стохастически независимы hk(x,y,z.,t)=Xxyk-1S(k-1)B(k)(t),

и одинаково распределены. Соответствующие сум- где B(k)(t) - k-кратная свертка распределений B(t),

марные процессы для совокупности N однолиней- k

ных СМО являются суммаминезависимых и одина- а S(k) в силу равенства (1) имеет вид S(к) = Т\г..

ково распределенных случайных процессов для од- =1

нолинейных СМ°. Поэтому производящую функ- Подставляя найденные для h(x,y,z,t) выражения

цию H(x,y,z,t) для процессов исходной системы в (9), получим

с неограниченным числом линий запишем в виде р (x y z t) =

H(x , y, z,t) = lim HN (x ,y,z,t,N) (6) и 0 ' ' ' t

tj nT » = _Xt + zXxY (1 _ r.)yk_1S(k _ 1) fB(r)dz,

Решение системы (3) найдем в виде ^ kJJ v Ч у

k =10

1 _2

H0(x,y,z,t,N) =1 + N p0(x,y,z,t) + oN), P1(x, y, z, u, t) = Ax J{B(u +1 _z) _ B(z)}dz,

Hk(x,y,z,u,t,N) = N р1( x,y,z,u,t) + N-2). (7) pk (x, y, z, u, t) = Xxyk_1S (k _ 1) X

Тогда для функций Fli(x,y,z,t) и Fk(x,y,z,u,t) зада- ' k1) (k}

чу (3,4) запишем следующим образом: XJ{B (z) B(u +t _z) _B (z)}dz. (10)

dF0( x, y, z, t) 0

-d-= _X + z у hk (x, y, z, t )(1 _ rk), Полученные равенства (10) решают поставлен-

dt k=1 ную задачу нахождения функций H0(x,y,z,t,N) и

dP1 (x, y, z,u, t) dP1(x, y, z, u, t) = Hk(x,y,z,u,t,N) из (8).

dt du

= _h1 (x, y, z, t) + xXB(u),

Исследование системы с неограниченным числом линий

dPk(x, y, z,u, t) dPk(x, y, z, u, t) = Применяя равенства (7) запишем функцию

dt du H(x,y,z,t,N) из (5) в виде

= _h (x, y, z, t) + yhk _1 (x, y, z, t)rk _B(u), (8) H (x, y, z, t, N) =

Fo(x,y,z,0)=0, Fk(x,y,z,u,0)=0. = H0(x, y, z, t, N) + lim У Hk (x, y, z, u, t, N) =

u k=1

, , ч dPk (x, y, z,0, t)

где hk(x,y,z,t) =—kV ' ' '. 1 2

du = 1 +—P(x, y, z, t) + o(N ). (11)

Нетрудно показать, что решение этой системы N

имеет вид Здесь функция F(x,y,z,t) имеет вид

P0 (x, y, z, t)=_ Xt + z У (1 _ rk) \hk (x, y, z, s)ds, P (x, у, z,t) = р0(x, у, z,t) +1™ У Pk(x, у, z,u,t).

k=1 0 k=1

t Подставляя сюда выражения (10) для функций

P1 (x, У, z, u, t) =Xx J{B(u +1 _ s)_B(s)}ds, (9) Fk(x,y,z,t), получим

^ (х, у, г, г) = (х -1) Хг + --) ад-1)+| у,л

Й1+(у-г)5(*) Г 0

Применяя равенства (6) и (11), найдем производящую функцию Н(х,у,£,0 трехмерного случайного процесса {|(0,«(0,т(0}, определяющего основные финансовые потоки процедуры пожизненной ренты.

Н (х, у, г, г) =

= 11т {1 +N ^(х, у, г, г) + о(N~2)| = е^

Таким образом, согласно (12), функция Н(х,у,г,/) имеет вид

Н (х, у, г, г) =

(х - 1)Хг +

= ехр

♦Ях^"-1) 5 (к-1) + [ у' -ГВ«' >(г)Л 1+(у-г)..(к) К 1

Системы с детерминированной продолжительностью фаз обслуживания

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

[О, если г < 1,

[1, если г > 1. Тогда для к-кратной свертки запишем

Ве)(,) = |0, если г < ', {1, если г > к.

Поэтому для любых значений / выполняется равенство

[0, если г < к,

|В(к\т)с1т =

г - к, если г > к.

При />1 функцию Н(х,у,г,1) из (9) можно записать в виде

Н (х, у, г, г)= (х - 1)Хг +

=ехр

+Хх£ -1)5 (к -1) + [+(у - г)5(к)

ук-1(г - к)

(13)

но поток, определяемый процессом «(/), - не пуас-соновский и определяется производящей функцией Н(1,у,1,0 числа фаз, завершенных за время

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

Процесс изменения прибыли

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

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

Обозначим Щ() - прибыль, накопленную плательщиком ренты за время / коммерческой деятельности по предоставлению услуг пожизненной ренты.

Так как прибыль равна разности доходов и расходов, то

следовательно, характеристическая функция величины прибыли ЩО имеет вид

в(а, г) = Ые-аШ (г) =

Г т({) п (г) у(1) 1

= Мехр+ а£д- + [.

{ ¡=1 ¡=1 ¡=1 ]

Полагая стохастическую независимость величин 4, 7, 4, принимающих свои значения в различные моменты времени наступления платежей для различных получателей ренты, а также их одинаковые по I распределения, для О (а,/) получим

в (а, г) = М {еа41(г )еащп( {)е -а4'т(г)} =

= Х М

¡еа41"( г) еа77'

>('),

хе

-а^т«)

|у(!) = V, п(г) = п, т(г) = т

где [/] - целая часть от числа

Найденное выражение (13) для производящей функции Н(х,у,£,0 трехмерного процесса {|(0,и(0,т(0} решает задачу анализа основных потоков математической модели процедуры пожизненной ренты.

Из вида (13) следует, что процессы и т(/) -пуассоновские, что вполне естественно для входящего и выходящего потоков системы M|GI|<x>,

хР(|, п, т, г) = = £ (МеаС1У (Меа771 )п(Ме^урп,т,г) =

т ,п V

= Н (х(-а), у(-а), г(а), г),

где х(а)=Меа1; у(а)=Меа; ,г(а)=Меа - характеристические функции случайных величин 4, 7 и 4 Таким образом, характеристическая функция величины прибыли плательщика ренты, накопленной за время имеет вид

О(а,0 = Н(х(-а),у(-а),г(а)Д (14) где функция Н(х,у,1,1) определяется равенством (13).

Исследование среднего значения величины прибыли

Из свойств характеристических функций следует, что

дв (а, г)

Ж1(г )=МЖ (г) =--

да

поэтому, в силу равенства (14), получим

№1(1)=М£Мт (0-М^МП(0-М£-МУ(0.

Пусть М£=1 - стоимость квартиры, например 2 млн р; Мг1=5 - доля от стоимости квартиры, которую следует выплатить клиенту за один месяц в качестве ренты; МС=51 - доля от стоимости квартиры, которую может потратить фирма на привлечение клиента. Следовательно, №() имеет вид W1(t)=Mm(t)- ЙЩО-^МЦО.

Средние значения для процессов у(0, п(0 и т(0 определим из вида (13) производящей функции Н(х,у,1,1), получим

Ж1 (г)=я£ [5* (к -1) - 5(к)](г - к) -

к=1 [ г ]

5(к)(г - к) -51Яг.

(15)

Из (15) очевидно следует, что среднее значение №() величины прибыли линейно зависит от интенсивности Я заключения новых договоров пожизненной ренты и существенно нелинейно №() зависит от t.

Определение функции дожития

Функцию $(к) в демографии называют функцией дожития, ее аналитическим выражением, наиболее адекватно соответствующим статистическим данным, является функция Гомперца-Мейкема

ЗД=ехр{-Ак-В (е-М)}, в которой параметры А, В и а определяются из статистических данных.

Для мужчин Российской Федерации известны следующие статистические данные [6] для безусловной функции дожития (табл. 1).

Таблица 1. Безусловная функции дожития для мужчин возраста х, лет

х 60 70 80 90

В(х) 0,651 0,434 0,188 0,003

к 0 120 240 360

S6o(k) 1 0,6667 0,2888 0,0046

ние условной функции дожития до к-го месяца мужчин старше 60 лет.

Значения параметров А, В и а для условной функции дожития

£60(АДа,к)=ехр{-Ак-В(еак-1)} найдем как решение системы трех нелинейных уравнений

£60(А,В,а,120)=0,6667, S60(A,B,а,240)=0,2888, £60(АД а,360)=0,0046 относительно А, В и а.

Это решение принимает следующие значения А=0,0028, В=0,0097, а=0,017, следовательно, для численных расчетов функцию £60(к) можно брать в виде

£60(к)=ехр{-0,0028к-0,0097(е°'017к-1)}.

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

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

Численный анализ математической модели

пожизненной ренты.

Экономические рекомендации

Подставляя найденную функцию дожития £60(к) в (15) и полагая 5=0 - затраты на привлечение новых получателей ренты учтены в остаточной стоимости недвижимости, определяющей величину ренты, выполним численный анализ среднего значения №¡(0 величины прибыли плательщика ренты при Я=1 (один новый клиент в месяц) и среднем значении остаточной стоимости квартиры в 2 млн р.

Изменение во времени величины прибыли ЩО для ряда значений 5 показано на рис. 3.

Из результатов численного анализа следует, что при 5<0,003 прибыль плательщика ренты положительна при всех t, при 5>0,0058 величина №[(0<0, то есть при таких значениях 5 плательщик ренты терпит только убытки от своей коммерческой деятельности. При 0,003<5<0,0058 функция №() имеет вид

Здесь х - возраст мужчины; £(х) - значения функции дожития для мужчин возраста х. Поэтому для условной функции дожития £60(к) для мужчин старше 60 лет эти значения будут следующими (табл. 2).

Таблица 2. Функция дожития до к-го месяца для мужчин старше 60 лет

Щ) =

< 0, при г < Т(5), > 0, при г > Т(5).

Здесь к определяется продолжительностью выплат рентных платежей в месяцах, а £60(к) - значе-

Продолжительность Т(5) первоначального периода издержек достаточно велика. Например, при 5=0,004 величина Т(5)=25 лет, а при 5=0,005 Т(5)=45 лет. При 5=0,004 максимальные убытки составляют 9,264 млн р. и достигаются на 18 году деятельности, а при 5=0,005 убытки достигают 31,728 млн р. и приходятся на 23 год коммерческой деятельности.

Из приведенных расчетов, очевидно, следует, что применение 5>0,003 для плательщика ренты экономически нецелесообразно.

а=0

к=1

При <5=0,003 и стоимости квартиры в 2 млн р. величина ренты составляет 6 тыс. р. в месяц, но в среднем эти выплаты покрывают лишь половину стоимости квартиры, так как средняя продолжительность дожития мужчины старше 60 лет составляет 172 месяца. Такая ситуация вряд ли может заинтересовать владельца недвижимости к заключению договора пожизненной ренты.

Выводы

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

СПИСОК ЛИТЕРАТУРЫ

1. Гнеденко Б.В., Коваленко И.Н. Введение в теорию массового обслуживания. - М.: КомКнига, 2005. - 400 с.

2. Хинчин А.Я. Работы по математической теории массового обслуживания. - М.: Физматлит, 1963. - 527 с.

3. Лившиц В.С., Пшеничников А.П., Харкевич А.Д. Теория телетрафика. - М.: Связь, 1979. - 224 с.

4. Галажинская О.Н. Бесконечно линейная бесконечно фазная система массового обслуживания со случайным прерыванием

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

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

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

обслуживания // Вестник Томского государственного университета. - 2006. - № 18. - С. 261-266.

5. Валентий Д.И., Кваша А.Я. Основы демографии. - М.: Мысль, 1989. - 285 с.

6. Четыркин Е.М. Пенсионные фонды. Зарубежный опыт для отечественных предприятий, актуарные расчеты. - М.: АО «Арго», 1993. - 100 с.

Поступила 17.01.2011 г.