[论文解读] The Laplace Transform Method for Linear Differential Equations of the Fractional Order

Igor Podlubný|ArXiv.org|Oct 30, 1997

Fractional Differential Equations Solutions参考文献 2被引用 173

一句话总结

本文提出了一种基于拉普拉斯变换的系统方法，用于求解具有常系数的线性分数阶微分方程，利用双参数马特格勒夫函数及其拉普拉斯变换。主要贡献在于推导出n项方程的分数阶格林函数的一般公式，通过级数反演实现解析解，并将适用范围扩展至标准和连续分数阶导数。

ABSTRACT

The Laplace transform method for solving of a wide class of initial value problems for fractional differential equations is introduced. The method is based on the Laplace transform of the Mittag-Leffler function in two parameters. To extend the proposed method for the case of so-called "sequential" fractional differential equations, the Laplace transform for the ''sequential'' fractional derivative is also obtained. Besides that, tools necessary for testing candidate solutions by direct substitution in corresponding equations are introduced: fractional derivatives of the Mittag-Leffler function and the rule for the fractional differentiation of integrals depending on a parameter. Definition of the fractional Green's function is given and some of its properties, necessary for constructing solutions of initial-value problems for fractional linear differential equations, are presented. Explicit expressions for the fractional Green's function for the special cases of one-, two-, three- and four-term equations are given, as well as the explicit expression for an arbitrary n-term fractional linear ordinary differential equation with constant coefficients. Several examples of solution of various types of fractional differential equations, including fractional diffusion-wave equation. The bibliography (54 items) covers many field of possible application of fractional derivatives.

研究动机与目标

• 开发一种统一且高效的方法，用于求解任意实数阶的线性分数阶微分方程。
• 克服现有方法（如级数展开或迭代法）在处理复杂方程时效率低下的局限性。
• 为初值问题提供一种基于拉普拉斯变换和马特格勒夫函数的系统化方法。
• 定义并推导出各种n项方程的分数阶格林函数的显式表达式。
• 将该方法扩展至标准和连续分数阶导数，确保解结构的一致性。

提出的方法

• 该方法采用双参数马特格勒夫函数 $ E_{\alpha,\beta}(z) $ 的拉普拉斯变换，其源自积分变换。
• 利用分数阶导数的拉普拉斯变换，特别是针对连续算子，推导出解核。
• 通过求解核 $ g_n(p) = \left(\sum_{k=0}^n a_k p^{\beta_k}\right)^{-1} $ 的逆拉普拉斯变换，建立一般解框架，其中 $ \beta_n > \cdots > \beta_0 $。
• 解以包含多项式系数和广义马特格勒夫函数的级数形式表示，支持逐项反演。
• 整合用于验证候选解的工具，包括马特格勒夫函数的分数阶导数及参数依赖积分的微分法则。
• 为一至四项方程推导出分数阶格林函数的显式公式，最终得出一般n项表达式。

实验结果

研究问题

• RQ1 如何系统地应用拉普拉斯变换求解具有任意实数阶导数的线性分数阶微分方程？
• RQ2 对于具有常系数的一般n项线性分数阶微分方程，其分数阶格林函数的结构是什么？
• RQ3 标准与连续分数阶导数的解有何不同？两者之间存在何种共同结构？
• RQ4 能否利用马特格勒夫函数及其拉普拉斯变换为初值问题推导出闭式解？
• RQ5 Wright函数和多项式展开在分数阶拉普拉斯变换反演过程中起什么作用？

主要发现

• 推导并应用了双参数马特格勒夫函数 $ E_{\alpha,\beta}(z) $ 的拉普拉斯变换，作为求解分数阶微分方程的基础工具。
• 为一至四项方程获得了分数阶格林函数的显式表达式，一般n项情形以包含多项式系数的级数形式表达。
• 对于n项方程 $ \sum_{k=0}^n a_k D^{\beta_k} y(t) = f(t) $，解表示为 $ y(t) = \int_0^t G_n(t-\tau) f(\tau) d\tau $，其中 $ G_n(t) $ 为分数阶格林函数。
• 一般n项方程的分数阶格林函数表达为 $ G_n(t) = \frac{1}{a_n} \sum_{m=0}^\infty \frac{(-1)^m}{m!} \sum_{\substack{k_0+\cdots+k_{n-2}=m \\ k_i \geq 0}} (m; k_0,\dots,k_{n-2}) \prod_{i=0}^{n-2} \left(\frac{a_i}{a_n}\right)^{k_i} t^{\gamma m + \delta} E_{\gamma, \epsilon}^{(m)}(\cdot) $，其中 $ \gamma = \beta_n - \beta_{n-1} $，且 $ E_{\gamma, \epsilon}^{(m)} $ 表示马特格勒夫函数的m阶导数。
• 该方法成功处理了标准与连续分数阶导数，两种形式下均得到一致的格林函数。
• 该方法为分数阶微积分中一大类初值问题提供了解析解，包括物理学、工程学和金融学中出现的问题。