[论文解读] Fractional variational calculus in terms of a combined Caputo derivative

Agnieszka B. Malinowska, Delfim F. M. Torres|arXiv (Cornell University)|Jul 5, 2010

Fractional Differential Equations Solutions参考文献 27被引用 11

一句话总结

本文提出了一种新的分数阶导数 ${{}^{C}D^{\beta,\beta}_{\beta}}$，通过凸参数 $\gamma$ 将左、右Caputo导数结合，从而扩展了变分问题的类别。该文推导了该广义导数的分数阶Euler–Lagrange方程、等周约束及横截条件，扩展了Agrawal（2007年）的前期工作，使物理与工程中非局部、具有记忆依赖性的系统建模更具灵活性。

ABSTRACT

We generalize the fractional Caputo derivative to the fractional derivative ${^CD^{α,β}_γ}$, which is a convex combination of the left Caputo fractional derivative of order $α$ and the right Caputo fractional derivative of order $β$. The fractional variational problems under our consideration are formulated in terms of ${^CD^{α,β}_γ}$. The Euler-Lagrange equations for the basic and isoperimetric problems, as well as transversality conditions, are proved.

研究动机与目标

• 通过引入以 $\gamma$ 为参数的左、右Caputo导数的凸组合，推广Caputo分数阶导数。
• 利用新的 ${{}^{C}D^{\beta,\beta}_{\beta}}$ 导数构建分数阶变分问题，以更灵活地模拟非局部与记忆依赖现象。
• 为该广义变分框架推导必要的最优性条件——Euler–Lagrange方程、等周约束及横截条件。
• 通过在 $\gamma = 1/2$ 时恢复为已知特例，扩展分数阶变分计算中现有成果，特别是Agrawal（2007年）的研究。
• 为多目标分数阶变分计算的未来研究提供基础，具有更广泛的物理适用性。

提出的方法

• 将组合Caputo导数 ${{}^{C}D^{\beta,\beta}_{\beta}}$ 定义为左、右Caputo导数（阶数分别为 $\alpha$ 与 $\beta$）的凸组合，其中 $\gamma \in [0,1]$。
• 通过泛函 $\mathcal{J}(y) = \int_a^b L(x, y(x), {{}^{C}D^{\beta,\beta}_{\beta}}y(x)) \, dx$ 构建分数阶变分问题。
• 通过标准变分原理推导分数阶Euler–Lagrange方程，得到如下系统：$D^{\beta,\alpha}_{1-\gamma} \partial_{N+1} F + \partial_y F = 0$。
• 通过引入拉格朗日乘子，将框架扩展至等周问题，引入约束 $\mathcal{G}^j(\mathbf{y}) = \int_a^b G^j[\mathbf{y}](x) \, dx = l_j$。
• 通过分析端点处的变分，推导自由边界问题的横截条件，确保在自然边界条件下满足最优性。
• 通过引入松弛变量 $\phi_j(x)$ 处理不等式约束，导出互补松弛条件 $\lambda_j \phi_j(x) = 0$。

实验结果

研究问题

• RQ1 如何将Caputo分数阶导数推广，使其以凸组合形式同时包含左、右导数？
• RQ2 以新 ${{}^{C}D^{\beta,\beta}_{\beta}}$ 导数构建的变分问题，其必要最优性条件（Euler–Lagrange方程）是什么？
• RQ3 在该广义分数阶框架下，等周约束如何影响Euler–Lagrange方程的结构？
• RQ4 当变分问题中边界值未固定时，会导出何种横截条件？
• RQ5 新导数能否作为特例恢复文献中已知结果，例如Agrawal（2007年）的研究？

主要发现

• 广义分数阶导数 ${{}^{C}D^{\beta,\beta}_{\beta}}$ 允许比标准Caputo导数更广泛的变分问题类别，尤其避免了 $y(a) = y(b) = 0$ 的严格限制。
• 基本问题的分数阶Euler–Lagrange方程推导为 $D^{\beta,\alpha}_{1-\gamma} \partial_{N+1} F + \partial_y F = 0$，其中 $F$ 为扩展空间中的拉格朗日量。
• 对于等周问题，当约束非退化时，$N$ 个Euler–Lagrange方程的系统为 $D^{\beta,\alpha}_{1-\gamma} \partial_{N+i} F + \partial_i F = 0$（$i = 2, \ldots, N$）。
• 为自由边界问题推导了横截条件，确保在端点未固定时仍满足最优性。
• 通过引入松弛变量 $\phi_j(x)$ 处理不等式约束，导出互补松弛条件 $\lambda_j \phi_j(x) = 0$，并构造了包含 $\lambda_j$ 的修正拉格朗日量。
• 当 $\gamma = 1/2$ 时，结果退化为Agrawal（2007年）的研究成果，证实了与前期工作的兼容性，并验证了该推广的合理性。