- •Оглавление
- •Предисловие
- •Основные сокращения и обозначения
- •Введение
- •Глава 1. Анализ моделей непрерывного времени
- •1. Математические и экономические предположения в моделях непрерывного времени
- •2. Процессы с непрерывными выборочными траекториями без редких событий
- •1. Финансовые производные
- •7. Определение стоимости отзываемого опциона
- •8. Разрывные стохастические процессы изменения цен акций
- •9. Определение стоимости опционов для разрывных стохастических процессов
- •10. Задачи определения стоимости опционов
- •Глава 3. Мартингалы и арбитраж на рынках ценных бумаг
- •1. Основные определения
- •2. Жизнеспособность и арбитраж
- •3. Модели рынка ценных бумаг
- •4. Конечная модель
- •6. Другие торговые стратегии
- •7. Обобщения
- •1. Постановки основных задач
- •2. Конечная теория
- •3. Непрерывная торговля
- •6. Иллюстративные примеры
- •Глава 5. Временная структура процентных ставок: мартингальный подход
- •1. Процесс дисконтированной цены облигации как мартингал
- •3. Случай, когда мгновенная процентная ставка является диффузионным процессом
- •1. Понятие о преобразовании Эсшера
- •2. Нейтральное к риску преобразование Эсшера
- •3. Формулы вычисления цен опционов
- •4. Опционы на несколько рисковых активов
- •5. Логарифмы цен акций как многомерный винеровский процесс
- •6. Цены активов, выплачивающих дивиденды
- •7. Определение цены бессрочного американского опциона
- •8. Логарифм цены акции как винеровский процесс
- •9. Русский опцион
- •10. Квазинепрерывные выборочные траектории
- •Литература
Теперь рассмотрим частный случай, когда акция не выплачивает
дивидендов (следовательно, θ 1 = 1 и θ 0 = − 2r/σ |
2). Покажем, что мак- |
симальное значение (6.67) получается при М → ∞ |
и что |
V(S, L, ∞ ) = S + (R − L) (L/S) − θ 0 = S + (R − L) (L/S) 2rσ 2 .
Этот результат можно также найти у Р. Мертона (1973).
Для доказательства этой формулы мы сначала заметим, что величина λ (S, L, М) является возрастающей функцией М и, следовательно, первое слагаемое правой части равенства (6.67) ограничено
величиной R λ (S, L, ∞ ) = R (L/S) − θ 0 . Второе слагаемое правой части равенства (6.67) можно оценить следующим образом
|
(М − K) µ (S, L, |
М) = (М − K) |
SLθ 0 |
− S θ 0 L |
= |
|
|
MLθ 0 |
− M θ 0 L |
||||
|
|
|
|
|
||
= |
M − K |
[S − L(L S)− θ 0 ]<S − L (L/S) 2r σ 2 = |
||||
|
M − L(L M )− θ 0 |
|
|
|
|
|
|
= |
lim |
(M − K) ( S, L, M) . |
|
||
|
|
M → ∞ |
|
|
|
|
§ 9. РУССКИЙ ОПЦИОН
Пусть М является числом таким, что М ≥ S. Пусть также
М(t) = mах [М, mах { S(и) | 0 ≤ и ≤ t} ],
что можно интерпретировать как максимум цены акции за время t. Заметим, что пара { S(t), М(t); t ≥ 0} является однородным марковским процессом. Термин «русский опцион» предложен Л. Шеппом и А. Ширяевым (Shepp, Shiryaev, 1993) для описания бессрочного американского опциона, платеж которого равен М(t), если он исполняется в момент времени t, t ≥ 0. То есть владелец русского опциона имеет привилегию получить максимум цены акции за время до того момента, в который он решил исполнить опцион. Цена опциона в момент
времени 0 является наибольшим значением по всем моментам оста- |
|
новки Т ≥ 0 величины Е [М(Т ) ехр {− rТ} ; h*]. |
~ |
|
|
Л. Шепп и А. Ширяев показали, что имеется число k , которое |
|
зависит только от r, ρ и σ , такое, что если S(0) > |
~ |
k М, оптимальной |
276
стратегией – это исполнение опциона в первый же момент времени t, |
||||
~ |
М(t). |
|
|
|
когда S(t) > k |
~ |
|
||
Покажем, |
как можно определить |
очень понятным образом. |
||
k |
Пусть k является числом 0 <k <1. Для текущей цены акции S = S(0) при kМ ≤ S рассмотрим стратегию исполнения опциона в момент ос-
тановки Тk = inf { t | S(t) = k М (t)} .
Стоимость этой стратегии обозначим R(S, М; k). Заметим, что
R(S, М; k) = М R(S /М, 1; k).
Отсюда и из определений λ и µ следует, что
R(S, М; k) = М λ (S, kМ, М) + R(М, М; k) µ (S, kМ, М) =
= М [λ (S, kМ, М) + R(1, 1; k) µ (S, kМ, М)].
Подставляя явные значения λ и µ ,
R(S, М; k) = М [λ (S/М, k, 1) + R(1, 1; k) µ (S/М, k, 1)] =
|
|
M |
|
S |
θ 0 |
||
= |
|
|
|
|
|||
|
|
|
|
|
|
||
kθ 0 |
− kθ 1 |
|
|
||||
|
|
|
M |
|
|||
|
|
|
|
|
|
|
− |
|
S θ 1 |
|
+ |
|
|
|
|
|
||
|
|||||
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
θ |
|
S |
θ 1 |
|
|
θ |
|
S |
θ 0 |
|
|
||
+ |
R(1, 1; k) |
k |
|
− |
k |
|
|
(6.68) |
||||||||
|
0 |
|
|
|
|
1 |
|
|
|
, |
||||||
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
M |
|
|
|
|
|
M |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где R(1, 1; k) определяется с помощью условия на границе при S = М. Условие исполнения русского опциона может быть получено путем следующих эвристических рассуждений. Если текущая цена акции S очень близка к М, можно быть уверенным «почти наверняка», что цена акции достигнет уровня М (и, следовательно, что максимум будет увеличиваться) до того, как опцион будет исполнен. Таким образом, если S близко к М, функция R (S, М; k) не зависит от точного
значения М и его производная по М
|
RМ (М, М; k) = 0. |
|
|
|
|
Отсюда и из (6.68) мы получаем условие |
|
|
|
||
[((1− θ 0 ) − (1− θ 1)) + R(1, 1; |
k)(kθ 0 (1− θ 1) − kθ 1 (1 |
− θ 0 ))] |
= 0, |
||
|
kθ 0 − |
kθ 1 |
|
|
|
|
|
|
|
277
что дает
R(1, 1; k) = |
|
θ 1 |
− |
θ 0 |
|
. |
|
kθ 1 (1− θ |
0 |
) − |
kθ 0(1− θ |
) |
|||
|
|
|
|
|
1 |
|
Подставим это выражение в формулу (6.68) и получим после упрощений следующий результат:
R(S, M ;k) = M |
(1− θ |
0 |
)(S M)θ |
1 − |
(1 |
− θ |
) (S M)θ 0 |
|
||
|
|
|
|
|
|
1 |
. |
|||
|
kθ 1 (1 |
− θ |
0 |
) − |
kθ |
0(1 |
− θ ) |
|||
|
|
|
||||||||
|
|
|
|
|
|
|
|
1 |
|
Теперь ясно, что оптимальное значение k является значением, которое минимизирует знаменатель, производная которого равна
(1 |
− θ |
0 |
) θ kθ |
1 − 1 + ( θ |
1 |
− |
1) |
θ |
0 |
kθ 0 − 1. |
|
|
||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||
Следовательно, оптимальное значение |
|
|
|
|
|
|||||||||||
|
|
|
~ |
|
|
θ 0 (1− θ 1) |
|
1 (θ 1 − θ |
0 ) |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
k |
|
|
|
|
|
|
|
|
|
||||
|
|
|
= |
|
θ 1(1− |
θ 0) |
|
|
|
|
|
, |
|
(6.69) |
||
и цена русского опциона равна |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
~ |
|
|
|
|
~ |
|
|
S ≤ |
M , |
|
|
|
R(S, M ; k ), если |
|
|
kM ≤ |
(6.70) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
||
|
M , |
|
|
если |
|
|
S ≤ kM . |
|
|
Формулы (6.69) и (6.70) являются эквивалентными формулам, полученным Л. Шеппом и А. Ширяевым.
§ 10. КВАЗИНЕПРЕРЫВНЫЕ ВЫБОРОЧНЫЕ ТРАЕКТОРИИ
Предположим, что выборочные траектории { S (t)} или, что эквивалентно, { Х(t)} не имеют скачков вниз (это предположение было использовано при получении цены (6.39)). Тогда имеет место следующая декомпозиция:
X(t) = Y(t) + v2W(t) − ct, t ≥ 0. |
(6.71) |
Здесь { Y(t)} является или составным пуассоновским процессом с положительными приращениями или пределом такого процесса;
278
{ W(t)} – стандартный винеровский процесс с нулевым средним и единичной дисперсией за единицу времени; последнее слагаемое ct представляет собой систематический дрейф. Производящая функция семиинвариантов случайной величины X(t) имеет вид
|
∞ |
|
zx |
|
|
2 |
|
2 |
|
|
|
|
|
|
∫ |
(e |
− |
1)[− dQ(x)] + v |
z |
/ 2 |
− |
|
, |
(6.72) |
|||
ln[M(z, t)] = t |
|
|
|
cz |
|||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
где Q (х) – неотрицательная и невозрастающая функция с Q(∞ ) = 0. Заметим, что для всякого положительного числа ε интеграл
∞
∫ (ezx − 1)[− dQ (x)],
ε
как функция от z, является производящей функцией семиинвариантов составного распределения Пуассона с пуассоновским параметром
λ(ε) = Q(ε)
ираспределением величины скачка
Р (х; ε) = |
Q(ε) − Q(x) |
, х ≥ ε. |
|
Q(ε) |
|
Для простоты обозначений предположим, что − Q (х) = q (x) dx для некоторой неотрицательной функции q(x). Пусть µ и σ 2 обозначают соответственно среднее и дисперсию { Х(t)} за единицу времени. Тогда
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
µ |
t = Е[Х(t)] = |
|
∫ xq(x)dx − |
|
t, |
|
|
|
(6.73) |
|||||||
|
|
c |
|
|
|
||||||||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
∞ |
|
2 |
|
|
|
|
|
|
2 |
|
|
|
σ |
|
|
|
∫ x |
q(x)dx + |
v |
|
t |
(6.74) |
||||||||
|
t = var[Х(t)] = |
|
|
|
|||||||||||||
и |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
∞ |
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
∫ x |
|
|
|
t. |
|
|
(6.75) |
|||||
|
Е [(Х(t) − µ t) |
] = |
|
q(x)dx |
|
|
|||||||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
279
Вообще для п ≥ 3 п-й семиинвариант случайной величины Х(t) равен
|
∞ |
n |
|
|
|
∫ x |
|
t. |
|
|
|
q(x)dx |
||
|
0 |
|
|
|
Из формул (6.4) и (6.72) следует, что
ln[M(z, t; h)] = ln[M(z + h, t)] − ln[M(h, t)] =
∞
=t ∫ (ezx − 1)ehx q(x)dx +
0
v2 z2 / 2 − (c − v2h)z . (6.76)
Таким образом, преобразование Эсшера с параметром h процесса, определяемого соотношением (6.71), имеет такой же вид со следующими модификациями:
q(x) → ehx q(x), |
(6.77) |
v2 → v2 (остается прежним), c → |
c − v2 h. |
Кроме того, из (6.76) следует, что (6.28) и (6.36) можно записать соответственно как
∞ |
(ex − 1)eh*xq(x)dx + v2h* = c + δ − ρ − |
v |
2 |
|
||||
∫ |
|
|
||||||
0 |
|
|
|
|
|
2 |
|
|
и |
|
|
|
|
|
|
|
|
∞ |
(eθ x − 1)eh*xq(x)dx + |
v |
2 |
θ |
2 |
|
|
|
∫ |
|
− (c − v2h*)θ = δ . |
(6.78) |
|||||
0 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Частный случай
Для модели, определяемой соотношениями (6.1) и (6.71), теперь предположим, что v = 0, т. е. S (t) = S (0) ехр{Y(t) − ct}, и что
q (x) = ахα − 1e− bx, x > 0,
где а > 0, α > − 1 и b > 0 являются параметрами.
В соответствии с модификацией (6.77) для h <b преобразование Эсшера такого процесса является членом этого же семейства с параметром b, замененным на
b(h) = b − h.
280
ПФМ случайной величины Y(t) имеет вид |
|
|
|
|
||||||||
|
∞ |
|
zx |
|
|
∞ |
zx |
|
α − 1 |
− bx |
|
|
|
t ∫ |
(e |
|
|
at ∫ (e |
− 1)x |
|
= |
||||
exp |
|
− 1)q(x)dx |
= exp |
|
e |
|
dx |
|||||
|
0 |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
b |
|
|
at |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
= |
b − |
z |
|
|
|
|
|
|
|
||
|
|
aΓ(α ) |
|
|
|
α |
|||||
|
|
|
|
b |
|
|
|||||
|
exp |
|
bα |
|
|
|
|
|
|||
|
|
|
|
|
|||||||
|
|
|
|
|
b − |
z |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
если |
α |
= 0, |
− |
|
|
, если |
α |
≠ 0. |
1 t |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, для α = 0 процесс { Y(t)} t ≥ 0 является гаммапроцессом; для α > 0 он является составным процессом Пуассона с пуассоновским параметром λ (а, α , b) = аГ(α ) ⁄ bα и гаммаплотностью вероятностей величины скачков
|
bα |
р (х; α , b) = |
Γ(α ) xα − 1e− bx , х > 0. |
Для − 1 <α <0 наиболее известным случаем является α = − ½, если { Y(t)} t ≥ 0 – обратный гауссов процесс с плотностью вероятностей случайной величины Y(t) вида
at |
|
− |
( |
bx − πat) |
2 |
|
|
|
|
|
|
|
, |
х > 0. |
|||||
|
3 2 |
exp |
|
|
|
|
|
||
x |
|
|
|
x |
|
|
|
||
|
|
|
|
|
|
|
|
Для b* = b(h*) = b − h* получаем уравнения
|
|
|
b * |
|
|
= |
|
c+ δ − ρ |
, если α = 0, |
|||
|
|
|
|
|
e |
|
a |
|||||
и |
|
b* − |
1 |
|
||||||||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
c + |
δ − ρ |
|
||
1 |
|
|
− |
|
|
1 |
|
= |
, если α ≠ 0. |
|||
|
(b * − |
1)α |
|
b*α |
aΓ(α ) |
|||||||
|
|
|
|
|
Решением уравнения (6.79) является
b* = |
1 |
, |
1− e− (c+ δ − ρ ) / a |
(6.79)
(6.80)
(6.81)
281
что при ρ = 0 превращается в формулу (6.15). В общем случае вид уравнения (6.79) не позволяет получить решение для b* в явной форме. Однако если α = 1 (величины скачков распределены экспоненциально), решение имеет вид
|
1 |
|
|
4a |
|
|
b* = |
|
|
1+ 1+ |
|
|
(6.82) |
2 |
|
c + δ − ρ |
. |
|||
|
|
|
|
|
||
Случай α = − ½ рассмотрен выше в § 3. |
|
|
||||
Для каждого фиксированного α |
можно определить параметры а, |
b и с методом моментов. Таким образом, можно считать известными µ , σ и третий центральный момент Х(1), который удобно записать как γσ3 (γ является коэффициентом асимметрии). Эти три момента можно
записать с помощью формул (6.73) − |
(6.75) в форме следующих ра- |
|||||||
венств |
|
|
|
|
|
|
|
|
∞ |
aΓ(α + 1) |
|
∞ |
x2q(x)dx = |
aΓ(α + 2) |
|
||
µ = ∫ xq(x)dx − c = |
− c , σ 2 = ∫ |
, |
||||||
bα + 1 |
bα + 2 |
|||||||
0 |
|
0 |
|
|
|
|||
и |
∞ |
|
|
|
|
|
|
|
|
|
aΓ(α + 3) |
|
|
|
|||
γσ3 = ∫ x3q(x)dx = |
. |
|
|
|||||
|
|
|
||||||
|
0 |
|
bα |
+ 3 |
|
|
|
Из этих равенств получим b = (α + 2)/γσ (при определении стоимости финансовой производной это b нужно заменить на b*)
а = |
(α |
+ 2)α |
+ |
2 |
(6.83) |
|||
Γ(α |
+ |
2)γα |
+ |
2σ α |
||||
и |
|
α |
+ 2 |
σ |
|
|
|
|
с = |
|
− |
µ . |
(6.84) |
||||
|
|
α |
+ 1 |
γ |
|
|
|
Формулы для отрицательных корней уравнения (6.78)
При предположениях v = 0 и q(x) (6.78) превращается в следующее:
∞
a ∫ (eθ x − 1)eh*x xα − 1e− bx
= ахα − 1e− bx, x > 0, уравнение
dx − cθ = δ
0
или
282
∞ |
δ + cθ |
|
|
|
∫ (eθ x − 1)xα − 1e− b*xdx = |
. |
(6.85) |
||
|
||||
0 |
a |
|
||
|
|
|
Явное выражение интеграла в левой части равенства (6.85) можно получить с помощью ПФМ случайной величины Y(t).
Если α = 0, то равенство (6.85) приобретает форму
b * |
|
|
|
cθ + |
δ |
|
|
|
|
= |
e |
a |
. |
(6.86) |
|||
b * − |
θ |
|||||||
|
|
|
|
|
|
Подстановка b*, полученного из формулы (6.86), в выражение (6.81) приводит к равенству
e− cθ /a + θ [e δ /a − e− (c−ρ )/a] = e δ /a.
С помощью формул (6.83) и (6.84) получаем
1 |
= |
γ2 |
|
c |
= |
γ |
σ − |
µγ |
||
|
|
и |
|
|
|
. |
||||
a |
4 |
a |
2 |
|||||||
|
|
|
|
|
2 |
Если α ≠ 0 и α > − 1, то уравнение (6.85) приобретает вид
1 |
− |
1 |
= |
δ + cθ |
, |
(6.87) |
|
(b * − θ )α |
b *α |
aΓ(α ) |
|||||
|
|
|
|
где b* определяется из уравнения (6.80). Уравнение (6.87) упростится при α = 1 к виду
1 |
− |
1 |
|
= |
δ + cθ |
, |
(6.88) |
|
b * − θ |
b * |
a |
||||||
|
|
|
|
что является квадратным уравнением относительно θ , а b* определя-
ется из (6.82),
а = |
27 |
и с = 3 |
σ |
− µ . |
3 |
|
|||
|
2 γ |
|||
|
2γσ |
В случае отсутствия дивидендов (ρ = 0) положительный корень
уравнения (6.88) равен θ 1 = 1, а отрицательный корень выражается в виде θ 0 = −δ b*/с.
283