[论文解读] Fractional discrete Laplacian versus discretized fractional Laplacian

Óscar Ciaurri, Luz Roncal+3 more|arXiv (Cornell University)|Jul 17, 2015

Nonlinear Partial Differential Equations参考文献 23被引用 12

一句话总结

本文通过半群理论在晶格 $\mathbb{Z}_h = h\mathbb{Z}$ 上引入了分数阶离散拉普拉斯算子 $(-\Delta_h)^s$，为连续分数阶拉普拉斯算子 $(-\Delta)^s$ 提供了严格逼近。对于 Hölder 连续函数，建立了 $\ell^\infty$ 误差估计，其阶为 $h^{\beta - 2s}$，在更高正则性下收敛性得到改善，并推导出具有精确衰减性质的显式点态核表示。

ABSTRACT

We define and study some properties of the fractional powers of the discrete Laplacian $$(-Δ_h)^s,\quad\hbox{on}~\mathbb{Z}_h = h\mathbb{Z},$$ for $h>0$ and $0

研究动机与目标

• 通过半群方法定义并研究晶格 $\mathbb{Z}_h$ 上的分数阶离散拉普拉斯算子 $(-\Delta_h)^s$。
• 当 $h \to 0$ 时，比较通过半群理论定义的离散分数阶拉普拉斯算子与连续分数阶拉普拉斯算子的离散化版本。
• 在函数正则性假设尽可能弱的条件下，建立该逼近的定量 $\ell^\infty$ 误差估计。
• 推导 $(-\Delta_h)^s$ 的显式点态公式，给出核表示，并证明 Hölder 正则性估计。
• 将分析扩展至负幂次，即离散分数阶积分 $(-\Delta_h)^{-s}$，并在高维空间中提供核的渐近估计。

提出的方法

• 通过半群公式定义 $(-\Delta_h)^s$：$(-\Delta_h)^s u_j = \frac{1}{\Gamma(-s)} \int_0^\infty \left( e^{t\Delta_h}u_j - u_j \right) \frac{dt}{t^{1+s}}$，其中 $e^{t\Delta_h}$ 是半离散热方程的解。
• 利用半离散热方程 $\partial_t w_j = \Delta_h w_j$ 的解（初始数据为 $u_j$）来表达分数阶幂。
• 通过修正贝塞尔函数和渐近展开，推导 $(-\Delta_h)^s$ 的显式点态核表示。
• 通过核差值与函数正则性的估计，建立 $(-\Delta_h)^s(r_h u)$ 与 $r_h(( -\Delta)^s u)$ 在 $\ell^\infty$ 范数下的误差界。
• 应用最陡下降法与渐近分析，推导一维和二维情形下核的衰减速率。
• 利用 Gamma 函数与修正贝塞尔函数 $I_k(t)$ 的性质，控制核的行为并证明一致有界性。

实验结果

研究问题

• RQ1 通过半群理论定义的分数阶离散拉普拉斯算子 $(-\Delta_h)^s$ 与 $h \to 0$ 时连续分数阶拉普拉斯算子 $(-\Delta)^s$ 的离散化版本相比如何？
• RQ2 在函数正则性假设尽可能弱的条件下，该离散逼近在 $\ell^\infty$ 范数下的收敛速率是多少？
• RQ3 能否为 $(-\Delta_h)^s$ 推导出具有核表示的显式点态公式？
• RQ4 $(-\Delta_h)^s$ 的 Hölder 正则性特征是什么？
• RQ5 在高维空间中，特别是 $\mathbb{Z}_h^2$ 中，核及其衰减特性如何表现？

主要发现

• 对于 $u \in C^{0,\beta}(\mathbb{R})$ 且 $2s < \beta$，$(-\Delta_h)^s(r_h u)$ 与 $r_h(( -\Delta)^s u)$ 之间的 $\ell^\infty$ 误差被控制在 $C[u]_{C^{0,\beta}} h^{\beta - 2s}$ 以内。
• 对于 $u \in C^{1,\beta}(\mathbb{R})$ 且 $2s < \beta$，一阶差分算子 $D_+( -\Delta_h)^s(r_h u)$ 的误差被控制在 $C[u]_{C^{1,\beta}} h^{\beta - 2s}$ 以内。
• 当 $2s > \beta$ 时，误差界改善为 $C[u]_{C^{1,\beta}} h^{\beta - 2s + 1}$，反映出在更高正则性下收敛性的提升。
• 对于 $u \in C^{k,\beta}(\mathbb{R})$ 且 $k + \beta - 2s$ 不为整数，同时满足 $2s < k + \beta$，$l$ 阶差分算子 $D_+^l$ 的 $\ell^\infty$ 误差被控制在 $C[u]_{C^{k,\beta}} h^{\beta - 2s + k - l}$ 以内，其中 $l = \lfloor k + \beta - 2s \rfloor$。
• 推导出 $(-\Delta_h)^s u_j$ 的显式点态公式，表示为与核 $K_s(j)$ 的卷积，且对大 $|j|$ 有 $|K_s(j)| \lesssim |j|^{-1-2s}$。
• 在二维情形下推导出核的渐近估计，表明其衰减行为与一维一致，并可推广至更高维空间。