[论文解读] On local analytic expansions of densities in the context of (micro-)hypoelliptic and classes of semi-elliptic equations

Joerg Kampen|arXiv (Cornell University)|Dec 2, 2010

Advanced Mathematical Modeling in Engineering参考文献 9被引用 3

一句话总结

本文为一类抛物型偏微分方程的转移密度发展了局部解析展开，通过利用基本解表示和迭代方法，将这些展开的应用从严格椭圆算子扩展到（微）次椭圆和半椭圆方程——例如Heston型模型和美式期权定价中出现的方程——从而在金融建模背景下实现了更高阶数值格式，显著提升了计算精度。

ABSTRACT

Explicit representations of densities for linear parabolic partial differential equations are useful in order to design computation schemes of high accuracy for a considerable class of diffusion models. Approximations of lower order based on the WKB-expansion have been used in order to compute Greeks in standard models of the interest rate market (cf. [2]). However, it turns out that for higher order approximations another related expansion leads to more accurate schemes. We compute a local explicit formula for a class of parabolic problems and determine a lower bound of the time horizon where it holds (given a certain bounded domain). Although the local analytic expansions hold only for strictly elliptic equations we show that the expansions can be used in order to design higher order schemes for various types of (micro)-hypoelliptic and semi-elliptic equations, e.g. the reduced market models considered in [7] or front fixing schemes for multivariate American derivatives [3].

研究动机与目标

• 将转移密度的局部解析展开从严格椭圆方程推广至（微）次椭圆和半椭圆抛物型PDE。
• 在有界区域内，为这类局部展开在有效时间范围上提供严格的下界。
• 将这些展开适配于数学金融中扩散模型的高精度数值格式。
• 建立自由边界问题（如美式期权定价）中迭代格式的收敛性与正则性。
• 在低维子空间上满足Hörmander条件的前提下，将严格椭圆算子的收敛结果推广至微次椭圆算子。

提出的方法

• 推导出基本解的显式局部解析展开形式 $ p(t,x;0,y) = \frac{1}{(4\pi t)^{n/2}} \exp\left(-\frac{d_R^2(x,y)}{4t}\right) \left(\sum_{k=0}^{\infty} d_k(x,y) t^k\right) $，适用于严格椭圆方程。
• 通过引入与时间相关的Riemann距离 $ d_R^2(t,x,y) $，将展开式推广至时变情形。
• 构建一个迭代格式 $ v^n $，通过求解一系列系数依赖于 $ F^{n-1} $ 的抛物型PDE，其中 $ F $ 为自由边界函数。
• 利用涉及变系数线性抛物型算子基本解 $ p_{v^n} $ 的卷积表示。
• 通过在超平面 $ \{x_1 = 1\} $ 上施加边界条件 $ v^n_{x_1} - v^n = -tK $，确保与自由边界条件的一致性。
• 通过包含基本解法向导数的积分方程求解边界源项 $ \phi_{v^n} $，其中强迫项 $ \Gamma_{v_n} $ 由初始数据和边界约束导出。

实验结果

研究问题

• RQ1 能否将严格椭圆方程密度的局部解析展开推广至金融模型中常见的（微）次椭圆和半椭圆PDE？
• RQ2 在有界区域内，此类局部展开的有效时间范围最大可达多长？
• RQ3 如何利用这些展开设计美式期权定价和简化市场模型的更高阶数值格式？
• RQ4 在微次椭圆算子背景下，自由边界问题的迭代格式在何种条件下收敛？
• RQ5 严格椭圆算子的正则性与收敛性结果在多大程度上可推广至满足低维子空间上Hörmander条件的算子？

主要发现

• 形式为 $ p(t,x;0,y) = \frac{1}{(4\pi t)^{n/2}} \exp\left(-\frac{d_R^2(x,y)}{4t}\right) \left(\sum_{k=0}^{\infty} d_k(x,y) t^k\right) $ 的局部解析展开适用于严格椭圆方程，并为更高阶近似提供了基础。
• 通过构造在维数 $ d \leq n $ 的子空间上满足Hörmander条件的迭代格式 $ v^n $，将展开推广至（微）次椭圆方程。
• 局部展开的有效时间范围存在下界，确保了在适配的Banach空间中迭代格式的收敛性。
• 对于形如 $ L = \frac{1}{2}\sum_{i=1}^m A_i^2 u + A_0 u $ 的微次椭圆算子，其迭代格式的收敛性与正则性已得证，其中向量场 $ A_i $ 在 $ \mathbb{R}^d \subset \mathbb{R}^n $ 上满足Hörmander条件。
• 该方法可实现Heston模型和美式期权前沿固定格式中希腊值与价值函数的高精度计算。
• 渐近分析表明，当 $ t \downarrow 0 $ 时，$ \sqrt{t} F_t $ 保持有界，支持了初始时刻附近迭代格式的稳定性。