ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2008 Управление, вычислительная техника и информатика № 1(2)
УДК 681.513.3
В.В. Домбровский, Д.В. Домбровский, Е.А. Ляшенко ДИНАМИЧЕСКАЯ ОПТИМИЗАЦИЯ ИНВЕСТИЦИОННОГО ПОРТФЕЛЯ ПРИ ОГРАНИЧЕНИЯХ НА ОБЪЕМЫ ВЛОЖЕНИЙ В ФИНАНСОВЫЕ АКТИВЫ
В работе представлена модель управления инвестиционным портфелем (ИП), учитывающая ограничения на объемы торговых операций. Цены рисковых финансовых активов подчиняются стохастическим разностным уравнениям со случайной волатильностью. Целью управления ИП является отслеживание гипотетического эталонного портфеля с заданной траекторией роста. Получены уравнения, определяющие оптимальные стратегии прогнозирующего управления ИП с обратной связью при ограничениях.
Ключевые слова: инвестиционный портфель, прогнозирующие стратегии управления, оптимизация при ограничениях.
Проблема выбора оптимальной структуры портфеля ценных бумаг (инвестиционного портфеля - ИП) - одна из основных в финансовом менеджменте. Существуют различные подходы к управлению ИП, которые отличаются способами описания эволюции цен финансовых активов, представления динамики ИП, целями управления (выбором функции риска) [1 - 6].
В работе рассматривается задача управления инвестиционным портфелем, состоящим из рисковых вложений (обыкновенных акций) и безрискового вклада (банковского счета), допускается также возможность займа по безрисковой ставке и участие в операциях продажи без покрытия. Предполагается, что цены рисковых финансовых активов описываются дискретными аналогами уравнений диффузионного типа со стохастической волатильностью. Целью управления является отслеживание желаемой траектории роста капитала, которая задается инвестором [5 - 6].
Для оптимизации инвестиционного портфеля предлагается использовать методологию управления с прогнозирующей моделью [6]. Привлекательной чертой такого подхода является возможность достаточно просто в явном виде учитывать ограничения на переменные состояния и управления. При этом получается стратегия управления с обратной связью, но удается избежать так называемого «проклятия размерности», которое препятствует синтезу управлений с обратной связью при ограничениях, если применять традиционные подходы с использованием метода динамического программирования. В случае оптимизации инвестиционного портфеля метод позволяет учесть ограничения на объемы торговых операций, например на размеры займов и операций продаж без покрытия.
1. Модель инвестиционного портфеля
Рассмотрим ИП, состоящий из п видов рисковых финансовых активов (обыкновенных акций) и одного безрискового финансового актива (банковского счета или надежных облигаций). Капитал, помещенный в рисковый актив г-го вида в
момент времени к, равен и,(к), г = 1, п, в безрисковый - и0(к) > 0, объем безриско-
вого займа равен ии+і(к) > 0. Тогда общий объем вложений (капитал портфеля) в момент к будет равен
П
V(к) = Х щ (к) + и0 (к)-ия+1 (к). (1)
І=1
Отметим, что если М;(к) < 0, і = 1, п, то это означает участие в операции «продажа без покрытия» [1] на сумму |иг(к)|. В момент времени к + 1 капитал портфеля станет равен (рассматривается самофинансируемый портфель)
П
V(к +1)= £ [1 + Пі (к + 1)]мг (к) + [1 + гх ]м0(к) -[1 + г2]мя+1 (к), (2)
І=1
где Пг(к + 1) - ставка доходности рисковых вложений на интервале [к; к+1] (случайная величина), т\ - неслучайная доходность безрисковых вложений, г2 - не-
случайная ставка по займу безрискового актива (т\ < г2).
Используя (1), уравнение (2) преобразуем к виду
П
V(к +1) = [1 + г]Г (к) +Е [п- (к)- г\щ (к) -[ -Г К+(к),
1=1
П
при этом и0 (к) = V(к) - X щ (к) + ии+1 (к).
І=1
При управлении портфелем учитываются следующие ограничения:
итіп (к) < щ (к) < итах (к), і = ~п ; (3)
0 < «я+1 (к) < <+? (к); (4)
п
0 < V(к) - х и (к) + Мп+1 (к) < «0тах (к). (5)
І=1
Если нижняя граница игтт (к) < 0 , і = 1, п, то для рискового актива і-го вида
допустимо участие в операции «продажа без покрытия» на сумму не более
|игшш (к)|; если игтт (к) > 0, і = 1, п, то операции «продажа без покрытия» для
рискового актива і-го вида запрещены; игшах (к), і = 1, п, определяют максимальный объем капитала, который можно вкладывать в рисковый актив і-го вида; м0тах (к) > 0 определяет максимальный объем капитала, который можно вкладывать в безрисковый актив, «тТ (к) > 0 определяет максимальный размер займа безрискового актива. Отметим, что величины игтт (к), і = 1, п , и и™* (к),
і = 1, п +1, на практике часто зависят от величины общего капитала ИП, что можно учесть, положив и™11 (к) = у• V(к), и“ж (к) = у" V(к), где у- и у" - постоянные коэффициенты.
Для описания эволюции цен рисковых финансовых активов используем модель вида
П
1 + цг- (к)+ X Оу [0(к),к]м>у (к)
і=і
Б, (к +1) = Б, (к)
где Б,(к) - цена г-й ценной бумаги в момент времени к (г = 1, п), р,г(к) - ожидаемая
доходность г-й ценной бумаги, ^(к) - вектор белых шумов размерности п с нулевым средним и единичной матрицей ковариации, а;у[0(к),к] - элементы матрицы волатильности ||агу[0(к),к]||, зависящие от 0(к) линейно, 0(к) - последовательность независимых д-мерных случайных векторов с известными первым и вторым моментами: М (0(к)} = 0(к), М{0(к) 0т(к)} = © 0 М{0т(к)^т(к) } = 0.
Необходимо определить стратегию управления инвестиционным портфелем путем перераспределения капитала между различными видами инвестиций так, чтобы капитал реального портфеля V (к) с наименьшими отклонениями следовал капиталу Р°(к) некоторого определяемого инвестором эталонного портфеля, эволюция которого описывается уравнением
V (к + 1) = [1 + ц° (к)] V °(к), где р°(к) - заданная желаемая доходность портфеля. В начальный момент V0 (0) = V (0) заданы.
2. Синтез прогнозирующих стратегий управления ИП
Для управления портфелем синтезируем стратегии с прогнозирующей моделью. Критерий качества управления со скользящим горизонтом (функция риска) имеет вид
!р-1г А -| 2 р—1 гг
ХГV(к + і/к)-V (к + і)] +Х и (к + і/к)Я(к + і)и(к + і/к) +
і=1 і=о
+ \у (к + р/к)-V0 (к + р)7V()}, (6)
где и(к + г/к) = [иі(к + г/к),...,и„(к + г'/к)], і = 0, р -1, - векторы прогнозирующих управлений, V(к + г'/к) - предсказанное состояние портфеля в момент времени к + г, Л(г) > 0 - матрица весовых коэффициентов управления.
Прогнозирующие управления определяются по следующему правилу: на каждом шаге к минимизируем функционал (6) по последовательности программных управлений и(к / к),...,и(к +р -1/к), зависящих от состояния системы в момент к. В качестве управления в момент к берем и(к) = и(к/к). Тем самым получаем управление и(к) как функцию состояния V (к), т. е. управление с обратной связью. Чтобы получить управления и(к + 1) на следующем шаге, процедура повторяется для текущего момента к + 1.
Теорема. Оптимальное управление на шаге к, минимизирующее критерий (6) при фиксированном горизонте прогнозирования р с учетом ограничений на управления (3) - (5), определяется уравнением
и(к) = [/„ 0И 0и]и(к/к),
где /„ - единичная матрица размерности п; 0„ - квадратная нулевая матрица размерности п; и(к/к) = [ит(к/к), ит(к+1/к), ., ит(к+р-1/к)]т - вектор прогнозирующих управлений, который определяется из решения задачи квадратичного программирования с критерием
У (к + р/к) = 2 2Т (к )0 (к )и (к) + ит (к) Н (к )и (к)
при ограничениях ишт (к) < 5(к)и(к) < итах (к),
где
<к) = \у (к), V0 (к)] Т ,
ит” (к) = [«“” (к), «г (к), 0, -V(к)]Т,
итах (к) = [мГ (к), «тг (к), «0тах (к)-V(к —,
И,, (к) ... И,р (к)
G(k), Н(к), 5 (к) - блочные матрицы:
а (к) = [^ (к) ... Ор (к)] ; И (к) =
Ир, (к) ... Ирр (к)
5(к) = [£(к) 0(я + 2)х(и+1)(р-1)] ,
блоки которых определяются следующими соотношениями:
С* (к) = {Г=2 Ат (к + у)} } } - 5);
В (к + г -1)<! П Ат (к + у)[ (р - я), г < я,
122 (Р - г) + К(к + г -1)>
НТ,
г = я,
г > я;
" 1 .. 0 0"
0 .. 1 0
0 .. 0 1
-1 .. . -1 1
я (к) =
0(я+2)х(я+1)О-1) - нулевая матрица размерности (п + 2)х(п +1)(р -1);
П=г Ат (к + ]) = 1; 2(з +1) = (з) + Ст • С ; б(0) = Ст • С ;
Ь11 (я) = Ат (к + р -1 - я)Q (я) А ( к + р -1 - я);
^2 (я) = Ат (к + р -1 - я)б(я)В0 (к + р-1 - я);
ь22 (я) = В0Т (к + Р-1- Я)£)Во (к + р -1- я) +
+м{х”=1 Б] [0(к + р -1 - 5), к + р -1 - 5] X
X2($)В у [0(к + р -1 - ^), к + р -1 - $]},
где А (к) = Ша§ (1 + тх (к), 1 + р0 (к));
Во (к) =
В у [0(к), к ] =
Ц (к) - Г (к) ... Цп (к) - г! (к) - [ (к) - г! (к)]'
0 ... 0 0
Ст1 у [0(к),к] ... апу [0(к),к] 0'
0 ... 0 0
; С = [1 -1]
ЛИТЕРАТУРА
1. Merton R.C. Continuous-time finance. Cambr.: Ma. Blackwell, 1990.
2. Bielecki T.R., Pliska S.R. Risk-sensitive dynamic asset management // Appl. Math. and Opt.
1999. No. 39. P. 337 - 360.
3. Young M.R. A Minimax portfolio selection rule with linear programming solution // Management Science. 1998. V. 44. No. 5. P. 673 - 683.
4. Browne S. Risk-constrained dynamic active portfolio management // Management Science.
2000. V. 46. No. 9. P. 1188 - 1199.
5. Домбровский В.В., Ляшенко Е.А. Линейно-квадратичное управление дискретными системами со случайными параметрами и мультипликативными шумами с применением к оптимизации инвестиционного портфеля // Автоматика и телемеханика. 2003. № 10. С. 50 - 65.
6. Домбровский В.В., Домбровский Д.В., Ляшенко Е.А. Управление с прогнозированием системами со случайными параметрами и мультипликативными шумами и применение к оптимизации инвестиционного портфеля // Автоматика и телемеханика. 2005. № 4. C. 84 - 97.
Статья представлена кафедрой математических методов и информационных технологий в экономике экономического факультета и кафедрой прикладной математики факультета прикладной математики и кибернетики Томского государственного университета, поступила в научную редакцию 8 июня 2007 г.