optByBatesFFT

Опционная цена по модели Бейтса с использованием FFT и FRFT

Синтаксис

[Price,StrikeOut] = optByBatesFFT(Rate,AssetPrice,Settle,Maturity,OptSpec,Strike,V0,ThetaV,Kappa,SigmaV,RhoSV,MeanJ,JumpVol,JumpFreq)

[Price,StrikeOut] = optByBatesFFT(___,Name,Value)

Описание

пример

[Price,StrikeOut] = optByBatesFFT(Rate,AssetPrice,Settle,Maturity,OptSpec,Strike,V0,ThetaV,Kappa,SigmaV,RhoSV,MeanJ,JumpVol,JumpFreq) вычисляет ванильную европейскую цену опции по модели Bates, используя методы Carr-Madan БПФ и Chourdakis FRFT.

пример

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

Примеры

свернуть все

Рабочий процесс для Графического изображения поверхности Опции цены с помощью модели Бейтса

Открыть Live Script

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

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

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

V0 = 0.04;
ThetaV = 0.05;
Kappa = 1.0;
SigmaV = 0.2;
RhoSV = -0.7;
MeanJ = 0.02;
JumpVol = 0.08;
JumpFreq = 2;

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

Расчет цен на опцию, а также вывод соответствующих забастовок. Если на 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 (will use the entire FFT strike grid)

% Compute option prices for the entire FFT strike grid
[Call, Kout] = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield);

% Show the lowest and highest strike values on the FFT strike grid
format
MinStrike = Kout(1) % Lowest possible strike in the current FFT strike grid

MinStrike = 2.9205e-135

MaxStrike = Kout(end) % Highest possible strike in the current FFT strike grid

MaxStrike = 1.8798e+138

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

ans = 7×2

   50.4929   29.4990
   58.8640   21.4545
   68.6231   12.8544
   80.0000    5.3484
   93.2631    1.2404
  108.7251    0.1648
  126.7505    0.0152

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

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

% Try a smaller number of FFT (or FRFT) 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)
[Call, Kout] = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT);

% Compare with numerical integration method
Range = (510:516);
Strike = Kout(Range);
CallFFT = Call(Range);
CallNI = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield);
Error = abs(CallFFT-CallNI);
table(Strike, CallFFT, CallNI, Error)

ans=7×4 table
    Strike     CallFFT       CallNI         Error  
    ______    _________    ___________    _________

    12.696       66.237         66.696      0.45912
    23.449        55.86         56.103      0.24239
    43.312       36.418         36.541      0.12246
        80       5.4029         5.3484     0.054469
    147.76     0.044921      0.0010864     0.043835
    272.93    0.0094655    -7.8249e-08    0.0094656
    504.11    0.0024986    -3.3873e-07    0.0024989

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

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

Strike = []; % Strike is not specified (will use the entire FFT or FRFT strike grid)
[Call, Kout] = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ...    
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001);

% Compare with numerical integration method
Strike = Kout(Range);
CallFFT = Call(Range);
CallNI = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, 'DividendYield', DividendYield);
Error = abs(CallFFT-CallNI);
table(Strike, CallFFT, CallNI, Error)

ans=7×4 table
    Strike    CallFFT    CallNI      Error   
    ______    _______    ______    __________

    79.76     5.4682     5.4682    1.5355e-08
    79.84     5.4281     5.4281    1.4833e-08
    79.92     5.3882     5.3882    1.4244e-08
       80     5.3484     5.3484     1.359e-08
    80.08     5.3088     5.3088    1.2875e-08
    80.16     5.2693     5.2693    1.2101e-08
    80.24       5.23       5.23    1.1272e-08

% 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; 

Call = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV,  MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001)

Call = 5.3484

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

Используйте Strike вход для определения ударов.

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

Call = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV,  MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001)

Вычислите Опцию цены для вектора забастовок и вектора дат той же длины

Используйте 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

Call = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV,  MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001) % Five values in vector output

Разверните область «Выходы поверхности»

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

[Call, Kout] = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV,  MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001, 'ExpandOutput', true) % (5 x 5) matrix output

Call = 5×5

    9.7516   11.4387   12.8395   14.0588   15.1361
    8.6554   10.3931   11.8344   13.0890   14.1980
    7.6432    9.4149   10.8865   12.1693   13.3046
    6.7153    8.5035    9.9952   11.2990   12.4553
    5.8705    7.6581    9.1594   10.4771   11.6491

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-by- NMaturities матрица может быть прямоугольной).

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

Call = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV,  MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001, 'ExpandOutput', true) % (5 x 6) matrix output

Call = 5×6

    7.5765    9.7516   11.4387   12.8395   14.0588   15.1361
    6.4020    8.6554   10.3931   11.8344   13.0890   14.1980
    5.3484    7.6432    9.4149   10.8865   12.1693   13.3046
    4.4173    6.7153    8.5035    9.9952   11.2990   12.4553
    3.6073    5.8705    7.6581    9.1594   10.4771   11.6491

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

Когда ExpandOutput является "true", выходы могут также быть NStrikes-by- 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

Call = optByBatesFFT(Rate, ManyAssetPrices, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV,  MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001, 'ExpandOutput', true) % (5 x 4) matrix output

Call = 5×4

    4.2033    6.6918    9.7516   13.2808
    3.5558    5.8112    8.6554   11.9993
    2.9906    5.0181    7.6432   10.7934
    2.5018    4.3096    6.7153    9.6652
    2.0825    3.6818    5.8705    8.6158

Постройте график поверхности цены опции

Используйте Strike вход для определения ударов. Увеличьте значение для NumFFT для поддержки более широкой области значений ударов. Кроме того, Maturity вход может быть вектором. Задайте ExpandOutput на "true" для вывода поверхности в виде NStrikes-by- 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;

Call = optByBatesFFT(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, ...
    V0, ThetaV, Kappa, SigmaV, RhoSV,  MeanJ, JumpVol, JumpFreq, ...
    'DividendYield', DividendYield, 'NumFFT', NumFFT, 'CharacteristicFcnStep', 0.065, ...
    'LogStrikeStep', 0.001, 'ExpandOutput', true);

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

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

Figure contains an axes. The axes with title Price contains an object of type surface.

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

свернуть все

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

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

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

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

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

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

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

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

Дата расчета опции, заданная как NINST-by- 1 или NColumns-by- 1 вектор с последовательными номерами дат, векторами символов даты, массивами datetime или строковыми массивами. The Settle дата должна быть перед Maturity дата.

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

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

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

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

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

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

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

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

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

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

`Strike` - значение цены опционной доставки
числовой

Опции цены значения, заданные как NINST-by- 1, NRows-by- 1, NRows-by- NColumns вектор страйк-цен.

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

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

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

`V0` - Начальное отклонение базового актива
числовой

Начальное отклонение нижнего актива, заданное как скалярное числовое значение.

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

`ThetaV` - Долгосрочное отклонение базовых активов
числовой

Долгосрочное отклонение базового актива, заданная в виде скалярного числового значения.

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

`Kappa` - Средняя скорость пересмотра для отклонения базового актива
числовой

Средняя скорость ревизии для нижнего актива, заданная в виде скалярного числового значения.

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

`SigmaV` - Волатильность отклонения базового актива
числовой

Волатильность отклонения базового актива, заданная в виде скалярного числового значения.

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

`RhoSV` - Корреляция между процессами Вайнера для базового актива и его отклонением
числовой

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

Типы данных: 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.

Пример:

[Price, StrikeOut] = optByBatesFFT(Rate,AssetPrice,Settle,Maturity,OptSpec,Strike,V0,ThetaV,Kappa,SigmaV,RhoSV,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 (SIA)
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 = BUS/252

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

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

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

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

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

`'VolRiskPremium'` - Премия за риск волатильности
`0` (по умолчанию) | число

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

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

`'LittleTrap'` - Флаг, указывающий на формулировку ловушки Little Heston
`true` (по умолчанию) | логическим со значениями `true` или `false`

Флаг, указывающий на формулировку Little Heston Trap Альбрехером и др., заданный как разделенная разделенными запятой парами, состоящая из 'LittleTrap' и логический:

true - Используйте композицию Albrecher et al.
false - Использовать исходную формацию Хестона.

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

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

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

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

`'CharacteristicFcnStep'` - Характеристическая функция, переменная интервала между сетками
`0.01` (по умолчанию) | число

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

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

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

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

Примечание

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

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

`'DampingFactor'` - Коэффициент затухания для композиции Карра-Мадана
`1.5` (по умолчанию) | число

Коэффициент затухания для композиции Карра-Мадана, заданная как разделенная запятой пара, состоящая из 'DampingFactor' и скалярное числовое значение.

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

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

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

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

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

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

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

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

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

свернуть все

`Price` - Опционные цены
числовой

Опционные цены, возвращенные как NINST-by- 1, или NRows-by- NColumns, в зависимости от ExpandOutput.

`StrikeOut` - Удары, соответствующие `Price`
числовой

Удары, соответствующие Price, возвращается как NINST-by- 1, или NRows-by- NColumns, в зависимости от ExpandOutput.

Подробнее о

свернуть все

Ванильные Опции

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

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

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

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

где:

St - цена базового актива на t времени.

K - цена доставки.

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

Модель диффузии стохастического перехода волатильности Бейтса

Модель Бейтса (Bates (1996)) является расширением модели Хестона, где, помимо стохастической волатильности, для моделирования внезапных движений цен на активы были добавлены и параметры диффузии перехода, подобные Мертону (1976).

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

$\begin{array}{l} d S_{t} = (r - q - λ_{p} μ_{J}) S_{t} d t + \sqrt{v_{t}} S_{t} d W_{t} + J S_{t} d P_{t} \\ d v_{t} = κ (θ - v_{t}) d t + σ_{v} \sqrt{v_{t}} d W_{t} \\ E [d W_{t} d W_{t}^{v}] = p d t \\ prob (d P_{t} = 1) = λ_{p} d t \end{array}$

где

r - непрерывная безрисковая ставка.

q - непрерывное дивидендное выражение.

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

v t является отклонением цены основного средства в момент 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}}$

v 0 является начальным отклонением цены актива при t = 0 (_v 0 > 0).

θ - долгосрочный уровень отклонений для (θ > 0).

κ - средняя скорость реверсии для (κ > 0).

σ v является волатильностью отклонения для (_σ v > 0 ).

p является корреляцией между процессами Вайнера W t и $W_{t}^{v}$ для (-1 ≤ p ≤ 1).

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

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

$λ_{p}$ - годовая частота (интенсивность) Пуассоновского процесса P t для ( $λ_{p}$ ≥ 0).

Функция характеристики $f_{B a t e s_{j} (ϕ)}$ для j = 1 (средняя мера цены актива) и j = 2 (нейтральная к риску мера) является:

$\begin{array}{l} f_{B a t e s (ϕ)} = \exp (C_{j} + D_{j} v_{0} + i ϕ \ln S_{t}) \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 ϕ) \\ m_{j} = {\begin{cases} m_{1} = \frac{1}{2} \\ m_{2} = - \frac{1}{2} \end{cases}} \\ C_{j} = (r - q) i ϕ τ + \frac{κ θ}{σ_{v}^{2}} [(b_{j} - p σ_{v} i ϕ + d_{j}) τ - 2 \ln (\frac{1 - g_{j} e^{d_{j} τ}}{1 - g_{j}})] \\ D j = \frac{b_{j} - p σ_{v} i ϕ + d_{j}}{σ_{v}^{2}} (\frac{1 - e^{d_{j} τ}}{1 - g_{j} e^{d_{j} τ}}) \\ g_{j} = \frac{b_{j} - p σ_{v} i ϕ + d_{j}}{b_{j} - p σ_{v} i ϕ - d_{j}} \\ d_{j} = \sqrt{{(b_{j} - p σ_{v} i ϕ)}^{2} - σ_{v}^{2} (2 u_{j} i ϕ - ϕ^{2})} \\ где j = 1, 2 : \\ u_{1} = \frac{1}{2}, u_{2} = - \frac{1}{2}, b_{1} = κ + λ_{V o l R i s k} - p σ_{v}, b_{2} = κ + λ_{V o l R i s k} \end{array}$

где

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

ƛ _VolRisk является премией за риск волатильности.

τ - время погашения (τ = T - t).

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

Определения для C j и _D j под «Ловушкой маленького Хестона» Albrecher et al. (2007):

$\begin{array}{l} C_{j} = (r - q) i ϕ τ + \frac{κ θ}{σ_{v}^{2}} [(b_{j} - p σ_{v} i ϕ - d_{j}) τ - 2 \ln (\frac{1 - ε_{j} e^{- d_{j} τ}}{1 - ε_{j}})] \\ D j = \frac{b_{j} - p σ_{v} i ϕ - d_{j}}{σ_{v}^{2}} (\frac{1 - e^{- d_{j} τ}}{1 - ε_{j} e^{- d_{j} τ}}) \\ ε_{j} = \frac{b_{j} - p σ_{v} i ϕ - d_{j}}{b_{j} - p σ_{v} i ϕ + d_{j}} \end{array}$

Формулировка Карра-Мадана

Формулировка Карра и Мадана (1999 год) является популярной модифицированной реализацией среды (1993 год).

Вместо того, чтобы вычислять вероятности P 1 и _P 2 как промежуточные шаги, Карр и Мадан разработали альтернативное выражение, так что взятие обратного преобразования Фурье дает саму опционную цену непосредственно.

$\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)

.r- характеристическая функциональная переменная.

α - коэффициент затухания.

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

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

P 2 - это вероятность _S t > K под нейтральной к риску мерой для модели.

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

Дискретизированная характерная переменная функции для интегрирования, <reservedrangesplaceholder5> <reservedrangesplaceholder4> (для j = 1,2,3..., N), имеет минимальное значение 0 и максимальное значение (N-1) (Δ <reservedrangesplaceholder0>), и это приближает непрерывный диапазон интегрирования от 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 - количество точек FFT/FRFT

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

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

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

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

Ссылки

[1] Albrecher, H., Mayer, P., Schoutens, W., and Tistaert, J. «The Little Heston Trap». Рабочий документ Linz and Graz University of Technology, K.U. Leuven, ING Financial Markets, 2006.

[2] Бейтс, Д. С. «Переходы и стохастическая волатильность: процессы обменного курса, неявные в опциях Deutsche Mark». Обзор финансовых исследований. Vol 9. № 1. 1996.

[3] Карр, П. и Д. Б. Мадан. «Оценка опций с использованием быстрого преобразования Фурье». Журнал вычислительных финансов. Vol 2. № 4. 1999.

[4] Chourdakis, K. «Опционное ценообразование с использованием дробного БПФ». Журнал вычислительных финансов. 2005.

[5] Heston, S. L. «A Closed-Form Решения for Опций with Stochastic Volatility with Applications to Bond and Currency Опций». Обзор финансовых исследований. Vol 6. № 2. 1993.

См. также

Темы

Введенный в R2018a

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

optByBatesFFT

Синтаксис

Описание

Примеры

Рабочий процесс для Графического изображения поверхности Опции цены с помощью модели Бейтса

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

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

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

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

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

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

Strike - значение цены опционной доставки числовой

V0 - Начальное отклонение базового актива числовой

ThetaV - Долгосрочное отклонение базовых активов числовой

Kappa - Средняя скорость пересмотра для отклонения базового актива числовой

SigmaV - Волатильность отклонения базового актива числовой

RhoSV - Корреляция между процессами Вайнера для базового актива и его отклонением числовой

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

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

JumpFreq - Годовая частота процесса перехода Пуассона числовой

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

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

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

'VolRiskPremium' - Премия за риск волатильности 0 (по умолчанию) | число

'LittleTrap' - Флаг, указывающий на формулировку ловушки Little Heston true (по умолчанию) | логическим со значениями true или false

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

'CharacteristicFcnStep' - Характеристическая функция, переменная интервала между сетками 0.01 (по умолчанию) | число

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

'DampingFactor' - Коэффициент затухания для композиции Карра-Мадана 1.5 (по умолчанию) | число

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

'ExpandOutput' - Флаг для расширения выходов false (выходы NINST-by- 1 векторы) (по умолчанию) | логическим со значением true или false

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

Price - Опционные цены числовой

StrikeOut - Удары, соответствующие Price числовой

Подробнее о

Ванильные Опции

Модель диффузии стохастического перехода волатильности Бейтса

Формулировка Карра-Мадана

Ссылки

См. также

Темы

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

Поддержка

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

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

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

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

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

`Strike` - значение цены опционной доставки
числовой

`V0` - Начальное отклонение базового актива
числовой

`ThetaV` - Долгосрочное отклонение базовых активов
числовой

`Kappa` - Средняя скорость пересмотра для отклонения базового актива
числовой

`SigmaV` - Волатильность отклонения базового актива
числовой

`RhoSV` - Корреляция между процессами Вайнера для базового актива и его отклонением
числовой

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

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

`JumpFreq` - Годовая частота процесса перехода Пуассона
числовой

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

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

`'VolRiskPremium'` - Премия за риск волатильности
`0` (по умолчанию) | число

`'LittleTrap'` - Флаг, указывающий на формулировку ловушки Little Heston
`true` (по умолчанию) | логическим со значениями `true` или `false`

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

`'CharacteristicFcnStep'` - Характеристическая функция, переменная интервала между сетками
`0.01` (по умолчанию) | число

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

`'DampingFactor'` - Коэффициент затухания для композиции Карра-Мадана
`1.5` (по умолчанию) | число

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

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

`Price` - Опционные цены
числовой

`StrikeOut` - Удары, соответствующие `Price`
числовой