Об одном семействе неаффинных моделей временной структуры доходности Текст научной статьи по специальности «Математика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА 2018 Управление, вычислительная техника и информатика № 42

УДК 336.51

DOI: 10.17223/19988605/42/6

Г.А. Медведев

ОБ ОДНОМ СЕМЕЙСТВЕ НЕАФФИННЫХ МОДЕЛЕЙ ВРЕМЕННОЙ СТРУКТУРЫ ДОХОДНОСТИ

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

Ключевые слова: уравнение временной структуры доходности; цена бескупонной облигации; модель СГЯ(1980); модель Ана-Гао; кривая доходности; форвардная кривая.

Предположим, что состояние финансового рынка описывается процентной ставкой г(0, которая следует однородному по времени марковскому процессу, порождаемому стохастическим дифференциальным уравнением [1]

Жг(0 = ц(г(0) ж + а(г(0) dw(t) с функцией дрейфа ц(х), функцией волатильности с(х) и стандартным винеровским процессом м>(?). Для удобства рассуждений будем обозначать функцию дрейфа т(г) = ц(г) - А,(г)а(г) и функцию диффузии 5(г) = 0,5 а2(г). Здесь А,(г) - так называемая рыночная цена риска [2]. Ранее [3] рассматривалась задача определения временной структуры доходности бескупонной облигации [4], когда функции т(г) и s(г) являются полиномами. Выяснялось, могут ли в этом случае кривые доходности быть полиномами или степенными рядами по переменной г. Оказалось, что это имеет место, если только т(г) и s(г) -полиномы не более чем первой степени. В этом случае модели временной структуры доходности являются аффинными.

В настоящей статье рассматривается подобная задача, но временная структура цены бескупонной облигации ищется в виде функционального ряда, отличающегося от степенного. Выяснено, что для некоторых случаев такие решения существуют. Получающаяся временная структура оказывается неаффинной и описывается вырожденными гипергеометрическими функциями. В это семейство входят такие известные модели процентных ставок, как модель СШ.(1980) [5] и модель Ана-Гао [6].

1. Общее уравнение для цены облигации и ее компонентов

Рассмотрим уравнение временной структуры для цены бескупонной облигации Р(г, т) [7]:

+ т(г) + КГ) -г?(г,т) = 0, Р(г, 0) = 1. (1)

дт дг дг

Здесь т(г) - функция дрейфа краткосрочной процентной ставки, а s(г) - квадрат ее волатильности.

Будем искать решение этого уравнения в форме

а(х)

Р(г, х) = £

и-0

(2)

где а(х), а и е„, п = 0, 1, 2, ..., - функция и коэффициенты, подлежащие определению. Соответствующие производные, использованные в уравнении (1), имеют вид:

дР(г, х) _ а'(х) ■

а(х)

-У (а + п)1

дх а(х) П=0 V г

д 2 Р(г, х)_ 1

\а+п

д Р(г, х)_ 1

дг

а(х)

У (а + п)

/ , ч\а+п+1

' а(х) х

п=0

V '

дг 2 а(х)2 п=о

У (а + п)(а + п +1)

/ , ч\а+п+2

а(х)

(3)

Предположим, что дрейф и волатильность краткосрочной процентной ставки таковы, что функции т(г) и s(г) являются полиномами порядкар и q соответственно:

Р ч

т(г) = Уткгк, *(г) = У^• (4)

к=0 к=0

Прежде чем подставлять явные выражения производных (3) и полиномов (4) в уравнение (1), преобразуем второе и третье слагаемые в левой части уравнения (1) к более удобной форме.

\ а+п-к+1 т ( Р

т (г) ^^ = -УУ(а+п) шка{хТп 11

д г „=п !г=п V г

с

= -У У (а + 1+к -1)тка(х)

а+ 1+к-1

V |Vа+1

1+к -1

1=1-Р \к =тах{0,1-;}

( ч

д2р(гх) ® Ч (а+п-к+2

•<г) —^ =УУ(а+п)(а+п+¡КФГЧI-I = У У (а+] +к - 2)(а+]+к - ¡Мх)^-2^

д г п=0 к=0 Vг / 1=2-ч \=пмх{0,2-1}

Запишем в аналогичной форме также первое и четвертое слагаемые левой части уравнения (1):

дР(г , х) дх

= -У((а +1) а'(х) а(х)а+1-1 с,- )Г1

1=0

а+1

- гР(г, х) = - У (

= -М а(х)

1=-1

а,+1+1

с1+1

)( 1

а+1

(5)

(6)

(7)

(8)

Теперь, подставляя выражения (5)-(8) в уравнение (1), получим

( р

У -1(1 |0)(а +1) а'(х) а(х)а+1-1 ^ -1(11 -1) а(х)а^41 с^ -1(111 -р) У (а +1 + к-1) тк а(х)а^+к-1 с^к-1 +

1 \ к=тах{0,1-1}

+10 |2 - Ч) У (а +1 + к - 2) (а +1 + к -1) ^ а(х)а+1+к-2 с ^-2

к=тах {0,2-1}

= 0.

(9)

Определенную сложность в выражении (9) вызывает тот факт, что суммирование по индексу ] для каждого слагаемого начинается по-разному: для первого слагаемого ] > 0, для второго слагаемого ] > -1, для третьего слагаемого ] > 1 - р, для четвертого слагаемого ] > 2 - q. Поэтому в выражениях слагаемых появились множители | к), представляющие индикаторные функции, равные единице, если ] > к, и нулю в противном случае. Равенство (9) должно выполняться равномерно по переменной г. При этом, поскольку функции г- (/ = 0, ±1, ±2, ...) являются линейно независимыми, коэффициенты перед этими функциями в выражении (9) должны быть равными нулю. Это приводит к системе уравнений для неизвестных параметров а, а(х) и Сп, п = 0, 1, 2, ..., в представлении (2) решения уравнения (1), если оно существует в таком виде. Заметим, что каждое слагаемое в каждом элементе суммы (9) имеет ненулевой множитель а(х)а, поэтому для упрощения на него можно сократить во всех элементах суммы.

с

п

г

с

п

п

г

\

г

2. Модель с нулевым дрейфом

Среди моделей процессов краткосрочной ставки г(0 с нулевым дрейфом широко известна модель СШ.(1980) [5], в которой ставка порождается в общем случае диффузионным процессом

¿г = с г1 dw. (10)

Для такого процесса плотность вероятностей / (х) выражается в виде [8]:

2

/ (х) = (у-1)(2у-1)-

X-1

V го )

\Г Л-2у X

х > г0 > 0.

V г0 ^

График этой плотности представлен на рис. 1 для случая у = 1,5, который чаще всего встречается в литературе. Именно в таком виде эту модель предлагали и авторы оригинальной статьи.

гсЛх)

0.30

0.25

0.20

0.15

0.10

0.05

^-,-,-,-,-,-,—,-,-,-,-,-,-,-,-,-,_ хг0

2 4 6 8 10

Рис. 1. Плотность вероятностей краткосрочной ставки для модели СШ.(1980)

Несмотря на то что модель известна давно, до сих пор не была описана временная структура ее бескупонной доходности. Оказывается, предлагаемый способ нахождения временной структуры позволяет это сделать. Примем в уравнении (10) у = 1,5 и 5 = 0,5 с2. Уравнение (1) для цены бескупонной облигации Р(г, т) приобретает вид

дР(г, т) 3 д2Р(г, т)

дт

- + sг

д г2

- гР(г, т) = о, Р(г, 0) = 1.

(11)

Будем искать решение этого уравнения в форме (2). Соответствующие производные имеют вид:

дР(г, т) а'(т)

дт

а(т) ±

X (а+п)\

п=0

а(т)

д 2 Р(г, т)

дг 2

1

—- X (а + п)(а+ п +1)

а(т) п=0

а(т)

а+п+2

После подстановки этих выражений в уравнение (11) получаем следующее равенство:

а'(т) а(т) Х

X (а + п)

а(т)

- а(т)Х

а(т)

а+п-1

а(т)Х (а + п)(а + п +1)|

п=0 ^ у п=0 ^ у п=0

Это равенство можно переписывается в более удобном виде:

а(т)

а+п-1

сп = 0 .

п=0

(а + п)а (т) Си + а(т) еп+1 - 5а(т)(а + п + 1)(а + п + 2)еп а(т)

а(т)

+ а(т)(1 -5а(а+1)) Сс

а(т)"

= 0.

Поскольку выражения (а(т)/г)к как функции переменной г для различных значений к являются линейно независимыми, а равенство должно выполняться равномерно по г, то коэффициенты перед этими выражениями для различных к должны быть равны нулю. И мы получаем систему уравнений относительно неизвестных а, а(т) и с„, п = 0, 1, 2, ...:

5а(а + 1) = 1, (12)

а' (т)

(а + п )

а(т)

2Сп + сп+1 -5 (а + п + 1)(а +п + 2)с„+1 = 0, п = С0, 1,2,„.

(13)

г

0

С

Сп.

п

г

г

Сп + 5

С

г

а-1

г

г

Из уравнения (12) определяется параметр а:

1

а = — 2

i+4" i

Л

Л

+ ст

-1

а

> 0.

(14)

Вообще говоря, уравнение (12) имеет два корня: положительный и отрицательный. Однако при отрицательном решении, как будет показано ниже, функция цены Р(г, х) приобретает свойства, которыми цена бескупонной облигации не обладает. Поэтому принимаем корень (14). Рассмотрим уравнение (13) для п = 0:

а'(х) ас0 + а(х)2с1 [1 -5 (а + 1)(а + 2)] = 0. С учетом равенства (12) его можно переписать как

а'(х) = а(х)2 2 ю/а2, (15)

где для краткости обозначено ю = С1/С0. Равенство (15) представляет собой дифференциальное уравнение относительно функции а(х). Решение уравнения имеет вид:

а(т) = -

а

(16)

2юх+^

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

а'(х) _2 ю а(х)2 а2

Теперь рассмотрим уравнение (13) для произвольного п > 1. Его можно записать как рекуррентное соотношение, определяющее коэффициент Сп+1 через коэффициент с„.

2(а + n) ю

[s (а + n + 1)(а + n + 2) -1] а2

-c„.

(17)

Заметим, что

2 (а + n)

2 (а + n)

2(а + 1)(n + а)

[s (а + n + 1)(а + n + 2) -1] а2

n + а

(а + n + 1)(а + n + 2) } 2 (n +1) а(n + 2(а +1) n +1 ^n + $

--1 | а '

а(а +1)

где для краткости обозначено = 2(а + 1), 9 = ^/а. Таким образом, последовательность коэффициентов {сп, п = 0, 1, 2, ... } выглядит следующим образом:

ю 9 (1 + а) „ _ „ (ю9)2 (1 + а)(2 + а)

а

Со, c1 = с0ю = c0 ю0 — , c2 = c0

$ 2 0 2 (1 + $) Тогда решение (2) уравнения (11) можно представить в форме:

'а(т)

____ (ю б)"-1 т-1 (k + а)

1х2х3 (1 + $)(2 + $)' С" = Со n! ¿¿(k + $)' -

P(r, т) = Со

, ю0 а(т) ^а ю 6 а(т) Y 1 T_r(k + а)

+ r Ji + ¿2.1 r J ~n\ I\(k + $)

(18)

n=2 ^ J ' k=0

Заметим, что среди специальных функций имеется так называемая вырожденная гипергеометрическая функция (функция Куммера) \F\(x,yz) (в обозначениях системы Wolfram Mathematica), которая определяется соотношением

х n n

i z

1F1(x, * z) = 1 + У О(x + k -1) = 1 + m Y zг(x +n) 1 У n\ ¿{(y + k-1) Г(x) ¿1 n\ Г(y + n)

Используя эти обозначения, цену Р(г, х) можно записать в виде:

V 1 а(х)Г Г(а + п)Л

P(r, т) = с

'а(т)

1 + -

= С

а(т)

1F а, ю 0

а(т)

У - |ю 9 Г(а) Уп! I г ) Г(^ + п)

V )

По своим экономическим свойствам цена облигации как функция срока до погашения х является непрерывной монотонно убывающей функцией, которая для всякого г > 0 имеет пределы [9]

V

/

0

lim P(r, т) = 1, lim P(r, x) = 0.

т—т—®

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

P(r, т) =

(1 + а) л/я

21+2а Г(а +1,5)

f Y \а

1

а, 2(1 + а), - ■

sr т

(19)

где 5 = 0,5 с2, а = 0,5 (л/1 + 4/5 -1) > 0, а Г(х) - гамма функция. Здесь предполагается, что а > 0. Когда

а < 0, используемая в формуле (19) гамма функция Г(х) может иметь нежелательные свойства. Например, при целочисленных отрицательных значениях аргумента она имеет неограниченные разрывы, на интервалах (2к, 2к + 1), к = 0, 1, 2, ..., она отрицательна и т.д., что не соответствует свойствам цены облигации. Поэтому отрицательные значения параметра а нежелательны.

Обычно временную структуру на практике представляют не через цену облигации, а с помощью доходности. По определению доходность до погашения бескупонной облигации (кривая доходности) у (г, т) и доходность форвардных ставок (форвардная кривая) Дг, т) определяются выражениями [10]:

У^ т) = -

ln P(r, т)

f (r, т) = -

д ln P(r, т)

дт

(20)

и, к сожалению, в компактном аналитическом виде не представляются и могут быть исследованы только численно. На рис. 2 представлены кривая доходностиу(г, т) и форвардная кривая_Дг, т) для следующих значений параметров: 5 = 0,8; г = 0,08. Для представления кривых «целиком» для всего интервала возможных значений сроков до погашения т е (0, да) использовано нелинейное преобразование сроков до погашения: и = 1 - е-рт, которое отображает положительную полуось (0, да) в единичный интервал (0, 1). Принятое при расчетах численное значение р = 1п10/30 = 0,07675 соответствует тому, что сроки до погашения от 0 до 30 отображаются в интервал (0, 0,9). В связи с таким преобразованием изменены и обозначения кривых: у(г, т) о Ди), Дг, т) о ^(и).

У(и),Р(и) 0.12

0.10

0.08

0.06

0.04

0.02

0.2 0.4 0.6 0.8 1.0

Рис. 2. Кривые доходности У(и) (сплошные) и форвардные кривые ^и) (пунктирные) для модели СШ.(1980) для значений ставки г е {0,05; 0,08; 0,12}

Как это обычно бывает, кривые стартуют от одного и того же значения доходности, равного процентной ставке г, и с увеличением срока до погашения стремятся при т ^ да (и ^ 1) к одному и тому же пределу, который в данном случае нулевого дрейфа равен нулю.

3. Модель с линейным дрейфом

Предположим теперь, что полиномы (4) определены так, что из их коэффициентов отличаются от нуля только два: ш\ = т > 0, 5з = 5 > 0, т.е. т(г) = тг, 5(г) = 5г3. В этом случае из равенства (9) можно получить соотношения:

1

v sr ^

т

- с0 (1 -а(а +1) 5) а(т)а = 0.

[—(а + /) а'(т) а(т)1 с1 — а(т)+1 еу+1 — (а + 1) ша(т)1 е^ +

+ (а +1 +1) (а +1 + 2) 5 а(т)1+1 с]+1 ] а(т)а = 0, 1 > 0.

Из этих соотношений получаем следующие уравнения для а, а(т) и Су, у > 1:

1 — а (а +1) 5 = 0.

— а а'(т) а(т)-— а шс0 — а(т) с + (а + 1)(а + 2) 5а(х) с1 = 0.

— (а +1) а'(т) а(т)—1 с^ — (а +1) ш с^ — а(т) Сj+1 +

+ (а + У + 1)(а + ] + 2) 5а(т) су+1 = 0, у > 1.

Для параметра а получается квадратное уравнение и, соответственно, два решения - положительное и отрицательное:

1

а, = — 1 2

4+1-1

1

а = —

2 2

4+1+1

Функция а(т) удовлетворяет дифференциальному уравнению первого порядка, которое имеет аналитическое решение

а(т) = -

- а m

- а2 m

(em т -1) (5(1 + а)(2 + а) -1) ю (em т -1) 2 ю

ю = ■

c0

Наконец, коэффициенты Су,у > 1, удовлетворяют следующему рекуррентному соотношению:

_ (а + ;') с1 а'(т) + ша(т)

7+1 125 + 75 (3 + 2а) + 5(2 + 3а + а 2) — 1 а(т)2 . Используя уже полученное аналитическое выражение для функции а(т), можно найти, что второй сомножитель этого выражения не зависит от т и имеет вид:

а'(т) + ша(т) 5 (1 + а)(2 + а) — 1

-2-=-ю ,

а(т)2 а

так что рекуррентное соотношение для Су можно представить в форме:

(а + Л

cj+1 =

(1 + j)(ß + j) 5 а

2 j '

2ю c,, ß = 1 + J4 +1 = 2(1 + а) > а.

Используя полученные результаты, находим, что п-е слагаемое ряда (2) можно записать так:

а(т)

sr (em х-1)

л П—1

^ п

а + к

п! к=0о ß+к

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

( V ( \

P(r, т) = с 0

m

s r(em т-1)

1F

а, 2(1 + а), -

m

s r(em т-1)

Наконец, выбирая параметр Со так, чтобы удовлетворить поведению функции Р(г, т), получим окончательное выражение для этой функции:

( V ( т \

а, 2(1 + а), —

р т) = +

21+2аГ(а + 3/2)

m

1F

m

s r(em т-1)

s r(e™т-1)

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

Функции y(r, т) и fr, т), определяемые формулами (20) через представление P(r, т), можно исследовать только численными методами. Правда, предельные значения этих функций могут быть найдены в аналитическом виде:

lim y(r, т) = lim f (r, т) = r, lim y(r, т) = lim f (r, т) = а m.

т^-0 т^-0 т^-го т^-го

s

s

c

s

- m

cn = c

0

r

\

У(и),Р(и)

Рис. 3. Функции 7(и) (сплошные) и Е(и) (пунктирные), вычисленные по формулам (20) для значений ставки г е {0,08; 0,36; 0,6}

Как видно, левый предел определяется только состоянием рынка и не зависит от параметров модели, а правый предел определяется только структурой модели и не зависит от состояния рынка в определенный момент времени. На рис. 3 эти функции в форме У(и) и ^(и) представлены для следующих значений параметров: 5 = 0,8; т = 0,5.

4. Модель с квадратичным дрейфом

Пусть теперь полиномы т(г) и 5(г) таковы, что р = 2, ц = 3, т.е. 1 - р = 2 - ц = - 1. Тогда компоненты суммы (9) будут отличаться от нуля только для у > -1, причем первое слагаемое отличается от нуля только дляу > 0. В этом случае получаем следующую систему уравнений: при у = -1

- с0 - а т2с0 + а (а+1)53с0 = 0 ; (21)

при у = 0

2 3

- а а'(т) а(т)-1 с0 - а(т) с - X (а + к -1) щ а(т)к-1 с^ + X (а + к - 2) (а + к -1) 5, а(т)к-2 с4_2 = 0 ;

к=1 к=2 (22)

при у = 1

2 3

- (а +1) а'(т) с1 - а(т)2 с2 (а + к) тк а(т)кск + Х (а + к -1)(а + к) 5ка(т)к-1Ск_1 = 0;

к=0 к=1 (23)

при у > 1

- (а + у) а'(т) а(т)у-1 с. - а(т)у+1 у -¿(а + у + к-1) т^С + £(а + у + к - 2)(а + у + к -1) ^(т)^-2 с+к_г = 0. (24)

к=0 к=0

Из уравнения (21), которое в предположении, что С0 Ф 0, имеет вид: а (а + 1) 53 = ат2 +1, определяется параметр а:

а1 =~^(т2 -53-д/453 + (т2 -53?)> а2 =-^(т2 -53 + д/453 + (т2 -53>2). (25)

253 253

Поскольку уравнение (21) - квадратное, оно имеет два корня, что означает, что решение уравнения (1) может иметь две составляющие вида (2), компромисс между которыми, а также начальное условие Р(г, 0) = 1 могут повлиять на выбор коэффициента С0.

Уравнение (22) является обыкновенным дифференциальным уравнением относительно функции а(т). Его решение имеет вид:

а(т) =---, (26)

ц + exp[(т + а£,е0)(щ -(а +1)52)]

где для компактности обозначено - = асо((1 + а^а - т\), ц = с\(\ + (а + 1)т2 - (а + 1)(а + 2)^э), а -постоянная интегрирования дифференциального уравнения, выбор которой делается в зависимости от свойств решения уравнения (1).

Из уравнения (23) определяется коэффициент С2, а уравнение (24) можно рассматривать как основу для конструирования рекуррентной формулы вычисления коэффициентов с„+\ через предыдущие коэффициенты с,] < п. Рассмотрим сначала уравнение (24). Оно позволяет выразить коэффициент с+\ через предыдущие коэффициенты с, с;-\, с-2 по формуле с _ а(т)(а + у)(а' (т) + а(х)(ш1 - (1 + а + j)s2 )с]. - [д(т)(-от0 + (а + уХ )с].-1 + (а + у - 2>0су_2 ](а + у -1) ^ 3+1 а(т)3[(1 + а + у) (2 + а + у)53 -(1 + а + у)т2-1]

Однако по определению коэффициентов Сп в выражении (2) они должны быть постоянными коэффициентами, не зависящими от переменной т. То есть что в формуле (27) правая часть равенства не должна зависеть от т. Это будет только тогда, когда то = 0, sо = 0, 5\ = 0, S2 = 0. Такое требование является необходимым условием существования нетривиального решения (2) и говорит о том, что нетривиальное решение имеет место не для любых полиномов т(г) и s(r) порядка 2 и 3 соответственно, а только для

т(г) = т\г + т2Г2, ^(г) = $эг3. (28)

Подстановка выявленных необходимых условий в формулу (27) для коэффициента Сп+\ приводит к рекуррентному соотношению

а + п а'(т) + а(т) т^ (ОС)}

(1 + а + п) (2 + а + п) 53 - (1 + а + п) т2 -1 а(Т)2

Заметим, что знаменатель первого сомножителя правой части равенства (29) может быть представлен в виде

(\ + а + п) (2 + а + п^з - (\ + а + п)т2- \ = 5з(! + п)(Р + п),

где

Р; = — (5 - д/4я3 + (т - 53)2) приа = а,

" (30)

Р =

Р2 = — (5 + д/4я3 + (т - 53 )2) при а = а2.

53

При выполнении выявленных необходимых условий функция а(т), определяемая по формуле (26), несколько упрощается:

а(т) =---, (3\)

ц + exp[(т + а % с0) тх]

где - = -т\ас0, ц = (\ + (а + \)т2 - (а + \)(а + 2)^з)с\. Подставляя в правую часть формулы (29) явное выражение функции о(т), определяемой по формуле (20), получим

а'(т) + а(т) т ц т 53 Рю а(т)2 - а

где ю = с\/с0. При этом зависимость от переменной т в правой части формулы (29) исчезает. Таким образом, рекуррентная формула (29) для коэффициента сп+\ преобразуется к окончательному виду:

„ _ р (а +п)ю сп (32)

п+' а(1 + п)(Р + п)

Теперь обратимся к решению последнего уравнения (23), из которого нужно определять с2. Поскольку среди необходимых условий есть равенство ^ = 0, то уравнение (23) будет совпадать с уравнением (24) при п = 1, а следовательно, коэффициент с2 вычисляется по формуле (31) при п = !.

Получается, что если полиномы т(г) и s(r) порядка 2 и 3 соответственно определяются выражениями (28), решение уравнения (1) может быть представлено в виде суммы двух рядов типа (2), каждый из которых имеет следующую структуру:

а(т)

г

/

1 +

а(т) юР | а г а ) Р

а(т) юР а

г

21 ( а (а +1) Л 2

Р (Р +1)

+ (а(т) юР I г а ) 3!

(

а (а + 1)(а + 2)

Л

а(т) юР г а

1 п—1

1 п

П! к=0

а + к Р + к

\

- +.

Р (Р + 1)(Р + 2)

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

а(т)

л а

V Г

1Ъ

а, Р,

а(т) юР г а)

(33)

Как уже было сказано, поскольку уравнение (21) имеет два решения (25), решение уравнения (1) может состоять из двух компонентов вида (2) с различными наборами параметров (а, Р), значения которых определяются формулами (25) и (30):

Р(г, т) = с

01

а(т)

1Ъ

аь Рь

а(т) юР1

а

+ с

02

а(т)

1Ъ

а2, ^

а(т) юР2

а

(34)

Прежде чем конкретизировать решение, сделаем некоторый предварительный анализ. Вначале рассмотрим свойства диффузионного процесса, задаваемого дрейфом и волатильностью, определяемыми функциями (4) и (28). Согласно принятым предположениям процесс краткосрочной процентной ставки г(0, соответствующий этим функциям, описывается уравнением:

<яГ(0 = (т1г(0 + т2г(02) & +^2^ г(03/2 Л. Маргинальная плотность вероятностей этого процесса имеет вид:

с2—у —5/г

/ (г) = * л, 5= ш > 0, у = < 2, 53 > 0, г > 0, г ' Г(2 — у) 53 53

где Г(х) - гамма функция. Принимая во внимание эти неравенства, обратим внимание, что параметры выражения (33) согласно формулам (25) и (30) принимают значения а1 < 0, а2 > 0, Р2 > 0, а Р1 может принимать положительные значения только тогда, когда параметр волатильности 53 > 4, что в реальных случаях практически не встречается.

Как известно, цена облигации при г > 0 является монотонно убывающей функцией по переменной т е (0, да) от Р(г, 0) = 1 до Р(г, да) = 0. Поэтому выражение (33) должно иметь такие же свойства. Функция 1-^\(х, у, г) имеет подходящие свойства только при х > 0, у > 0, г е (—да, 0). Поэтому первое слагаемое в представлении (34) должно отсутствовать. Кроме того, чтобы аргумент г у ^ принимал значения в интервале (—да, 0) при изменении т в интервале (0, да), необходимо постоянную интегрирования в выражении (20) определить равенством = 1п(Рю53)/ас0ть Тогда

а(т) =

X

ц + exp[(т + а^ с0) ш1] (еш1т — 1) 53

Наконец, для того чтобы выполнялось требование Р(г, 0) = 1, необходимо, чтобы не определенный до сих пор параметр С0 определялся равенством С0 = Г(Р — а)/Г(Р). Таким образом, окончательный вид решения (2) уравнения (1) приобретает в рассматриваемом случае окончательный вид:

Р(г, т) =

Г(Р —а)

(

Г(Р)

ш,

V ^(е^—1)

1 ъ

а, Р,

ш

г53(ещт — 1)

где параметры а и Р определяются с помощью формул (25) и (30):

а = ^(ш2 — 5 +44^+Ш2

— 53)2) > 0,

Р = —(53 +л/

453 + (Ш2 — 53)2) > 0.

а

с

0

с

0

V

а

а

2

г

г

г

г

— ш

а

3

Заметим, что это решение полностью совпадает с решением, полученным другим путем в [6], где в обозначениях авторов т\ = к9 - А,\ > 0, тг = - к - А,г < 0, S3 = с2/2. В принципе, используя выражение (23), можно по формулам (20) найти аналитические выражения для кривой доходности y(r, т) и форвардной кривой fr, т). Однако эти выражения являются очень громоздкими, и практичнее привлекать численные методы для выражения этих функций для необходимых числовых параметров. На рис. 4 представлены функции Y(u) и F(u), которые взаимно однозначно связаны с кривыми y(r, т) и fr, т), как объяснено выше, для модели с квадратичным дрейфом при следующих значениях параметров: s = 0,8; m\ = 0,2; тг = 1. Предельные значения этих функций, как и в предыдущей модели, могут быть найдены в аналитическом виде:

lim Y(u) = lim F(u) = r, lim Y(u) = lim F(u) = a m1.

u—-0 u—0 u—u—

0.0 0.2 0.4 0.6 0.8 1.0

Рис. 4. Кривые доходности У(ы) (сплошные) и форвардные кривые ^и) (пунктирные) для модели с квадратичным дрейфом для значений ставки г е {0,08; 0,25; 0,6}

Заключение

В статье представлены модели, для которых могут быть найдены кривые доходности бескупонных облигаций и соответствующие форвардные кривые, не относящиеся к классу аффинных моделей. К сожалению, модели, допускающие такие решения, немногочисленны и, в частности, включают в себя некоторые известные модели: модель CIR(1980) [5] и модель Ана-Гао [6]. Сформулируем требования к структуре модели краткосрочной процентной ставки, которая допускала бы получение временной структуры цены облигации в форме (2).

Параметры ряда (2) определяются уравнением (9), из которого, собственно, получается система уравнений относительно неизвестных а, а(х) и сп, п = 0, 1, 2, ...

1. Для получения нетривиального решения (т.е. для наличия Сп Ф 0) необходимо, чтобы степени p и q полиномов т(г) и s(r), определяющих дрейф и волатильность краткосрочной процентной ставки, удовлетворяли одному из следующих условий: {p < 2, q = 3}, {p = 2, q < 3}, {p > 2, q = p + 1}. В этих случаях находятся уравнения, из которых определяется положительный параметр а.

2. Другое необходимое условие связано с существованием функции а(х), не зависящей от индекса суммирования ряда (2).

3. Кроме того, необходимо, чтобы коэффициенты {сп} не зависели от переменной х.

Одновременное выполнение этих необходимых условий существенно сужает семейство моделей, для которых решение уравнения временной структуры (1) имеет вид (2).

Вообще говоря, удается найти только один основной случай {p = 2, q = 3} - модель с квадратичным дрейфом Ана-Гао [6], из которого при m2 ^ 0 получается модель с линейным дрейфом. И далее при mi ^ 0 имеем модель CIR(1980) [5] с нулевым дрейфом.

ЛИТЕРАТУРА

1. Merton R. Continuous-Time Finance. Oxford : Blackwell Publ., 1990. 732 p.

2. Duffee G.R. Term Premia and Interest Rate Forecasts in Affine Models // J. of Finance. 2002. V. 57, No. 1. P. 405-443.

3. Медведев Г.А. Полиномиальные модели временной структуры доходности // Вестник Томского государственного универ-

ситета. Управление, вычислительная техника и информатика. 2017. № 2 (39). С. 39-48.

4. The Handbooks of fixed income securities / F. Fabozzi, T. Fabozzi, eds. New York : IRWIN Prof. Publ., 1995. 1402 p.

5. 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.

6. 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.

7. Vasicek O.A. An equilibrium characterization of the term structure // J. of Financial Economics. 1977. V. 5. P. 177-188.

8. Медведев Г.А. Плотности вероятностей процессов процентных ставок доходности // Вестник Томского государственного

университета. Управление, вычислительная техника и информатика. 2016. № 3 (36). C. 35-48.

9. Hull J.C. Options, futures, and other derivatives. Boston : Prentice Hall, 1998. 841 p.

10. Медведев Г.А. Стохастические процессы финансовой математики. Минск : Белорус. гос. ун-т, 2005. 243 c.

Медведев Геннадий Алексеевич, д-р физ.-мат. наук, профессор. E-mail: MedvedevGA@bsu.by Белорусский государственный университет (г. Минск, Беларусь).

Поступила в редакцию 24 ноября 2017 г.

Medvedev Gennady A. (Belarusian State University, Minsk, Republic Belarus). On a family of non-affine models of the term structure of yield.

Keywords: equation of term structure of yield; price of zero-coupon bond; CIR(1980) model; Ahn-Gao model; yield curve; forward curve.

DOI: 10.17223/19988605/42/6

The equation of term structure for the price of a zero-coupon bond is considered, the solution of which in analytical form is known, basically, for the simplest models and has an affine structure with respect to the short-term rate. In this paper, we construct solutions of this equation for a family of term structure models that are based on short-term rate processes, in which the square of volatility is proportional to the third power of the short-term rate in stochastic differential equations. The solution of the equation is sought in the form of a definite functional series and, as a result, is reduced to a confluent hypergeometric function. Three versions of the underlying stochastic differential equations for short-term rate processes are considered: with zero drift, linear drift, and quadratic drift. Numerical examples are given for the yield curve and the forward rate curve for these versions. Some conditions for the existence of nontrivial solutions of the equation of term structure in the family of processes under consideration are formulated.

Unfortunately, models that admit such solutions are few and, in particular, include some well-known models: the CIR (1980) model and the Ahn-Gao model. The requirements for the structure of the short-term interest rate model, which would allow the receipt of a term structure of the bond price in the form considered in the article, are reduced to the following.

1. To obtain a non-trivial solution, it is necessary that the degrees of polynomials determining the drift and volatility of the short-term interest rate satisfied certain constraints.

2. Another necessary condition is connected with the fact that the functional series is power-law with respect to a function that does not depend on the index of summation of the series.

3. In addition, it is necessary that the coefficients of the functional series do not depend on the maturity of the bond. Simultaneous fulfillment of these necessary conditions significantly narrows the family of models for which the solution of the

time-structure equation has the form considered in the article.

REFERENCES

1. Merton, R. (1990) Continuous-Time Finance. Oxford: Blackwell Publ.

2. Duffee, G.R. (2002) Term Premia and Interest Rate Forecasts in Affine Models. Journal of Finance. 57(1). pp. 405-443. DOI:

10.1111/1540-6261.00426

3. Medvedev, G.A. (2017) Polynomial models of yield term structure. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie,

vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(39). pp. 39-48. (In Russian). DOI: 10.17223/19988605/39/6

4. Fabozzi, F. &, Fabozzi, T. (eds) (1995) The Handbooks of Fixed Income Securities. New York: IRWIN Prof.

5. Cox, J.C., Ingersoll, J.E. & Ross, S.A. (1980) An analysis of variable rate loan contracts. Journal of Finance. 35. pp. 389-403. DOI:

10.1111/j.1540-6261.1980.tb02169.x

6. Ahn, D.-H. & Gao, B. (1999) A parametric nonlinear model of term structure dynamics. The Review of Financial Studies. 12(4). pp.

721-762. DOI: 10.1093/rfs/12.4.721

7. Vasicek, O.A. (1977) An equilibrium characterization of the term structure. Journal of Financial Economics. 5. pp. 177-188. DOI:

10.1016/0304-405X(77)90016-2

8. Medvedev, G.A. (2016) The probability density of the processes of yield interest rates. Vestnik Tomskogo gosudarstvennogo uni-

versiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 3(36). pp. 35-48. (In Russian). DOI: 10.17223/19988605/36/4

9. Hull, J.C. (1998) Options, futures, and other derivatives. Boston: Prentice Hall.

10. Medvedev, G.A. (2005) Stokhasticheskie protsessy finansovoy matematiki [Stochastic Processes of Financial Mathematics]. Minsk: BSU.