optSensByMertonFFT

Цена опции и чувствительность моделью Merton76 с помощью БПФ и FRFT

Синтаксис

[PriceSens,StrikeOut] = optSensByMertonFFT(Rate,AssetPrice,Settle,Maturity,OptSpec,Strike,Sigma,MeanJ,JumpVol,JumpFreq)

[PriceSens,StrikeOut] = optSensByMertonFFT(___,Name,Value)

Описание

пример

[PriceSens,StrikeOut] = optSensByMertonFFT(Rate,AssetPrice,Settle,Maturity,OptSpec,Strike,Sigma,MeanJ,JumpVol,JumpFreq) вычисляет европейскую цену опции ванили и чувствительность моделью Merton76, с помощью БПФ Топкого-места-Madan и методов Chourdakis FRFT.

пример

[PriceSens,StrikeOut] = optSensByMertonFFT(___,Name,Value) добавляют дополнительные аргументы пары "имя-значение".

Примеры

свернуть все

Рабочий процесс для графического вывода поверхности чувствительности опции Используя модель Merton76

Скрипт Open Live Script

Используйте optSensByMertonFFT чтобы калибровать сетку забастовки БПФ для чувствительности, вычислите чувствительность опции и постройте поверхности чувствительности опции.

Задайте переменные опции и параметры модели Merton76

AssetPrice = 80;
Rate = 0.03;
DividendYield = 0.02;
OptSpec = 'call';

Sigma = 0.16;
MeanJ = 0.02;
JumpVol = 0.08;
JumpFreq = 2;

Вычислите цены опции за целый БПФ (или FRFT) сетка забастовки, не задавая "забастовку"

Вычислите чувствительность опции и также выведите соответствующие забастовки. Если Strike вход пуст ([] ), чувствительность опции будет вычислена на целом БПФ (или FRFT) сетка забастовки. БПФ (или FRFT) сетка забастовки определяется как exp(log-strike grid), где каждый столбец сетки логарифмической забастовки имеет NumFFT точки с LogStrikeStep разрядка, которые примерно сосредоточены вокруг каждого элемента log(AssetPrice). Значение по умолчанию для NumFFT 2^12. В дополнение к чувствительности в первом выходе дополнительный последний выход содержит соответствующие забастовки.

Settle = datenum('29-Jun-2017');
Maturity = datemnth(Settle, 6);
Strike = []; % Strike is not specified

[Delta, Kout] = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta");

% Show the lowest and highest strike values on the FFT strike grid
format
[Kout(1) Kout(end)]

ans = 1×2
10¹³⁸ ×

    0.0000    1.8798

% Show a subset of the strikes and corresponding option sensitivities
Range = (2046:2052);
[Kout(Range) Delta(Range)]

ans = 7×2

   50.4929    0.9895
   58.8640    0.9801
   68.6231    0.8816
   80.0000    0.5283
   93.2631    0.1551
  108.7251    0.0241
  126.7505    0.0025

Измените количество БПФ (или FRFT) точки и сравните с optSensByMertonNI

Попробуйте различное количество БПФ (или FRFT) точки и сравните результаты с численным интегрированием. В отличие от optSensByMertonFFT, который использует БПФ (или FRFT) методы для быстрого расчета через целую область значений забастовок, optSensByMertonNI функционируйте использует прямое численное интегрирование, и это обычно медленнее, специально для нескольких забастовок. Однако значения вычисляются optSensByMertonNI может служить сравнительным тестом для корректировки настроек для optSensByMertonFFT.

% Try a smaller number of FFT points 
% (e.g. for faster performance or smaller memory footprint)
NumFFT = 2^10; % Smaller than the default value of 2^12
Strike = []; % Strike is not specified (will use the entire FFT strike grid)
[Delta, Kout] = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT);

% Compare with numerical integration method
Range = (510:516);
Strike = Kout(Range);
DeltaFFT = Delta(Range);
DeltaNI = optSensByMertonNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta");
Error = abs(DeltaFFT-DeltaNI);
table(Strike, DeltaFFT, DeltaNI, Error)

ans=7×4 table
    Strike     DeltaFFT      DeltaNI        Error   
    ______    __________    __________    __________

    12.696       0.89726       0.99002      0.092766
    23.449       0.93421       0.99002       0.05581
    43.312       0.94691       0.99001      0.043093
        80       0.50983       0.52827      0.018446
    147.76      0.004147    0.00019101      0.003956
    272.93      0.001071     1.547e-09      0.001071
    504.11    0.00030521    5.7578e-10    0.00030521

Внесите дальнейшие корректировки в БПФ (или FRFT)

Если значения в выходе DeltaFFT существенно отличаются от тех в DeltaNI, попытайтесь внести изменения в optSensByMertonFFT настройки, такие как CharacteristicFcnStep, LogStrikeStep, NumFFT, DampingFactor, и так далее. Обратите внимание на то, что, если (LogStrikeStep * CharacteristicFcnStep) 2*pi/NumFFT, БПФ используется. В противном случае FRFT используется.

Strike = []; % Strike is not specified (will use the entire FFT or FRFT strike grid)
[Delta, Kout] = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, 'LogStrikeStep', 0.001);

% Compare with numerical integration method
Strike = Kout(Range);
DeltaFFT = Delta(Range);
DeltaNI = optSensByMertonNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta");
Error = abs(DeltaFFT-DeltaNI);
table(Strike, DeltaFFT, DeltaNI, Error)

ans=7×4 table
    Strike    DeltaFFT    DeltaNI      Error   
    ______    ________    _______    __________

    79.76     0.53701     0.53701    5.6407e-12
    79.84      0.5341      0.5341    5.3257e-12
    79.92     0.53119     0.53119    5.0099e-12
       80     0.52827     0.52827    4.6956e-12
    80.08     0.52536     0.52536    4.3811e-12
    80.16     0.52245     0.52245    4.0653e-12
    80.24     0.51953     0.51953    3.7503e-12

% Save the final FFT (or FRFT) strike grid for future reference. For
% example, it provides information about the range of Strike inputs for
% which the FFT (or FRFT) operation is valid.
FFTStrikeGrid = Kout;
MinStrike = FFTStrikeGrid(1) % Strike cannot be less than MinStrike

MinStrike = 47.9437

MaxStrike = FFTStrikeGrid(end) % Strike cannot be greater than MaxStrike

MaxStrike = 133.3566

Вычислите чувствительность опции для одной забастовки

Однажды желаемый БПФ (или FRFT) настройки определяются, используют Strike введите, чтобы задать забастовки вместо того, чтобы обеспечить пустой массив. Если заданные забастовки не совпадают со значением на БПФ (или FRFT) сетка забастовки, выходные параметры интерполированы на заданных забастовках.

Settle = datenum('29-Jun-2017');
Maturity = datemnth(Settle, 6);
Strike = 80;

Delta = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, 'LogStrikeStep', 0.001)

Delta = 0.5283

Вычислите чувствительность опции для вектора забастовок

Используйте Strike введите, чтобы задать забастовки.

Settle = datenum('29-Jun-2017');
Maturity = datemnth(Settle, 6);
Strike = (76:2:84)';

Delta = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, 'LogStrikeStep', 0.001)

Delta = 5×1

    0.6727
    0.6013
    0.5283
    0.4565
    0.3883

Вычислите чувствительность опции для вектора забастовок и вектора дат тех же длин

Используйте Strike введите, чтобы задать забастовки. Кроме того, Maturity введите может быть вектор, но он должен совпадать с длиной Strike вектор, если ExpandOutput аргумент пары "имя-значение" не установлен в "true".

Settle = datenum('29-Jun-2017');
Maturity = datemnth(Settle, [12 18 24 30 36]); % Five maturities
Strike = [76 78 80 82 84]'; % Five strikes

Delta = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001) % Five values in vector output

Delta = 5×1

    0.6419
    0.5907
    0.5565
    0.5311
    0.5110

Расширьте Выходные параметры для поверхности

Установите ExpandOutput аргумент пары "имя-значение" "true" расширять выходные параметры в NStrikes- NMaturities матрицы. В этом случае они - квадратные матрицы.

[Delta, Kout] = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001, 'ExpandOutput', true) % (5 x 5) matrix output

Delta = 5×5

    0.6419    0.6305    0.6245    0.6204    0.6173
    0.5922    0.5907    0.5905    0.5905    0.5905
    0.5422    0.5507    0.5565    0.5607    0.5637
    0.4927    0.5112    0.5229    0.5311    0.5372
    0.4447    0.4725    0.4898    0.5020    0.5110

Kout = 5×5

    76    76    76    76    76
    78    78    78    78    78
    80    80    80    80    80
    82    82    82    82    82
    84    84    84    84    84

Вычислите чувствительность опции для вектора забастовок и вектора дат различных длин

Когда ExpandOutput "true", NStrikes не должны совпадать с NMaturities. Таким образом, выход NStrikes- NMaturities матрица может быть прямоугольной.

Settle = datenum('29-Jun-2017');
Maturity = datemnth(Settle, 12*(0.5:0.5:3)'); % Six maturities
Strike = (76:2:84)'; % Five strikes

Delta = optSensByMertonFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001, 'ExpandOutput', true) % (5 x 6) matrix output

Delta = 5×6

    0.6727    0.6419    0.6305    0.6245    0.6204    0.6173
    0.6013    0.5922    0.5907    0.5905    0.5905    0.5905
    0.5283    0.5422    0.5507    0.5565    0.5607    0.5637
    0.4565    0.4927    0.5112    0.5229    0.5311    0.5372
    0.3883    0.4447    0.4725    0.4898    0.5020    0.5110

Вычислите чувствительность опции для вектора забастовок и вектора цен активов

Когда ExpandOutput "true", выходом может также быть NStrikes- NAssetPrices прямоугольная матрица путем принятия вектора цен активов.

Settle = datenum('29-Jun-2017');
Maturity = datemnth(Settle, 12); % Single maturity
ManyAssetPrices = [70 75 80 85]; % Four asset prices
Strike = (76:2:84)'; % Five strikes

Delta = optSensByMertonFFT(Rate, ManyAssetPrices, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, 'OutSpec', "delta", ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001, 'ExpandOutput', true) % (5 x 4) matrix output

Delta = 5×4

    0.3796    0.5157    0.6419    0.7472
    0.3315    0.4637    0.5922    0.7043
    0.2874    0.4137    0.5422    0.6592
    0.2474    0.3664    0.4927    0.6128
    0.2117    0.3224    0.4447    0.5657

Постройте поверхности чувствительности опции

Используйте Strike введите, чтобы задать забастовки. Увеличьте значение для NumFFT поддерживать более широкую область значений забастовок. Кроме того, Maturity введите может быть вектор. Установите ExpandOutput к "true" выводить поверхности как NStrikes- NMaturities матрицы.

Settle = datenum('29-Jun-2017');
Maturity = datemnth(Settle, 12*[1/12 0.25 (0.5:0.5:3)]');
Times = yearfrac(Settle, Maturity);
Strike = (2:2:200)';

% Increase 'NumFFT' to support a wider range of strikes
NumFFT = 2^13;

[Delta, Gamma, Rho, Theta, Vega] = optSensByMertonFFT(...
    Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    Sigma, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield, ...
    'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, 'LogStrikeStep', 0.001, ...
    'OutSpec', ["delta", "gamma", "rho", "theta", "vega"], ...
    'ExpandOutput', true);

[X,Y] = meshgrid(Times,Strike);

figure;
surf(X,Y,Delta);
title('Delta');
xlabel('Years to Option Expiry');
ylabel('Strike');
view(-112,34);
xlim([0 Times(end)]);

figure;
surf(X,Y,Gamma)
title('Gamma')
xlabel('Years to Option Expiry')
ylabel('Strike')
view(-112,34);
xlim([0 Times(end)]);

figure;
surf(X,Y,Rho)
title('Rho')
xlabel('Years to Option Expiry')
ylabel('Strike')
view(-112,34);
xlim([0 Times(end)]);

figure;
surf(X,Y,Theta)
title('Theta')
xlabel('Years to Option Expiry')
ylabel('Strike')
view(-112,34);
xlim([0 Times(end)]);

figure;
surf(X,Y,Vega)
title('Vega')
xlabel('Years to Option Expiry')
ylabel('Strike')
view(-112,34);
xlim([0 Times(end)]);

Входные параметры

свернуть все

`Rate` — Постоянно составляемая безрисковая процентная ставка
десятичное число

Постоянно составляемая безрисковая процентная ставка в виде скалярного десятичного значения.

Типы данных: double

`AssetPrice` — Текущая цена базового актива
числовой

Текущая цена базового актива в виде числового значения с помощью скаляра или NINST- 1 или NColumns- 1 вектор.

Для получения дополнительной информации о соответствующих размерностях для AssetPrice, смотрите аргумент пары "имя-значение" ExpandOutput.

Типы данных: double

`Settle` — Расчетный день опции
последовательный номер даты | вектор символов даты | datetime | массив строк

Расчетный день опции в виде NINST- 1 или NColumns- 1 вектор с помощью последовательных чисел даты, векторов символов даты, массивов datetime или строковых массивов. Settle дата должна быть перед Maturity дата.

Для получения дополнительной информации о соответствующих размерностях для Settle, смотрите аргумент пары "имя-значение" ExpandOutput.

Типы данных: double | char | datetime | string

`Maturity` — Дата погашения опции
последовательный номер даты | вектор символов даты | datetime | массив строк

Дата погашения опции в виде NINST- 1 или NColumns- 1 вектор с помощью последовательных чисел даты, векторов символов даты, массивов datetime или строковых массивов.

Для получения дополнительной информации о соответствующих размерностях для Maturity, смотрите аргумент пары "имя-значение" ExpandOutput.

Типы данных: double | char | datetime | string

`OptSpec` — Определение опции
массив ячеек вектора символов со значениями `'call'` или `'put'` | массив строк со значениями `"call"` или `"put"`

Определение опции в виде NINST- 1 или NColumns- 1 вектор с помощью массива ячеек из символьных векторов или строковых массивов со значениями 'call' или 'put'.

Для получения дополнительной информации о соответствующих размерностях для OptSpec, смотрите аргумент пары "имя-значение" ExpandOutput.

Типы данных: cell | string

`Strike` — Значение цены исполнения опциона опции
числовой

Значение цены исполнения опциона опции в виде NINST- 1, NRows- 1, NRows- NColumns вектор цен исполнения опциона.

Если этот вход является пустым массивом ([]), цены опции вычисляются на целом БПФ (или FRFT) сетка забастовки, которая определяется как exp(log-strike grid). Каждый столбец сетки логарифмической забастовки has'NumFFT' точки с 'LogStrikeStep' разрядка, которые примерно сосредоточены вокруг каждого элемента log(AssetPrice).

Для получения дополнительной информации о соответствующих размерностях для Strike, смотрите аргумент пары "имя-значение" ExpandOutput.

Типы данных: double

`Sigma` — Энергозависимость базового актива
числовой

Энергозависимость актива подчиненного в виде скалярного числового значения.

Типы данных: double

`MeanJ` — Среднее значение случайного размера скачка процента
десятичное число

Среднее значение случайного размера скачка процента (J) в виде скалярного десятичного значения, где log(1+J) нормально распределено со средним значением (log(1+MeanJ)-0.5*JumpVol^2) и стандартное отклонение JumpVol.

Типы данных: double

`JumpVol` — Стандартное отклонение `log`(1+J)
десятичное число

Стандартное отклонение log(1+J), где J случайный размер скачка процента в виде скалярного десятичного значения.

Типы данных: double

`JumpFreq` — Ежегодная частота процесса скачка Пуассона
числовой

Ежегодная частота процесса скачка Пуассона в виде скалярного числового значения.

Типы данных: double

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример:

[PriceSens,StrikeOut] = optSensByMertonFFT(Rate,AssetPrice,Settle,Maturity,OptSpec,Strike,Sigma,MeanJ,JumpVol,JumpFreq,'Basis',7)

`'Basis'` — Основание дневного количества инструмента
0 (значение по умолчанию) | числовые значения: `0`,1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13

Дневное количество инструмента в виде разделенной запятой пары, состоящей из 'Basis' и скаляр с помощью поддерживаемого значения:

0 = фактический/фактический
1 = 30/360 (СИА)
2 = Фактический/360
3 = Фактический/365
4 = 30/360 (PSA)
5 = 30/360 (ISDA)
6 = 30/360 (европеец)
7 = Фактический/365 (японский язык)
8 = фактический/фактический (ICMA)
9 = Фактический/360 (ICMA)
10 = Фактический/365 (ICMA)
11 = 30/360E (ICMA)
12 = Фактический/365 (ISDA)
13 = ШИНА/252

Для получения дополнительной информации смотрите Основание.

Типы данных: double

`'DividendYield'` — Постоянно составляемая доходность базовых активов
0 (значение по умолчанию) | числовой

Постоянно составляемый базовый актив уступает в виде разделенной запятой пары, состоящей из 'DividendYield' и скалярное числовое значение.

Типы данных: double

`'OutSpec'` — Задайте выходные параметры
`["price"]` (значение по умолчанию) | массив строк со значениями `"price"\delta\Gamma`, `"vega"\rho`, и `"theta"` | массив ячеек из символьных векторов со значениями `'price'\delta\Gamma`, `'vega'\rho`, и `'theta'`

Задайте выходные параметры в виде разделенной запятой пары, состоящей из 'OutSpec' и NOUT- 1 или 1- NOUT массив строк или массив ячеек из символьных векторов с поддерживаемыми значениями.

Примечание

"vega" чувствительность с уважением начальная энергозависимость sqrt (V0).

Пример: OutSpec = ["price","delta","gamma","vega","rho","theta"]

Типы данных: string | cell

`'NumFFT'` — Количество узлов решетки в переменной характеристической функции
4096 (значение по умолчанию) | числовой

Количество узлов решетки в переменной характеристической функции и в каждом столбце сетки логарифмической забастовки в виде разделенной запятой пары, состоящей из 'NumFFT' и скалярное числовое значение.

Типы данных: double

`'CharacteristicFcnStep'` — Интервал сетки переменной характеристической функции
0.01 (значение по умолчанию) | числовой

Интервал сетки переменной характеристической функции в виде разделенной запятой пары, состоящей из 'CharacteristicFcnStep' и скалярное числовое значение.

Типы данных: double

`'LogStrikeStep'` — Интервал сетки логарифмической забастовки
`2*pi/NumFFT/CharacteristicFcnStep` (значение по умолчанию) | числовой

Интервал сетки логарифмической забастовки в виде разделенной запятой пары, состоящей из 'LogStrikeStep' и скалярное числовое значение.

Примечание

Если (LogStrikeStep*CharacteristicFcnStep) 2*pi/NumFFT, БПФ используется. В противном случае FRFT используется.

Типы данных: double

`'DampingFactor'` — Затухание фактора для формулировки Топкого-места-Madan
1.5 (значение по умолчанию) | числовой

Затухание фактора для формулировки Топкого-места-Madan в виде разделенной запятой пары, состоящей из 'DampingFactor' и скалярное числовое значение.

Типы данных: double

`'Quadrature'` — Тип квадратуры
`"simpson"` (значение по умолчанию) | вектор символов с values:`'simpson'` или `'trapezoidal'` | массив строк со значениями: `"simpson"` или `"trapezoidal"`

Тип квадратуры в виде разделенной запятой пары, состоящей из 'Quadrature' и односимвольный векторный массив или массив строк со значением 'simpson' или 'trapezoidal'.

Типы данных: char | string

`'ExpandOutput'` — Отметьте, чтобы расширить выходные параметры
`false` (выходными параметрами является `NINST`- `1` векторы) (значение по умолчанию) | логический со значением `true` или `false`

Отметьте, чтобы расширить выходные параметры в виде разделенной запятой пары, состоящей из 'ExpandOutput' и логическое:

true — Если true, выходными параметрами является NRows- NColumns матрицы. NRows количество борьбы за каждый столбец, и это определяется Strike входной параметр. Например, Strike может быть NRows- 1 вектор или NRows- NColumns матрица. Если Strike пусто, NRows равно NumFFT. NColumns определяется размерами AssetPrice, Settle, Maturity, и OptSpec, который должен все быть или скаляром или NColumns- 1 векторы.
false — Если false, выходными параметрами является NINST- 1 векторы. Кроме того, входные параметры Strike, AssetPrice, Settle, Maturity, и OptSpec должен все быть или скаляр или NINST- 1 векторы.

Типы данных: логический

Выходные аргументы

свернуть все

`PriceSens` — Цены опции или чувствительность
числовой

Цены опции или чувствительность, возвращенная как NINST- 1, или NRows- NColumns, В зависимости от ExpandOutput. Аргумент пары "имя-значение" OutSpec определяет типы и порядок выходных параметров.

`StrikeOut` — Забастовки, соответствующие `Price`
числовой

Забастовки, соответствующие Price, возвращенный как NINST- 1, или NRows- NColumns, В зависимости от ExpandOutput.

Больше о

свернуть все

Опция ванили

vanilla option является категорией опций, которая включает только самые стандартные компоненты.

Опция ванили имеет дату истечения срока и прямую цену исполнения опциона. Американские параметры стиля и европейские параметры стиля оба категоризированы как опции ванили.

Выплата для опции ванили следующие:

Для вызова: $\max (S t - K, 0)$
Для помещенного: $\max (K - S t, 0)$

где:

St является ценой базового актива во время t.

K является ценой исполнения опциона.

Для получения дополнительной информации см. Опцию Ванили.

Модель диффузии скачка Мертона

Модель диффузии скачка Мертона (Мертон (1976)) является различным расширением модели Black-Scholes, где внезапные перемещения цен активов (оба вверх и вниз) моделируются путем добавления параметров диффузии скачка с Пуассоновским процессом.

Стохастическое дифференциальное уравнение:

$\begin{array}{l} d S_{t} = (r - q - λ_{p} μ_{j}) S_{t} d t + σ S_{t} d W_{t} + J S_{t} d P_{t} \\ prob (d P_{t} = 1) = λ_{p} d t \end{array}$

где

r является непрерывным безрисковым уровнем.

q является непрерывной дивидендной доходностью.

W _t является процессом Вайнера.

J является случайным условным выражением размера скачка процента на появлении скачка, где ln(1+J) нормально распределено со средним значением $\ln (1 + μ_{J}) - \frac{δ^{2}}{2}$ и стандартное отклонение δ, и (1+J) имеет логарифмически нормальное распределение:

$\frac{1}{(1 + J) δ \sqrt{2 π}} \exp {{\frac{- [\ln (1 + J) - (\ln (1 + μ_{J}) - \frac{δ^{2}}{2}]}{2 δ^{2}}}^{2}}$

μ _J является средним значением J для (μ _J>-1).

δ является стандартным отклонением ln(1+J) для (δ ≥ 0).

ƛ _p является ежегодной частотой (интенсивность) Пуассоновского процесса P _t для (ƛ _p ≥ 0).

σ является энергозависимостью цены активов на (σ> 0).

Характеристическая функция $f_{M e r t o n 76_{j}} (ϕ)$ для j = 1 (мера цен активов) и j = 2 (нейтральная к риску мера):

$\begin{array}{l} f_{M e r t o n 76_{j}} = f_{B S_{j}} \exp (λ_{p} τ {(1 + μ_{j})}^{m_{j} + \frac{1}{2}} [{(1 + μ_{j})}^{i ϕ} e^{δ^{2} (m_{j} i ϕ + \frac{{(i ϕ)}^{2}}{2})} - 1] - λ_{p} τ μ_{j} i ϕ) \\ где для j = 1, 2 : \\ f_{B S_{1}} (ϕ) = \frac{f_{B S_{2}} (ϕ - i)}{f_{B S_{2}} (- i)} \\ f_{B S 2} (ϕ) = \exp (i ϕ [\ln S_{t} + (r - q - \frac{σ^{2}}{2}) τ] - \frac{ϕ^{2} σ^{2}}{2} τ) \\ m_{1} = \frac{1}{2}, m_{2} = - \frac{1}{2} \end{array}$

где

ϕ является переменной характеристической функции.

τ является временем к зрелости (τ = T - t).

i является модульным мнимым числом (i ² =-1).

Формулировка топкого-места-Madan

Топкое место и Мадан (1999) формулировка являются популярной модифицированной реализацией Хестона (1993) среда.

Вместо того, чтобы вычислять вероятности, P ₁ и P ₂ как промежуточное звено продвигается, Топкое место и Мадан разработали альтернативное выражение так, чтобы взятие его обратного преобразования Фурье дало саму цену опции непосредственно.

$\begin{array}{l} C a l l (k) = \frac{e^{- α k}}{π} \int_{0}^{\infty} Re [e^{- i u k} ψ (u)] d u \\ ψ (u) = \frac{e^{- r τ} f_{2} (ϕ = (u - (α + 1) i))}{α^{2} + α - u^{2} + i u (2 α + 1)} \\ P u t (K) = C a l l (K) + K e^{- r τ} - S_{t} e^{- q τ} \end{array}$

где

r является непрерывным безрисковым уровнем.

q является непрерывной дивидендной доходностью.

S _t является ценой активов во время t.

τ время к зрелости (τ = T-t).

Call (K) является досрочной ценой в забастовке K.

Put (K) является помещенной ценой в забастовке K.

i является модульным мнимым числом (i ² =-1).

ϕ переменная характеристической функции.

α фактор затухания.

u является переменной характеристической функции для интегрирования, где ϕ = (u - (α + 1) i).

f ₂ (ϕ) является характеристической функцией для P ₂.

P ₂ является вероятностью S _t> K под нейтральной к риску мерой для модели.

Чтобы применить БПФ или FRFT к этой формулировке, переменная характеристической функции для интегрирования, u, дискретизируется в NumFFT(N) точка с размером шага CharacteristicFcnStep (Δu) и логарифмическая забастовка k дискретизируются в точки N с размером шага LogStrikeStep(Δk).

Дискретизированная переменная характеристической функции для интегрирования, u _j (для j = 1,2,3, …, N), имеет минимальное значение 0 и максимальное значение (N-1) (Δu), и это аппроксимирует непрерывный диапазон интегрирования от 0 до бесконечности.

Дискретизированная сетка логарифмической забастовки, k _n (для n = 1, 2, 3, N) приблизительно сосредоточена вокруг ln(S _t), с минимальным значением

$\ln (S_{t}) - \frac{N}{2} Δ k$

и максимальное значение

$\ln (S_{t}) + (\frac{N}{2} - 1) Δ k$

Где минимальная допустимая забастовка

$S_{t} \exp (- \frac{N}{2} Δ k)$

и максимальная допустимая забастовка

$S_{t} \exp [(\frac{N}{2} - 1) Δ k]$

В результате дискретизации выражение для колл-опциона становится

$C a l l (k_{n}) = Δ u \frac{e^{- α k_{n}}}{π} \sum_{j = 1}^{N} Re [e^{- i Δ k Δ u (j - 1) (n - 1) e^{i u_{j}} [\frac{N Δ k}{2} - \ln (S_{t})]} ψ (u_{j})] w_{j}$

где

Δu является размером шага дискретизированной переменной характеристической функции для интегрирования.

Δk является размером шага дискретизированной логарифмической забастовки.

N является количеством точек FRFT или БПФ.

w _j является весами для квадратуры, используемой в аппроксимации интеграла.

БПФ используется, чтобы выполнить вышеупомянутое выражение, если Δk и Δu подвергаются следующему ограничению:

$Δ k Δ u = (\frac{2 π}{N})$

в противном случае функции используют метод FRFT, описанный в Chourdakis (2005).

Ссылки

[1] Убавляет, D. S. “Скачки и стохастическая энергозависимость: процессы обменного курса, неявные в опциях немецкой марки”. Анализ финансовых исследований. Vol 9. № 1. 1996.

[2] Топкое место, P. и Д.Б. Мадан. “Оценка опции Используя быстрое преобразование Фурье”. Журнал вычислительных финансов. Vol 2. № 4. 1999.

[3] Продолжение следует, R. и П. Танков. Финансовое моделирование с процессами скачка. Chapman & Hall/CRC Press, 2004.

[4] Chourdakis, K. “Оценка опции Используя дробный БПФ”. Журнал вычислительных финансов. 2005.

[5] Мертон, R. “Оценка опции, Когда Базовые Возвраты Запаса Прерывисты”. Журнал Финансовой Экономики. Vol 3. 1976.

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

optSensByMertonFFT

Синтаксис

Описание

Примеры

Рабочий процесс для графического вывода поверхности чувствительности опции Используя модель Merton76

Входные параметры

`Rate` — Постоянно составляемая безрисковая процентная ставка
десятичное число

`AssetPrice` — Текущая цена базового актива
числовой

`Settle` — Расчетный день опции
последовательный номер даты | вектор символов даты | datetime | массив строк

`Maturity` — Дата погашения опции
последовательный номер даты | вектор символов даты | datetime | массив строк

`OptSpec` — Определение опции
массив ячеек вектора символов со значениями `'call'` или `'put'` | массив строк со значениями `"call"` или `"put"`

`Strike` — Значение цены исполнения опциона опции
числовой

`Sigma` — Энергозависимость базового актива
числовой

`MeanJ` — Среднее значение случайного размера скачка процента
десятичное число

`JumpVol` — Стандартное отклонение `log`(1+J)
десятичное число

`JumpFreq` — Ежегодная частота процесса скачка Пуассона
числовой

Аргументы в виде пар имя-значение

`'Basis'` — Основание дневного количества инструмента
0 (значение по умолчанию) | числовые значения: `0`,1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13

`'DividendYield'` — Постоянно составляемая доходность базовых активов
0 (значение по умолчанию) | числовой

Примечание

`'NumFFT'` — Количество узлов решетки в переменной характеристической функции
4096 (значение по умолчанию) | числовой

`'CharacteristicFcnStep'` — Интервал сетки переменной характеристической функции
0.01 (значение по умолчанию) | числовой

`'LogStrikeStep'` — Интервал сетки логарифмической забастовки
`2*pi/NumFFT/CharacteristicFcnStep` (значение по умолчанию) | числовой

Примечание

`'DampingFactor'` — Затухание фактора для формулировки Топкого-места-Madan
1.5 (значение по умолчанию) | числовой

`'Quadrature'` — Тип квадратуры
`"simpson"` (значение по умолчанию) | вектор символов с values:`'simpson'` или `'trapezoidal'` | массив строк со значениями: `"simpson"` или `"trapezoidal"`

Выходные аргументы

`PriceSens` — Цены опции или чувствительность
числовой

`StrikeOut` — Забастовки, соответствующие `Price`
числовой

Больше о

Опция ванили

Модель диффузии скачка Мертона

Формулировка топкого-места-Madan

Ссылки

Смотрите также

Темы

Введенный в R2018a

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

Поддержка

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

optSensByMertonFFT

Синтаксис

Описание

Примеры

Рабочий процесс для графического вывода поверхности чувствительности опции Используя модель Merton76

Входные параметры

Rate — Постоянно составляемая безрисковая процентная ставка десятичное число

AssetPrice — Текущая цена базового актива числовой

Settle — Расчетный день опции последовательный номер даты | вектор символов даты | datetime | массив строк

Maturity — Дата погашения опции последовательный номер даты | вектор символов даты | datetime | массив строк

OptSpec — Определение опции массив ячеек вектора символов со значениями 'call' или 'put' | массив строк со значениями "call" или "put"

Strike — Значение цены исполнения опциона опции числовой

Sigma — Энергозависимость базового актива числовой

MeanJ — Среднее значение случайного размера скачка процента десятичное число

JumpVol — Стандартное отклонение log(1+J) десятичное число

JumpFreq — Ежегодная частота процесса скачка Пуассона числовой

Аргументы в виде пар имя-значение

'Basis' — Основание дневного количества инструмента0 (значение по умолчанию) | числовые значения: 0,1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13

'DividendYield' — Постоянно составляемая доходность базовых активов0 (значение по умолчанию) | числовой

Примечание

'NumFFT' — Количество узлов решетки в переменной характеристической функции4096 (значение по умолчанию) | числовой

'CharacteristicFcnStep' — Интервал сетки переменной характеристической функции0.01 (значение по умолчанию) | числовой

'LogStrikeStep' — Интервал сетки логарифмической забастовки 2*pi/NumFFT/CharacteristicFcnStep (значение по умолчанию) | числовой

Примечание

'DampingFactor' — Затухание фактора для формулировки Топкого-места-Madan1.5 (значение по умолчанию) | числовой

'Quadrature' — Тип квадратуры "simpson" (значение по умолчанию) | вектор символов с values:'simpson' или 'trapezoidal' | массив строк со значениями: "simpson" или "trapezoidal"

Выходные аргументы

PriceSens — Цены опции или чувствительность числовой

StrikeOut — Забастовки, соответствующие Price числовой

Больше о

Опция ванили

Модель диффузии скачка Мертона

Формулировка топкого-места-Madan

Ссылки

Смотрите также

Темы

Введенный в R2018a

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

Поддержка

`Rate` — Постоянно составляемая безрисковая процентная ставка
десятичное число

`AssetPrice` — Текущая цена базового актива
числовой

`Settle` — Расчетный день опции
последовательный номер даты | вектор символов даты | datetime | массив строк

`Maturity` — Дата погашения опции
последовательный номер даты | вектор символов даты | datetime | массив строк

`OptSpec` — Определение опции
массив ячеек вектора символов со значениями `'call'` или `'put'` | массив строк со значениями `"call"` или `"put"`

`Strike` — Значение цены исполнения опциона опции
числовой

`Sigma` — Энергозависимость базового актива
числовой

`MeanJ` — Среднее значение случайного размера скачка процента
десятичное число

`JumpVol` — Стандартное отклонение `log`(1+J)
десятичное число

`JumpFreq` — Ежегодная частота процесса скачка Пуассона
числовой

`'Basis'` — Основание дневного количества инструмента
0 (значение по умолчанию) | числовые значения: `0`,1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13

`'DividendYield'` — Постоянно составляемая доходность базовых активов
0 (значение по умолчанию) | числовой

`'NumFFT'` — Количество узлов решетки в переменной характеристической функции
4096 (значение по умолчанию) | числовой

`'CharacteristicFcnStep'` — Интервал сетки переменной характеристической функции
0.01 (значение по умолчанию) | числовой

`'LogStrikeStep'` — Интервал сетки логарифмической забастовки
`2*pi/NumFFT/CharacteristicFcnStep` (значение по умолчанию) | числовой

`'DampingFactor'` — Затухание фактора для формулировки Топкого-места-Madan
1.5 (значение по умолчанию) | числовой

`'Quadrature'` — Тип квадратуры
`"simpson"` (значение по умолчанию) | вектор символов с values:`'simpson'` или `'trapezoidal'` | массив строк со значениями: `"simpson"` или `"trapezoidal"`

`PriceSens` — Цены опции или чувствительность
числовой

`StrikeOut` — Забастовки, соответствующие `Price`
числовой