УДК 336.76
ПОИСК ОПТИМАЛЬНОГО КРЕДИТНОГО РЫЧАГА 0 УСЛОВИЯХ МАКСИМИЗАЦИИ ОЖИДАЕМОЙ СКОРОСТИ РОСТА СТОИМОСТИ ПОРТФЕЛЯ
О.И. Кривошеев
Рассмотрена задача управления портфелем на базе портфеля ценных бумаг Марковица, испытывающим независимые случайные приращения логарифма стоимости, описываемые винеровским процессом с трендом и выступающей в качестве безрискового актива наличности. Получены выражение для оптимального кредитного рычага, характеризующего структуру портфеля, оптимальное время или порог коррекции рычага в условиях ненулевых издержек операций с активами. Найдены равновесные доходность актива и рычаг в условиях выравнивания доходностей собственного капитала на всех рынках. Описанный подход адаптирован для предприятий реальной экономики, оптимизирующих свое финансовое состояние. Выведены формулы, описывающие отток капитала.
Ключевые слова: САРМ, винеровский процесс, инвестирование, кредитный рычаг, портфельное инвестирование, акции, случайные блуждания, доходность.
введение
Для моделирования инвестиционного процесса крайне важно описать поведение инвесторов по управлению своим финансовым состоянием. В настоящей статье рассматривается класс стратегий [1, 2] инвестирования на бирже, состоящих в поддержании фиксированного соотношения собственных и инвестированных средств (т. е. фиксированного кредитного рычага). Можно также говорить о доле долга в располагаемых средствах, которая взаимнооднозначно связана с кредитным рычагом. Целевым функционалом избирается ожидаемая средняя скорость прироста собственного капитала. Пример такого подхода к инвестированию в индекс Standard & Poor 500 за 2004—2014 гг. (рис. 1) приводит к доходности как функции кредитного рычага, изображенной на рис. 2: оптимум в 7 % годовых достигнут при кредитном рычаге 1,8, т. е. близком к 2.
Видно, что в результате вариации уровня кредитного рычага получается кривая доходности, похожая на параболу (далее будет продемонстрировано, что в некоторых естественных предположениях это должна быть именно парабола). Важно, что оптимальный рычаг или оптимальное финансовое состояние предприятия в данном подходе никак не зависит от рисковых предпочтений инвестора. Более традиционный подход состоит в
формировании оптимального портфеля на основе рассмотрения одновременно множества допустимых альтернатив и предпочтений инвестора, задаваемых, например, функцией полезности. Ряд
k, S&P 500
1800,00
j^
1400,00 IV"v4.Ji .л/
1000,00
600,00 ту
200,00
2004-01-14 2006-10-10 2009-07-06 2012-04-01 t
Рис. 1. Изменение индекса Standard & Poor 500 за 2004—2014 гг.
Рис. 2. Зависимость доходности на собственный капитал от кредитного рычага I, который поддерживался постоянным (с учетом изменения котировок) на каждом шаге ряда
Рис. 3. Зонт Марковица и линия САРМ: углы «зонта» соответствуют рыночным активам, сам зонт — их положительным комбинациям, линия САРМ — Парето-граница комбинаций рыночных активов, смешанных с безрисковым активом — деньгами
положений этого подхода был сформулирован в статье Г.М. Марковица (1952), где впервые предложена модель, ныне известная под его именем. Множество доступных инвестору альтернатив описывалось хорошо известным «зонтом Марковица» (рис. 3), представляющим собой криволинейную область в координатах «волатильность — доходность», своеобразными «шипами» излома границы которой являются отдельные биржевые активы. Внутренние точки этого зонта несколько нетривиально получаются как линейные (выпуклые) комбинации чистых активов, при которых доходности обычным образом усредняются с положительными весами, а волатильности или среднеквадратичные отклонения ст складываются с теми же весами на основе сложения нормально распределенных случайных величин с неединичной матрицей корреляции. Парето-оптимальные портфели образуют «юго-восточную» границу этого множества.
В дальнейшем первоначальная модель была доработана Д. Трейнером и У. Шарпом путем введения важного дополнительного безрискового актива с доходностью i = b и нулевым отклонением ст = 0. В роли актива выступали американские гособлигации (деньги). Модель получила название CAPM — Capital Asset Pricing Model [2—11], а множество Парето-оптимальных альтернатив «развернулось» в прямую касательную линию, проходящую через одну — крайнюю — точку «зонта Мар-ковица» и точку (b, 0) — безрискового актива; ее также называют линией САРМ. Точкой касания является так называемый портфель Марковица. И теперь можно считать, что любой Парето-оп-тимальный портфель состоит из гособлигаций (т. е. попросту денег) и этого портфеля Марко-вица. В большинстве работ и учебников также часто предполагают, что портфель Марковица —
это просто вся совокупность торгуемых на бирже бумаг, взятых с их весом, именуемая «индекс».
В настоящей статье предлагается предельно рафинированный подход — мы используем широко популяризированную Блейком и Шоулзом биномиальную модель ценовых блужданий (аппроксимирующую винеровский процесс) для анализа поведения инвестора, максимизирующего математическое ожидание скорости роста собственного капитала (существенно, что нами усреднению подвергается не сам доход, а фактически его нелинейное преобразование — IRR. средний IRR не равен IRR усредненного денежного потока).
Любой исследуемый портфель лежит на линии САРМ и состоит из денег и портфеля Марковица, взятых с весами, в сумме равными единице. Нулевой кредитный рычаг — когда все вложено в гособлигации/деньги, единичный — когда все вложено в портфель Марковица. Портфели между этими точками имеют положительный рычаг, меньший единицы, портфели на линии САРМ выше портфеля Марковица имеют рычаг, больший единицы (часть денег взята в кредит).
Отличие предлагаемого в данной статье подхода от классической модели Марковица в том, что рассматривается «естественная» целевая функция поведения инвестора, что позволяет сделать более конкретные выводы и вывести определенные лаконичные формулы, полезные для математического моделирования, в частности, формулы доходности и оптимального кредитного рычага. Удается показать, что эти формулы пригодны не только для управления биржевыми активами, но и для управления финансовым состоянием реальных предприятий, что позволяет в дальнейшем (уже вне рамок данного изложения) моделировать динамику оптимального инвестиционного поведения как поддержание кредитного рычага (а также как релаксационное приведение структуры активов к найденному оптимуму) и органично вводить учет инвестиций в моделях экономики типа моделей общего равновесия, добавив к ним новый экзогенный параметр или, в зависимости от модификации, внутреннюю переменную — волатильность.
Структура изложения. Сначала выводится формула оптимального кредитного рычага — оптимального уровня долга. Далее вводятся транзакци-онные издержки. Это позволяет получить оптимальное время коррекции кредитного рычага, при котором находится компромисс между потерями от отклонения от оптимума уровня долга (кредитного рычага) и потерями от потенциально бесконечных на фрактально-бесконечной броуновской траектории транзакционных издержек.
Затем рассматриваются 100 %-но коррелированные рынки с разной волатильностью и перетоком капиталов и выводятся формулы для равно-
весного уровня доходности на этих рынках — как уровня, при котором доходность собственного капитала сравнивается со среднемировой.
Тот же подход при известной скорости роста каждого локального рынка позволяет вычислить скорость оттока капитала, требуемую для поддержания местной (локальной) доходности физического капитала и кредитного рычага на равновесном уровне.
Потом обсуждается возможность переноса результатов для финансовых рынков на реальные производственные активы (математически это переход от интеграла к производной). Наконец, выводится несколько вариантов оптимального инвестиционного правила, позволяющего рационально моделировать инвестиционный процесс.
1. дискретная последовательность цен
Принятые обозначения:
¡1 — средняя (ожидаемая) доходность физического капитала актива;
¡с — средняя (ожидаемая) доходность собственного капитала;
¡с(1)/1 — удельная рентабельность;
К — основной капитал;
1
I =
— кредитный рычаг (денежный
1 - 9
мультипликатор);
М — математическое ожидание;
г — риск;
V — приращение (случайное) собственного капитала за период;
8 — конечное отклонение портфеля за время т, в гауссовом приближении 8 = Ттф;
А = \1 — 1*\/1* — отклонение рычага от оптимума;
9 — доля заемных средств — уровень долга;
— долгосрочный темп роста экономики; т — размер шага дискретной аппроксимации стохастического процесса по времени;
т* — оптимальное время пересмотра портфеля в условиях издержек;
ф — уровень генерации дисперсии изменения логарифма цен актива в винеровском процессе, приближающем их блуждание — квадрат волатильности;
* — индекс оптимального показателя; „п — индекс равновесного показателя;
х означает случайную реализацию величины х, например, используются /.
Малый инвестор на рынке последовательно видит цены актива рт = (р0, р1, ..., рп - 1, рп). Можно,
но не обязательно, считать, что эти значения равно отстоят на временной оси за суммарное время Т = пт, где т — интервал между последовательными значениями, п — максимальный индекс последовательности. При стратегии постоянного кредитного рычага действия инвестора таковы. На каждом дискретном шаге ? весь собственный капитал С t плюс С(1 — 1) кредитных средств помещается в актив, актив испытывает непрогнозируемое приращение — изменение цены вверх или вниз, далее все продается (долг погашается), и полученный новый капитал С { + 1 претерпевает такое же вложение. На практике это не оптимально: при ненулевых издержках взаимно аннулирующиеся транзакции купли-продажи естественно заменить одной транзакцией коррекции рычага: покупкой, если собственный капитал возрос, и продажей в противном случае. Проследим за ростом собственного капитала.
При нулевой ставке процента1 случайный коэффициент его роста:
V = р + 1/Р - (1 - 1) или V = 1((р, + 1/р) - 1) + 1.
Целевая функция инвестора состоит в максимизации доходности собственного капитала:
/с(/) = 1п(^(1 Ж2(0 ... Vn(l ))/Т ^ тах.
Считается, что У() > 0. Если последнее не выполнено даже в одном случае происходит банкротство инвестора, в связи с чем формально положим его доходность ¡с ^ -да. Величины ¡с и Vk являются функциями кредитного рычага, они определены и конечны для каждой положительной т
последовательности р = (р0, р1, ..., рп - 1, рп), хотя, в силу условия У к, Vk(l) > 0, не всюду.
Если итоговая цена выше начальной, а сам ряд немонотонен (см. рис. 1 и 4), легко показать, что функция ¡с(1) в области конечных значений унимодальна (см. рис. 2).
Возьмем минимальное отношение последовательных цен р { + х/р{ < 1 компонент ¡с
^ = 1п[1(р + 1/р - 1) + 1], как и ¡с в целом, убывает до -да при
I ^
1 - р t + 1/Р/
1 При ненулевой ставке процента следует писать несколько сложнее: V = - (I - 1)(1 + Ь).
о 1п р, к о о О о ° о О о о Шаг, £ К Р, о о о о о Шаг, г
а о о о о б
Рис. 4. Произвольная немонотонная последовательность прологарифмированных цен активов (а) и их приращений (б)
производная /с(/) при / = 0 в значительной мере определяется итоговым приращением логарифма цены актива 1п(рп/р0), особенно если мы рассматриваем все более мелкие разбиения с |рк + 1 — рк| П рк, как это имеет место в пределе для непрерывной зависимости р(/). Верно утверждение, что, если 1п(рп/р0) > 0 и существует t: (р{ + 1/р{ < 1), то имеется
конечный оптимальный кредитный рычаг /*, максимизирующий доходность стратегии постоянного рычага для данного ряда. Примером такой функции доходности служит изображенная на рис. 2 доходность /с(/). (Ряд на рис. 1, для которого она посчитана, не монотонен, а конечная цена выше начальной.)
На практике ряд заранее не известен — есть лишь предположения о вероятностях траекторий, поэтому целевая функция должна записываться как математическое ожидание доходности /с(/) по распределению рядов
1С = [И1п(У1У2У3
Ю]/(пт),
(1)
где V — случайный коэффициент роста собственного капитала, который зависит как от случайных событий, так и от размера рычага /, т — шаг по времени, п — число шагов.
2. непрерывный ценовой ряд
Если мы имеем непрерывную функцию, например, такую, как изображена на рис. 5, то можно рассмотреть ряд дискретных последовательностей
Р = (р0, Р1, ..., Рп - 1, Рп) с увеличивающимся п и уменьшающимся расстоянием т между последовательными значениями рк и рк + 1. Вызывает интерес существование предела при мелкости таких разбиений, стремящейся к нулю: т ^ 0.
Легко показать, что для рядов, описываемых гладкими функциями, такая предельная функция
I (/) будет прямолинейной, с производной 1п— /Т, с р 0 где Т — длина ряда; т. е. практически полезной содержательной теории построить нельзя, так как максимум либо в нуле, либо на бесконечности. Из практики же известно, что биржевые ряды всюду фрактальны. Наиболее часто предполагают, что ряд описывается винеровским процессом блуждания логарифма цены (броуновским движением с трендом). Для дальнейшего изложения зафиксируем это предположение.
Пусть все компоненты вектора р = (р0, р1, ..., рп - 1, рп) отстоят на одинаковое время т. В предположении, что каждое относительное увеличение цены р? + 1/р( актива распределено одинаково и не зависит от предыстории, в выражении (1) можно рассматривать только первый сомножитель:
¿с(/) = [М 1п( ^1(/ ))]/т,
(2)
Рис. 5. Логарифм приращения цены актива от времени (в шагах). Винеровский процесс с трендом (случайная траектория)
Рис. 6. Экспериментальные зависимости доходности от кредитного рычага (положительная часть) на шести последовательных интервалах винеровской траектории типа изображенной на рис. 5
Величина т1- = р ^ + 1/р^ _ 1 представляет собой фактическую доходность за один шаг т, она состоит из тренда ехр(т/1) и зависящего от промежутка т
т + 5(т)
случайного отклонения ехр(8п): ть = е от-
т/1 + 5(т)
куда V = I - 1 + 1е
В данной работе не предполагается возможность потерь кредитора или, по крайней мере, ее учета. Такая ситуация возникает на бирже, когда спекулянт рискует своим резервом: если оценка активов минус долг равна нулю, то активы распродаются, а счет закрывается.
Предваряя формальные выкладки, проведем численный эксперимент. Возьмем несколько реализаций ряда винеровских блужданий логарифма цены (с трендом) типа изображенных на рис. 5, и применим к ним — к каждой реализации по отдельности набор стратегий с поддержанием фиксированного кредитного рычага I. Результаты представлены на рис. 6: доходность отложена по вертикали, кредитный рычаг I по горизонтали.
Как видим, кривые доходности от рычага образуют зависимости, похожие на параболы. Они различаются высотой, но максимизирующий рычаг I* довольно стабилен. Попробуем теоретически рассчитать значение максимальной доходности г* и *.
3.теоретическое вычисление оптимального кредитного рычага
При броуновском блуждании нужно максимизировать по I целевую функцию:
где Ф(Х) — несмещенная функция распределения Гаусса с единичной дисперсией фт — дисперсия распределения для стоящей под дифференциалом сложной функции. Так как 1пх ^ -да,
х — +0
быстрое затухание гауссовой функции плотности распределения существенно для корректности интеграла, который берется по области, где
т/1 + 5
I - 1 + е I > е > 0 — по области без слишком больших по модулю отрицательных 8, что в линеаризованном случае сводится к условию, в рамках которого мы пренебрежем областью за пределами 2—3 ст: 8 > (е - 1)/1 + Ц, е > 0.
Исходя из допущения, что управлять портфелем можно сколь угодно быстро, подстраиваясь под изменения в рыночной ситуации, воспользуемся простой и наиболее удобной биномиальной моделью процесса, приводящего к броуновскому блужданию, использованной Блейком и Шоулзом. Основываясь на центральной предельной теореме теории вероятности, рассмотрим (порождаемый подбрасыванием симметричной «монеты») дрейф с
трендом ¡1 и случайным отклонением ±8, 8 = л/фг , за шаг длины т. При этом цена актива гарантированно возрастет на т11, а цена портфеля с рычагом I на т/^1, отклонение портфеля составит ±81. Чтобы рассчитать среднюю рентабельность собственного капитала ¡с в соответствии с формулой (3), достаточно усреднения коэффициентов приращения при каждом выпадении монеты, генерирующей на микроуровне броуновский процесс.
От процесса блуждания доходности при V- =
т/1 ± 5
= е ^ 1 + тг^ ± 8 (рис. 7, а) можно перейти к блужданиям при постоянном рычаге (рис. 7, б), где
т/1 ± 5
V = 1 + 1е = 1 + 1тг1 ± 81 (в линейном приближении) для любого I. Посчитаем доходность:
1п(1 + гст) = 11п{(1ехр(тг1 - 8) + 1 - I) х
х(1ехр(тг! - 8) + 1-1)} = 11п(1ехр(тг\ - 8) + 1 - I) +
+ 21п(1 ехр(тг1 - 8) + 1 - I)
или
1 г т и + 5
с = Ит т 11п(1 - 1 + е 1 1)йФ
т — 0 т
^ тах ,
.Тфт 1 е[0;+»)
/ = Ит 4- 1п{(1 ехр(т/, + 8) + 1 - I) х
с т — +0 т2 1
X (I ехр(тг1 - 8) + 1 - I)}.
ш
\ + ц +8
l + ii +8
l + Zij-
1 + /г, +
1 + li1 - 18
1 + Иг - 18
Рис. 7. Блуждания доходности: а — при единичном рычаге; б — при рычаге, поддерживающемся на уровне ,
Для удобства от логарифмов можно перейти к степеням, записав
= lim
т —^ +0
J(lexp (т/j + 8) + 1 - l)(lexp (t/j - 8) + 1 - l) — - 1
Вычислениям предварим численный расчет. На рис. 8. приведен численный предел функций
/с = -21п{(/ ехр(т/1 + 8) + 1 - /) х
х (l ехр(т/1 - 8) + 1 - l)}
и
'c =
J(lexp(т+ 8) + 1 - l)(lexp(т/^ - 8) + 1 - l) — 1
при ^ = 1, ф = 1 и достаточно малом т, когда зависимость от т не прослеживается, и график параболы /с = //"1 — (/ — /)ф/2.
Рис. 8. Кривая средней доходности собственного капитала от рычага I при данных доходности и темпе генерации дисперсии ф (¿! = 1, ф = 1); I* = 1,5 — оптимальный кредитный рычаг
Теоретический расчет сводится к работе с начальными членами рядов Тейлора: раскладывая выражение под корнем до первого порядка по т, имея в виду, что первый порядок по т — это второй
порядок по 8, ибо 8 = тф, получим последовательно 1 + /т = ^(/ехр(т/^ + 8) + 1 - /)(/ехр(т^ - 8) + 1 - /) = = 7(/(1 + т/^ + 8 + 82/2) + 1 - /) х
х (l( 1 + т - 8 + 82/2) + 1 - l) + 0(т2) =
= 1 + /(т^ + фт/2))2 - /2фт + 0(т2) = = 1 + /т/1 — (/2 — /)фт + 0(т2)/2.
Окончательно
2
/ст = /т/1 — (/ — / )фт + о(т)/2.
Здесь учтено, что в винеровском процессе 82 = тф
и ех = 1 + х + х2/2 + о(х2), „ДТ* = 1 + 0,5х + о(х),
2
/ст = /т/0 — (/ — /)ф/2.
Приближенно получим зависимость доходности от кредитного рычага (см. рис. 8):
/С = l/j - (l2 - l)ф/2,
(3)
откуда максимизирующий доходность рычаг и доходность есть l * = /j/ф + 1/2, и /'С = l *2ф/2 или /* = (/j/ф + 0,5)2ф/2.
Определим некоторую величину /0 = lim (lc/l)
0 l — +0 c
как предел доходности собственного капитала, отнесенной к соответствующему ей кредитному рычагу l:
/0 = lim -С = lim
0 l — +0 l l — +0
/j -1 (l -1 )ф"
= /1 + 1 ф.
/c =
>
Практически все формулы с использованием величины ¡0 приобретают более лаконичный вид. Например, оптимальный рычаг
I* = ¡0/Ф, (4)
доходность собственного капитала в оптимуме
¡с* = ¡0 /(2ф).
Обсудим несколько неожиданный результат: при ¡1 = 0 оптимальная доходность положительна
1
¡с = (2) 2 > 0, т. е. даже в отсутствии тренда можно заработать на шумовых блужданиях логарифма цены, если поместить в актив половину собственного капитала, а другую половину в безрисковую наличность.
Рассмотрим биномиальные блуждания с отклонениями ех и е х без тренда, когда ¡1 = 0. При вычисленном выше оптимальном кредитном рычаге I * = 1/2 доходность
; _ 11„ (1 + ех 1 + е х
¡с = 2-1п (-2--2-
-х---- > 0
— действительно, оказывается положительна.
12
Отрицательное слагаемое - ^ (I - I)ф в формуле
(3) можно интерпретировать как плату за риск собственного капитала г = ф1(1 - 1)/2. С такой точки зрения риск физического капитала, определенный как гк = г/1, составит гк = ф(1 - 1)/2. Отметим также что, при I . 1 формулы с использованием как ¡1, так и ¡0 асимптотически совпадают.
4. случаи многих активов или регионов с синхронным изменением цен
Теперь рассмотрим разные активы, синхронно меняющие свою стоимость, причем амплитуда колебаний цен на разных рынках разная. В интерпретации реального сектора, это могут быть также разные регионы, в которых синхронно меняются цены, в том числе цены на экспортный товар. Рассмотренные формулы позволяют для каждого региона с его собственными параметрами фг и ¡1г (здесь г — индекс региона) вычислить оптимальный рычаг I* и соответствующую ему рентабельность собственного капитала ¡ :
Рис. 9. Долгосрочная перспектива — семейство кривых доходности собственного капитала при разных уровнях ф в предположении (после) выравнивания доходности собственного капитала: линия максимумов — прямая
Очевидно, инвесторы, решая, куда вложить свои денежные средства, смогут выбирать регион или рынок с максимальной доходностью, что приведет к выравниванию доходности в разных секторах экономики и регионах, таким образом, возникнет новая «мировая константа» ¡с. При этом пики всех кривых доходностей для разных фг образуют одну горизонтальную прямую линию (рис. 9).
При всюду одинаковой доходности собственного капитала ¡с доходность актива и кредитный рычаг должны быть
'0 eqr
= л/2фг ¡с
¡2 ¡0
1*г = 2ф или 1*г
= (т + 0,5) ф .
¡1вЧг = ¡0eqr - фг/2 = л/2^ - фг/2,
^г = V2 ¡с/ ф•.
5. частота коррекции оптимального кредитного рычага
Рассчитаем оптимальное время управления портфелем для случая ненулевых издержек изменения портфеля. Поскольку длина траектории броуновского движения между любыми двумя точками бесконечна, то при сколь угодно малом, но ненулевом уровне комиссии у попытка корректировать любые сколь угодно малые отклонения от оптимального рычага будет приводить к бесконечным издержкам. Очевидно, их корректировать надо либо по достижении определенного отклонения А1 = \1 - 1*\, либо, что почти то же самое, через определенное время т. Перепишем доходность собственного капитала как разложение по степе-
Поправка
Рис. 10. Отрицательная поправка к доходности собственного капитала, обусловленная издержками у в зависимости от порога пересмотра портфеля (в зависимости от характерного времени реакции т)
ням (l - l *), где l * — значение кредитного рычага в максимуме доходности:
ic = li - ф/2/2 = const - ф(/ - l*)2/2 =
= ico - ф(/ - l*)2/2,
где ic0 = i0l*/2.
Разовые издержки коррекции y|l - l *| надо разделить на время между коррекциями, чтобы учесть со знаком минус в формуле для доходности:
ic = ic - Ф(/ - l*)2/2 - y|l - l*|/т.
Новый кредитный рычаг — это отношение размера активов /*(л/Фг + 1) за время т к приросту
собственного капитала за этот же период 1 + Тфг l*. Значит, изменение кредитного рычага
|l - l*| * |l* - /*(ДТ + i)/(i + l*)| * * |l* - г2|7ф^ = l*|l* - 1|ТФ^,
откуда ic = ц0/2 - Ф/*2|/* - 1|2фт/2 - yl*|l* - 1|ТФМ — если акцентироваться на зависимости от управляемого нами времени т, это функционал вида
ic = C - Вт - АД/т (рис. 10), где А, В и Cне зависят от т.
Для нахождения максимума доходности продифференцируем по т полученное выражение для доходности и приравняем производную к нулю:
|2Ф , 1,
2''r ^тт3/2
d i = - ф l*2|l* - 1|2ф + 1 yl*|l* - 1| ^ 1 = о,
йт c 2
3/2
откуда т =
JL
и т*
J.
т
2/3
11* - 1 ГДф~' 11* - 12/31*2/3ф
Если оптимальное правило должно быть указано в терминах относительного отклонения, то пор-
тфель должен пересматриваться при достижении относительным отклонением от оптимума \1 - 1*\/1* величины А:
А* = |l - l*|/l* = |l* - 1|ТФ
2/3
.1/2
|l* -12/31*2/V
2/3
= Y1/3Il* - 1| Y l* 1/3
С учетом соотношений (см. § 6) фГ2 = i0l* = ц0
и выравнивания доходности leq = , i0eq = Т^оФ
2/3
т* =
^2/3l *2/3
|l* - 12/3l*2/3ф |l* - 12/3ф1/3;
I,* J2/3
А* = |l - l*|/l* = у1/31 l - 1 =
l *1/3
r-o - 1
Ф
2/3 ,1/6 l/3 Ф_ 1/6 , Vo
|l - l *| = l *2/3|l * - 1|2/3y1/3 =
- 1
2/3 1/3
/ y1/3
ф
1/3
Если пренебречь единицей при большом рычаге I * . 1, то оптимальное время пересмотра
т* = у2/3/(ф1/3 ц0/3), соответствующий порог коррекции для относительного отклонения
|М = т1/3 ф 1/6 I* ц1/6 .
г-0
6. моделирование инвестиционного поведения в реальной экономике
В дальнейшем модель может быть использована для обобщения статической модели Вальраса [1] в формате введения в последнюю инвестиций. На основе уравнения изменения логарифма
рычага — 1п1 = -¡I + ц (где ц = — 1пК, см. § 10) может быть сформулировано оптимальное инвестиционное правило ц = ¡1, если кредитный рычаг уже
в равновесии, и ц = ¡1 + — 1п1, если кредитный ры-
dt
чаг почему-либо подвергся отклонению.
2 — 1 Заменой — 1п1 ^ - (I - I*) (релаксационный dt т
член, моделирующий стремление экономических
Отметим, что в этой формуле релаксационный член можно поделить на /, получив й 1п/ = 1 (/ — /*) или / = 1 (/ — /*).
й? /т й? т
агентов привести рычаг к оптимальному) и подстановкой ц = 1пК получаем
- 1пК = // + 1 (/ — /*)
— т '
(
т — время релаксации
отождествляемое нами с
2/3
рассчитанным выше временем т*
|/* - 12/3/*2/3ф/ Этой скорости изменения капитала соответствует объем инвестиций I, пропорциональный ц + где - — скорость выбытия фондов.
Имеем (в расчете на единицу капитала):
I = ^- + // + т (/ — /*). В отсутствии перетока
капитала также может быть использовано выражение (4) и формулы
d lnK = \ I + 1 f I - ^ dt 1 т v ф.
или
d inK = /i i + 1 f i -1 -iY
dt 1 т v 2
7. время оценки изменения оптимального рычага
Ранее мы анализировали отклонение рычага / от точно известного рычага / *. Для измерения величин /1 и ф при достаточном объеме данных могли бы быть задействованы статистические тесты, однако реальные ряды не подчиняются предположениям модели винеровского процесса. Но кривая /с(/) определена для любой реализации ряда — независимо от модели, и оптимальный рычаг может быть измерен по достаточно большой предыстории, как это показано на рис. 2 для ряда 8 & Р 500, приведенного на рис. 1. На реальном рынке тенденции постоянно изменяются, как видно на том же рис. 1. Если подвергается отклонению оптимальный рычаг, все должно определяться временем, на котором можно детектировать это отклонение.
Для простоты считаем, что нам нужно детектировать изменение от начального значения /1н, которое скачком превращается в конечное значение /1к. По существу нам надо сравнить отклонение тренда (/1н — /1к)Тоткл и среднее отклонение,
обусловленное шумом л/Тотклф (которое доминирует на малых временах).
Первое из них превзойдет второе после того, как они впервые сравняются при
T = Y
откл 2 '
(г1н - г1к)
Для простоты полагая одну из этих величин равной нулю, получим время, на котором метод позволяет детектировать изменение ^ на 100 % (т. е. на
величину порядка самой себя), Тоткл = ф/ г'2 = 1/(2ic), т. е. длина ряда, по которому нашим методом может быть измерен оптимальный рычаг не может быть меньше т2 = 1/г*. Для ситуации, представленной
на рис. 2, это т2 = -1-- « 7 лет, что лишь
2 • 0,07 год-1
немного меньше 10 лет — длины ряда на рис. 1.
Это второе естественное характерное время при принятии инвестиционных решений. В формуле d lnK = г1 + 1 (I — I *) теперь либо рычаг dt т
I * = arg max ic(l), измеренный по предыстории при-
i > о
мерной длины /с , либо можно предложить идейно аналогичную формуле (2) схему со скользящим пересчетом кривой доходности
/Д) := /* lnVT(0 + z'c(1 )(1 - Т /*),
где /* = max / (/) — оптимальная скорость роста
i > о с
собственного капитала, т — шаг расчета (он может совпадать с т*), V.(1) — приращение собственного капитала за интервал длины т, т /С — параметр скользящего среднего, отвечающий характерному времени равному /с .
Как показано в § 9 параметр ф также можно интерпретировать через отклонение рентабельности
реального актива ф = т0А/ (рис. 11), где А/ — отклонение рентабельности реального актива, т0 — характерное время отклонения (его можно назвать временем корреляции).
Рис. 11. Рентабельность в зависимости от порога реакции (в зависимости от характерного времени реакции)
8. расчет перетока капитала
При выравнивании доходности ввиду возможности перетока в среднем d lnK = 9, где 9 — сред-
dt
няя скорость роста рыночной ниши данной отрасли. Поскольку не весь капитал инвестируется в данном секторе, то появляется уравнение, описывающее межсекторный, межрегиональный и меж-страновый переток инвестиций. Из условия 2гс = ц0, где ц0 — темп роста глобальной экономики, оптимальный рычаг определяется как l*q = /ф .
Доходность на капитал г = Ф (точная формула i1eq = - ф/2, см. § 10, формула 16).
Тогда усредненный отток капитала с единицы физического капитала K измеряется как разность средней доходности и скорости роста: cK = г* — 9 s
= У^ОФ — 9. (Отрицательный знак этой величины
будет означать приток капитала от внешних инвесторов в экономику.)
Если умножим его на стоимость активов pKK, то получим полный отток капитала:
C = PkKck = (г* — 9)PkK s (Т^Ф — 9)PkK.
Отметим, содержательное заключение состоит в том, что для борьбы с оттоком капитала правительство должно проводить стабилизационную политику — снижать параметр ф = т0 Аг2 или обеспечивать экономический рост, в частности, посредством защиты рынка.
В более точном варианте отток капитала
C = PkKCk = (^ — Ф/2 — 9)РкK.
9. применение к реальной экономике
Ключом перехода к реальной экономике служит осознание факта, что рентабельности биржевого и реального активов соотносятся как интеграл и производная. Воспользуемся ступенчатой моделью ряда рентабельности. С равной вероятностью 0,5 IRR принимает значение г^ + Аг или г'1 — Аг (см. рис. 11).
Рассмотрим завод без долга. Реинвестируя всю прибыль, создадим фиктивный «биржевой» актив, растущий как K(t) = K(t0)exp[(t — t0)/j при посто-
t
янной рентабельности г и как K(t) = K(t0)exp J/(t)dt
при переменной При случайных и независимых пиках высоты /^ + Аг или /^ — Аг и длины т0 рост (логарифма) цен этого актива описывается броуновским процессом с ф = т0А/2, после чего могут быть применены все выведенные формулы.
10. инвестиционное правило и долгосрочный темп роста экономики
Рассмотрим размер удельного долга 9 = Б/(ркК), где Б — долг, в зависимости от номинальных или
реальных инвестиций ц = 1 ^: — 9 = — г1 + (1 — 9)ц
К ш dt
или — 1п/ = — г'1/ + ц. Считается что цена фондов
рк неизменна.
Вывод требует расчета производной удельного долга — Б и некоторых преобразований, в рк Л К
основе которых лежат соотношения 4 Б = I — F,
Л
F = (г + d)pKK, d K = — — dK. Поскольку l =
1 - e:
то — 1п/ = —!— 9, где — — скорость выбытия, I — Л 1 - 9 Л
инвестиции.
Если ц0 — трендовая скорость роста экономики, то в равновесии 0 = —/'1/ + ц0, где также /^ — средняя рентабельность, получим ^Г = ц0 и ц0 = 2/с,
откуда = Тц07ф , /0ег = ТЦ0ф . С учетом соотношения /'1 = /0 — ф/2 получим
= — ф/2.
(5)
Впрочем, для удобства может оказаться целесообразно пользоваться приближенной формулой
г1е® ~ г0е®
С = л/^Ф .
(6)
Для моделирования реальной экономики бывает очень полезно применять поддерживающее фиксированный рычаг инвестиционное правило
ц(^ = /'(^/* или, что то же,
dlnK dt
= /(t)l *.
Оно задает мгновенный инвестиционный спрос вашей модели, который начнет явно зависеть от вектора цен после представления текущей рентабельности функцией этого вектора.
о
заключение и выводы
Для наиболее примитивного, но универсального предположения о характере изменения (логарифма) цен на рынке (винеровский процесс с трендом) рассчитан оптимальный в смысле максимизации скорости роста (рентабельности) алгоритм управления портфелем. Рассчитан оптимальный рычаг I * = ¡0/ф и соответствующая ему доходность ¡* = ¡2 /(2ф).
В долгосрочном равновесии после выравнивания доходности с учетом равенств ¡1 = ц0 = 2¡c
(см. § 10) получили I = >0/ф , ¡ = л/цоФ .
В условиях транзакционных издержек на коррекцию рычага при «сверхдлинном рычаге» коррекцию оптимально осуществлять через время
т = у2/3/(ф1/3ц02/3) или т = у2/3/(/*4/3ф)
или при переходе порога относительного отклонения
1/3,1/6
1М = У^Ь. или А = \/ - 1*\/1* = у1/3/*1/3.
1 ц0/
Эти формулы (ассимптотические, I* ^ да) оптимального порога отклонения (времени коррекции) очевидно не логичны при I * = 1 — потому при малых рычагах в диапазоне 0 < I * < 2, может быть и чуть выше, надлежит пользоваться более точными формулами:
J.
2/3
у2/3 ¡ *2/3
А = \i - i*\/i * = у1/3
\i* -12/3г2/3ф \i* -12/3ц2/3ф1/3 |/* -12/3 =
i *1/3
-1
2/3 Y1/3 ф16 ' 1/6'
\i - i*\ = i*2/3\i* - 1\2/3у1/3 = |7^7ф - 1\2/3у1/3 Ц
1/3
ф
1/3
Чтобы применить вторые варианты формул, следовало бы проверить предположения I* = ^/ц0/ф,
¡* = >0ф , невыполнение которых может указывать на неадекватность модели, пере- или недооценен-ность активов, неравновесие, в том числе надувание ценового пузыря. Также надлежит проконтролировать и винеровский характер приращения (логарифмов) цен активов.
При моделировании инвестиционного поведения в реальной экономике могут понадобиться инвестиционные правило, сохраняющее кредитный
рычаг вида — 1пК = ¡(йГ, и его обобщение на случай dt
отклонения от оптимума: d lnK = i(t)l* + - (l — l*),
dt т
где т = т*. Также следует иметь в виду, что для измерения трендового показателя ix или оптимального рычага l путем максимизации доходности ic(l) := V(l)/T + i'c(l)(l — т/T) с погрешностью менее 100 % потребуется время T > 1/i*.
Замечание по методологии исследования. Существенный момент заключается в том, что операция усреднения применялась не к сумме дохода с последующим расчетом IRR — рентабельности, а к распределению уже рассчитанных значений IRR. Для броуновского движения это означает максимизацию медианного дохода, который (в этой постановке) совпал со средней доходностью. Перестановка расчета IRR и операции усреднения (когда усреднение сперва применяется к доходу) приведет к радикально новым выводам.
литература
1. Кривошеев О.И. Построение обобщенной модели Вальраса с инвестициями на базе оптимизации рациональными инвесторами кредитного рычага // Тез. докл. седьмой между-нар. конф. «Управление развитием крупномасштабных систем» MLSD'2013. - М., 2013. - С 195-198.
2. Шарп У.Ф., Александер Г.Дж., Бейли Дж. — Инвестиции: пер. с англ. — М.: Инфра-М, 1997.
3. Шапкин А.С. Экономические и финансовые риски. Оценка, управление, портфель инвестиций. — М.: Юнити, 2003.
4. Black Fischer, Jensen M.C., and Scholes M. The Capital Asset Pricing Model: Some Empirical Tests, P. 79—121 in M. Jensen ed. — URL: http://www.efalken.com/LowVolClassics/black-jensenscholes.pdf (дата обращения 9.07.2015).
5. French C. W. The Treynor Capital Asset Pricing Model // Journal of Investment Management. — 2003. — Vol. 1, N 2. — P. 60—72.
6. Markowitz, H.M. The early history of portfolio theory: 1600—1960 // Financial Analysts Journal. — 1999. — Vol. 55, N 4.
7. Sharpe W.F. Capital asset prices: A theory of market equilibrium under conditions of risk // Journal of Finance. — 1964. — Vol. XIX, N 3. — P. 425—442.
8. Tobin J. Liquidity preference as behavior towards risk. — URL: http://dido.wss.yale.edu/P/cm/m19/m19-01.pdf (дата обращения 9.07.2015).
9. Treynor J.L. Toward a Theory of Market Value of Risky Assets. Unpublished manuscript. A final version was published in 1999, in Asset Pricing and Portfolio Performance: Models, Strategy and Performance Metrics. R.A. Korajczyk (ed.). — London: Risk Books. — P. 15—22.
10. Triana P. Lecturing Birds on Flying: Can Mathematical Theories Destroy the Financial Markets? — N.-Y.: Wiley, 2009. — 400 c.
11. Altman E. Corporate Financial Distress and Bankruptcy: 3rd edition. — N.-Y.: John Wiley and Sons, 2005.
Статья представлена к публикации членом редколлегии
В.В. Клочковым.
Кривошеев Олег Игоревич — преподаватель,
Московский государственный университет экономики
статистики и информатики; науч. сотрудник,
Институт проблем управления им. В.А. Трапезникова РАН,
Н [email protected], [email protected].
=
т
ф