Предлагается модель изменения страхового капитала Пенсионного фонда Российской Федерации в зависимости от суммы поступивших страховых взносов и суммы выплаченной страховой части пенсии застрахованным лицам. Изучаются основные характеристики капитала фонда.
Investigation of Russian Federation retirement fund insurance capital modification processmathematical model.pdf Проводимая с начала 2002 г. пенсионная реформа вне-сла кардинальные изменения в принципы работы пенсион-ной системы. В соответствии с реформой трудовые пенсиибудут складываться из трех составляющих: базовой, страхо-вой и накопительной. Базовая пенсия будет являться мини-мальным размером трудовой пенсии, неким гарантирован-ным минимумом и будет выплачиваться из федеральногобюджета за счет средств от уплаты единого социальногоналога (ЕСН). Страховая часть пенсий зависит от результа-тов работы конкретного работника и будет равна суммевзносов, которые будут перечисляться работодателем в егопользу. Все взносы, поступающие на индивидуальный счетзастрахованного, будут суммироваться и к достижению импенсионного возраста составят определенный страховой ка-питал. Таким образом, размер страховой части трудовойпенсии определяется исключительно величиной страховогокапитала, образуемого из взносов по обязательному пенси-онному страхованию, перечисленных за застрахованныхлиц. Данная часть пенсионных выплат будет осуществлять-ся из Пенсионного фонда(t) = r 1−e− t−t + e− t−t,1( 2( 0) )( )2 0 22(t) =12be− t−t −1r+ +1 ( )1 0 02ra 12b r t t⎛ ⎞+⎜⎝ − ⎟⎠ − +, (3)где a1 , b1 - первые начальные моменты случайныхвеличин, определяемых функциями распределенияA(x) и B(x) ; (t0)=0, (t0)= 0.Доказательство. Обозначим 1 = и выполним в(2) замену j =y, S =z, 21 P(j,S,t) = (y,z,t,),тогда уравнение (2) примет вид2(y,z,t, ) (r ( 12)y) ( y, z,t, )t + + + =0( , ,,) ()zr y z ut dAu= − − +012 y (y,z u,t, )dB(u)+ + ++2(y+ )(y+ ,z,t,).Раскладывая функции (y ,z u,t,) в ряд поприращениям аргументов с точностью до o() и вы-полнив элементарные преобразования, получим урав-нение{(r 2y)} {(ra112b1y)}t y z = − − − − ,которое является вырожденным уравнением Фоккера- Планка для плотности (y,z,t) распределения ве-роятностей значений двумерного диффузионногопроцесса с коэффициентами переноса (r− 2y) и(ra1−12b1y) и коэффициентами диффузии, равныминулю. Следовательно, рассматриваемый процесс яв-ляется детерминированным.В силу полученных результатов можно утвер-ждать, что процессы (t) и (t) удовлетворяют сис-теме обыкновенных дифференциальных уравнений{ 21 1( ) ( ),( ) 12 ( ),t r tt ra b t = − = − решение которой имеет вид (3). Теорема доказана.Отметим, что в стационарном режиме2 = r, апервые начальные моменты связаны отношением21 12 1b a= .Далее проведем исследование процесса отклоне-ния процессов 1 j(t)и 1 S(t)от найденных средних(t) и (t) . Для этого рассмотрим предельный, при , процесс для последовательности⎧ j(t) − (t) ,S(t) − (t)⎫⎨ ⎬⎩ ⎭. (4)Теорема 2. При предельный процесс{y(t),z(t)} для последовательности (4) является дву-мерным гауссовским диффузионным процессом с ко-эффициентами переносаA1(y,z,t) = −2y; A2(y,z,t) = −12yb1 (5)и диффузииB11(t) =r+ 2(t), B22(t) =ra2 +12(t)b2,B12 =ra1, (6)где a1 , b1; a2 , b2 - первые и вторые начальные мо-менты соответствующих случайных величин.Доказательство. Обозначим 1 = 2и выполнимв (2) заменуj2 = (t) + y, S2 = (t) + z,21 P(j,S,t)=H(y,x,t,).Получим уравнениеH(y,z,t, ) (t) H(y,z,t, ) (t) H(y,z,t, )t y z − − + +(r+ (2 +12)((t) + y)H( y, z,t,) =2( )0( , ,, ) ()t zr Hy z ut dAu += − − +012 ( (t) y) H(y,z u,t, )dB(u)+ + + ++2((t)+ (y+ ))H(y+ ,z,t,).Раскладывая функции H(y ,z u,t,) в ряд поприращениям аргументов с точностью до o(2) и вы-полнив элементарные преобразования, получим при 0 :H { 2yH} { 12b1yH}t y z =− − − − + { } { }2 22 2 2 2 21 ( ) 1 12 ( )2 2r t H ra b t Hy z + + + − + 2ra1 H .y z+ Данное уравнение является уравнением Фоккера -Планка для плотности H(y,z,t) распределения веро-ятностей значений двумерного диффузионного про-цесса {y(t), z(t)} с коэффициентами переноса (5) идиффузии (6).В силу того, что коэффициенты диффузии не зави-сят от пространственных переменных y и z , а коэф-фициенты переноса зависят от них линейно, диффу-зионный процесс {y(t), z(t)} является гауссовским.Покажем это.Двумерный диффузионный процесс {y(t), z(t)}определяется системой двух стохастических диффе-ренциальных уравнений{ 1 11 1 12 22 21 1 22 2( ) ( ) ( ),( ) ( ) ( ),dy t A dt dw t dw tdz t A dt dw t dw t= + += + + (7)где w1(t) и w2(t) - независимые стандартные вине-ровские процессы.Исходя из полученных значений коэффициентовдиффузии (6), определим параметры ij системы (7).Пусть 21 = 0 , тогда1222 22 1222( ) ( ), ( ) ( ) ,( )t B t t B tB t = =21211 1122(t) B B .B = − (8)В этом случае система (7) принимает вид{ 2 11 1 12 21 22 2( ) ( ) ( ) ( ) ( ),( ) 12 ( ) ( ),dy t ydt t dw t t dw tdz t b ydt t dw t= − + + = − + (9)где коэффициенты ij определяются равенствами (8)и не зависят от значений процесса {y(t), z(t)}. Тео-рема доказана.Используя формулу Ито, найдем решение системы(7) при нулевых начальных условиях в виде:2 2 20 00 011 1 12 21 22 2( ) ( ) ( ) ( ) ( ) ,( ) 12 ( ) ( ) ( ) (10)t tt s st tt tt ty t e s e dw s s e dw sz t b y s ds s dw s− ⎧ ⎧ ⎫ ⎪ ⎪ ⎪ = + ⎨ ⎬ ⎪⎪⎪ ⎪⎭ ⎩ ⎨⎪=− + ⎪⎪⎩ Таким образом, мы получили явное выражениепроцесса {y(t), z(t)}, являющегося двумерным гаус-совским процессом.Найдем корреляционную функцию процесса z(t){ } 1 2 1 2 ( , ) ( ), ( ) z R t t =M z t z t . Используя формулы (10) исвойства стандартных винеровских процессов, полу-чим( 2 1 2 0 2 1 2 021 (max( , ) ) (min( , ) )1 22( , ) t t t t t tzR t t = a r e− − −e− − −2 1 2 ) 2 ( )2 12 02e t t 1 ra 12b r min(t ,t ) t − − ⎛ ⎞− + +⎜⎝ − ⎟⎠ −. (11)В силу замены (4) процесс S(t) изменения страхо-вого капитала Пенсионного фонда Российской Феде-рации (процесс изменения) имеет видS(t) = (t) + z(t), (12)где детерминированная функция (t) определяетсяформулой (3), а z(t) - гауссовский случайный про-цесс, определяемый системой (9) стохастическихдифференциальных уравнений, решение которой за-писано в виде (10). Очевидно, что математическоеожидание процесса S(t) составляет (t) , а его кова-риационная функция определяется равенством (11).Так как изменение гауссовского процесса полностьюопределяется его математическим ожиданием и кова-риационной функцией, то для процесса изменения ве-личины страхового капитала S(t) проведено его пол-ное исследование, из которого достаточно просто по-лучаются все его характеристики.Таким образом, в работе показано, что процессизменения страхового капитала Пенсионного фондаможно аппроксимировать гауссовским процессом (12)с математическим ожиданием (t) и ковариацион-ной функцией Rz(t1,t2), где Rz(t1,t2) - корреляци-онная функция процесса z(t) , определяемая равенст-вом (11).
Гарайшина Ирина Рашитовна | Томский государственный университет | аспирантка кафедры теории вероятностей и математической статистики факультета прикладной математики и кибернетики | irina_g@asf.ru |
Назаров Анатолий Андреевич | Томский государственный университет | профессор, доктор технических наук, заведующий кафедрой теории вероятностей и математической статистики факультета прикладной математики и кибернетики | nazarov@fpmk.tsu.ru |
Гихман И.И., Скороход А.В. Стохастические дифференциальные уравнения. Киев: Наукова думка, 1968. 354 с.
Гнеденко Б.В., Коваленко И.Н. Введение в теорию массового обслуживания. М.: Наука, 1987. 336 с.
Федеральный закон от 17.12.01 г. № 173-ФЗ «О трудовых пенсиях в Российской Федерации».
Федеральный закон от 15.12.01 г. № 167-ФЗ «Об обязательном пенсионном страховании в Российской Федерации».