Документация

Квадратичная модель

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

Позвольте обозначают, что ковариационная матрица уровней актива возвращается.

Классическая модель среднего отклонения состоит из минимизации портфельного риска, как измерено

подвергните набору ограничений.

Ожидаемый доход должен быть не меньше, чем минимальной ставкой портфеля, возвращают это желания инвестора,

сумма инвестиционных частей должна составить в целом в общей сложности один,

и, будучи частями (или проценты), они должны быть числами между нулем и один,

Поскольку цель минимизировать портфельный риск квадратична, и ограничения линейны, получившаяся задача оптимизации является квадратичной программой или QP.

Проблема С 225 активами

Давайте теперь решим QP с 225 активами. Набором данных является от ИЛИ-БИБЛИОТЕКИ [Чанг, T.-J., Мид, N., Бизли, J.E. и Sharaiha, Y.M., "Эвристика для кардинальности ограничила оптимизацию портфеля" Computers & Operations Research 27 (2000) 1271-1302].

Мы загружаем набор данных и затем настраиваем ограничения для основанного на проблеме подхода. В этом наборе данных нормы прибыли располагаются между-0.008489 и 0.003971; мы выбираем желаемый промежуточный возврат, например, 0.002 (0,2 процента).

Загрузите набор данных, сохраненный в MAT-файле.

load('port5.mat','Correlation','stdDev_return','mean_return')

Вычислите ковариационную матрицу от корреляционной матрицы.

Covariance = Correlation .* (stdDev_return * stdDev_return');
nAssets = numel(mean_return); r = 0.002;     % number of assets and desired return

Создайте задачу оптимизации, цель и ограничения

Создайте задачу оптимизации для минимизации.

portprob = optimproblem;

Создайте переменную 'x' вектора оптимизации с элементами nAssets. Эта переменная представляет часть богатства, которое инвестируют в каждый актив, так должен находиться между 0 и 1.

x = optimvar('x',nAssets,'LowerBound',0,'UpperBound',1);

Целевой функцией является 1/2*x'*Covariance*x. Включайте эту цель в проблему.

objective = 1/2*x'*Covariance*x;
portprob.Objective = objective;

Сумма переменных равняется 1, означая, что целый портфель инвестируют. Выразите это как ограничение и поместите его в проблему.

sumcons = sum(x) == 1;
portprob.Constraints.sumcons = sumcons;

Средняя доходность должна быть больше, чем r. Выразите это как ограничение и поместите его в проблему.

averagereturn = dot(mean_return,x) >= r;
portprob.Constraints.averagereturn = averagereturn;

Решите проблему с 225 активами

Установите некоторые опции и вызовите решатель.

Установите опции включать итеративное отображение и устанавливать более трудный допуск завершения оптимальности.

options = optimoptions('quadprog','Display','iter','TolFun',1e-10);

Вызовите решатель и измерьте тактовое стеной время.

tic
[x1,fval1] = solve(portprob,'Options',options);
toc

 Iter            Fval  Primal Infeas    Dual Infeas  Complementarity  
    0    7.212813e+00   1.227500e+02   1.195948e+00     2.217295e-03  
    1    8.160874e-04   3.615084e-01   3.522160e-03     2.250524e-05  
    2    7.220766e-04   3.592574e-01   3.500229e-03     3.378157e-05  
    3    4.309434e-04   9.991108e-02   9.734292e-04     2.790551e-05  
    4    4.734300e-04   5.551115e-16   7.771561e-16     4.242216e-06  
    5    4.719034e-04   6.661338e-16   3.122502e-16     8.002618e-07  
    6    3.587475e-04   4.440892e-16   3.035766e-18     3.677066e-07  
    7    3.131814e-04   8.881784e-16   3.686287e-18     9.586695e-08  
    8    2.760174e-04   7.771561e-16   1.463673e-18     1.521063e-08  
    9    2.345751e-04   1.110223e-15   1.138412e-18     4.109608e-09  
   10    2.042487e-04   1.221245e-15   1.084202e-18     6.423267e-09  
   11    1.961775e-04   1.110223e-16   9.757820e-19     6.068329e-10  
   12    1.949281e-04   4.440892e-16   9.215718e-19     4.279951e-12  

Minimum found that satisfies the constraints.

Optimization completed because the objective function is non-decreasing in 
feasible directions, to within the value of the optimality tolerance,
and constraints are satisfied to within the value of the constraint tolerance.

Elapsed time is 0.292806 seconds.

Постройте результаты.

plotPortfDemoStandardModel(x1.x)

Проблема с 225 активами с ограничениями группы

Мы теперь добавляем к образцовым ограничениям группы, которые требуют, чтобы 30% денег инвестора инвестировали в активы 1 - 75, 30% в активах 76 - 150, и 30% в активах 151 - 225. Каждая из этих групп активов могла быть, например, различными отраслями промышленности, такими как технология, автомобильная, и фармацевтическая. Ограничения, которые получают это новое требование,

Добавьте ограничения группы в существующие равенства.

grp1 = sum(x(1:75)) >= 0.3;
grp2 = sum(x(76:150)) >= 0.3;
grp3 = sum(x(151:225)) >= 0.3;
portprob.Constraints.grp1 = grp1;
portprob.Constraints.grp2 = grp2;
portprob.Constraints.grp3 = grp3;

Вызовите решатель и измерьте тактовое стеной время.

tic
[x2,fval2] = solve(portprob,'Options',options);
toc

 Iter            Fval  Primal Infeas    Dual Infeas  Complementarity  
    0    7.212813e+00   1.227500e+02   3.539920e-01     5.253824e-03  
    1    7.004556e-03   2.901399e+00   8.367185e-03     2.207460e-03  
    2    9.181962e-04   4.095630e-01   1.181116e-03     3.749424e-04  
    3    7.515047e-04   3.567918e-01   1.028932e-03     3.486333e-04  
    4    4.238346e-04   9.005778e-02   2.597127e-04     1.607718e-04  
    5    3.695008e-04   1.909891e-04   5.507829e-07     1.341881e-05  
    6    3.691407e-04   6.146337e-07   1.772508e-09     6.817457e-08  
    7    3.010636e-04   7.691892e-08   2.218223e-10     1.837302e-08  
    8    2.669065e-04   1.088252e-08   3.138350e-11     5.474712e-09  
    9    2.195767e-04   8.122574e-10   2.342425e-12     2.814320e-08  
   10    2.102910e-04   2.839773e-10   8.189470e-13     1.037476e-08  
   11    2.060985e-04   6.713696e-11   1.936133e-13     2.876950e-09  
   12    2.015107e-04   0.000000e+00   8.131516e-19     1.522226e-10  
   13    2.009670e-04   4.440892e-16   8.673617e-19     5.264375e-13  

Minimum found that satisfies the constraints.

Optimization completed because the objective function is non-decreasing in 
feasible directions, to within the value of the optimality tolerance,
and constraints are satisfied to within the value of the constraint tolerance.

Elapsed time is 0.162406 seconds.

Постройте результаты, наложенные на результаты предыдущей проблемы.

plotPortfDemoGroupModel(x1.x,x2.x);

Сводные данные результатов до сих пор

Мы видим из второй столбиковой диаграммы, что в результате дополнительных ограничений группы портфель теперь более равномерно распределяется через три группы актива, чем первый портфель. Эта наложенная диверсификация также привела к небольшому увеличению риска, как измерено целевой функцией (см., что столбец маркировал "f(x)" для последней итерации в итеративном отображении для обоих выполнений).

Проблема с 1000 активами Используя случайные данные

В порядке показать, как решатель ведет себя на большей проблеме, мы будем использовать случайным образом сгенерированный набор данных с 1000 активами. Мы генерируем случайную корреляционную матрицу (симметричный, положительно-полуопределенный с единицами на диагонали) использование функции gallery в MATLAB®.

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

rng(0,'twister');
nAssets = 1000; % desired number of assets

Создайте случайные данные

Сгенерируйте средние значения возвратов между-0.1 и 0.4.

a = -0.1; b = 0.4;
mean_return = a + (b-a).*rand(nAssets,1);
r = 0.15;                                     % desired return

Сгенерируйте стандартные отклонения возвратов между 0,08 и 0.6.

a = 0.08; b = 0.6;
stdDev_return = a + (b-a).*rand(nAssets,1);

Загрузите корреляционную матрицу, которая была сгенерирована с помощью Correlation = gallery('randcorr',nAssets). (Генерация корреляционной матрицы этого размера требует времени, так загрузите предсгенерированный вместо этого.)

load('correlationMatrixDemo.mat','Correlation');

Вычислите ковариационную матрицу от корреляционной матрицы.

Covariance = Correlation .* (stdDev_return * stdDev_return');

Создайте задачу оптимизации, цель и ограничения

Создайте задачу оптимизации для минимизации.

portprob2 = optimproblem;

Создайте переменную 'x' вектора оптимизации с элементами nAssets.

x = optimvar('x',nAssets,'LowerBound',0,'UpperBound',1);

Включайте целевую функцию в проблему.

objective = 1/2*x'*Covariance*x;
portprob2.Objective = objective;

Включайте ограничения, что сумма переменных равняется 1, и средняя доходность больше, чем r.

sumcons = sum(x) == 1;
portprob2.Constraints.sumcons = sumcons;
averagereturn = dot(mean_return,x) >= r;
portprob2.Constraints.averagereturn = averagereturn;

Решите проблему с 1000 активами

Вызовите решатель и измерьте тактовое стеной время.

tic
x3 = solve(portprob2,'Options',options);
toc

 Iter            Fval  Primal Infeas    Dual Infeas  Complementarity  
    0    2.142849e+01   5.490000e+02   3.031839e+00     5.210929e-03  
    1    9.378552e-03   6.439102e+00   3.555978e-02     6.331676e-04  
    2    1.128129e-04   3.705915e-03   2.046582e-05     1.802721e-05  
    3    1.118804e-04   1.852958e-06   1.023291e-08     1.170562e-07  
    4    8.490176e-05   7.650016e-08   4.224702e-10     7.048637e-09  
    5    3.364597e-05   4.440892e-16   3.062871e-18     1.037370e-09  
    6    1.980189e-05   2.220446e-16   8.876905e-19     8.465558e-11  

Minimum found that satisfies the constraints.

Optimization completed because the objective function is non-decreasing in 
feasible directions, to within the value of the optimality tolerance,
and constraints are satisfied to within the value of the constraint tolerance.

Elapsed time is 0.913357 seconds.

Сводные данные

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

Более тщательно продуманные исследования возможны при помощи функций, специально предназначенных для оптимизации портфеля в Financial Toolbox™.