Смешано-целочисленная оптимизация портфеля квадратичного программирования: основанный на решателе

В этом примере показано, как решить задачу оптимизации портфеля Смешано-целочисленного квадратичного программирования (MIQP) с помощью intlinprog Решатель Смешано-целочисленного линейного программирования (MILP). Идея состоит в том, чтобы итеративно решить последовательность проблем MILP, которые локально аппроксимируют проблему MIQP. Для подхода, основанного на проблеме смотрите Смешано-целочисленную Оптимизацию Портфеля Квадратичного программирования: основанный на проблеме.

Проблемная схема

Когда Markowitz показал ("Выбор Портфеля", J. Финансируйте Объем 7, Выпуск 1, стр 77-91, март 1952), можно описать много задач оптимизации портфеля как проблемы квадратичного программирования. Предположим, что у вас есть набор N активы и хотят выбрать портфель, с x(i) будучи частью ваших инвестиций, которые находятся в активе i. Если вы знаете вектор r из среднего значения возвращается из каждого актива и ковариационной матрицы Q из возвратов, затем для данного уровня нерасположенности к риску λ вы максимизируете настроенный риском ожидаемый доход:

maxx(rTx-λxTQx).

quadprog решатель решает эту проблему квадратичного программирования. Однако в дополнение к простой проблеме квадратичного программирования, вы можете хотеть ограничить портфель во множестве путей, таких как:

  • Наличие не больше, чем M активы в портфеле, где M <= N.

  • Наличие, по крайней мере, m активы в портфеле, где 0 < m <= M.

  • Наличие полунепрерывных ограничений, значение также x(i)=0, или fminx(i)fmax для некоторых фиксированных частей fmin>0 и fmaxfmin.

Вы не можете включать эти ограничения в quadprog. Трудность является дискретной природой ограничений. Кроме того, в то время как решатель смешано-целочисленного линейного программирования intlinprog действительно обрабатывает дискретные ограничения, это не обращается к квадратичным целевым функциям.

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

Начните путем моделирования ограничений.

Моделирование дискретных ограничений

x вектор из частей распределения активов, с 0x(i)1 для каждого i. Чтобы смоделировать количество активов в портфеле, вам нужны переменные индикатора v таким образом, что v(i)=0 когда x(i)=0, и v(i)=1 когда x(i)>0. Чтобы получить переменные, которые удовлетворяют этому ограничению, установите v вектор, чтобы быть бинарной переменной и наложить линейные ограничения

v(i)fminx(i)v(i)fmax.

Эти неравенства оба осуществляют это x(i) и v(i) нуль в точно то же время, и они также осуществляют это fminx(i)fmax каждый раз, когда x(i)>0.

Кроме того, чтобы осуществить ограничения на количество активов в портфеле, наложите линейные ограничения

miv(i)M.

Объективные и последовательные линейные аппроксимации

Как сначала сформулировано, вы пытаетесь максимизировать целевую функцию. Однако все решатели Optimization Toolbox™ минимизируют. Поэтому сформулируйте проблему как минимизацию отрицания цели:

minxλxTQx-rTx.

Эта целевая функция нелинейна. intlinprog Решатель MILP требует линейной целевой функции. Существует стандартный метод, чтобы переформулировать эту проблему в одну с линейными объективными и нелинейными ограничениями. Введите слабую переменную z представлять квадратичный термин.

minx,zλz-rTx таким образом , что xTQx-z0,z0.

Когда вы итеративно решаете приближения MILP, вы включаете новые линейные ограничения, каждое из которых аппроксимирует нелинейное ограничение локально около текущей точки. В частности, для x=x0+δ где x0 постоянный вектор и δ переменный вектор, первый порядок, который приближение Тейлора к ограничению

xTQx-z=x0TQx0+2x0TQδ-z+O(|δ|2).

Заменяя δ x-x0 дает

xTQx-z=-x0TQx0+2x0TQx-z+O(|x-x0|2).

Для каждого промежуточного решения xk вы вводите новое линейное ограничение в x и z как линейная часть выражения выше:

-xkTQxk+2xkTQx-z0.

Это имеет форму Axb, где A=2xkTQ, существует a -1 множитель для z назовите, и b=xkTQxk.

Этот метод добавления новых линейных ограничений к проблеме называется сокращающим плоским методом. Для получения дополнительной информации смотрите Дж. Э. Келли младшего "Плоский Сокращением Метод для Решения Выпуклых Программ". Дж. Сок. Indust. Прикладная Математика. Издание 8, № 4, стр 703-712, декабрь 1960.

MATLAB® Problem Formulation

Описывать проблемы для intlinprog решатель, необходимо сделать следующее:

  • Решите то, что представляют ваши переменные

  • Опишите нижние и верхние границы в терминах этих переменных

  • Дайте линейные матрицы равенства и неравенства

Имейте первое N переменные представляют x вектор, следующее N переменные представляют двоичный файл v вектор и последняя переменная представляют z слабая переменная. Существуют 2N+1 переменные в проблеме.

Загрузите данные для проблемы. Эти данные имеют 225 ожидаемых доходов в векторном r и ковариация возвратов в 225 225 матричном Q. Данные эквивалентны в Использовании Квадратичного программирования на примере Задач оптимизации Портфеля.

load port5
r = mean_return;
Q = Correlation .* (stdDev_return * stdDev_return');

Определите номер активов как N.

N = length(r);

Установите индексы для переменных

xvars = 1:N;
vvars = N+1:2*N;
zvar = 2*N+1;

Нижние границы всего 2N+1 переменные в проблеме являются нулем. Верхние границы первого 2N переменные один, и последняя переменная не имеет никакой верхней границы.

lb = zeros(2*N+1,1);
ub = ones(2*N+1,1);
ub(zvar) = Inf;

Определите номер активов в решении быть между 100 и 150. Включите это ограничение в проблему в форме, а именно,

miv(i)M,

путем записи двух линейных ограничений формы Axb:

iv(i)M

i-v(i)-m.

M = 150;
m = 100;
A = zeros(1,2*N+1); % Allocate A matrix
A(vvars) = 1; % A*x represents the sum of the v(i)
A = [A;-A];
b = zeros(2,1); % Allocate b vector
b(1) = M;
b(2) = -m;

Включайте полунепрерывные ограничения. Возьмите минимальную ненулевую часть активов, чтобы быть 0.001 для каждого типа актива и максимальной части, чтобы быть 0.05.

fmin = 0.001;
fmax = 0.05;

Включайте неравенства x(i)fmax(i)*v(i) и fmin(i)*v(i)x(i) как линейные неравенства.

Atemp = eye(N);
Amax = horzcat(Atemp,-Atemp*fmax,zeros(N,1));
A = [A;Amax];
b = [b;zeros(N,1)];
Amin = horzcat(-Atemp,Atemp*fmin,zeros(N,1));
A = [A;Amin];
b = [b;zeros(N,1)];

Включайте ограничение, которое портфель составляет 100%, которые инвестируют, означая xi=1.

Aeq = zeros(1,2*N+1); % Allocate Aeq matrix
Aeq(xvars) = 1;
beq = 1;

Установите коэффициент нерасположенности к риску λ к 100.

lambda = 100;

Задайте целевую функцию λz-rTx как вектор. Включайте нули для множителей v переменные.

f = [-r;zeros(N,1);lambda];

Решите задачу

Чтобы решить задачу итеративно, начните путем решения задачи с текущими ограничениями, которые еще не отражают линеаризации. Целочисленные ограничения находятся в vvars вектор.

options = optimoptions(@intlinprog,'Display','off'); % Suppress iterative display
[xLinInt,fval,exitFlagInt,output] = intlinprog(f,vvars,A,b,Aeq,beq,lb,ub,options);

Подготовьте останавливающееся условие к итерациям: остановитесь когда слабая переменная z в 0,01% истинного квадратичного значения. Установите более трудные погрешности, чем значение по умолчанию, чтобы помочь гарантировать, что проблема остается строго выполнимой, когда ограничения накапливаются.

thediff = 1e-4;
iter = 1; % iteration counter
assets = xLinInt(xvars); % the x variables
truequadratic = assets'*Q*assets;
zslack = xLinInt(zvar); % slack variable value
options = optimoptions(options,'LPOptimalityTolerance',1e-10,'RelativeGapTolerance',1e-8,...
                      'ConstraintTolerance',1e-9,'IntegerTolerance',1e-6);

Сохраните историю вычисленных истинных квадратичных и слабых переменных для графического вывода.

history = [truequadratic,zslack];

Вычислите квадратичные и слабые значения. Если они отличаются, то добавьте другое линейное ограничение и решите снова.

В синтаксисе тулбокса, каждом новом линейном ограничении Axb прибывает из линейной аппроксимации

-xkTQxk+2xkTQx-z0.

Вы видите что новая строка A=2xkTQ и новый элемент в b=xkTQxk, с z термин, представленный-1 коэффициентом в A.

После того, как вы найдете новое решение, используйте линейное ограничение на полпути между старыми и новыми решениями. Этот эвристический способ включать линейные ограничения может быть быстрее, чем простое взятие нового решения. Чтобы использовать решение вместо промежуточной эвристики, прокомментируйте "На полпути" линия ниже и не прокомментируйте следующую.

while abs((zslack - truequadratic)/truequadratic) > thediff % relative error
    newArow = horzcat(2*assets'*Q,zeros(1,N),-1); % Linearized constraint
    rhs = assets'*Q*assets;                       % right hand side of the linearized constraint
    A = [A;newArow];
    b = [b;rhs];
    % Solve the problem with the new constraints
    [xLinInt,fval,exitFlagInt,output] = intlinprog(f,vvars,A,b,Aeq,beq,lb,ub,options);
    assets = (assets+xLinInt(xvars))/2; % Midway from the previous to the current
%     assets = xLinInt(xvars); % Use the previous line or this one
    truequadratic = xLinInt(xvars)'*Q* xLinInt(xvars);
    zslack = xLinInt(zvar);
    history = [history;truequadratic,zslack];
    iter = iter + 1;
end

Исследуйте уровень решения и сходимости

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

plot(history)
legend('Quadratic','Slack')
xlabel('Iteration number')
title('Quadratic and linear approximation (slack)')

Figure contains an axes. The axes with title Quadratic and linear approximation (slack) contains 2 objects of type line. These objects represent Quadratic, Slack.

Каково качество решения MILP? output структура содержит ту информацию. Исследуйте абсолютный разрыв между внутренне вычисленными границами на цели в решении.

disp(output.absolutegap)
     0

Абсолютный разрыв является нулем, указывая, что решение MILP точно.

Постройте оптимальное выделение. Используйте xLinInt(xvars), не assets, потому что assets не может удовлетворить ограничениям при использовании на полпути обновление.

bar(xLinInt(xvars))
grid on
xlabel('Asset index')
ylabel('Proportion of investment')
title('Optimal Asset Allocation')

Figure contains an axes. The axes with title Optimal Asset Allocation contains an object of type bar.

Можно легко видеть, что все ненулевое распределение активов между полунепрерывными границами fmin=0.001 и fmax=0.05.

Сколько ненулевые активы там? Ограничение состоит в том, что существует между 100 и 150 ненулевыми активами.

sum(xLinInt(vvars))
ans = 100

Каковы ожидаемый доход для этого выделения и значение доходности с учетом риска?

fprintf('The expected return is %g, and the risk-adjusted return is %g.\n',...
    r'*xLinInt(xvars),-fval)
The expected return is 0.000595107, and the risk-adjusted return is -0.0360382.

Более тщательно продуманные исследования возможны при помощи функций, специально предназначенных для оптимизации портфеля в Financial Toolbox™. Для примера, который показывает, как использовать класс Портфеля, чтобы непосредственно обработать полунепрерывный и ограничения кардинальности, смотрите Оптимизацию Портфеля с Полунепрерывным и Ограничения Кардинальности (Financial Toolbox).

Похожие темы

Для просмотра документации необходимо авторизоваться на сайте