Скачиваний:
94
Добавлен:
01.05.2014
Размер:
3.45 Mб
Скачать

Теперь рассмотрим частный случай, когда акция не выплачивает

дивидендов (следовательно, θ 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 величины Е [М(Т ) ехр {− } ; h*].

~

 

Л. Шепп и А. Ширяев показали, что имеется число k , которое

зависит только от r, ρ и σ , такое, что если S(0) >

~

k М, оптимальной

276

стратегией – это исполнение опциона в первый же момент времени t,

~

М(t).

 

 

когда S(t) > k

~

 

Покажем,

как можно определить

очень понятным образом.

k

Пусть k является числом 0 <k <1. Для текущей цены акции S = S(0) при S рассмотрим стратегию исполнения опциона в момент ос-

тановки Тk = inf { t | S(t) = k М (t)} .

Стоимость этой стратегии обозначим R(S, М; k). Заметим, что

R(S, М; k) = М R(S /М, 1; k).

Отсюда и из определений λ и µ следует, что

R(S, М; k) = М λ (S, , М) + R(М, М; k) µ (S, , М) =

= М [λ (S, , М) + R(1, 1; k) µ (S, , М)].

Подставляя явные значения λ и µ ,

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) = ахα 1ebx, 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α − 1ebx , х > 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α − 1ebx

= ахα − 1ebx, x > 0, уравнение

dx cθ = δ

0

или

282

δ + cθ

 

 

(eθ x 1)xα − 1eb*xdx =

.

(6.85)

 

0

a

 

 

 

 

Явное выражение интеграла в левой части равенства (6.85) можно получить с помощью ПФМ случайной величины Y(t).

Если α = 0, то равенство (6.85) приобретает форму

b *

 

 

 

cθ +

δ

 

 

 

=

e

a

.

(6.86)

b *

θ

 

 

 

 

 

 

Подстановка b*, полученного из формулы (6.86), в выражение (6.81) приводит к равенству

ecθ /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