Метод Галёркина и проекционные методы: решение УЧП финансовой математики

Математический базис метода Галёркина строится на слабой (вариационной) формулировке УЧП. Вместо того чтобы требовать выполнения дифференциального уравнения Блэка-Шоулза $mathcal{L}V = 0$ в каждой точке пространства, мы умножаем его на гладкую "пробную" (тестовую) функцию $v(x)$ из подходящего пространства Соболева (обычно $H^1$) и интегрируем по всей области: $int (mathcal{L}V)v dx = 0$. Применяя интегрирование по частям (формулу Грина), мы понижаем порядок производных, что снижает требования к гладкости решения и позволяет находить цены для опционов с разрывными выплатами (например, бинарных).

Суть метода состоит в проецировании бесконечномерного пространства Соболева на конечномерное подпространство. Мы ищем приближенное решение в виде линейной комбинации базисных функций: $V_h(x, t) = sum c_i(t) phi_i(x)$. Базисные функции $phi_i$ обычно выбираются как кусочно-полиномиальные функции (сплайны или "шляпки" на треугольной сетке). Подставляя это разложение в вариационную формулировку, УЧП элегантно превращается в систему обыкновенных дифференциальных уравнений: $M dot{c}(t) + S c(t) = 0$, где $M$ — матрица масс (Mass Matrix), а $S$ — матрица жесткости (Stiffness Matrix), элементы которых вычисляются интегралами от базисных функций.

Метод конечных элементов обладает колоссальным преимуществом при оценке барьерных опционов со сложными криволинейными барьерами. Алгоритмы триангуляции Делоне позволяют сгущать расчетную сетку (Mesh) исключительно в тех областях, где кривизна цены (Гамма опциона) максимальна — около страйка и барьеров, экономя вычислительные ресурсы. Для интегрирования по времени полученной системы ОДУ традиционно используется $ heta$-схема (включая неявную схему Эйлера и Кранка-Николсона), гарантирующая безусловную устойчивость расчетов.

Развитие проекционных методов привело к созданию Разрывного метода Галёркина (Discontinuous Galerkin, DG), который снимает ограничение на непрерывность базисных функций между элементами сетки. Это критически важно в финансовой математике при моделировании скачкообразной диффузии (где интегро-дифференциальный оператор порождает сильные градиенты) и в кредитном риск-менеджменте для моделирования волн дефолтов. Программная реализация методов Галёркина на C++ с использованием библиотек линейной алгебры (PETSc) позволяет квантам инвестиционных банков решать сложнейшие задачи ценообразования с аналитической точностью.

Список литературы:
1. Achdou Y., Pironneau O. Computational Methods for Option Pricing. — SIAM, 2005.
2. Toivanen J. Numerical Methods for Financial Engineering. — Springer, 2010.
3. Марчук Г.И. Методы вычислительной математики. — М.: Наука, 1989.
Эксперты: Ив Ашду, Оливье Пиронно, Гурий Марчук.