График в этом примере показывает трехмерную поверхность. Для каждой точки на поверхности высота (значение z) представляет собой сумму гамм для каждого варианта в портфеле, взвешенную на сумму каждого варианта. Ось X представляет изменяющуюся цену, а ось Y - время. График добавляет четвертый размер, показывая дельту в качестве цвета поверхности. Этот пример имеет приложения в хеджировании. Сначала настройте портфель с произвольными данными. Текущие цены варьируются от $20 до $90 за каждый вариант. Затем установите соответствующие цены упражнений для каждого варианта.

Range = 20:90;
PLen = length(Range);
ExPrice = [75 70 50 55 75 50 40 75 60 35];

Установите все безрисковые процентные ставки на уровне 10% и установите сроки погашения в днях. Установите для всех летучих веществ значение 0,35. Задайте количество опций каждого прибора и выделите пространство для матриц.

Rate = 0.1*ones(10,1);
Time = [36  36  36  27  18  18  18  9  9  9];
Sigma = 0.35*ones(10,1);
NumOpt = 1000*[4  8  3  5  5.5  2  4.8  3  4.8  2.5];
ZVal = zeros(36, PLen);
Color = zeros(36, PLen);

Для каждого инструмента создайте матрицу (размера Time около PLen) цен на каждый период.

for i = 1:10

    Pad = ones(Time(i),PLen);
    NewR = Range(ones(Time(i),1),:);

Создание вектора периодов времени 1 кому Time и матрица раз, по одному столбцу для каждой цены.

  T = (1:Time(i))';
  NewT = T(:,ones(PLen,1));

Используйте функции гамма-чувствительности Блэка-Шоулза и дельта-чувствительности blsgamma и blsdelta для вычисления гамма и дельта.

    ZVal(36-Time(i)+1:36,:) = ZVal(36-Time(i)+1:36,:) ...
        + NumOpt(i) * blsgamma(NewR, ExPrice(i)*Pad, ...
        Rate(i)*Pad, NewT/36, Sigma(i)*Pad);

    Color(36-Time(i)+1:36,:) = Color(36-Time(i)+1:36,:) ...
        + NumOpt(i) * blsdelta(NewR, ExPrice(i)*Pad, ...
        Rate(i)*Pad, NewT/36, Sigma(i)*Pad);

end

Нарисуйте поверхность как сеть, задайте точку обзора и измените ось X на противоположную из-за точки обзора. Оси находятся в диапазоне от 20 до 90, от 0 до 36 и от - ∞ до ∞.

mesh(Range, 1:36, ZVal, Color);
view(60,60);
set(gca, 'xdir','reverse', 'tag', 'mesh_axes_3');
axis([20 90  0 36  -inf inf]);

Добавьте заголовок и метки осей и нарисуйте рамку вокруг графика. Аннотируйте цвета полосой и помечайте ее меткой.

title('Call Option Portfolio Sensitivity');
xlabel('Stock Price ($)');
ylabel('Time (months)');
zlabel('Gamma');
set(gca, 'box', 'on');
colorbar('horiz');