[论文解读] Fractional calculus for power functions

Bartłomiej Dyda|arXiv (Cornell University)|Mar 17, 2011

Mathematical functions and polynomials参考文献 22被引用 4

一句话总结

本文推导了在 $\mathbb{R}^d$ 中形如 $(1 - |x|^2)_+^p$ 和 $x_d(1 - |x|^2)_+^p$ 的幂函数的分数阶拉普拉斯算子的显式公式，结果通过超几何函数表示。关键贡献在于以超几何函数 ${}_2F_1$ 表达 $\Delta^{eta/2}u_p^{(d)}(x)$ 的闭式表达，从而通过变分法实现对单位球内分数阶拉普拉斯算子第一特征值的精确估计。

ABSTRACT

We calculate the fractional Laplacian for functions of the form $u(x)=(1-|x|^2)_+^p$ and $v(x)=x_d u(x)$. As an application, we estimate the first eigenvalues of the fractional Laplacian in a ball.

研究动机与目标

• 计算 $\Delta^{\alpha/2}$ 在 $\mathbb{R}^d$ 中径向与反对称幂函数 $u_p^{(d)}(x) = (1 - |x|^2)_+^p$ 和 $v_p^{(d)}(x) = x_d(1 - |x|^2)_+^p$ 上的作用。
• 建立这些函数的分数阶拉普拉斯算子与超几何函数之间的联系，从而实现对特征值问题的解析处理。
• 将推导出的公式应用于估计单位球 $B \subset \mathbb{R}^d$ 中分数阶拉普拉斯算子的第一特征值 $\lambda_1$ 与 $\lambda_*$，特别关注径向与反对称特征函数。
• 通过由幂函数构造的测试函数，提供 $\lambda_1$ 与 $\lambda_*$ 的变分界，改进已知估计。

提出的方法

• 利用分数阶拉普拉斯算子的奇异积分定义及贝塔与伽马函数的性质，推导出 $\Delta^{\alpha/2}u_p^{(d)}(x)$ 的显式公式。
• 将结果表示为 $\Phi_{p,\alpha}^{(d)}(|x|^2)$，其中包含超几何函数 ${}_2F_1\left(\frac{\alpha + d}{2}, -p + \frac{\alpha}{2}; \frac{d}{2}; |x|^2\right)$。
• 利用恒等式 $u_p^{(d)}(x) = {}_2F_1\left(\frac{d}{2}, -p; \frac{d}{2}; |x|^2\right)$（当 $|x| < 1$ 时）将分数阶拉普拉斯算子重写为紧凑的超几何形式。
• 应用变分原理 $\lambda_1 \leq \mathcal{E}(u)/\|u\|_{L^2}^2$，选取测试函数 $u = (1 - |x|^2)_+^{\alpha/2} + \eta(1 - |x|^2)_+^{1 + \alpha/2}$，以推导上界。
• 通过涉及 $d$、$\alpha$ 与常数的闭式表达，推导出使瑞利商最小的最优 $\eta_{\text{min}}$。
• 建立 $d+2$ 维中径向特征函数与 $d$ 维中 $x_d$-反对称特征函数之间的一一对应关系，证明 $\lambda_1^{(d+2)} = \lambda_*^{(d)}$。

实验结果

研究问题

• RQ1 当 $|x| < 1$ 时，分数阶拉普拉斯算子 $\Delta^{\alpha/2}$ 作用于幂函数 $u_p^{(d)}(x) = (1 - |x|^2)_+^p$ 的显式形式是什么？
• RQ2 分数阶拉普拉斯算子如何作用于反对称函数 $v_p^{(d)}(x) = x_d(1 - |x|^2)_+^p$？
• RQ3 能否使用幂函数作为测试函数来估计单位球中分数阶拉普拉斯算子的特征值问题？
• RQ4 在维度 $d+2$ 中的径向特征值与在维度 $d$ 中的反对称特征值之间存在何种关系？
• RQ5 单位球中分数阶拉普拉斯算子的第一特征值 $\lambda_1$ 与第一反对称特征值 $\lambda_*$ 的精确下界与上界是什么？

主要发现

• 当 $|x| < 1$ 时，$u_p^{(d)}(x)$ 的分数阶拉普拉斯算子为 $\Phi_{p,\alpha}^{(d)}(|x|^2) = \mathcal{A}_{d,-\alpha} B(-\alpha/2, p+1) \pi^{d/2} \Gamma(d/2)^{-1} \cdot {}_2F_1\left(\frac{\alpha + d}{2}, -p + \frac{\alpha}{2}; \frac{d}{2}; |x|^2\right)$。
• 对于反对称函数 $v_p^{(d)}(x)$，分数阶拉普拉斯算子满足 $\Delta^{\alpha/2}v_p^{(d)}(x) = x_d \Phi_{p,\alpha}^{(d+2)}(|x|^2)$。
• 单位球中分数阶拉普拉斯算子的第一特征值 $\lambda_1$ 满足 $\lambda_1 \geq \mu_{d,\alpha}$，其中 $\mu_{d,\alpha}$ 由比较函数 $\psi$ 导出。
• 第一反对称特征值 $\lambda_*$ 满足 $\lambda_* = \lambda_1^{(d+2)}$，确立了维度转移对应关系。
• 通过 $u_{j + \alpha/2}^{(d)}$ 的线性组合计算 $\lambda_1$ 与 $\lambda_*$ 的数值上界，最优 $\eta$ 由涉及 $d$ 与 $\alpha$ 的闭式表达给出，从而在所有维度与 $\alpha \in (0,2)$ 范围内实现紧密估计。
• 当 $\alpha = 1.5$，$d = 3$ 时，$\lambda_1$ 的上界为 $4.61848$，$\lambda_*$ 的上界为 $5.05977$，数值最小化过程表现出接近最优的收敛性。