Формула цены и плотность Бридена-Литценбергера

От цены колла как интеграла выплаты к риск-нейтральной плотности через опционную бабочку.

optionsbreeden-litzenberger

На этой странице

Колл и пут

Вспомним графики опционов на экспирации. Страйк K=100K = 100.

Колл после страйка растёт вместе с акцией. Пут растёт, когда акция падает ниже страйка. Это зеркальные ломаные.

Колл и пут на экспирации, K = 100.

Бабочка

Чтобы вывести формулу Бридена-Литценбергера, нужна опционная бабочка: покупка и продажа коллов с разными страйками.

Бабочка собирается так: купить колл со страйком K−dKK - dK, продать два колла со страйком KK, купить колл со страйком K+dKK + dK. На графике dK=20dK = 20: шатёр от 80 до 120, в страйке выплата 20.

Бабочка на экспирации, K = 100, dK = 20.

Плотность

Функция плотности - это функция, с которой интегрируют функцию выплат, чтобы получить цену опциона в момент времени tt.

Выплата сама от времени не зависит: для колла это всегда max⁡(S−K,0)\max(S - K, 0). Плотность задаёт веса возможных цен в экспирацию. Произведение интегрируют, и выходит цена контракта в момент tt. Для этого момента плотность берётся от цены StS_t до экспирации.

C(St,t)=∫−∞∞f(x) P(x,T∣St,t) dxC(S_t, t) = \int_{-\infty}^{\infty} f(x)\, P(x, T \mid S_t, t)\, dx

Здесь f(x)f(x) - выплата, если в экспирацию цена равна xx. P(x,T∣St,t)P(x, T \mid S_t, t) - плотность: вес этого xx при условии, что в момент tt акция стоит StS_t.

Дальше из этих выплат выводится формула плотности. Для колла выплата (ST−K)(S_T - K) ненулевая только выше страйка:

C(S0,T)=∫K∞(ST−K)⋅f(ST) dSTC(S_0, T) = \int_{K}^{\infty} (S_T - K) \cdot f(S_T) \, dS_T

Цена бабочки

Цена бабочки - разность трёх коллов:

Ctbf=Ct(K−dK)−2Ct(K)+Ct(K+dK)C_{t}^{bf} = C_{t}(K-dK) - 2C_{t}(K)+ C_{t}(K + dK)

Нормировка на dK

Нормируем на dKdK:

CtbfdK=Ct(k−dK)−2Ct(K)+Ct(K+dK)dK≈dK∂2Ct(K)∂K2\frac{C_{t}^{bf}}{dK} = \frac{C_{t}(k-dK) - 2C_{t}(K)+ C_{t}(K + dK)}{dK} \approx dK\frac{\partial^2C_{t}(K)}{\partial K^2}

Приближённо это вторая производная цены по страйку, умноженная на dKdK.

Подстановка выплаты

Перепишем выражение, подставив в интеграл выплату, не цену. В знаменателе один dKdK:

CtbfdK=∫K−dKK+dK(S−(K−dK))+−2(S−K)++(S−(K+dK))+dK⋅P(S,T∣S0,0) dS≈≈P(K,T∣S0,0)dK\begin{aligned} \frac{C_{t}^{bf}}{dK} &= \int_{K-dK}^{K+dK} \frac{(S-(K-dK))^{+} - 2(S-K)^{+} + (S-(K+dK))^{+}}{dK} \cdot P(S,T|S_0, 0)\, dS \approx \\ &\approx P(K,T|S_0, 0)dK \end{aligned}

Откуда берётся множитель dK

Разберём, как получилось последнее соотношение. Дисконтирование не учитываем: цена бабочки - это интеграл её выплаты по плотности.

Ctbf=∫0∞g(S) P(S,T∣S0,0) dS.C_{t}^{bf} = \int_{0}^{\infty} g(S)\, P(S, T \mid S_0, 0)\, dS.

Ниже разберём, как отсюда получается CtbfdK≈P(K,T∣S0,0) dK\dfrac{C_{t}^{bf}}{dK} \approx P(K, T \mid S_0, 0)\, dK.

Площадь шатра

Интеграл выплаты без плотности - это площадь под графиком. У треугольника площадь равна половине произведения основания на высоту:

12⋅2⋅dK⋅dK=dK2\frac{1}{2} \cdot 2 \cdot dK \cdot dK = dK^2 g(S)=(S−K+dK)+−2(S−K)++(S−K−dK)+g(S) = (S-K+dK)^+ - 2(S-K)^+ + (S-K-dK)^+ g(S)={0,S≤K−dKS−K+dK,K−dK≤S≤KdK−(S−K),K≤S≤K+dK0,S≥K+dKg(S) = \begin{cases} 0, & S \le K-dK \\ S-K+dK, & K-dK \le S \le K \\ dK - (S-K), & K \le S \le K+dK \\ 0, & S \ge K+dK \end{cases}
  • Треугольник (“шатёр”) с пиком g(K)=dKg(K) = dK в точке S=KS=K
  • Носитель: [K−dK, K+dK][K-dK,\ K+dK]
  • Площадь: ∫0∞g(S) dS=dK2\displaystyle\int_0^\infty g(S)\, dS = dK^2

Теорема о среднем

Если ff непрерывна на [a,b][a, b], то найдётся точка c∈[a,b]c \in [a, b], для которой

∫abf(x) dx=f(c) (b−a).\int_a^b f(x)\, dx = f(c)\,(b - a).

Интеграл равен значению функции в одной точке, умноженному на длину отрезка. f(c)f(c) и есть среднее значение ff на [a,b][a, b].

Если рядом стоит неотрицательный вес gg, теорема звучит так: найдётся c∈[a,b]c \in [a, b], для которой

∫abf(x) g(x) dx=f(c)∫abg(x) dx.\int_a^b f(x)\, g(x)\, dx = f(c) \int_a^b g(x)\, dx.

При g=1g = 1 это то же самое утверждение. Интеграл произведения равен значению ff в одной точке, умноженному на суммарный вес ∫g\int g.

На примере

Шатёр gg занимает весь отрезок [K−dK,K+dK][K - dK, K + dK], поэтому под интегралом плотность во многих точках. PP непрерывна и g≥0g \ge 0, значит по теореме о среднем с весом найдётся точка ξ\xi внутри шатра, ξ∈[K−dK,K+dK]\xi \in [K - dK, K + dK]:

Ctbf=P(ξ,T∣S0,0)∫0∞g(S) dS=P(ξ,T∣S0,0) dK2.C_{t}^{bf} = P(\xi, T \mid S_0, 0) \int_0^\infty g(S)\, dS = P(\xi, T \mid S_0, 0)\, dK^2.

Площадь ∫g=dK2\int g = dK^2 уже посчитана. Делим на dKdK:

CtbfdK=P(ξ,T∣S0,0) dK.\frac{C_{t}^{bf}}{dK} = P(\xi, T \mid S_0, 0)\, dK.

Справа плотность в некоторой точке шатра, умноженная на dKdK. Страйк KK появится, когда шатёр сожмётся.

На графике K=100K = 100, dK=20dK = 20. Шатёр стоит на [80,120][80, 120], площадь 202=40020^2 = 400. Цена такой бабочки равна 400⋅P(ξ)400 \cdot P(\xi), и ξ\xi лежит между 8080 и 120120. Деление цены на dK=20dK = 20 даёт 20⋅P(ξ)20 \cdot P(\xi).

Теорема не указывает, где именно лежит ξ\xi. Положение ξ\xi зависит от PP. Точки, где высота самого шатра равна средней dK/2dK/2, лежат на K±dK/2K \pm dK/2. Это точки функции gg.

Чтобы получить плотность в страйке, шатёр сужают. Любая точка отрезка [K−dK,K+dK][K - dK, K + dK] удалена от KK не больше чем на dKdK, поэтому при dK→0dK \to 0 получается ξ→K\xi \to K и P(ξ)→P(K)P(\xi) \to P(K).

Делить удобнее сразу на площадь:

CtbfdK2=Ct(K−dK)−2Ct(K)+Ct(K+dK)dK2=P(ξ,T∣S0,0).\frac{C_{t}^{bf}}{dK^2} = \frac{C_{t}(K - dK) - 2C_{t}(K) + C_{t}(K + dK)}{dK^2} = P(\xi, T \mid S_0, 0).

В пределе dK→0dK \to 0 вторая разность становится второй производной:

P(K,T∣S0,0)=∂2Ct∂K2.P(K, T \mid S_0, 0) = \frac{\partial^2 C_t}{\partial K^2}.

Итого: формула плотности

P(K,T∣S0,0)=∂2Ct∂K2P(K,T\mid S_0,0) = \frac{\partial^2 C_t}{\partial K^2}

Смысл: риск-нейтральная плотность цены базового актива в точке KK равна второй производной цены колла по страйку, взятой в этой же точке KK.

Как получили:

  1. Записали C(K)C(K) через интеграл выплаты по плотности PP.
  2. Разложили функцию выплат бабочки. Вторая разность C(K−dK)−2C(K)+C(K+dK)C(K-dK)-2C(K)+C(K+dK) равна интегралу этой треугольной выплаты g(S)g(S) по плотности P(S)P(S).
  3. По теореме о среднем этот интеграл =P(ξ,T)⋅dK2= P(\xi,T)\cdot dK^2, где ξ→K\xi\to K при dK→0dK\to0.
  4. Делим на dK2dK^2, устремляем dK→0dK\to0 → получаем именно значение плотности в точке KK.

Практический смысл: имея на рынке цены коллов по непрерывному спектру страйков, можно восстановить всю риск-нейтральную плотность будущей цены актива. Для этого берётся вторая производная цены опциона по страйку, а на практике вторая разность.