Рассматривается задача управления инвестиционным портфелем с учетом ограничений на объемы торговых операций. Целью управления является превышение в среднем капитала индексного портфеля. Предлагается использовать методологию управления с прогнозирующей моделью для синтеза управления с обратной связью при ограничениях.
Adaptive feedback strategies of active investmentmanagement under portfolio constraints .pdf Проблема управления инвестиционным портфелем (ИП) является одной из основных в управлении финансами и представляет как теоретический, так и практический интерес. Выделяют два основных подхода к ее решению. Классический подход, предложенный в [1, 2], и последующие его модификации [3 - 5], исходят из предположения о том, что при формировании своего портфеля инвестор, с одной стороны, хотел бы минимизировать риск портфеля (обычно дисперсию портфеля или связанные с ней меры риска), с другой - получать желаемую доходность (либо в двойственной постановке - максимизировать доходность при ограниченном риске). При этом задача оптимизации структуры портфеля (определения оптимальных долей вложений в различные виды активов) решается в статической постановке (однопериодные модели) и в зависимости от выбора функции риска и способов учета неопределенности сводится к решению задач квадратичного, стохастического или линейного программирования [1 - 5].Второй подход основан на построении динамических моделей ИП и использовании для выбора оптимальной структуры портфеля методов теории стохастического управления [6 - 17]. Классическая оптимизационная проблема в динамической постановке заключается в определении стратегии управления ИП, максимизирующей некоторую интегральную функцию полезности, зависящую от уровня текущего потребления и конечного богатства, которая имеет достаточно условный характер и обычно подбирается из класса функций, для которых можно получить приемлемое аналитическое решение [6]. За исключением весьма ограниченного набора функций полезности, такой подход приводит к трудной проблеме численного решения уравнений динамического программирования [7]. Используются также и другие критерии оптимизации портфеля, в частности в [9] критерий, чувствительный к риску. В работах [10 - 14] задача управления ИП формулируется как динамическая задача слежения за капиталом некоторого гипотетического эталонного портфеля, имеющего задаваемую инвестором желаемую доходность.В работе [14] рассматривается так называемая задача активного управления, целью которой является превышение в среднем капитала некоторого индексного (базового) портфеля. Индексный портфель состоит из различных видов акций в заданных постоянных долях, пропорциональных объемам этих активов на рынке, и отражает состояние и структуру рынка в целом [15, 16]. Достаточно полный обзор методов оптимизации ИП в динамической постановке с использованием различных критериев дан в [17].Известно, что реальные модели должны учитывать ограничения на объемы торговых операций, таких как покупка-продажа отдельных видов финансовых активов, продажа без покрытия [16], заем средств по безрисковой ставке. В динамических моделях учет подобных ограничений приводит к проблеме, известной как проклятие размерности, если решать задачу с помощью классических методов, например динамического программирования [6, 7].В данной работе для решения задачи управления ИП при ограничениях предлагается использовать метод прогнозирующего управления, получивший широкое признание и применение в практике управления сложными технологическими процессами при наличии ограничений [12, 13, 18]. Решается задача активного управления. Получены оптимальные стратегии управления с обратной связью, расчет которых включает решение последовательности задач квадратичного программирования. Представлены результаты численного моделирования, подтверждающие эффективность и работоспособность предложенного подхода.1. Модель инвестиционного портфеляРассмотрим инвестиционный портфель (ИП), состоящий из n рисковых ценных бумаг (обыкновенные акции) и безрискового финансового актива (банковский счет или надежные облигации). Динамика капитала инвестиционного портфеля V(k) описывается уравнением [13]:V (k +1) = [1 + rl (k )]V (k) + +1К (k) -1 (k)]ut (k) -[r2 (k)-п (k)K+1 (k) ,(11где r\t (k) - доходность рискового актива i-го вида, rt(k)- доходность безрискового актива, r2(k) - ставка по займу безрискового актива, и,(к) - объем капитала, вложенного в i-й вид акций, ии+1 (k) > 0 - объем заемных средств,ио (k) = V(k) - Xщ (k) + un+l (k) > 0- капитал, помещенный в безрисковый актив. Начальное значение V(0) известно. Управляющими являются переменные и,(к), i = 1, n +1, на которые наложены следующие ограничения:и™11 (k) < ui (k) < Mmax (k), i = 1,n ;(2) 0 < un+1 (k) < M„m+T (k) ;(3)0 < V(k) - it ut (k) + un+l (k) < omax (k) ,(4)i=1где ишах (k) > 0, i = 0, n +1. Если нижняя граница ишт (к) < 0, то для акций i-го вида допустимо участие в операции продажи без покрытия [16] на сумму не больше, чем |иШ1П (к)|; если ишт (к) > 0 , то операции продажи без покрытия дляакций i-го вида запрещены; ишах (к) определяют максимальный объем капитала, который можно вкладывать в акции i-го вида; и™* (к) - максимальный размерзайма безрискового актива; иЦ13* (к) - максимальный объем капитала, которыйможно вкладывать в безрисковый актив.Отметим, что величины ишт (к), i = 1,n , и и™* (k), j = 0,n +1, на практикечасто зависят от величины общего капитала ИП, что можно учесть, положиви™11 (k) = у) V (к), и™* (k) = yj V (k), где yj и у" - постоянные коэффициенты.Для описания эволюции цен рисковых финансов^гх активов используем модель вида1 + Иг. (к) + aj (k) Wj (k)| ,(5)St (k +1) = St (k)где Si (k) > 0 - цена ценной бумаги i-го вида в момент времени k, i = 1, n, цг- (k) -ожидаемая доходность на интервале k +1, агу (k) - элементы матрицы волатиль-ности ||стгу (k)|| размерности n х n, Wj (k) - белый шум с нулевым средним и единичной дисперсией (j = 1, n). Модель (5) является дискретной аппроксимацией (сточностью до второго порядка малости) непрерывной модели цен, так называемой модели геометрического броуновского движения [19]. Доходность рискового актива i-го вида для модели цен (5) имеет видП (к) = Ц, (к) + £ cjj (k) wj (k).2. Управление инвестиционным портфелемРассмотрим задачу активного управления. Предположим, что задан некоторый базовый портфель, который содержит набор финансовых активов в заданных долях в соответствии со структурой и состоянием рынка в целом [15, 16]. При пассивном управлении [16] инвестор формирует ИП, распределяя капитал в тех же пропорциях. Активное управление заключается в том, чтобы превосходить доходность индексного портфеля путем перераспределения инвестиций.Динамика капитала индексного портфеля V(k) описывается следующимуравнением [14]:V (к + 1) =1 + X с.1ц, (к) + X CTj(k ) Wj (к )i=1 V j=11" Z ai I r1 (k)i=1V? (k) ,(6)где аг- > 0 , i = 1, n, - доля капитала ИП, вложенная в i-й рисковый актив,n1 - Z ai > 0 - доля капитала ИП, вложенная в безрисковый актив. В начальныймомент капитал V(0) = V (0) задан.В общем случае индексный и инвестиционный портфели могут содержать различные виды акций. Если i-й вид акций включается в индексный портфель, но невключается в ИП, то ишт (k) = ишах (k) = 0 . Если i-й вид акций включен в ИП, ноне включен в индексный портфель, то аг- = 0 .Для управления портфелем синтезируем стратегии с прогнозирующей моделью. Критерий качества управления со скользящим горизонтом (функция риска) имеет видJ(k + p/k) = EJ £_V(k + i)-(1 + ß)V(k + i)]2 ++ uJ (k + i/k)R(k + i)u(k + ilk)jV(k), V(k)) ,(7 где E{..y...} - оператор условного математического ожидания,u(k + i/ k) = [щ(к + i/ k), .„ un(k + i/ k)] (i = 0,p-1) - векторы прогнозирующих управлений, p - горизонт прогноза, R (k + i) > 0 - весовые матрицы, параметр ß> 0 характеризует склонность инвестора к риску.Прогнозирующие управления определяются по следующему правилу: на каждом шаге k минимизируем функционал (7) по последовательности программных управлений u(k/k),...,u(k+_>-1/k), зависящих от состояния системы в момент k. В качестве управления в момент k берем u(k) = u(k/k). Тем самым получаем управление u(k) как функцию состояния V(k), т. е. управление с обратной связью. Чтобы получить управления u(k+1) на следующем шаге, процедура повторяется для текущего момента k+1.Теорема. Оптимальное управление на шаге k, минимизирующее критерий (7) на траекториях системы (1), (6) при фиксированном горизонте прогнозирования p с учетом ограничений на управления (2) - (4), определяется уравнениеми(к) = [/я+1,оя+1, ...,оя+1 ]и(k),где 1„ - единичная матрица размерности n, 0n - квадратная нулевая матрица размерности n, U(к) = [^ит (к/к), uт (к +1/k),..., uт (k + p-I/к)J - вектор прогнозирующих управлений, который определяется из решения задачи квадратичного программирования с критериемY (к + p/k) = 2 zT (к )G (к )U (к) + UT (к) H (к )U (к) ,(8) при ограниченияхumin (k) блоки которых определяются следующими соотношениями: Gs (k) = {ПГ=0 A (k + j)} А2 (i> - s);(10) (11)(12)(13)(14)BqT (k + i- 1)Щ A (k + уп42 (/?-s), i < s,His (k) 4 (p - i) + R(k + i -1),i = s, i > s;гдеП£ A (k + j ) = 1,Q(s +1) = Lu (s) + CT C , ß(0) = CT C ; Li (s) = £ Aj (k + _> -1 -s)Q(s) Aj (k + _> -1 -s) ;Li(s) = £ Aj (k + _>-1 -s)ß(s)Bj (k + p-1 -s);L22 (s) = £ BT (k + _> -1 -s)Q(s)Bj (k + _> -1 -s),j=0(k) = diag| 1 + ri (k), 1 + ri (k) + 2ai ((i (k)- i (k)) j;i=1(15) (16)(17)(18)(19)Ц1 (к) - ГХ (k) ... Ц„ (к) - r! (к) Г! (к) - Г2 (к)_0 ... 0 0Aj (k) = diag | 0, 2 aiCTj (k) j;j (k) ... a „j (k) 0 _ 0 ... 0 0 C = [1, -(1 + ß)].Доказательство. Введя вектор z(k + i) = [V(k + i), F(k + i)] T, можно переписать систему уравнений (1), (6) в следующем виде:z (k + i +1) = А0 (k) z (k + i) + B0 (k) u (k + i/h) +-£ A (k) wj (k)z(k + i) + £ B (k) wj (k)и(k + i/k),(23)с начальным условием z(k) = [V(k), F(k)]T, где A0(k), B0(k), A7(k), B7(k), j = 1,определяются уравнениями (19) - (22) соответственно. Критерий (7) перепишется следующим образом:J (k + p/k) = Е { zт (к + i) C TCz (к + i) +ni=l+ 2 и1 (k + i/k)R(k + i)u(k + i/k )/V(k),V(k)(24)i=0Выражая последовательно все z (k + i), i = 1, p, через z(k) с использованием уравнения системы (23) и подставляя результат в критерий (24), получимJ (к + p/k) = zT (к) Lj (J? -1) z (k) ++2zT(k)I (Ii A (k + s )\ Ll2 (p - r - l)u (k + r/k)r=1 ls=0+2Р1 I uT(k + r/k)BT(k + r)( ft Ao (k + s)^L2(p-1-j)u(k + jlk)-r=0 j=r+l ls=r+l J+ 2 uT (k + r/k)[L22 (/7-1 -r) + R(k + r)]u(k + r/k),где L11(s), L12(s), L22(s) определяются по формулам (15) - (18). Выражение (25) нетрудно записать в матричном виде:J (k + p/k ) = zT (k )[Q (/7 )-R (k )]z (k ) ++2 zT (k )G (k )U (k) + UT (k) H (k )U (k),где матрицы G(k) и H(k) имеют вид (12), (13). Ограничения вида (2) - (4) в матричном виде выглядят как (9), где матрицы wmm(k), wmax(k), S(k) имеют вид (10), (11) и (14) соответственно.Итак, имеем задачу минимизации критерия (26) при ограничениях (9), которая эквивалентна задаче квадратичного программирования с критерием (8) при ограничениях (9).3. Численное моделированиеРассмотрим рынок, на котором присутствуют акции четырех видов с ожидаемыми доходностями \iv = 0,011, \i2 = 0,009, ц3 = 0,007 , ц4 = 0,006 и постоянными волатильностями о{ = 0,033 , ст2 = 0,029, ст3 = 0,026 , ст4 = 0,023 ; wt (k) -независимые последовательности белых гауссовских шумов с нулевыми средними и единичными дисперсиями. На рис. 1 приведена динамика цен акций, по оси абсцисс отложены временные интервалы, а по оси ординат - цены акций.0,6 I 1 1 1 1 1 1 1 1 1 Есть возможность вкладывать капитал в банк по ставке тх = 0,001 и брать кредит в банке по ставке тх = 0,002. Инвестор включает в свой ИП акции только первого, второго и третьего типов, поэтому и™11 (к) = (к) = 0 . На объемы торговых операций наложены следующие ограничения: операции продажи без покрытия для инвестора запрещены и™11 (к) = 0, i = 1, n ; объем займа в банке ограниченвеличиной, равной половине начального капитала ИП и™^ (к) = 2V (0). В начальный момент капитал ИП V (0) = 1. Коэффициент превышения ß = 0,1 (10%).Индексный портфель содержит акции только второго, третьего и четвертого типов и характеризуется следующими долями вложений: а! = 0, а2 = 0,3, а3 = 0,3, а4 = 0,4. Матрица весовых коэффициентов в функции риска R = diag(0,01; 0,01; 0,01; 0,01; 0,01). Задача решается на промежутке k е [1; 100].Результаты решения представлены на рис. 2 и 3, где на оси абсцисс указаны номера временных интервалов, а по оси ординат - суммы вложений.3и 1 1 1 1 1 1 1 Г"2\1,510,5041,6ЗаключениеВ данной работе рассмотрена задача управления инвестиционным портфелем при ограничениях на объемы торговых операций, таких как инвестирование в отдельные виды финансовых активов, продажа без покрытия, заем средств по безрисковой ставке. Целью управления портфелем является превышение капитала индексного портфеля, отражающего состояние рынка (активное управление).Для синтеза оптимальных стратегий адаптивоного управления применялся метод прогнозирующего управления, который позволяет получить управление с обратной связью при ограничениях на управляющие переменные.Результаты численного моделирования подтверждают работоспособность и эффективность предложенной модели.
Tsay R.S. Analysis of Financial Time Series. John Wiley & Sons, Inc., 2002.
Rawlings J.B. Tutorial: Model Predictive Control Technology // Proceed. Amer. Control Conf. June. 1999. San Diego. P. 662 - 676.
Runggaldier W.J. On stochastic control in finance. Mathematical Systems Theory in Biology, Communication, Computation and Finance (D. Gilliam and J. Rosental, eds). IMA Book Series (MTNS-2002), Springer Verlag. 2002.
Шарп У., Александер Г., Бейли Дж. Инвестиции. М.: ИНФРА-М, 1997.
Browne S. Risk-Constrained Dynamic Active Portfolio Management // Management Sci. 2000. V. 46. No. 9. P. 1188 - 1199.
Домбровский В.В., Ляшенко Е.А. Динамическая модель управления инвестиционным портфелем на финансовом рынке со стохастической волатильностью // Автоматика и вычислительная техника. 2003. № 5. С. 12 - 21.
Домбровский В.В., Домбровский Д.В., Ляшенко Е.А. Управление с прогнозирующей моделью системами со случайными зависимыми параметрами при ограничениях и применение к оптимизации инвестиционного портфеля // Автоматика и телемеханика. 2006. № 12. C. 71 - 85.
Домбровский В.В., Домбровский Д.В., Ляшенко Е.А. Управление с прогнозированием системами со случайными параметрами и мультипликативными шумами и применение к оптимизации инвестиционного портфеля // Автоматика и телемеханика. 2005. № 4. C. 84 - 97.
Домбровский В.В., Федосов Е.Н. Модель управления инвестиционным портфелем в пространстве состояний на нестационарном диффузионно-скачкообразном финансовом рынке // Автоматика и вычислительная техника. 2002. № 6. С. 13 - 24.
Герасимов Е.С., Домбровский В.В. Динамическая сетевая модель управления инвестиционным портфелем при квадратической функции риска // Автоматика и телемеханика. 2002. № 2. C. 119 - 128.
Bielecki J.R., Pliska S.R. Risk-Sensitive Dynamic Asset Management // Appl. Math. Optim. 1999. V. 39. P. 337 - 360.
Li D., Ng W.-L. Optimal Dynamic Portfolio Selection: Multi-Period Mean-Variance Formulation // Math. Finance. 2000. No. 10. P. 387 - 406.
Merton R.C. Continuous-time finance. Cambr. Ma. Blackwell, 1990.
Kushner H.J. Consistency Issues for Numerical Methods for Variance Control with Applications to Optimization in Finance // IEEE Trans. on Automatic Contr. 1999. V. 44. No. 12. P. 2283 - 2296.
Dupakova J. Portfolio Optimisation via Stochastic Programming: Methods of Output Analysis // Math. Methods Operat. Res. 1999. No. 50. P. 245 - 270.
Young M.R. A Minimax Portfolio Selection Rule with Linear Programming Solution // Management Science. 1998. V. 44. No. 5. P. 673 - 683.
Golub B., Holmer M., McKendall R., Pohiman I., Zenios S.A. A Stochastic Programming Model for Money Management // Europ. J. Operation. Res. 1995. V. 85. P. 282 - 296.
Markowitz H. Portfolio Selection // J. Finance. 1952. V. 7. № 1. P. 77 - 91.
Tobin J. Liquidity Preference as Behavior Towards Risk // Review of Economic Studies. 1958. V. 26. No. 1. P. 65 - 86.