[论文解读] Stability Analysis of Fractional Differential Equations with Unknown Parameters

Mehmet Emir Köksal|arXiv (Cornell University)|Sep 15, 2017

Fractional Differential Equations Solutions参考文献 43被引用 18

一句话总结

本文提出了一种用于带未知参数的分数阶微分方程（FDEs）参数稳定性分析的图形化D-分解方法，避免了复杂的解析计算。通过在参数空间中推导实根边界、复根边界和无穷远根边界，该方法可直观识别稳定性区域，实现对大范围参数区间的鲁棒稳定性评估。该方法通过基准FDEs得到验证，展示了在具有分数阶动力学的工程系统中可靠且实用的稳定性分析能力。

ABSTRACT

In this paper, the stability of fractional differential equations (FDEs) with unknown parameters is studied. FDEs bring many advantages to model the physical systems in the nature or man-made systems in the industry. Because this representation has a property between linear differential equations and nonlinear differential equations. Therefore, the designer may use the FDEs to model complex systems instead of nonlinear differential equations which have hard mathematical background. Using the graphical based D-decomposition method, we investigate the parametric stability analysis of FDEs without complicated mathematical analysis. To achieve this, stability boundaries are obtained firstly, and then the stability region set depending on the unknown parameters is found. The applicability of the presented method is shown considering some benchmark equations which are often used to verify the results of a new method. Simulation examples shown that the method is simple and give reliable stability results.

研究动机与目标

• 解决分数阶微分方程（FDEs）在未知或可变参数下的稳定性分析挑战，其复杂性高于经典稳定性分析。
• 通过实现从零到无穷的完整参数区间分析，克服现有方法仅限于小范围参数不确定性的局限。
• 为工程师和研究人员开发一种实用、可视化且计算高效的稳定性评估方法，无需依赖复杂的解析解。
• 提供系统性框架，用于识别一致阶与非一致阶FDEs在参数空间中的稳定性区域。
• 通过可视化参数变化时稳定性如何演变，实现鲁棒性分析，支持分数阶系统的系统设计与控制。

提出的方法

• 将D-分解方法应用于将s平面中的稳定性边界（虚轴）映射到参数域中的三类边界：实根边界、复根边界和无穷远根边界。
• 基于Caputo或Grünwald-Letnikov分数阶导数定义，利用FDE传递函数的特征广义多项式，推导出稳定性边界的解析表达式。
• 以Matignon稳定性准则为基础，要求所有根满足 |arg(−λi)| > απ/2，从而在参数空间中定义稳定性条件。
• 在可能的情况下，将FDE传递函数转换为一致阶形式，利用三角恒等式简化边界曲线的推导。
• 通过绘制推导出的边界方程（例如，a = −bω−α sin(0.5πα)/sin(πα)，c = −aω²α cos(πα) − bωα cos(0.5πα)）随参数变化的图像，生成稳定性区域。
• 通过二维或三维图像（必要时用颜色表示第四维）可视化结果，展示稳定性区域随参数变化的演化过程。

实验结果

研究问题

• RQ1 如何在不依赖复杂解析方法的前提下，分析具有未知参数的分数阶微分方程的稳定性？
• RQ2 在系数变化的FDEs中，将稳定区域与不稳定区域分隔开的参数边界是什么？
• RQ3 当参数如阻尼或刚度（a, b, c）在大范围区间内变化，包括正负值时，稳定性区域如何变化？
• RQ4 所提出的方法能否有效验证并扩展已知基准FDEs（如Basset方程或工业加热炉模型）的稳定性结果？
• RQ5 一致阶与非一致阶分数阶对稳定性区域的形状和鲁棒性有何影响？

主要发现

• 所提出的D-分解方法通过图形化可视化成功识别出具有未知参数的FDEs的稳定性区域，消除了对复杂解析解的依赖。
• 对于α = 0.5的一致阶FDE，复根边界表现出对称性：若(a, c)位于曲线上，则(c, a)也位于曲线上，而这一性质在一般非一致阶情况下不成立。
• 在示例2中，图8展示了α = 0.2且b ∈ [−5, −1]时的稳定性区域，表明随着α增大，区域缩小，并在α → 0时填满右上象限。
• 当α = 0.8时，若将b替换为−b，稳定性区域关于原点与α = 0.2情况对称，验证了方法的一致性。
• 在示例3中，工业加热炉的标称参数（a = 14994, b = 6009.5, c = 1.69）位于计算出的稳定性区域内，无需额外计算即验证了Sondhi与Hote的稳定性结论。
• 图12中的三维稳定性区域表明，当b为正时，所有正a与c值，以及所有负a与c值（在椭圆区域内）均保持稳定性，表明在宽广参数范围内具有鲁棒性。