О приближённом представлении дробного броуновского движения

Основная информация

• ID статьи: 2503.04575
• Название: On Approximate Representation of Fractional Brownian Motion
• Автор: Konstantin A. Rybakov
• Классификация: math.PR (теория вероятностей)
• Дата публикации: март 2025 г. (препринт arXiv, версия v2 от 14 октября 2025 г.)
• Ссылка на статью: https://arxiv.org/abs/2503.04575

Аннотация

В данной работе исследуется ортогональное разложение дробного броуновского движения относительно полиномов Лежандра. Такое разложение имеет как теоретическое, так и практическое значение, поскольку может быть использовано для приближения и моделирования дробного броуновского движения в непрерывном времени. В статье приводятся соотношения для среднеквадратической ошибки приближения и проводится сравнение с ранее полученными результатами.

Исследовательский контекст и мотивация

1. Основная проблема: как эффективно приближать и моделировать дробное броуновское движение (fractional Brownian motion, FBM) в непрерывном времени
2. Значимость: дробное броуновское движение находит широкое применение в финансовом моделировании, физике, инженерии и других областях; его моделирование критично для проверки теоретических результатов и решения практических задач
3. Ограничения существующих методов:
   • Традиционные методы на дискретной временной сетке обеспечивают точность только в узлах сетки
   • Существующие универсальные методы моделирования гауссовских процессов имеют высокую вычислительную сложность
   • Отсутствуют методы непрерывного времени, позволяющие точно вычислять ошибку приближения
4. Исследовательская мотивация: на основе предыдущей работы автора 2 предложить более простой метод вычисления коэффициентов разложения, снизить вычислительную сложность и повысить точность

Основные вклады

1. Теоретический вклад: вывод явных соотношений для коэффициентов разложения ядра, связанного с интегральным представлением дробного броуновского движения (относительно полиномов Лежандра)
2. Улучшение метода: предоставление более простого метода вычисления коэффициентов разложения по сравнению с предыдущей работой 2, снижение вычислительной сложности
3. Анализ ошибок: установление уравнений для точного вычисления среднеквадратической ошибки полиномиального приближения дробного броуновского движения
4. Приближение в непрерывном времени: реализация приближённого представления дробного броуновского движения в непрерывном времени, а не только в дискретных временных точках

Подробное описание метода

Постановка задачи

Для дробного броуновского движения $B_H(t)$ (показатель Хёрста $H \in (0,1)$) найти коэффициенты его разложения относительно ортогонального базиса полиномов Лежандра и реализовать приближение в непрерывном времени:

$B_H(t) = \sum_{i=0}^{\infty} B_H^i \hat{P}(i,t)$

где $\hat{P}(i,t)$ — нормализованные полиномы Лежандра.

Основная математическая база

Определение дробного броуновского движения:

• Функция ковариации: $R_H(t,\tau) = \frac{t^{2H} + \tau^{2H} - |t-\tau|^{2H}}{2}$
• Интегральное представление: $B_H(t) = \int_0^t k_H(t,\tau) dB(\tau)$

Ядро: $k_H(t,\tau) = a_H(t-\tau)^{H-1/2} \,_2F_1\left(\frac{1}{2}-H, H-\frac{1}{2}, H+\frac{1}{2}, 1-\frac{t}{\tau}\right)\mathbf{1}(t-\tau)$

где $a_H = \sqrt{\frac{\pi H(1-2H)}{\Gamma(2-2H)\cos\pi H}}$ (при $H \neq 1/2$).

Основные теоретические результаты

Теорема 1 (явное выражение коэффициентов разложения): Для $j$-го столбца матрицы $K^H$: $K^H_{*j} = a_H\Gamma(3/2-H)\sqrt{\frac{2j+1}{T}} \sum_{k=0}^j \frac{l_{jk}(3/2-H)_k}{T^k(H+1/2+k)k!} F^{H+1/2+k}$

Для $H \neq 1/2$ коэффициенты разложения удовлетворяют: $K^H_{ij} = a_H\Gamma(1/2-H)\sqrt{\frac{2j+1}{T}} F^{H+1/2}_i \sum_{k=0}^j (-1)^{j-k} \frac{1/2-H+k}{H+1/2+k} \Pi^{(ij)}_k$

Реализация алгоритма

1. Вычисление коэффициентов: использование рекуррентных соотношений (12)-(13) для вычисления коэффициентов разложения степенных функций
2. Построение матрицы: конструирование матрицы коэффициентов разложения $K^H$ согласно теореме 1
3. Генерация случайных величин: генерирование независимых стандартных нормальных случайных величин $V_i$
4. Построение приближения: $\tilde{B}_H(t) = \sum_{i=0}^{L-1} \tilde{B}^H_i \hat{P}(i,t)$

Экспериментальная установка

Параметры

• Временной интервал: $T = [0,1]$
• Показатель Хёрста: $H = 0.1, 0.2, \ldots, 0.9$
• Длина усечения: $L = 4, 8, 16, 32, 64, 128$

Метрики оценки

Среднеквадратическая ошибка приближения: $\varepsilon = E\int_T (B_H(t) - \tilde{B}_H(t))^2 dt$

Теоретическая формула: $\varepsilon = \|k_H(\cdot)\|^2_{L^2(T^2)} - \|\bar{K}^H\|^2$

где $\|k_H(\cdot)\|^2_{L^2(T^2)} = \frac{T^{2H+1}}{2H+1}$

Результаты экспериментов

Основные результаты

Таблица 1: Среднеквадратическая ошибка приближения $\varepsilon$

Ключевые выводы

1. Анализ сходимости:
   • При $H = 0.5$ сходимость наиболее быстрая
   • При приближении $H$ к граничным значениям (0 или 1) скорость сходимости значительно снижается
   • Относительная ошибка составляет $\varepsilon/(2H+1)^{-1}$
2. Сравнение методов:
   • Новый метод качественно близок к результатам работы 2
   • Даже при малых значениях $L$ ошибки обоих методов очень близки
   • Это указывает на то, что представление матричного произведения оказывает незначительное влияние на точность, но обеспечивает вычислительную устойчивость
3. Анализ скорости сходимости:
   • Скорость сходимости разложения степенной функции $t^\alpha$ составляет $O(L^{-(2\alpha+1)})$
   • При $\alpha \to -1/2$ скорость сходимости стремится к нулю, что объясняет низкую сходимость при граничных значениях показателя Хёрста

Связанные работы

1. Обзор классических методов:
   • Представление Мандельброта-Ван Несса 25
   • Интегральное представление Декрюсфонда-Юстюнеля 15
   • Обзоры различных методов моделирования 4-7
2. Связь данной работы с предыдущими исследованиями:
   • Основана на предыдущем спектральном методе автора 2
   • Использует математическое описание спектральной формы систем управления 16,17
   • Расширяет представление винеровского процесса полиномами Лежандра 29,30
3. Расширения приложений:
   • Применимо к дробному процессу Орнштейна-Уленбека
   • Подходит для моделирования дробного броуновского моста
   • Может быть расширено на итерированные представления случайных интегралов

Заключение и обсуждение

Основные выводы

1. Успешно выведены явные формулы для коэффициентов разложения ядра дробного броуновского движения
2. Установлен метод точного вычисления среднеквадратической ошибки приближения
3. Реализовано приближённое представление дробного броуновского движения в непрерывном времени
4. Подтверждена эффективность и вычислительные преимущества метода

Ограничения

1. Требования к точности вычислений: новый метод требует высокоточной арифметики с плавающей точкой или символических вычислений
2. Граничные эффекты: при приближении показателя Хёрста к 0 или 1 сходимость замедляется
3. Вычислительная сложность: хотя вычисление коэффициентов упрощено, всё ещё требуются матричные операции

Направления будущих исследований

1. Расширение на многомерное дробное броуновское движение
2. Применение к более сложным случайным процессам (например, дробным стохастическим дифференциальным уравнениям)
3. Разработка более эффективных численных алгоритмов

Глубокая оценка

Преимущества

1. Теоретическая строгость: обеспечивает полный математический вывод и строгий анализ ошибок
2. Методологическая инновация: упрощает вычисление коэффициентов разложения, повышает практическую применимость метода
3. Достаточная экспериментальная база: метод проверен на различных значениях показателя Хёрста и длинах усечения
4. Практическая ценность: предоставляет новый инструмент для моделирования дробного броуновского движения в непрерывном времени

Недостатки

1. Требования к вычислениям: высокие требования к численной точности могут ограничить практическое применение
2. Производительность на границах: сходимость при экстремальных значениях показателя Хёрста требует улучшения
3. Анализ сравнения: недостаточно подробное сравнение производительности с другими существующими методами

Влияние

1. Научный вклад: вносит значительный вклад в теорию дробного броуновского движения и методы вычислений
2. Практическая ценность: имеет потенциальное применение в финансовой математике, физическом моделировании и других областях
3. Методологическое значение: предоставляет справочную базу для ортогональных разложений других случайных процессов

Области применения

1. Приложения, требующие моделирования дробного броуновского движения в непрерывном времени
2. Сценарии с требованиями точного контроля ошибки приближения
3. Численное решение дробных стохастических дифференциальных уравнений
4. Ценообразование финансовых производных и управление рисками

Список литературы

В статье цитируется 37 соответствующих источников, включая:

• Основную теорию дробного броуновского движения 22-25
• Методы моделирования случайных процессов 1,3,4-7
• Спектральные методы и ортогональные разложения 16,17,26
• Теорию дробного исчисления 32