Об одном семействе неаффинных моделей временной структуры доходности | Вестн. Том. гос. ун-та. Управление, вычислительная техника и информатика. 2018. № 42. DOI: 10.17223/19988605/42/6

Об одном семействе неаффинных моделей временной структуры доходности

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

On a family of non-affine models of the term structure of yield.pdf Предположим, что состояние финансового рынка описывается процентной ставкой r(t), которая следует однородному по времени марковскому процессу, порождаемому стохастическим дифференциальным уравнением [1] dr(t) = ц(г(0) dt + a(r(0) dw(t) с функцией дрейфа ц(х), функцией волатильности с(х) и стандартным винеровским процессом w(t). Для удобства рассуждений будем обозначать функцию дрейфа m(r) = ц(г) - A,(r)a(r) и функцию диффузии s(r) = 0,5 a2(r). Здесь X(r) - так называемая рыночная цена риска [2]. Ранее [3] рассматривалась задача определения временной структуры доходности бескупонной облигации [4], когда функции m(r) и s(r) являются полиномами. Выяснялось, могут ли в этом случае кривые доходности быть полиномами или степенными рядами по переменной r. Оказалось, что это имеет место, если только m(r) и s(r) -полиномы не более чем первой степени. В этом случае модели временной структуры доходности являются аффинными. В настоящей статье рассматривается подобная задача, но временная структура цены бескупонной облигации ищется в виде функционального ряда, отличающегося от степенного. Выяснено, что для некоторых случаев такие решения существуют. Получающаяся временная структура оказывается неаффинной и описывается вырожденными гипергеометрическими функциями. В это семейство входят такие известные модели процентных ставок, как модель CIR(1980) [5] и модель Ана-Гао [6]. 1. Общее уравнение для цены облигации и ее компонентов Рассмотрим уравнение временной структуры для цены бескупонной облигации P(r, т) [7]: +m(r)^ + s(r) д2рМ-rP(r,т) = о, P(r, 0) = 1. (1) дт дr дr Здесь m(r) - функция дрейфа краткосрочной процентной ставки, а s(r) - квадрат ее волатильности. Будем искать решение этого уравнения в форме а(т) (2) с n P(r, х) = У r n=0 где а(х), а и cn, n = 0, 1, 2, ..., - функция и коэффициенты, подлежащие определению. Соответствующие производные, использованные в уравнении (1), имеют вид: л а+n dP(r, х) _ a (х) ■ с n -У (а + n)\ Si a(x) to V r / , ч\а+n+1 ' a(x) a(x) 8 P(r, х)_ 1 У (а + n) n 8r a(x) V r n=0 / , N\а+n+2 a(x) 8 2 P(r, х)_ 1 (3) r к У (а + n)(а + n +1) 8r 2 a(x)2 n=o Предположим, что дрейф и волатильность краткосрочной процентной ставки таковы, что функции m(r) и s(r) являются полиномами порядкаp и q соответственно: p q m(r) = Уmkrk, s(r) = Уskrk. (4) k=0 k=0 Прежде чем подставлять явные выражения производных (3) и полиномов (4) в уравнение (1), преобразуем второе и третье слагаемые в левой части уравнения (1) к более удобной форме. \ a.+n-k+1 т ( p = -У У (а + j + k -1) mka(i) а+j+k-1 (5) (6) (7) (8) с j+k -1 \ / j \ а+j m (r) -Pj1 = -УУ(а+n) mka{xT+n 11 8 r „=n t=n V r j=1-p Vk =max{0,1-j} ( q Vj^+J 82P(rx) °° q f1\ +2 s(r) -^-j-1 =УУ(а+n)(а+n+1)sta(x)a+ncn 1-1 = У У (а+j +k - 2)(а+j+k - ШхГ+1-2 с^ 8 r n=0 k=0 Vr) j=2-q V t=max{0,2- j} 8P(r, х) 8х = -V( а(х) j=-1 а++1 , - rP(r, х) = - у ( Cj+1 Запишем в аналогичной форме также первое и четвертое слагаемые левой части уравнения (1): а+j = -УУ((а + j) а'(х) а(х)а+-1 Cj )f1 j=0 а+j )f 1 Теперь, подставляя выражения (5)-(8) в уравнение (1), получим ( p У -1 (j |0)(а + j) a' (х) а(х)а+j-1 с} -1 (j\-1) а(х)а++1 j -1 (j |1 - p) У (а + j + k-1) mka(zT+j+k-1 cJ+k_, + j V k=max{0,1-j} A /j^+j +1 (j \2 - q) £ (а + j + k - 2) (а + j + k -1) sk a(zT+J+k-2 с j+k-2 k=max{0,2-j} ) (9) = 0. r Определенную сложность в выражении (9) вызывает тот факт, что суммирование по индексу j для каждого слагаемого начинается по-разному: для первого слагаемого j > 0, для второго слагаемого j > -1, для третьего слагаемого j > 1 - p, для четвертого слагаемого j > 2 - q. Поэтому в выражениях слагаемых появились множители I(j | к), представляющие индикаторные функции, равные единице, если j > к, и нулю в противном случае. Равенство (9) должно выполняться равномерно по переменной r. При этом, поскольку функции r-j (j = 0, ±1, ±2, ...) являются линейно независимыми, коэффициенты перед этими функциями в выражении (9) должны быть равными нулю. Это приводит к системе уравнений для неизвестных параметров а, а(х) и Cn, n = 0, 1, 2, ..., в представлении (2) решения уравнения (1), если оно существует в таком виде. Заметим, что каждое слагаемое в каждом элементе суммы (9) имеет ненулевой множитель а(х)а, поэтому для упрощения на него можно сократить во всех элементах суммы. 2. Модель с нулевым дрейфом Среди моделей процессов краткосрочной ставки r(t) с нулевым дрейфом широко известна модель CIR(1980) [5], в которой ставка порождается в общем случае диффузионным процессом dr = с r7dw. (10) Для такого процесса плотность вероятностей f (х) выражается в виде [8]: Л-2у X 2 x-1 V ro J f (х) = (у-1)(2у-1) - х > Г0 > 0. r V ro J 0 График этой плотности представлен на рис. 1 для случая у = 1,5, который чаще всего встречается в литературе. Именно в таком виде эту модель предлагали и авторы оригинальной статьи. 0.30 0.25 0.20 0.15 0.10 0.05 ^-,-,-,-,-,-,-,-,-,-,-,-,-,-,-,-,_ xr0 2 4 6 8 10 Рис. 1. Плотность вероятностей краткосрочной ставки для модели CIR(1980) Несмотря на то что модель известна давно, до сих пор не была описана временная структура ее бескупонной доходности. Оказывается, предлагаемый способ нахождения временной структуры позволяет это сделать. Примем в уравнении (10) у = 1,5 и s = 0,5 с2. Уравнение (1) для цены бескупонной облигации P(r, т) приобретает вид дP(r, т) 3 д2P(r, т) (11) - + sr 5т - rP(r, т) = 0, P(r, 0) = 1. д r1 Будем искать решение этого уравнения в форме (2). Соответствующие производные имеют вид: а+п+2 д 2 P(r, т) dr 2 5P(r, т) а ' (т) а(т) а(т) 1 -- Z (а + п)(а+ п +1) Z (а+п)\ c Сп • п а(т) п= дт а(т) п=0 r r п=0 После подстановки этих выражений в уравнение (11) получаем следующее равенство: а+п-1 a(т)Z (а + п)(а + п +1)| п=0 V J п=0V J п=0 Это равенство можно переписывается в более удобном виде: а+п-1 а' (т) а(т) Z а(т) а(т) а(т) - а(т)Z Z (а + п) сп = 0 . Сn + S c r (а + п)а (т) cn + а(т) cn+1 - sa(т)(а + п + 1)(а + п + 2)cn а(т) а-1 а(т)" а(т) + а(т)(1 -sа(а+1)) cc = 0. r r п=0 Поскольку выражения (а(т)/г)к как функции переменной r для различных значений к являются линейно независимыми, а равенство должно выполняться равномерно по r, то коэффициенты перед этими выражениями для различных к должны быть равны нулю. И мы получаем систему уравнений относительно неизвестных а, а(т) и с„, n = 0, 1, 2, ...: sa(a + 1) = 1, (12) а' (т) (13) (а + п ) гСп + c+ -s (а + п + ^)(а +п + 2)c+1 = ^^ n = ^ iU^... Из уравнения (12) определяется параметр а: 1+4 -1 V8 Л 1 а = - 2 + ст (14) -1 > 0. а У Вообще говоря, уравнение (12) имеет два корня: положительный и отрицательный. Однако при отрицательном решении, как будет показано ниже, функция цены P(r, т) приобретает свойства, которыми цена бескупонной облигации не обладает. Поэтому принимаем корень (14). Рассмотрим уравнение (13) для n = 0: а'(т) аc0 + а(т)2с1 [1 -s (а + 1)(а + 2)] = 0. С учетом равенства (12) его можно переписать как а'(т) = а(т)2 2 ю/а2, (15) где для краткости обозначено ю = cJcq. Равенство (15) представляет собой дифференциальное уравнение относительно функции о(т). Решение уравнения имеет вид: (16) а а(т) = - 2ют+^ с точностью до константы которая при необходимости определяется из свойств цены облигации. Заметим, что из равенства (15) следует, что а'(т) _2 ю а(т)2 а2 Теперь рассмотрим уравнение (13) для произвольного n > 1. Его можно записать как рекуррентное соотношение, определяющее коэффициент Cn+1 через коэффициент с„. V 2(а + n) ю (17) -c„. [s (а + n + 1)(а + n + 2) -1] а2 Заметим, что 2(а + 1)(n + а) 2 (а + n) 2 (а + n) 9 n + а (а + n + 1)(а + n + 2) } 2 (n +1) а(n + 2(а +1) n +1 ^n + Z --1 | а ' а(а +1) _ (ю 9)"-1 у-1 (к + а) 1х2х3 (1 + Z)(2 + ^ С" = С0 n! П(к + Z)' . где для краткости обозначено = 2(а + 1), 9 = ^/а. Таким образом, последовательность коэффициентов {cn, n = 0, 1, 2, . } выглядит следующим образом: ю 9 (1 + а) „ _ „ (ю9)2 (1 + а)(2 + а) а С0, c1 = с0ю = c0 ю 9 - , c2 = c0 Z 2 0 2 (1 + Z) Тогда решение (2) уравнения (11) можно представить в форме: [s (а + n + 1)(а + n + 2) -1] а2 а(т) (18) P(r, т) = C0 , ю9 а(т) ^а ю 9 а(т) Y 1 ""T^ + а) + r Ji + Ч r J ~п\ П(к + Z) n=2 ^ J ' к=0 Заметим, что среди специальных функций имеется так называемая вырожденная гипергеометрическая функция (функция Куммера) 1F1(x,j,z) (в обозначениях системы Wolfram Mathematica), которая определяется соотношением х n n i z 1F1(x, y, z) = 1 + Z - п(x +k -1) = 1 + I(y) Z - r(x +n) Z n\ пП(У + к -1) ВД Z n\ r(y + n) Используя эти обозначения, цену P(r, т) можно записать в виде: Г(3) V 1 L* а(т) Г Г(а + п)Л Чт) а(тУ а(т) P(r, т) = с. 1+ =с 1F1 а, Z, ю 9 Z - |ю 9 Г(а) Zn\I r J Г(3 + n) \ j По своим экономическим свойствам цена облигации как функция срока до погашения т является непрерывной монотонно убывающей функцией, которая для всякого r > 0 имеет пределы [9] lim P(r, т) = 1, lim P(r, x) = 0. Эти требования можно удовлетворить, определяя соответствующим образом до сих пор неопределенные константы со и Окончательное выражение для цены бескупонной облигации приобретает вид: (1 + а) л/я 1 (19) P(r, т) = а, 2(1 + а), - ■ 1 v sr sr т где 5 = 0,5 с1, а = 0,5 (л/1 + 4/s -1) > 0, а Г(х) - гамма функция. Здесь предполагается, что а > 0. Когда а < 0, используемая в формуле (19) гамма функция Г(х) может иметь нежелательные свойства. Например, при целочисленных отрицательных значениях аргумента она имеет неограниченные разрывы, на интервалах (2к, 2к + 1), к = 0, 1, 2, ..., она отрицательна и т.д., что не соответствует свойствам цены облигации. Поэтому отрицательные значения параметра а нежелательны. Обычно временную структуру на практике представляют не через цену облигации, а с помощью доходности. По определению доходность до погашения бескупонной облигации (кривая доходности) y(r, т) и доходность форвардных ставок (форвардная кривая) fr, т) определяются выражениями [10]: д ln P(r, т) дГ ln P(r, т) У(П т) = - (20) f (r, т) = т и, к сожалению, в компактном аналитическом виде не представляются и могут быть исследованы только численно. На рис. 2 представлены кривая доходности y(r, т) и форвардная кривая fr, т) для следующих значений параметров: s = 0,8; r = 0,08. Для представления кривых «целиком» для всего интервала возможных значений сроков до погашения т е (0, да) использовано нелинейное преобразование сроков до погашения: u = 1 - е-рт, которое отображает положительную полуось (0, да) в единичный интервал (0, 1). Принятое при расчетах численное значение р = ln10/30 = 0,07675 соответствует тому, что сроки до погашения от 0 до 30 отображаются в интервал (0, 0,9). В связи с таким преобразованием изменены и обозначения кривых: y(r, т) о Y(u), fr, т) о F(u). Y(u),F(u) 0.12 0.10 0.08 0.06 0.04 0.2 0.4 0.6 0.8 1.0 Рис. 2. Кривые доходности Y(u) (сплошные) и форвардные кривые F(u) (пунктирные) для модели CIR(1980) для значений ставки r е {0,05; 0,08; 0,12} 0.02 Как это обычно бывает, кривые стартуют от одного и того же значения доходности, равного процентной ставке r, и с увеличением срока до погашения стремятся при т ^ да (u ^ 1) к одному и тому же пределу, который в данном случае нулевого дрейфа равен нулю. 3. Модель с линейным дрейфом Предположим теперь, что полиномы (4) определены так, что из их коэффициентов отличаются от нуля только два: m\ = m > 0, S3 = s > 0, т.е. m(r) = mr, s(r) = sr3. В этом случае из равенства (9) можно получить соотношения: - c0 (1 -a(a +1) s) a(x)a = 0. [-(a + j) a'(x) a(x)j-1 Cj - a(x)j+1 cj+1 - (a + j) ma(x)j Cj + + (a + j +1) (a + j + 2) s a(x)j+1 cj+1 ] a(x)a = 0, j > 0. Из этих соотношений получаем следующие уравнения для a, а(х) и Cj, j > 1: 1 - a (a +1) s = 0. - a a'(x) a(x)-- a mC0 - a(x) c + (a + 1)(a + 2) sa(x) c1 = 0. - (a + j) a'(x) a(x)-1 Cj - (a + j) m Cj - a(x) cj+1 + + (a + j + 1)(a + j + 2) sa(x) j = 0, j > 1. Для параметра a получается квадратное уравнение и, соответственно, два решения - положительное и отрицательное: 1 a, = - 1 2 4+1-1 1 a9 = - 2 2 4+1+1 s s Функция а(х) удовлетворяет дифференциальному уравнению первого порядка, которое имеет аналитическое решение - a2 m - a m C ю = ■ c0 a(x) = - (em х -1) (s(1 + a) (2 + a) -1) w (em х - 1) 2 w Наконец, коэффициенты Cj, j > 1, удовлетворяют следующему рекуррентному соотношению: _ (a + j) Cj a'(x) + ma(x) j+1 j 2 s + js (3 + 2a) + s(2 + 3a + a 2) - 1 a(x)2 . Используя уже полученное аналитическое выражение для функции а(х), можно найти, что второй сомножитель этого выражения не зависит от х и имеет вид: a'(x) + ma(x) s (1 + a)(2 + a) - 1 -2-=-ю , a(x)2 a так что рекуррентное соотношение для Cj можно представить в форме: (a + j) CJ+1 = 2j s 2ю c,, P = 1 + J4 +1 = 2(1 + a) > a. (1 + j)(P + j) s a Используя полученные результаты, находим, что n-е слагаемое ряда (2) можно записать так: a + к Л П-1 ^ п a(x) - m Cn = 0 sr (em х-1) r п! к=0 Р+к Таким образом, цена бескупонной облигации P(r, х) может быть представлена в аналитическом виде с помощью вырожденной гипергеометрической функции: ( V ( \ m P(r, х) = с 0 1F s r(em х-1) s r(em х-1) m a, 2(1 + a), - \ Наконец, выбирая параметр Co так, чтобы удовлетворить поведению функции P(r, х), получим окончательное выражение для этой функции: ( V ( т \ a, 2(1 + a), - s r(em х- 1) m 1F 21+2ar(a + 3/2) s r(e™х- 1) Как и в предыдущем случае здесь используется положительный корень a. Функции y(r, х) и fr, х), определяемые формулами (20) через представление P(r, х), можно исследовать только численными методами. Правда, предельные значения этих функций могут быть найдены в аналитическом виде: lim y(r, х) = lim f (r, х) = r, lim y(r, х) = lim f (r, х) = a m. х^-0 х^-0 х^-да х^-да p r х) = (1 + m Y(U),F(U) Рис. 3. Функции Y(u) (сплошные) и F(u) (пунктирные), вычисленные по формулам (20) для значений ставки r е {0,08; 0,36; 0,6} Как видно, левый предел определяется только состоянием рынка и не зависит от параметров модели, а правый предел определяется только структурой модели и не зависит от состояния рынка в определенный момент времени. На рис. 3 эти функции в форме Y(u) и F(u) представлены для следующих значений параметров: 5 - 0,8; m - 0,5. 4. Модель с квадратичным дрейфом Пусть теперь полиномы m(r) и s(r) таковы, что p - 2, q - 3, т.е. 1 - p - 2 - q - - 1. Тогда компоненты суммы (9) будут отличаться от нуля только для j > -1, причем первое слагаемое отличается от нуля только для j > 0. В этом случае получаем следующую систему уравнений: при j - -1 -c0 - a m2c0 +a (а+1)s3c0 - 0 ; (21) при j - 0 2 3 -а a'(x) a(x)-1 c0 - a(x) c - Е (a + k - 1) щ a(x)k-1 ck_x + Е (a + k - 2) (a + k - 1) sk a(x)k-2 - 0 ; k-1 k-2 (22) при j - 1 2 3 - (a +1) a'(x) c1 - a(x)2 c2 - Е(a + k) mk a(x)kck + Е (a + k - 1)(a + k) Ska(x)k-1ck-1 - 0; k=0 k=1 (23) при j > 1 -(a + j) a'(x) a(x)j-1 cj - a(x)j+1 j - £(a + j + k - 1) mka(x)J+kAcJ+k_ 1 + £(a + j + k - 2)(a + j + k - 1) ska(x)+k-1 c^ - 0. (24) k-0 k-0 Из уравнения (21), которое в предположении, что С0 Ф 0, имеет вид: a (a + 1) S3 - am2 +1, определяется параметр a: a1 -T^(m2 - S3 - V4s3 + (m2 - s3)2)' a2 - T1 (m2 - S3 +V4s3 + (m2 - S3>2). (25) 2s3 2s3 Поскольку уравнение (21) - квадратное, оно имеет два корня, что означает, что решение уравнения (1) может иметь две составляющие вида (2), компромисс между которыми, а также начальное условие P(r, 0) - 1 могут повлиять на выбор коэффициента С0. Уравнение (22) является обыкновенным дифференциальным уравнением относительно функции a(x). Его решение имеет вид: а(т) =-X-, (26) ц + ехр[(т + а£,c0)(m -(а +1)s2)] где для компактности обозначено X = асо((1 + а^а - m\), ц = ci(i + (а + 1)ш2 - (а + 1)(а + 2)s3), а -постоянная интегрирования дифференциального уравнения, выбор которой делается в зависимости от свойств решения уравнения (1). Из уравнения (23) определяется коэффициент C2, а уравнение (24) можно рассматривать как основу для конструирования рекуррентной формулы вычисления коэффициентов c„+\ через предыдущие коэффициенты Cj, j < n. Рассмотрим сначала уравнение (24). Оно позволяет выразить коэффициент Cj+\ через предыдущие коэффициенты Cj, c;-\, c;-2 по формуле с _ а(т)(а + j)(a' (т) + a(x)(mi - (1 + а + j)^ )c]. - [д(т)(-от0 + (а + j)si )су- + (а + j - ](а + j -1) ^ 3+1 а(т)3[(1 + а + j) (2 + а + j) s3 - (1 + а + j) m2-1] Однако по определению коэффициентов Cn в выражении (2) они должны быть постоянными коэффициентами, не зависящими от переменной т. То есть что в формуле (27) правая часть равенства не должна зависеть от т. Это будет только тогда, когда то = 0, S0 = 0, s\ = 0, S2 = 0. Такое требование является необходимым условием существования нетривиального решения (2) и говорит о том, что нетривиальное решение имеет место не для любых полиномов m(r) и s(r) порядка 2 и 3 соответственно, а только для m(r) = m\r + m2r2, s(r) = S3r3. (28) Подстановка выявленных необходимых условий в формулу (27) для коэффициента Cn+\ приводит к рекуррентному соотношению а + n а'(т) + а(т) m^ (OQ} (1 + а + n) (2 + а + n) s3 - (1 + а + n) m2 -1 а(Т)2 Заметим, что знаменатель первого сомножителя правой части равенства (29) может быть представлен в виде (\ + а + n) (2 + а + n)s3 - (\ + а + n)m2 - \ = S3(\ + n)(P + n), где Pj = - (s - V4s3 + (m2 - s3 )2) при а = а, " (30) P = P2 = - (s3 +J4s3 + (m2 - s )2) при а = а • S3 При выполнении выявленных необходимых условий функция а(т), определяемая по формуле (26), несколько упрощается: а(т) =-X-, (3\) ц + ехр[(т + а с0) mj] где X = -m^C0, ц = (\ + (а + \)m2 - (а + \)(а + 2)s3)C\. Подставляя в правую часть формулы (29) явное выражение функции о(т), определяемой по формуле (20), получим а'(т) + а(т) m ц m s3 Рю а(т)2 X а где ю = C\/C0. При этом зависимость от переменной т в правой части формулы (29) исчезает. Таким образом, рекуррентная формула (29) для коэффициента Cn+\ преобразуется к окончательному виду: „ _ Р (а + п)ю cn (32) n+' а(1 + n)(P + n) Теперь обратимся к решению последнего уравнения (23), из которого нужно определять C2. Поскольку среди необходимых условий есть равенство S0 = 0, то уравнение (23) будет совпадать с уравнением (24) при n = 1, а следовательно, коэффициент C2 вычисляется по формуле (31) при n = L Получается, что если полиномы m(r) и s(r) порядка 2 и 3 соответственно определяются выражениями (28), решение уравнения (1) может быть представлено в виде суммы двух рядов типа (2), каждый из которых имеет следующую структуру: 21 f а (а +1) Л 2 а / й(т) a(x) юР | а r а у Р a(x) юр а 1 + с 0 Р (Р +1) r r \ + f a(x) юр V r а У 3! Р (Р + 1)(Р + 2) Используя вновь вырожденную гипергеометрическую функцию, полученный результат можно компактно записать в аналитическом виде Л 1 и-1 П "! к=0 а (а + 1)(а + 2) a(x) юр r а а + к Р + к - +. a(x) юР r а у л а a(x) а, Р, 1F (33) с 0 V r У V Как уже было сказано, поскольку уравнение (21) имеет два решения (25), решение уравнения (1) может состоять из двух компонентов вида (2) с различными наборами параметров (а, Р), значения которых определяются формулами (25) и (30): P(r, x) = с. а а 2 a(x) a(x) юР1 a(x) a(x) юР2 (34) 1F аь Рь 1 F1 а2, Р2х + с 01 02 r r r r а а Прежде чем конкретизировать решение, сделаем некоторый предварительный анализ. Вначале рассмотрим свойства диффузионного процесса, задаваемого дрейфом и волатильностью, определяемыми функциями (4) и (28). Согласно принятым предположениям процесс краткосрочной процентной ставки r(t), соответствующий этим функциям, описывается уравнением: dr(t) = (m1r(t) + m2r(t)2) dt +J2^ r(t)3/2 dt. Маргинальная плотность вероятностей этого процесса имеет вид: ^2-у -5/r f (r) = 5-y г* л, 5 = Ш > 0, у = < 2, s3 > 0, r > 0, r ' Г(2-у) s3 s3 где Г(х) - гамма функция. Принимая во внимание эти неравенства, обратим внимание, что параметры выражения (33) согласно формулам (25) и (30) принимают значения а1 < 0, а2 > 0, Р2 > 0, а Р1 может принимать положительные значения только тогда, когда параметр волатильности s3 > 4, что в реальных случаях практически не встречается. Как известно, цена облигации при r > 0 является монотонно убывающей функцией по переменной x е (0, да) от P(r, 0) = 1 до P(r, да) = 0. Поэтому выражение (33) должно иметь такие же свойства. Функция 1F1(x, y, z) имеет подходящие свойства только при x > 0, y > 0, z е (-да, 0). Поэтому первое слагаемое в представлении (34) должно отсутствовать. Кроме того, чтобы аргумент z у 1F1 принимал значения в интервале (-да, 0) при изменении x в интервале (0, да), необходимо постоянную интегрирования в выражении (20) определить равенством = Ы^ю^/а^га! Тогда - ш a(x) = ц + exp [(x + а ^ с0) m1] (em1x-1) s3 Наконец, для того чтобы выполнялось требование P(r, 0) = 1, необходимо, чтобы не определенный до сих пор параметр c0 определялся равенством c0 = Г(Р - а)/Г(Р). Таким образом, окончательный вид решения (2) уравнения (1) приобретает в рассматриваемом случае окончательный вид: Г(Р-а) ( а Ш, m а, Р, 1F Г(Р) rs3(em1% -1) P(r, x) = V r^e^-1) где параметры а и Р определяются с помощью формул (25) и (30): 4S3 + (ш2 - S3)2) > 0. Р = - (S3 +J - S3)2) > 0, а = (m2 -S Заметим, что это решение полностью совпадает с решением, полученным другим путем в [6], где в обозначениях авторов т\ = к9 - А,\ > 0, тг = - к - Л,2 < 0, S3 = с2/2. В принципе, используя выражение (23), можно по формулам (20) найти аналитические выражения для кривой доходности y(r, т) и форвардной кривой fir, т). Однако эти выражения являются очень громоздкими, и практичнее привлекать численные методы для выражения этих функций для необходимых числовых параметров. На рис. 4 представлены функции Y(u) и F(u), которые взаимно однозначно связаны с кривыми y(r, т) и fr, т), как объяснено выше, для модели с квадратичным дрейфом при следующих значениях параметров: s = 0,8; m\ = 0,2; тг = 1. Предельные значения этих функций, как и в предыдущей модели, могут быть найдены в аналитическом виде: limY(u) = limF(u) = r, limY(u) = limF(u) = a u-0 u-0 u-u- 0.0 0.2 0.4 0.6 0.8 1.0 Рис. 4. Кривые доходности Y(u) (сплошные) и форвардные кривые F(u) (пунктирные) для модели с квадратичным дрейфом для значений ставки r е {0,08; 0,25; 0,6} Заключение В статье представлены модели, для которых могут быть найдены кривые доходности бескупонных облигаций и соответствующие форвардные кривые, не относящиеся к классу аффинных моделей. К сожалению, модели, допускающие такие решения, немногочисленны и, в частности, включают в себя некоторые известные модели: модель CIR(\980) [5] и модель Ана-Гао [6]. Сформулируем требования к структуре модели краткосрочной процентной ставки, которая допускала бы получение временной структуры цены облигации в форме (2). Параметры ряда (2) определяются уравнением (9), из которого, собственно, получается система уравнений относительно неизвестных a, а(т) и сп, п = 0, 1, 2, ... 1. Для получения нетривиального решения (т.е. для наличия Сп Ф 0) необходимо, чтобы степени p и q полиномов m(r) и s(r), определяющих дрейф и волатильность краткосрочной процентной ставки, удовлетворяли одному из следующих условий: {p < 2, q = 3}, {p = 2, q < 3}, {p > 2, q = p + 1}. В этих случаях находятся уравнения, из которых определяется положительный параметр a. 2. Другое необходимое условие связано с существованием функции о(т), не зависящей от индекса суммирования ряда (2). 3. Кроме того, необходимо, чтобы коэффициенты {сп} не зависели от переменной т. Одновременное выполнение этих необходимых условий существенно сужает семейство моделей, для которых решение уравнения временной структуры (1) имеет вид (2). Вообще говоря, удается найти только один основной случай {p = 2, q = 3} - модель с квадратичным дрейфом Ана-Гао [6], из которого при т2 - 0 получается модель с линейным дрейфом. И далее при т\ - 0 имеем модель CIR(\980) [5] с нулевым дрейфом.

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

уравнение временной структуры доходности, цена бескупонной облигации, модель CIR(1980), модель Ана-Гао, кривая доходности, форвардная кривая, equation of term structure of yield, price of zero-coupon bond, CIR(1980) model, Ahn-Gao model, yield curve, forward curve

Авторы

ФИООрганизацияДополнительноE-mail
Медведев Геннадий АлексеевичБелорусский государственный университетпрофессор, доктор физико-математических наук, профессор факультета прикладной математики и информатикиMedvedevGA@bsu.by
Всего: 1

Ссылки

Merton R. Continuous-Time Finance. Oxford : Blackwell Publ., 1990. 732 p.
Duffee G.R. Term Premia and Interest Rate Forecasts in Affine Models // J. of Finance. 2002. V. 57, No. 1. P. 405-443.
Медведев Г.А. Полиномиальные модели временной структуры доходности // Вестник Томского государственного универ ситета. Управление, вычислительная техника и информатика. 2017. № 2 (39). С. 39-48.
The Handbooks of fixed income securities / F. Fabozzi, T. Fabozzi, eds. New York : IRWIN Prof. Publ., 1995. 1402 p.
Cox J.C., Ingersoll J.E., Ross S.A. An analysis of variable rate loan contracts // J. of Finance. 1980. V. 35. P. 389-403.
Ahn D.-H., Gao B. A parametric nonlinear model of term structure dynamics // The Review of Financial Studies. 1999. V. 12(4). P. 721-762.
Vasicek O.A. An equilibrium characterization of the term structure // J. of Financial Economics. 1977. V. 5. P. 177-188.
Медведев Г.А. Плотности вероятностей процессов процентных ставок доходности // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2016. № 3 (36). C. 35-48.
Hull J.C. Options, futures, and other derivatives. Boston : Prentice Hall, 1998. 841 p.
Медведев Г.А. Стохастические процессы финансовой математики. Минск : Белорус. гос. ун-т, 2005. 243 c.
 Об одном семействе неаффинных моделей временной структуры доходности | Вестн. Том. гос. ун-та. Управление, вычислительная техника и информатика. 2018. № 42. DOI: 10.17223/19988605/42/6

Об одном семействе неаффинных моделей временной структуры доходности | Вестн. Том. гос. ун-та. Управление, вычислительная техника и информатика. 2018. № 42. DOI: 10.17223/19988605/42/6