[论文解读] Integral equations, fractional calculus and shift operator

D. Babusci, G. Dattoli+1 more|arXiv (Cornell University)|Jul 29, 2010

Fractional Differential Equations Solutions参考文献 4被引用 4

一句话总结

本文提出了一种基于分数阶导数和广义平移算子（特别是以 $x\partial_x$ 为基础的伸缩算子）的新型运算演算框架，用于求解一类积分方程，包括广义 Lamb-Bateman 方程和贝塞尔型方程。其关键贡献在于提出了一种统一方法，通过 $x\partial_x$ 的分数阶幂表示解，从而利用积分表示和特殊函数（如贝塞尔-赖特函数与广义超几何形式）获得解析解。

ABSTRACT

We present an extension of a previously developed method employing the formalism of the fractional derivatives to solve new classes of integral equations. This method uses different forms of integral operators that generalizes the exponential shift operator.

研究动机与目标

• 将分数阶微积分与平移算子形式化方法扩展至经典 Lamb-Bateman 方程以外的新类积分方程求解。
• 为涉及非整数次幂算子 $x\partial_x$ 的积分方程开发系统性求解方法。
• 将指数平移算子推广至伸缩与广义函数算子，以求解具有非线性缩放核的积分方程。
• 建立 $\Gamma(x\partial_x + 1)$ 及其逆的运算恒等式，以求解具有指数核和幂律核的方程。
• 推测并数值验证非整数次幂 $x\partial_x$ 的分数阶斯特林数恒等式。

提出的方法

• 该方法利用拉普拉斯变换恒等式 $a^{-\nu} = \frac{1}{\Gamma(\nu)}\int_0^\infty e^{-sa} s^{\nu-1} ds$，将 $ (x\partial_x)^{-\nu} $ 表示为指数伸缩算子 $ \hat{E}(-s) = \exp(-s x\partial_x) $ 的积分形式。
• 伸缩算子 $ \hat{E}(\lambda) = \exp(\lambda x\partial_x) $ 的作用为 $ \hat{E}(\lambda)g(x) = g(e^\lambda x) $，从而可将积分方程转化为运算形式。
• 对于具有高斯或指数核的方程，该方法计算如 $ \int_0^\infty e^{-y^2} \hat{E}(-y^2) dy $ 的积分，导出分数阶算子如 $ \sqrt{x\partial_x} $。
• 方程 $ \int_0^\infty dy\, u(e^{-y^2}x) = f(x) $ 的解为 $ u(x) = \frac{2}{\sqrt{\pi}} \sqrt{x\partial_x} f(x) $，随后通过 $ (x\partial_x)^{-1/2} $ 表示为积分形式。
• 利用广义平移算子恒等式 $ \exp(\lambda q(x)\partial_x)g(x) = g\left[F^{-1}(\lambda + F(x))\right] $，其中 $ F(x) = \int^x d\xi / q(\xi) $，求解具有非线性缩放的方程。
• 该方法进一步推广至形如 $ \int_0^\infty dt\, f(x g(t)) $ 的积分，重写为 $ \hat{O}(x\partial_x)f(x) $，其中 $ \hat{O}(\mu) = \int_0^\infty [g(t)]^\mu dt $。

实验结果

研究问题

• RQ1 如何利用分数阶导数与广义平移算子求解涉及非整数次幂 $ x\partial_x $ 的积分方程？
• RQ2 对于核为 $ e^{-y^2} $、$ e^{-y} $ 或 $ (1+y^2)^{-\nu} $ 的积分方程，其解的运算形式是什么？
• RQ3 广义平移算子形式化方法能否应用于具有非线性缩放变换（如 $ F^{-1}(-y^2 + F(x)) $）的方程？
• RQ4 如何从运算角度表达 $ \Gamma(x\partial_x + 1) $ 的逆，以求解具有指数核的方程？
• RQ5 对于非整数 $ \nu \in (0,1) $，$ (x\partial_x)^\nu $ 的结构是什么？能否通过广义斯特林数将其展开为 $ x^k \partial_x^k $ 的形式？

主要发现

• 修正的 Lamb-Bateman 方程 $ \int_0^\infty dy\, u(e^{-y^2}x) = f(x) $ 的解为 $ u(x) = \frac{2}{\sqrt{\pi}} \sqrt{x\partial_x} f(x) $，可简化为 $ u(x) = \frac{2}{\pi} \int_0^x \frac{\partial_\xi f(\xi)}{\sqrt{\ln(x/\xi)}} d\xi $。
• 对于 $ \int_0^\infty dy\, e^{-y} u(y x) = f(x) $，解为 $ u(x) = \left\{ \Gamma(x\partial_x + 1) \right\}^{-1} f(x) $；当 $ f(x) = e^{-x} $ 时，结果为 $ u(x) = J_0(2\sqrt{x}) $，即零阶贝塞尔函数。
• 当 $ f(x) = e^{-x} $ 且核为 $ u(y^m x) $ 时，解为 $ u(x) = W_0(x|m) = \sum_{k=0}^\infty \frac{(-1)^k}{k! (m k)!} x^k $，即零阶贝塞尔-赖特函数。
• 方程 $ \int_0^\infty dy\, u(\sqrt{x^2 + 2y^2}) = f(x) $ 的解为 $ u(x) = -\frac{2}{\sqrt{\pi}} \sqrt{ -\frac{1}{x}\partial_x } f(x) $，其推导基于恒等式 $ \exp(\lambda \frac{1}{x}\partial_x) g(x) = g(\sqrt{x^2 + 2\lambda}) $。
• 对于 $ \int_0^\infty dt\, \frac{1}{(1+t^2)^{\nu x\partial_x}} \sin x $，解为 $ F(x;\nu) = \frac{\sqrt{\pi}}{2} \sum_{n=0}^\infty \frac{(-1)^n}{(2n+1)!} \frac{\Gamma((2n+1)\nu - 1/2)}{\Gamma((2n+1)\nu)} x^{2n+1} $，并经数值验证。
• 本文推测 $ (x\partial_x)^\nu = \sum_{k=0}^\infty S(k,\nu) x^k \partial_x^k $ 对于 $ 0 < \nu < 1 $ 成立，数值证据支持该恒等式，但严格证明仍待完成。