[论文解读] Calculus of Variations with Classical and Fractional Derivatives

Tatiana Odzijewicz, Delfim F. M. Torres|arXiv (Cornell University)|Jul 4, 2010

Fractional Differential Equations Solutions参考文献 30被引用 3

一句话总结

本文提出了一种新颖的变分框架，通过在拉格朗日函数中同时包含 $ y' $ 和 $ {}_aD_t^eta y $，统一了经典导数与分数阶导数在变分法中的应用。该框架推导了适用于任意边界条件和等周约束问题的欧拉-拉普拉斯型最优性条件，并确保在 $ k=0 $ 时精确恢复经典变分法，而不同于以往方法仅在 $ \alpha \to 1 $ 的极限下实现恢复。其主要贡献在于提出了一种一致的扩展方法，避免了标准黎曼-刘维尔分数阶导数所施加的限制性条件 $ y(a)=y(b)=0 $。

ABSTRACT

We give a proper fractional extension of the classical calculus of variations. Necessary optimality conditions of Euler-Lagrange type for variational problems containing both classical and fractional derivatives are proved. The fundamental problem of the calculus of variations with mixed integer and fractional order derivatives as well as isoperimetric problems are considered.

研究动机与目标

• 将经典变分法扩展至单一变分框架中，同时包含整数阶与分数阶导数。
• 通过确保在 $ k=0 $ 时精确恢复经典情况，解决以往分数阶变分法中仅在 $ \alpha \to 1 $ 极限下可恢复经典理论的局限性。
• 为涉及混合导数 $ y' + k\,{}_aD_t^\alpha y $ 的变分问题建立必要的最优性条件，支持任意边界条件 $ y(a)=y_a $、$ y(b)=y_b $。
• 在分数阶设定下处理等周问题，提供超越以往工作的系统性扩展。

提出的方法

• 引入依赖于 $ t $、$ y $ 和组合导数 $ v = y' + k\,{}_aD_t^\alpha y $ 的广义拉格朗日函数，其中 $ k \in \mathbb{R} $ 且 $ \alpha \in (0,1) $。
• 对黎曼-刘维尔分数阶导数应用分部积分法，推导扩展泛函的欧拉-拉普拉斯方程。
• 推导出包含左、右黎曼-刘维尔导数的分数阶微分方程形式的必要最优性条件。
• 利用Mittag-Leffler函数表达变分问题中出现的分数阶微分方程的解。
• 应用拉格朗日乘子法处理等周约束，导出增强型拉格朗日函数形式。
• 通过解析示例验证该方法，包括经典情形（$ k=0 $）、极限情形 $ \alpha \to 1 $，以及 $ \alpha = 1/2 $ 的分数阶等周问题。

实验结果

研究问题

• RQ1 如何在单一变分公式中一致地扩展经典变分法，以同时包含经典导数与分数阶导数？
• RQ2 对于涉及混合导数 $ y' + k\,{}_aD_t^\alpha y $ 且具有任意边界条件的变分问题，其必要最优性条件是什么？
• RQ3 如何在分数阶变分问题中整合等周约束，同时保持与经典情形的一致性？
• RQ4 为何以往的分数阶变分方法无法精确恢复经典理论，这一局限性如何克服？
• RQ5 当 $ \alpha = 1/2 $、$ k=1 $ 且给定边界值时，分数阶等周问题的极值解的解析形式是什么？

主要发现

• 所提出的框架在 $ k=0 $ 时能精确恢复经典变分法，而不同于以往方法仅在 $ \alpha \to 1 $ 极限下实现恢复。
• 对于基本问题，必要最优性条件被推导为包含左、右黎曼-刘维尔导数的分数阶微分方程。
• 对于等周问题，极值解为 $ y(t) = \int_0^t E_{1-\alpha,1}(-k(t-\tau)^{1-\alpha}) \xi \, d\tau $，该解满足约束条件与边界条件。
• 当 $ k=0 $ 时，极值解退化为 $ y(t) = \xi t $，与标准问题的经典极小化器一致。
• 在 $ \alpha \to 1 $ 的极限下，分数阶问题退化为经典变分问题，且极小化器被重新缩放为 $ y(t) = \frac{\xi}{k+1}t $，验证了其一致性。
• 当 $ k=1 $、$ \alpha = 1/2 $ 且 $ \xi=1 $ 时，极值解为 $ y(t) = -\left(1 - e^t \mathrm{erfc}(\sqrt{t}) + \frac{2\sqrt{t}}{\sqrt{\pi}} \right) $，该解满足分数阶等周约束与边界条件。