Оптимальная торговля при мгновенном и персистентном price impact, предсказуемых доходностях и многомасштабной стохастической волатильности

6.5/10

Patrick Chan · бывш. Program in Applied and Computational Mathematics, Princeton University · Ronnie Sircar, Iosif Zimbidis · Department of Operations Research & Financial Engineering, Princeton University · 24 июля 2025

Оригинал: Chan, P., Sircar, R. and Zimbidis, I. «Optimal Trading under Instantaneous and Persistent Price Impact, Predictable Returns and Multiscale Stochastic Volatility», v1, 23 июля 2025 — arxiv.org/abs/2507.17162 (PDF).

Рисунки воспроизведены из оригинальной публикации.

Ключевые слова: optimal trading, price impact, return predictability, multiscale stochastic volatility.

Аннотация

Мы рассматриваем задачу динамической оптимизации портфеля, включающую предсказуемые доходности, мгновенные транзакционные издержки, price impact и стохастическую волатильность, — расширяя классические результаты Gârleanu and Pedersen (2013), где волатильность постоянна. Построение оптимальной портфельной стратегии в этой общей постановке затруднено нелинейностью возникающих уравнений Гамильтона–Якоби–Беллмана (ГЯБ). Для решения предлагается многомасштабное разложение по волатильности, захватывающее динамику стохастической волатильности на разных временных масштабах: сингулярное возмущение для быстрого mean-reverting фактора волатильности и регулярное — для медленного. Также введена аппроксимация для малого price impact с демонстрацией численной точности. Формально выведены асимптотические приближения до второго порядка; симуляции Монте-Карло показывают, как включение этих поправок улучшает Profit and Loss (PnL) результирующей портфельной стратегии.

1. Введение

Tractable-каркас для динамической оптимизации портфеля при рыночных трениях предложили Gârleanu and Pedersen (2013, 2016): они вывели торговые стратегии в замкнутой форме с учётом предсказуемости доходностей, мгновенных и персистентных транзакционных издержек. Модель учитывает мгновенные издержки от механики лимитного стакана и персистентные — price impact в литературе, — отражающие длительное влияние сделок на цену актива. Опираясь на литературу по price impact, Webster (2023) даёт обзор современной теории с обширными ссылками на optimal execution и market impact и практическими соображениями по реализации.

В этой работе мы расширяем эту линию исследований, вводя стохастическую волатильность в задачу динамической торговли. Конкретно — многомасштабную модель стохастической волатильности с медленным и быстрым факторами по Fouque et al. (2011). Работа также даёт дополнительные инсайты в структуру и решение системы уравнений, возникающей даже при постоянной волатильности в наиболее общей постановке с персистентными и мгновенными издержками. Gârleanu and Pedersen (2013) не выводят явные решения в этом случае, лишь утверждая существование единственного решения при подходящих условиях. Мы показываем, что при относительно малом price impact система допускает численно tractable аппроксимации. Сходный асимптотический подход использован в Ekren and Muhle-Karbe (2019), где оба вида трения — мгновенное и персистентное — малы, что приводит к решению через алгебраическое уравнение Riccati.

В отличие от случая постоянной волатильности у Gârleanu and Pedersen (2013), при стохастической волатильности уравнение ГЯБ становится нелинейным. Однако метод асимптотической аппроксимации Fouque et al. (2011) показывает, что задача остаётся tractable как возмущение вокруг полностью решаемого решения при постоянной волатильности. Ниже — ключевые аспекты построения нашей динамической оптимизационной постановки.

Предсказуемость доходностей. Одна из главных задач трейдера — прогнозировать движение актива и извлекать прибыль из прогнозов. Прогнозы не ограничиваются простыми краткосрочными проекциями, а часто включают сложные модели с mean-reversion и momentum.

Транзакционные издержки. При исполнении стратегий через изменение позиций возникает трение: bid-ask спред, комиссии. Инвесторы хотели бы следовать теоретически «оптимальному» портфелю, но издержки вызывают отклонение от портфеля, вычисленного без трений. В литературе их называют мгновенными: после сделки они не влияют на цену актива.

Price impact. Ключевой фактор при выборе стратегии: влияние торговых решений на цену. Первые модели появились в 1980-х; значимость быстро признали в академии и индустрии. В литературе price impact иногда называют персистентными транзакционными издержками.

Стохастическая волатильность. Играет ключевую роль в ценообразовании активов, объясняя volatility smile и другие эмпирические факты. Chacko and Viceira (2005) подчёркивают важность многофакторного моделирования — мы опираемся на эту идею, поскольку анализ high-frequency и low-frequency данных даёт существенно разные результаты.

Мы объединяем все перечисленные элементы в одной постановке динамической оптимизации портфеля. Задача остаётся столь же tractable, как у Gârleanu and Pedersen, хотя модель волатильности расширена от простого постоянного случая до двухфакторной многомасштабной стохастической волатильности. Явные поправочные члены вызывают отклонение оптимальной скорости торговли от случая постоянной волатильности; эти поправки существенно улучшают PnL портфеля, что мы количественно оцениваем симуляциями Монте-Карlo в численном разделе. Вклад работы также в том, что многомасштабные модели стохастической волатильности ранее применялись в безтрictionных задачах вроде оптимизации портфеля Merton (Fouque et al., 2017b), а здесь — в постановке с рыночными трениями.

Кроме того, для широкого класса practically relevant моделей стохастической волатильности (Heston, exponential OU, 3/2-model) поправки к стратегиям постоянной волатильности выводятся в замкнутой форме и соответствуют рациональному экономическому поведению.

В частности, при быстрой динамике трейдерам рекомендуется снижать плечо, когда текущая волатильность превышает долгосрочное среднее, независимо от корреляции доходность–волатильность. При медленной динамике корреляция играет более существенную роль: при положительной корреляции инвестор оптимально снижает скорость торговли, ожидая, что более высокая ожидаемая доходность сопровождается более высокой волатильностью. Мы также показываем, что влияние медленной стохастической волатильности сильнее, чем быстрой, в контексте нашей бесконечно-горизонтной задачи optimal trading: ведущий член разложения по быстрой волатильности тождественно обращается в нуль, и для захвата основных эффектов нужно вычислять второй порядок.

Более свежие работы фокусируются на динамике price impact. Carmona and Webster (2019) показали локальную вогнутость propagator в sophisticated моделях для общего стакана, подтвердив численные результаты Bouchaud et al. (2009). Muhle-Karbe et al. (2024b) демонстрируют, что случайность уровня impact эффективно захватывает эффект вогнутости, давая простую линейную time-varying модель. Brokmann et al. (2025) предлагают аппроксимировать нелинейный price impact линейными торговыми стратегиями.

Недавние исследования включают machine и deep learning в оптимизацию портфеля: Buehler et al. (2019) — hedging path-dependent payoffs нейросетями при торговом трении; Han et al. (2018) — deep learning для high-dimensional forward-backward SDE в нелинейном стохастическом управлении; Heaton et al. (2017) — deep portfolios для сложной структуры взаимодействий активов. Deep learning и robust control применяются к model uncertainty: Pham et al. (2022) — robust dynamic mean-variance при drift ambiguity; Lin et al. (2022) — multiagent deep RL для portfolio management; Cai and Yu (2025) — Bayesian learning под minimax rule. Muhle-Karbe et al. (2024a) сравнивают решения Gârleanu and Pedersen с neural network методами: при оптимальной калибровке линейных impact-параметров к underlying nonlinear impact модели политика конкурентна с NN, кроме крайне illiquid рынков. Abi Jaber and Neuman (2022) и Neuman and Voß (2022) включают predictive finite-variation signals в propagator-модели с transient и temporary impact.

Асимптотическая аппроксимация при многомасштабном моделировании стохастической волатильности применялась в optimal investment и derivative pricing. По Fouque et al. (2011), singular и regular perturbation methods эффективно аппроксимируют линейные задачи ценообразования опционов. Наш анализ сходен с этой линией, но отличается в ключевых аспектах: мы явно включаем транзакционные издержки и price impact; рассматриваем mean-variance оптимизацию на бесконечном торговом горизонте; даём явное вычисление членов до второго порядка поправки по быстрой стохастической волатильности.

Обзор. В таблице 1 суммированы модели динамической торговли в литературе. Type — continuous/discrete; Cont./Disc.; (g)BM — (geometric) Brownian motion; SV — stochastic volatility; pred. — predictability; Lrr — linear rebalancing rules.

Таблица 1. Модели динамической торговли: постановки, модели и подходы к решению.
РаботаTypePrice dynamicsTrading frictionPrice ImpactObjectiveSolution
Merton (1971)Cont.gBMNoneNoUtilityAnalytic
Liu and Loewenstein (2002)Cont.gBMProportionalNoUtilityAnalytic
Kraft (2005)Cont.gBM+HestonNoneNoUtilityAnalytic
Chacko and Viceira (2005)Cont.gBM+3/2NoneNoUtilityAnalytic
Moallemi and Saglam (2012)Disc.BM+pred.QuadraticNoMean-Variancelrr
Obizhaeva and Wang (2013)Cont.BMNoneLinearExecution CostAnalytic
Bichuch and Sircar (2014)Cont.gBM+SVProportionalNoUtilityAsymptotic
Fouque et al. (2014)Cont.gBM+pred.NoneNoUtilityAsymptotic
Fouque et al. (2017b)Cont.gBM+SVNoneNoUtilityAsymptotic
Gârleanu and Pedersen (2013, 2016)Cont./Disc.BM+pred.QuadraticLinearMean-VarianceAnalytic
Passerini and Vazquez (2015)Cont.BM+pred.Linear-quadraticNoMean-VarianceApproximate
Curato et al. (2017)Cont.BMLocally concaveLocally concaveExecution CostNumerical
Collin-Dufresne et al. (2020)Cont./Disc.Markovian (switching)QuadraticNoneExecution CostAnalytic
Neuman and Voß (2022)Cont.BM+pred.NoneLinearExecution CostAnalytic
Abi Jaber and Neuman (2022)Cont.BM+pred.QuadraticGeneral propagatorExecution CostAnalytic
Muhle-Karbe et al. (2023)Cont.BM+partially pred.QuadraticNoneUtilityAnalytic
Настоящая работаCont.BM+pred.+SVQuadraticLinearMean-VarianceAsymptotic

1.1. Структура работы

В разделе 2 представлена непрерывно-временная модель с многомасштабной стохастической волатильностью; выведено уравнение ГЯБ и аналитическое решение при постоянной волатильности. В подразделе 2.5 — случай постоянной волатильности Gârleanu and Pedersen; в 2.6 — аппроксимация для относительно малого price impact. Далее, для управляемости изложения, анализируем два фактора отдельно. Раздел 3 — быстрый фактор волатильности (singular perturbation для ГЯБ). Раздел 4 — медленный фактор (regular perturbation). Раздел 5 — объединение быстрых и медленных результатов при многомасштабной стохастической волатильности. Раздел 6 — численные примеры. Раздел 7 — заключение и направления расширения.

2. Модель рыночной торговли

2.1. Рыночное трение

При исполнении рыночных заявок возникает издержка от влияния торговых решений на ликвидность и, следовательно, на цену. В литературе эта издержка обычно моделируется квадратичной по скорости торговли (Gârleanu and Pedersen, 2016). Обозначим $q_t$ — позицию (число акций) в момент $t$:

\[ dq_t = u_t\,dt, \]

где $u_t$ — скорость торговли. Издержка вводится при допущении, что после исполнения ликвидность возвращается к исходному уровню и не оказывает lasting impact на цену. Её называют мгновенной (instantaneous) или temporary transaction cost:

\[ C_{\mathrm{instant}}(u_t) = \frac{K}{2}u_t^2, \]

где $K > 0$ — константа уровня транзакционных издержек. При $u_t > 0$ цена исполнения превышает unaffected price, разница пропорциональна скорости торговли.

Кроме temporary costs при исполнении, учитываем персистентную (transient) транзакционную издержку — price impact: покупка или продажа сдвигает цену вверх или вниз, эффект сохраняется после сделки, поскольку ликвидность может не восстанавливаться мгновенно.

Вводим unperturbed price $P_t$ и perturbed price $S_t$ (Webster, 2023):

\[ S_t = P_t + l_t, \]

где $(l_t)_{t\ge 0}$ — накопленный price impact. По Gârleanu and Pedersen (2016) и Webster (2023), impact эволюционирует линейно по изменению позиции и экспоненциально затухает:

\[ dl_t = -\beta l_t\,dt + \lambda\,dq_t = (\lambda u_t - \beta l_t)\,dt, \]

где $\beta, \lambda > 0$. Webster (2023) называет $\beta$ decay rate impact и $\lambda$ — мерой рыночной ликвидности. Для простоты параметры постоянны во времени.

2.2. Функция богатства трейдера

Функция богатства состоит из трёх компонент: позиция в акции, штраф за риск и транзакционные издержки. Процентная ставка нулевая. Трейдер следует self-financing стратегии для максимизации финансового gain от позиции:

\[ \Pi_s = \int_0^s q_t\,dS_t. \]

Из (3): $q_t\,dS_t = q_t(dP_t + dl_t) = q_t\,dP_t - (\beta q_t l_t - \lambda q_t u_t)\,dt$.

Далее — член risk aversion: трейдер максимизирует прибыль и минимизирует risk exposure. Динамика цены одного актива с одним return predictor $x_t$:

\[ dP_t = \alpha_t\,dt + \sigma(Y_t, Z_t)\,dB_t, \]

где $B_t$ — стандартное броуновское движение, $Y_t, Z_t$ — факторы стохастической волатильности. Без потери общности drift $\alpha_t$ (alpha signal) разлагаем на константу $\bar\alpha$ и zero-mean intraday компонент $x_t$: $\alpha_t := \bar\alpha + x_t$. Сигнал $x_t$ — процесс Ornstein–Uhlenbeck:

\[ dx_t = -\kappa x_t\,dt + \sqrt{\eta}\,dW_t^{(0)}, \]

где $W_t^{(0)}$ независимо от $B_t$. Risk aversion — штраф в objective:

\[ -\frac{\gamma}{2}\sigma^2(Y_t, Z_t) q_t^2\,dt, \]

где $\gamma$ — параметр risk aversion. С учётом temporary cost из (2) трейдер выбирает $(u_t)_{t\ge 0}$ для максимизации discounted expected excess returns с штрафом за риск и издержки:

\[ \max_{u}\ \mathbb{E}\left[\int_0^\infty e^{-\rho(s-t)}\left(q_s x_s - \beta q_s l_s - \lambda q_s u_s - \frac{\gamma}{2}\sigma^2(Y_s,Z_s)q_s^2 - \frac{K}{2}u_s^2\right)ds\right], \]

где $\rho > 0$ — discount factor. При совместной марковости $(x_t, Y_t, Z_t)$ определяем функцию ценности

\[ v(q,l,x,y,z) = \sup_u \mathbb{E}_{q,l,x,y,z}\left[\int_0^\infty e^{-\rho t}\left(q_t x_t - \beta q_t l_t - \lambda q_t u_t - \frac{\gamma}{2}\sigma^2(Y_t,Z_t)q_t^2 - \frac{K}{2}u_t^2\right)dt\right], \]

где $\mathbb{E}_{q,l,x,y,z}[\cdot] = \mathbb{E}[\cdot \mid q_0=q, l_0=l, x_0=x, y_0=y, z_0=z]$, $(q,l,x)\in\mathbb{R}^3$.

2.3. Многомасштабная стохастическая волатильность

Работаем в каркасе multiscale stochastic volatility (Fouque et al., 2011, 2017b): один быстрый и один медленный фактор волатильности. Волатильность — функция $\sigma(Y_t, Z_t)$; факторы:

\[ \frac{1}{\varepsilon}\,dY_t = b(Y_t)\,dt + \frac{1}{\sqrt{\varepsilon}}\,a(Y_t)\,dW_t^{(1)}, \] \[ dZ_t = \delta c(Y_t)\,dt + \delta g(Y_t)\,dW_t^{(2)}, \]

где $W_t^{(0)}, W_t^{(1)}, W_t^{(2)}$ — стандартные броуновские движения на filtered probability space с мгновенными корреляциями

\[ d\langle W^{(0)}, W^{(i)}\rangle_t = \rho_i\,dt,\quad i=1,2,\qquad d\langle W^{(1)}, W^{(2)}\rangle_t = \rho_{12}\,dt, \]

при $|ρ_i| < 1$, $|ρ_{12}| < 1$ и положительной определённости ковариационной матрицы трёх BM. Brownian motions независимы от $B_t$. Модель задаётся коэффициентами $\bar\alpha, \kappa, \eta, \sigma, a, b, c, g$; малые $\varepsilon, \delta$ характеризуют быструю и медленную вариацию факторов $Y$ и $Z$.

Предполагаем $Y_t \stackrel{d}{=} Y_{t/\varepsilon}^{(1)}$, где $Y^{(1)}$ — ergodic процесс с инвариантным распределением $\Phi$, не зависящим от $\varepsilon$. Аналогично $Z_t \stackrel{d}{=} Z_{\delta t}^{(1)}$ — diffusion с drift $c$ и diffusion $g$; для slow-scale asymptotics при $\delta\downarrow 0$ ergodicity $Z^{(1)}$ не требуется.

2.4. Уравнение Гамильтона–Якоби–Беллмана

Без потери общности полагаем $\bar\alpha = 0$. Supremum — по admissible стратегиям: $F_t$-progressively measurable, square-integrable, с единственным strong solution для (6), (8), (9). Dynamic programming даёт уравнение ГЯБ:

\[ \rho v = \sup_u\left\{\frac{1}{\varepsilon}\mathcal{L}_0 v + \frac{1}{\sqrt{\varepsilon}}\mathcal{L}_1 v + \delta\mathcal{M}_1 v + \delta\mathcal{M}_2 v + \frac{1}{\delta}\mathcal{M}_3 v + \mathcal{L}_2(\sigma(y,z))v + \sup_u\left[-\lambda qu - \frac{K}{2}u^2 + uv_q + \lambda uv_l\right]\right\} = 0, \]

с операторами

\[ \mathcal{L}_2(\sigma^2(y,z)) = \frac{\eta}{2}\frac{\partial^2}{\partial x^2} - \kappa x\frac{\partial}{\partial x} - \beta l\frac{\partial}{\partial l} + qx - \beta ql - \frac{\gamma}{2}\sigma^2 q^2, \] \[ \mathcal{L}_0 = \frac{a(y)^2}{2}\frac{\partial^2}{\partial y^2} + b(y)\frac{\partial}{\partial y},\quad \mathcal{L}_1 = \sqrt{\eta}\,\rho_1 a(y)\frac{\partial^2}{\partial x\partial y}, \] \[ \mathcal{M}_1 = \sqrt{\eta}\,\rho_2 g(z)\frac{\partial^2}{\partial x\partial z},\quad \mathcal{M}_2 = \frac{g(z)^2}{2}\frac{\partial^2}{\partial z^2} + c(z)\frac{\partial}{\partial z},\quad \mathcal{M}_3 = \rho_{12} a(y)g(z)\frac{\partial^2}{\partial y\partial z}. \]

$\mathcal{L}_0$ и $\mathcal{M}_2$ — infinitesimal generators процессов $Y^{(1)}$ и $Z^{(1)}$. Максимизация по $u$ квадратична; оптимальная скорость торговли

\[ u^* = \frac{1}{K}(\lambda q + v_q + \lambda v_l), \]

подстановка даёт нелинейное PDE (10), не решаемое напрямую аналитически или численно. Подход — perturbation вокруг случая постоянной волатильности (Gârleanu and Pedersen, 2013).

2.5. Решение при постоянной волатильности

При постоянной $\sigma$ функция ценности $v(q,l,x)$ не зависит от $y, z$. Уравнение ГЯБ упрощается:

\[ \rho v = \mathcal{L}_2(\sigma)v + \frac{1}{2K}(\lambda q + v_q + \lambda v_l)^2. \]

По Gârleanu and Pedersen (2013), ansatz:

\[ v(q,l,x) = -\tfrac{1}{2}A_{qq}q^2 - \tfrac{1}{2}A_{ll}l^2 - \tfrac{1}{2}A_{xx}x^2 - A_{ql}ql - A_{qx}qx - A_{xl}xl - A_0. \]

Подстановка в (11) даёт алгебраическую систему для коэффициентов $A$ (уравнения (13)–(19) оригинала). Система нелинейна и сложнее случая без price impact ($\lambda=\beta=0$). В общем случае замкнутого решения нет, но при достаточно большом $\gamma$ существует единственное положительное решение (Gârleanu and Pedersen, 2013). При $\lambda=\beta=0$ — замкнутая форма для $A_{qq}, A_{xx}, A_{qx}, A_0$, остальные нулевые.

Решая подсистему (13), (14), (16) для $A_{qq}, A_{ql}, A_{ll}$, получаем квартическое уравнение для $A_{ql}$; корни дают кандидаты, остальные уравнения линейны и решаются последовательно.

Оптимальная торговая стратегия. Оптимальная скорость:

\[ u^*(q,x,l) = \frac{1}{K}\left[(\lambda - A_{qq} + \lambda A_{ql})q + (A_{ql} + \lambda A_{ll})l + (A_{qx} + \lambda A_{xl})x\right]. \]

В терминах target portfolio и tracking speed (Gârleanu and Pedersen, 2013, 2016):

\[ u^* = r_c(\sigma^2)\left(\mathrm{aim}_t^c(\sigma^2) - q\right), \]

где

\[ r_c(\sigma^2) = \frac{1}{K}(A_{qq} - \lambda - \lambda A_{ql}),\qquad \mathrm{aim}_t^c(\sigma^2) = \frac{(A_{ql}+\lambda A_{ll})l_t + (A_{qx}+\lambda A_{xl})x_t}{A_{qq}-\lambda-\lambda A_{ql}}, \]

при $r_c \ne 0$. Трейдер отслеживает time-varying target $\mathrm{aim}_t^c$ со скоростью $r_c$. Зависимость от $\sigma$ будет расширена в случае стохастической волатильности.

2.6. Аппроксимация малого price impact

Система (13)–(19) для leading order нелинейна. Предлагаем аналитический метод упрощения. Эмпирически (Webster, 2023) price impact мал относительно размера позиции. Вводим scaling parameter $\theta$ в dynamics impact:

\[ dl_t = -\theta\beta l_t\,dt + \theta\lambda\,dq_t, \]

$\theta$ мало. Обозначим $A^\theta$ — решение системы с заменами $\beta \mapsto \theta\beta$, $\lambda \mapsto \theta\lambda$. Разложение:

\[ A_{qq}^\theta = A_{qq}^{(0)} + \theta A_{qq}^{(1)} + \theta^2 A_{qq}^{(2)} + \cdots \]

Подстановка даёт полином по $\theta$; коэффициенты при каждой степени обращаются в нуль, что позволяет решать рекурсивно.

2.6.1. Член нулевого порядка $\theta^0$

При $\theta=0$ — результат Gârleanu and Pedersen (2013) без price impact:

\[ \frac{1}{2K}(A_{qq}^{(0)})^2 + \rho A_{qq}^{(0)} - \frac{\gamma\sigma^2}{2} = 0,\quad \frac{1}{K}A_{qq}^{(0)}A_{qx}^{(0)} - \kappa A_{qx}^{(0)} - 1 - \rho A_{qx}^{(0)} = 0, \] \[ \frac{1}{2K}(A_{qx}^{(0)})^2 - \kappa A_{xx}^{(0)} - \frac{\rho}{2}A_{xx}^{(0)} = 0,\quad \frac{\eta}{2}A_{xx}^{(0)} - \rho A_0^{(0)} = 0, \]

с $A_{ql}^{(0)} = A_{ll}^{(0)} = A_{xl}^{(0)} = 0$. Решение:

\[ A_{qq}^{(0)} = \frac{K}{\gamma\sigma^2}\left(\sqrt{\rho^2 + \frac{4\gamma\sigma^2}{K}} - \rho\right),\quad A_{qx}^{(0)} = \frac{1}{\kappa+\rho + A_{qq}^{(0)}/K},\quad A_{xx}^{(0)} = \frac{(A_{qx}^{(0)})^2 K}{2\kappa+\rho},\quad A_0^{(0)} = \frac{\eta}{2\rho}A_{xx}^{(0)}. \]

2.6.2. Член первого порядка $\theta^1$

Линейная система для $A^{(1)}$; ключевые результаты:

\[ A_{qq}^{(1)} = \frac{2\lambda A_{qq}^{(0)}}{2A_{qq}^{(0)} + K\rho},\quad A_{ql}^{(1)} = -\frac{\beta K}{2A_{qq}^{(0)} + K\rho},\quad A_{ll}^{(1)} = 0, \] \[ A_{qx}^{(1)} = A_{qx}^{(0)}\frac{A_{qx}^{(0)}A_{ql}^{(1)}}{\rho K + A_{qq}^{(0)}},\quad A_{xl}^{(1)} = \frac{A_{qx}^{(0)}A_{ql}^{(1)}}{K(\kappa+\rho)},\quad A_{xx}^{(1)} = \frac{2A_{qx}^{(0)}A_{qx}^{(1)}}{K(2\kappa+\rho)},\quad A_0^{(1)} = \frac{\eta}{2\rho}A_{xx}^{(1)}. \]

Член второго порядка — в приложении A. Точность первого порядка демонстрируется численно в разделе 6.

3. Один быстрый фактор

Сначала анализируем задачу при одном быстром факторе волатильности $\sigma = \sigma(Y_t)$. PDE функции ценности $v(q,l,x,y)$:

\[ \rho v = \mathcal{L}_2(\sigma^2(y))v + \frac{1}{2K}(\lambda q + v_q + \lambda v_l)^2 + \frac{1}{\varepsilon}\mathcal{L}_0 v + \frac{1}{\sqrt{\varepsilon}}\mathcal{L}_1 v. \]

Singular perturbation при $\varepsilon\downarrow 0$; асимптотическое разложение:

\[ v^\varepsilon(q,l,x,y) = v^{(0)}(q,l,x) + \sqrt{\varepsilon}\,v^{(1)}(q,l,x) + \varepsilon\,v^{(2)}(q,l,x,y) + \varepsilon^{3/2}v^{(3)}(q,l,x,y) + \cdots \]

3.1. Член нулевого порядка $v^{(0)}$

Порядок $\varepsilon^{-1}$: $\mathcal{L}_0 v^{(0)} = 0$ — $v^{(0)}$ не зависит от $y$. Порядок $\varepsilon^{-1/2}$: $v^{(1)}$ также не зависит от $y$. Порядок $\varepsilon^0$:

\[ \mathcal{L}_0 v^{(2)} + \mathcal{L}_2(\sigma(y))v^{(0)} + \frac{1}{2K}(\lambda q + v_q^{(0)} + \lambda v_l^{(0)})^2 = 0. \]

Условие разрешимости (Fredholm alternative):

\[ \mathcal{L}_2(\langle\sigma^2\rangle)v^{(0)} + \frac{1}{2K}(\lambda q + v_q^{(0)} + \lambda v_l^{(0)})^2 = 0, \]

где $\langle h\rangle = \int h(y)\,\Phi(dy)$ — expectation под инвариантным распределением $Y^{(1)}$. Это PDE случая постоянной волатильности с $\sigma$ заменённой на $\langle\sigma\rangle$; $v^{(0)}$ — квадратичная форма как в (12).

Вычитая (24) из (23):

\[ \mathcal{L}_0 v^{(2)} = \frac{\gamma}{2}\left(\sigma^2(y) - \langle\sigma^2\rangle\right)q^2 =: h(y). \]

Poisson equation с решением

\[ v^{(2)}(q,y) = -\int_0^\infty P_t h(y)\,dt + C(q,l,x) = -\tfrac{1}{2}q^2\varphi(y),\quad \varphi(y) := \frac{\gamma}{2}\int_0^\infty P_t\left[\sigma^2(y)-\langle\sigma^2\rangle\right]dt. \]

Пример 1. $\sigma^2(y)=y$, CIR: $b(y)=\chi(\mu-y)$, $a(y)=\psi\sqrt{y}$. Тогда $v^{(2)} = -\frac{\gamma}{2\chi}(y-\mu)q^2$.

Пример 2. Exponential OU (Masoliver and Perelló, 2006): $\sigma(y)=me^y$, $b(y)=-\theta y$, $a(y)=\hat\sigma$. Интеграл $\int_0^\infty P_t h(y)\,dt$ не имеет простой замкнутой формы; для малых флуктуаций вокруг среднего — аппроксимация через incomplete gamma functions.

3.2. Член первого порядка $v^{(1)}$

Порядок $\sqrt{\varepsilon}$ даёт Poisson equation для $v^{(3)}$; условие разрешимости — линейное однородное PDE для $v^{(1)}$, допускающее тривиальное решение $v^{(1)}\equiv 0$. Следовательно, $v^{(3)}$ также не зависит от $y$.

3.3. Член второго порядка $v^{(2)}$

На порядке $\varepsilon$, с $v^{(1)}=0$, получаем линейное однородное уравнение для $C$; без потери общности $C\equiv 0$. Ведущая поправка к функции ценности — (26).

3.4. Оптимальная торговая стратегия

С поправкой управление:

\[ u^*(q,l,x,y) = \frac{1}{K}\left[(\lambda-A_{qq}+\lambda A_{ql})q + (A_{ql}+\lambda A_{ll})l + (A_{qx}+\lambda A_{xl})x - \varepsilon q\,\varphi(y)\right]. \]

В терминах tracking:

\[ u_t^* = r_f\left(\mathrm{aim}_t^f - q\right),\quad r_f := r_c(\langle\sigma^2\rangle) + \frac{1}{K}\varepsilon\varphi(Y_t), \] \[ \mathrm{aim}_t^f := \mathrm{aim}_t^c(\langle\sigma^2\rangle)\left(1 - \varepsilon\frac{\varphi(Y_t)}{A_{qq}-\lambda-\lambda A_{ql}} + \cdots\right). \]

Главная поправка к оптимальной скорости: $u^{(2)}(q,y) = -\frac{q}{K}\varphi(y)$.

Пример 1 (продолжение). CIR: $u^{(2)}(q,y) = -\frac{\gamma}{\theta K}q(y-\mu)$ — снижать exposure при волатильности выше long-term average. Скорректированные параметры: $r_f = r_c + \varepsilon\frac{\gamma}{\theta K}(y-\mu)$; $\mathrm{aim}_t^f = \mathrm{aim}_t^c\left(1 - \varepsilon\frac{\gamma}{\theta(A_{qq}-\lambda-\lambda A_{ql})}(y-\mu)\right)$. Скорость $r_f$ растёт с $y$; target portfolio уменьшается при $y$ выше среднего. Эффекты ослабляются большим $K$ или меньшим $\gamma$.

4. Один медленный фактор

Рассматриваем случай одного slow-scale фактора. PDE $v(q,l,x,z)$:

\[ \rho v = \mathcal{L}_2(\sigma^2(z))v + \frac{1}{2K}(\lambda q + v_q^\delta + \lambda v_l^\delta)^2 + \sqrt{\delta}\,\mathcal{M}_2 v^\delta + \sqrt{\delta}\,\mathcal{M}_1 v^\delta. \]

Разложение:

\[ v^\delta(q,l,x,z) = v^{(0)}(q,l,x,z) + \sqrt{\delta}\,v^{(1)}(q,l,x,z) + \delta\,v^{(2)}(q,l,x,z) + \cdots \]

4.1. Член нулевого порядка $\delta^0$

При $\delta=0$ — PDE как (11) с явной $z$-зависимостью $\sigma$. Решение — параметризованный случай постоянной волатильности:

\[ v^{(0)}(q,l,x,z) = -\tfrac{1}{2}A_{qq}(z)q^2 - \tfrac{1}{2}A_{ll}(z)l^2 - \tfrac{1}{2}A_{xx}(z)x^2 - A_{ql}(z)ql - A_{qx}(z)qx - A_{xl}(z)xl - A_0(z). \]

4.2. Член первого порядка

На порядке $\sqrt{\delta}$:

\[ -\rho v^{(1)} - \beta l v_l^{(1)} - \kappa x v_x^{(1)} + \frac{\eta}{2}v_{xx}^{(1)} + (v_q^{(1)}+\lambda v_l^{(1)})(\lambda q + v_q^{(0)}+\lambda v_l^{(0)}) + \mathcal{M}_1 v^{(0)} = 0. \]

Ansatz $v^{(1)}(q,l,x,z) = B_q(z)q + B_l(z)l + B_x(z)x$. Приравнивание коэффициентов при $q,l,x$ даёт линейную систему для $B_q, B_l, B_x$ с

\[ \Gamma(z) = \beta\lambda - \rho A_{qq}(z) - \beta A_{qq}(z) + \lambda\rho - K\rho^2 + \beta\lambda A_{ql}(z) + 2\lambda\rho A_{ql}(z) - K\beta\rho + \lambda^2\rho A_{ll}(z). \]

При $\lambda=\beta=0$: $\Gamma(z)>0$ для $\sigma(z)>0$. При малом price impact предполагаем $\Gamma(z)>0$. Решение $B_q(z), B_l(z), B_x(z)$ — формулы (28) оригинала (длинные выражения через производные $A$ и $g(z)$).

4.3. Оптимальная торговая стратегия

С первой поправкой:

\[ u^* = \frac{1}{K}\left[(\lambda-A_{qq}(z)+\lambda A_{ql}(z))q + (A_{ql}(z)+\lambda A_{ll}(z))l + (A_{qx}(z)+\lambda A_{xl}(z))x + \sqrt{\delta}(B_q(z)+\lambda B_l(z))\right]. \]

В терминах tracking: $u_t^* = r_s(Z_t)(\mathrm{aim}_t^s - q)$ с

\[ r_s := r_c(\sigma^2(Z_t)),\qquad \mathrm{aim}_t^s := \mathrm{aim}_t^c(\sigma^2(Z_t)) + \sqrt{\delta}\,\frac{B_q(Z_t)+\lambda B_l(Z_t)}{A_{qq}-\lambda-\lambda A_{ql}}. \]

При малом price impact $B_q$ и $\rho_2$ имеют противоположные знаки при монотонно возрастающей $\sigma$ по $Z_t$. При положительной корреляции $\rho_2$ между волатильностью и доходностями инвестор оптимально снижает скорость торговли; при отрицательной — увеличивает. Производные для $B$-членов — приложение B; второй порядок для одного медленного фактора — приложение C.

5. Многомасштабная двухфакторная модель волатильности

Обогащённая модель с быстрым и медленным факторами. PDE — (10) с $v^{\varepsilon,\delta}$. Разложение по обоим масштабам:

\[ v^{\varepsilon,\delta} = v^{\varepsilon,0} + \sqrt{\delta}\,v^{\varepsilon,1} + \delta\,v^{\varepsilon,2} + \cdots;\quad v^{\varepsilon,0} = v^{(0)} + \sqrt{\varepsilon}\,v^{(1,0)} + \varepsilon\,v^{(2,0)} + \cdots \]

$v^{(0)}$ — решение случая постоянной волатильности с $\bar\sigma^2(z) := \langle\sigma^2(\cdot,z)\rangle$. Поправки из раздела 3:

\[ v^{(1,0)} = 0,\qquad v^{(2,0)}(q,z,y) = -\tfrac{1}{2}q^2\varphi(y,z),\quad \varphi(y,z) := \frac{\gamma}{2}\int_0^\infty P_t\left[\sigma^2(y,z)-\bar\sigma^2(z)\right]dt. \]

Первая поправка по медленному масштабу на порядке $\sqrt{\delta}$: $v^{(0,1)}(q,l,x,z) = B_q(z)q + B_l(z)l + B_x(z)x$ с коэффициентами из (28), где $A$ определены при $\bar\sigma^2(z)$.

Оптимальная торговая стратегия. С поправками до порядков $\varepsilon$ и $\delta$:

\[ u^{\varepsilon,\delta} = \frac{1}{K}\left[(\lambda-A_{qq}+\lambda A_{ql})q + (A_{ql}+\lambda A_{ll})l + (A_{qx}+\lambda A_{xl})x + \sqrt{\delta}(B_q+\lambda B_l) - \varepsilon q\,\varphi(y,z)\right]. \]

Или $u^* = r_m(\mathrm{aim}_t^m - q)$ с

\[ r_m := r_c(\bar\sigma^2(Z_t)) + \frac{1}{K}\varepsilon\varphi(Y_t,Z_t), \] \[ \mathrm{aim}_t^m := \mathrm{aim}_t^c(\bar\sigma^2(Z_t))\left(1 - \varepsilon\frac{\varphi(Y_t,Z_t)}{A_{qq}(Z_t)-\lambda-\lambda A_{ql}(Z_t)} + \sqrt{\delta}\,\frac{B_q(Z_t)+\lambda B_l(Z_t)}{A_{qq}(Z_t)-\lambda-\lambda A_{ql}(Z_t)+\varepsilon\varphi(Y_t,Z_t)}\right). \]

Замечание 5.1. И target, и tracking speed затронуты multiscale stochastic volatility в представлении выше.

6. Примеры и численные решения

Численные примеры демонстрируют эффективность вычисления асимптотической аппроксимации в practically relevant моделях и улучшение PnL симуляциями Мonте-Carlo. Цель — оценить различие PnL при исполнении с постоянной скоростью и с corrected optimal rate.

6.1. Численная аппроксимация price impact

Для Monte Carlo в slow-scale режиме используем метод приложения B. Полное решение нелинейной системы на каждом шаге и траектории computationally prohibitive; применяем first-order аппроксимацию из §2.6 для всех коэффициентов $A$. Сравниваем с exact solutions и вычисляем normalized error по широкому диапазону $\sigma$ (рисунок 1).

Рисунок 1
Рисунок 1. Нормированная ошибка каждого коэффициента при first-order аппроксимации малого price impact (§2.6). Параметры: $\rho=0.2$, $\gamma=5$, $K=1$, $\lambda=\beta=1$, $\theta=0.1$, $\kappa=1$, $\eta=1$.

Также вычисляем normalized error для четырёх главных производных коэффициентов, необходимых для $B_q$ и $B_l$ (рисунок 2). На обоих рисунках ошибки порядка $O(\theta^2)$, что подтверждает использование first-order аппроксимации в Monte Carlo. Normalized error быстро убывает с ростом $\sigma$ для всех коэффициентов и их производных.

Рисунок 2
Рисунок 2. Нормированная ошибка ключевых производных для вычисления $B$-членов. Параметры как на рисунке 1.

6.2. Пример быстрого фактора

Динамика цены и волатильности — модель Heston: $\sigma(y)=\sqrt{y}$, $b(y)=\chi(\mu-y)$, $a(y)=\psi\sqrt{y}$ при Feller condition $\psi^2 < 2\mu$. Оптимальная стратегия с second-order поправкой (§3):

\[ u^*(q,l,x,y) = \frac{1}{K}\left((\lambda-A_{qq}+\lambda A_{ql})q + (A_{ql}+\lambda A_{ll})l + (A_{qx}+\lambda A_{xl})x\right) - \frac{1}{K}\frac{\gamma}{\theta}\varepsilon(y-\mu)q. \]

На рисунке 3 — trading speed и aim portfolio как функции $K$. Сплошная линия ($y=\mu$) — скорость без поправки $r_c$. При волатильности выше $\mu$ поправка увеличивает adjusted speed $r_f$; ниже $\mu$ — уменьшает. Для target portfolio $\mathrm{aim}^f$ — противоположный эффект.

Рисунок 3
Рисунок 3. Оптимальная скорость торговли и aim portfolio при разных $K$, с second-order поправкой и без. Параметры: $\gamma=5$, $\lambda=\beta=0.1$, $\kappa=\eta=1$, $\rho_1=0.5$, $K=1$, $x=1$, $\mu=0.2$, $\chi=1$, $\psi=0.25$, $\varepsilon=0.25$.

Горизонт $T=2$ года; Monte Carlo сравнивает PnL (i) constant-rate optimal strategy и (ii) corrected strategy с fast mean-reverting SV (рисунок 4). Second-order поправка существенно улучшает performance: средний прирост 53.27 bps относительно constant volatility strategy при $Y_0=0.4$. Дисперсия PnL ниже в corrected case: 0.6286 против 0.6515.

Рисунок 4
Рисунок 4. Monte Carlo: сравнение PnL портфеля при constant volatility strategy и corrected strategy. Параметры как на рис. 3 с $\gamma=2$, $q_0=0$, $l_0=0$, $x_0=5$.

6.3. Пример медленного фактора

Kraft (2005): $\sigma(z)=\sqrt{z}$, CIR для $Z_t$: $c(z)=m-z$, $g(z)=\beta\sqrt{z}$ при Feller $\beta^2 < 2m$. Trading speed в этом режиме не меняется; сравниваем aim portfolio (рисунок 5). При положительной корреляции BM доходностей и фактора $Z$ target $\mathrm{aim}_t^s$ больше, чем при constant volatility; при отрицательной — меньше.

Рисунок 5
Рисунок 5. Сравнение aim portfolio при constant и corrected strategies по диапазону $\sigma(z)$. Слева: $x=0.2$, $l=0$; справа: $x=0.2$, $l=0.05$. Параметры: $\gamma=5$, $\lambda=\beta=0.1$, $\kappa=\eta=1$, $\rho_2=0.5$, $K=1$, $m=0.2$, $\beta=0.25$, $\delta=0.5$.

Monte Carlo при $T=2$ года (рисунок 6): empirical PDFs PnL gains при разных начальных $Z_0$. Статистика — таблица 2. Хотя $\delta$ в slow-factor setting меньше $\varepsilon$ в fast-scale примере, PnL gains сопоставимы по величине — slow factor сильнее корректирует constant volatility strategy.

Рисунок 6
Рисунок 6. Monte Carlo: эмпирические PDF PnL gains под corrected strategy при разных начальных $Z_0$. Параметры как на рис. 5, $q_0=0$, $l_0=0$, $x_0=2$, $\delta=0.252$.
Таблица 2. Статистика PnL gains (bps) под corrected strategy при разных начальных $Z_0$.
$Z_0$Mean (bps)Std (bps)95% CI
0.134.2532235.7638[29.6322, 38.8741]
0.214.726679.2333[13.1736, 16.2796]
0.38.312739.4801[7.5389, 9.0865]
0.45.016524.3310[4.5396, 5.4934]
0.52.382515.0484[2.0875, 2.6774]
0.60.680410.7493[0.4698, 0.8911]

7. Заключение

Эмпирические исследования показывают, что стохастическую волатильность эффективно моделируют два фактора — быстрый и медленный (Fouque et al., 2011). Мы расширяем классическую задачу optimal execution при постоянной волатильности (Gârleanu and Pedersen, 2013) этим multiscale volatility framework.

Влияние стохастической волатильности на optimal execution систематически изучается асимптотическими аппроксимациями, сохраняющими tractability оптимальной стратегии. Выведены первые два члена асимптотического разложения функции ценности Gârleanu and Pedersen (2013) отдельно для одного быстрого и одного медленного фактора; в §5 показано объединение при general multiscale SV. Rigorous accuracy analysis для option pricing — Fouque et al. (2011); для Merton problem — Fouque et al. (2022); для нашей постановки — future work.

Предложен метод small parameter approximations для системы Gârleanu and Pedersen (2016) при relatively small price impact. Численные примеры и Monte Carlo показывают, что resulting strategies ведут к более эффективному portfolio management и улучшению PnL при volatility fluctuations на distinct time scales.

Направления расширения:

  1. Stochastic liquidity. Almgren (2012): $\lambda$ варьируется в течение дня (пики на open/close). Реалистичное расширение — time-varying $\lambda$ с joint dynamics.
  2. Nonlinear price impact. Более sophisticated propagator models vs tractability optimal strategy.
  3. Partial information and expert opinions. Return signal полностью observable как OU; на практике — partial observation (Fouque et al., 2017a; Muhle-Karbe et al., 2023).

Приложения A, B и C

Приложение A содержит полные формулы second-order small price impact approximation (§2.6). Приложение B — метод вычисления производных коэффициентов $A$ для slow-scale correction. Приложение C — explicit formulas second-order approximation для одного медленного фактора. Технические доказательства и длинные алгебраические выкладки опущены в этом переводе. См. оригинал: arxiv.org/pdf/2507.17162.

Литература

Оригинал статьи: Chan, Sircar and Zimbidis, «Optimal Trading under Instantaneous and Persistent Price Impact, Predictable Returns and Multiscale Stochastic Volatility», arXiv:2507.17162