Опционная цена по модели Бейтса с помощью численного интегрирования
добавляет необязательные аргументы пары "имя-значение". Price
= optByBatesNI(___,Name,Value
)
optByBatesNI
использует численное интегрирование, чтобы вычислить опцию цены и затем построить опцию поверхность цены.
Задайте переменные опции и Параметры модели 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;
Вычислите цену опции для одиночного забастовки
Settle = datenum('29-Jun-2017'); Maturity = datemnth(Settle, 6); Strike = 80; Call = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, V0, ... ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ... 'DividendYield', DividendYield)
Call = 5.3484
Вычислите опционные цены для вектора забастовок
The Strike
вход может быть вектором.
Settle = datenum('29-Jun-2017'); Maturity = datemnth(Settle, 6); Strike = (76:2:84)'; Call = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, V0, ... ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ... 'DividendYield', DividendYield)
Call = 5×1
7.5765
6.4020
5.3484
4.4173
3.6073
Вычислите Опцию цены для вектора забастовок и вектора дат той же длины
Используйте 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 = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, V0, ... ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ... 'DividendYield', DividendYield) % Five values in vector output
Call = 5×1
9.7516
10.3931
10.8865
11.2990
11.6491
Разверните выход для поверхности
Установите ExpandOutput
аргумент пары "имя-значение" в "true"
чтобы развернуть выход в NStrikes
-by- NMaturities
матрица. В этом случае это квадратная матрица.
Call = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, V0, ... ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ... 'DividendYield', DividendYield, '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.9951 11.2990 12.4553
5.8705 7.6581 9.1594 10.4771 11.6491
Вычислите Опцию цены для вектора забастовок и вектора дат разной длины
Когда 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 = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, V0, ... ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ... 'DividendYield', DividendYield, '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.9951 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 = optByBatesNI(Rate, ManyAssetPrices, Settle, Maturity, OptSpec, Strike, V0, ... ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ... 'DividendYield', DividendYield, 'ExpandOutput', true) % (5 x 4) matrix output
Call = 5×4
4.2033 6.6918 9.7516 13.2808
3.5558 5.8111 8.6554 11.9993
2.9906 5.0181 7.6432 10.7934
2.5018 4.3096 6.7153 9.6651
2.0825 3.6818 5.8705 8.6158
Постройте график поверхности цены опции
The Strike
и 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)'; Call = optByBatesNI(Rate, AssetPrice, Settle, Maturity, OptSpec, Strike, V0, ... ThetaV, Kappa, SigmaV, RhoSV, MeanJ, JumpVol, JumpFreq, ... 'DividendYield', DividendYield, '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]);
Rate
- Постоянно сложная процентная ставка без рискаПостоянно сложенная процентная ставка без риска, заданная как скалярное десятичное значение.
Типы данных: double
AssetPrice
- Текущая базовая цена активаТекущая базовая цена актива, заданная в виде числа значения с использованием скаляра или NINST
-by- 1
или NColumns
-by- 1
вектор.
Для получения дополнительной информации о соответствующих размерностях для AssetPrice
, см. аргумент пары "имя-значение" ExpandOutput
.
Типы данных: double
Settle
- Дата расчета опции Дата расчета опции, заданная как NINST
-by- 1
или NColumns
-by- 1
вектор с последовательными номерами дат, векторами символов даты, массивами datetime или строковыми массивами. The Settle
дата должна быть перед Maturity
дата.
Для получения дополнительной информации о соответствующих размерностях для Settle
, см. аргумент пары "имя-значение" ExpandOutput
.
Типы данных: double
| char
| datetime
| string
Maturity
- Дата погашения опцииДата погашения опции, заданная как 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
вектор страйк-цен.
Для получения дополнительной информации о соответствующих размерностях для 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 = optByBatesNI(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'
- Флаг, указывающий на формулировку «ловушка маленького Хестона»true
(по умолчанию) | логическим со значениями true
или false
Флаг, указывающий на формулировку «Little Heston Trap» Альбрехером и
др., заданный как разделенная разделенными запятой парами, состоящая из 'LittleTrap' и логический:
true
- Используйте
композицию Albrecher et al.
false
- Использовать исходную формацию Хестона.
Типы данных: logical
'AbsTol'
- Абсолютный допуск ошибок для численного интегрирования1e-10
(по умолчанию) | числоАбсолютный допуск ошибки для численного интегрирования, заданный как разделенная разделенными запятой парами, состоящая из 'AbsTol'
и скалярное числовое значение.
Типы данных: double
'RelTol'
- Относительный допуск ошибок для численного интегрирования1e-6
(по умолчанию) | числоОтносительная погрешность для численного интегрирования, заданный как разделенная разделенными запятой парами, состоящая из 'RelTol'
и скалярное числовое значение.
Типы данных: double
'IntegrationRange'
- Численная область значений интегрирования, используемый для аппроксимации непрерывного интеграла по [0 Inf]
[1e-9 Inf]
(по умолчанию) | векторЧисленная область значений интегрирования, используемый для аппроксимации непрерывного интеграла по [0 Inf]
, заданная как разделенная разделенными запятой парами, состоящая из 'IntegrationRange'
и a 1
-by- 2
вектор, представляющий [LowerLimit UpperLimit]
.
Типы данных: double
'Framework'
- Среды для вычисления цен на опции и чувствительности с помощью численного интегрирования моделей"heston1993"
(по умолчанию) | строку со значениями "heston1993"
или "lewis2001"
| вектор символов со значениями 'heston1993'
или 'lewis2001'
Среда для вычисления цен на опции и чувствительности с помощью численного интегрирования моделей, заданная как разделенная разделенными запятой парами, состоящая из 'Framework'
и скалярный строковый или символьный вектор со следующими значениями:
"heston1993"
или 'heston1993'
- Метод, используемый в Heston (1993)
"lewis2001"
или 'lewis2001'
- Метод, используемый в Lewis (2001)
Типы данных: char
| string
'ExpandOutput'
- Флаг для расширения выходовfalse
(выходы NINST
-by- 1
векторы) (по умолчанию) | логическим со значением true
или false
Флаг для расширения выходов, заданный как разделенная разделенными запятой парами, состоящая из 'ExpandOutput'
и логический:
true
- Если true
, выходы NRows
-by- NColumns
матрицы. NRows
количество ударов для каждого столбца и определяется Strike
вход. Для примера, Strike
может быть NRows
-by- 1
вектор, или NRows
-by- NColumns
матрица. 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
.
A vanilla option - это категория опций, которая включает только самые стандартные компоненты.
Ванильная опция имеет срок годности и прямолинейную цену доставки. Опции в американском стиле и опции в европейском стиле классифицируются как опции ванили.
Выплата для ванильной опции следующая:
Для вызова:
Для размещения:
где:
St - цена базового актива на t времени.
K - цена доставки.
Для получения дополнительной информации смотрите Опцию Vanilla.
Модель Бейтса (Bates (1996)) является расширением модели Хестона, где, помимо стохастической волатильности, для моделирования внезапных движений цен на активы были добавлены и параметры диффузии перехода, подобные Мертону (1976).
Стохастическое дифференциальное уравнение является:
где
r - непрерывная безрисковая ставка.
q - непрерывное дивидендное выражение.
S t является ценой актива в момент t.
v t является отклонением цены основного средства в момент t.
J - размер случайного процентного перехода, обусловленный происходящим переходом, где ln
(1 + J) обычно распределяется со средним и стандартное отклонение, и (1 + J) имеет логнормальное распределение:
v 0 является начальным отклонением цены актива при t = 0 (v 0 > 0).
θ - долгосрочный уровень отклонений для (θ > 0).
κ - средняя скорость реверсии для (κ > 0).
σ v является волатильностью отклонения для (σ v > 0 ).
p является корреляцией между процессами Вайнера W t и для (-1 ≤ p ≤ 1).
μ J является средним значением J для (μ J > -1 ).
δ - стандартное отклонение ln
(1 + J) для (δ ≥ 0 ).
- годовая частота (интенсивность) Пуассоновского процесса P t для ( ≥ 0).
Функция характеристики для j = 1 (средняя мера цены актива) и j = 2 (нейтральная к риску мера) является:
где
ϕ - переменная функции характеристики.
ƛ VolRisk является премией за риск волатильности.
τ - время погашения (τ = T - t).
i является единичным мнимым числом для (i2= -1).
Определения для C j и D j под «Ловушкой маленького Хестона» Albrecher et al. (2007):
Численное интегрирование используется для вычисления непрерывного интеграла для обратного преобразования Фурье.
Метод численного интегрирования в среде Heston (1993) основан на следующих выражениях:
где
r - непрерывная безрисковая ставка.
q - непрерывное дивидендное выражение.
S t является ценой актива в момент t.
K - это удар.
τ время к зрелости (τ = T - t).
Call (K) - цена вызова при K забастовки.
Put (K) - положительная цена при забастовке K
i является модулем мнимым числом (i2= -1)
.r- характеристическая функциональная переменная.
f j (.r) является характеристической функцией для P j (j = 1,2).
P 1 - это вероятность S t > K под измерением цены актива для модели.
P 2 - это вероятность S t > K под нейтральной к риску мерой для модели.
Где j = 1,2 так, что f 1 (в) и f 2 (в) являются характеристическими функциями для вероятностей P 1 и P 2, соответственно.
Эта среда выбирается со значением по умолчанию “Heston1993”
для Framework
аргумент пары "имя-значение".
Численное интегрирование используется для вычисления непрерывного интеграла для обратного преобразования Фурье.
Метод численного интегрирования в среде Льюиса (2001) основан на следующих выражениях:
где
r - непрерывная безрисковая ставка.
q - непрерывное дивидендное выражение.
S t является ценой актива в момент t.
K - это удар.
τ время к зрелости (τ = T - t).
Call (K) - цена вызова при K забастовки.
Put (K) - положительная цена при забастовке K
i является модулем мнимым числом (i2= -1)
.r- характеристическая функциональная переменная.
u - переменная функции характеристики для интегрирования, где .
f 2 (.r) является характеристической функцией для P 2.
P 2 - это вероятность S t > K под нейтральной к риску мерой для модели.
Эта среда выбирается со значением “Lewis2001”
для Framework
аргумент пары "имя-значение".
[1] Бейтс, Д. С. «Переходы и стохастическая волатильность: процессы обменного курса, неявные в опциях Deutsche Mark». Обзор финансовых исследований. Vol 9. № 1. 1996.
[2] Heston, S. L. «A Closed-Form Решения for Опций with Stochastic Volatility with Applications to Bond and Currency Опций». Обзор финансовых исследований. Vol 6. № 2. 1993.
[3] Lewis, A. L. «A Simple Option Formula for General Jump-Diffusion and Other Exponential Levy Processes». Envision Financial Systems and OptionCity.net, 2001.
optByBatesFFT
| optByHestonFFT
| optByHestonNI
| optByMertonFFT
| optByMertonNI
| optSensByBatesFFT
| optSensByBatesNI
| optSensByHestonFFT
| optSensByHestonNI
| optSensByMertonFFT
| optSensByMertonNI
У вас есть измененная версия этого примера. Вы хотите открыть этот пример с вашими правками?
1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.
2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.
3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.
4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.
5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.