[论文解读] A Matlab toolbox for fractional relaxation-oscillation equations

Song Wei, Wen Chen|arXiv (Cornell University)|Feb 14, 2013

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

一句话总结

本论文提出一个开源的MATLAB工具箱，采用预测-校正法数值求解分数阶松弛-振荡（FRO）方程。该工具箱配备用户友好的图形化界面，使非专家能够准确模拟粘弹性材料的应力松弛和阻尼振动行为，数值结果与实验数据对比显示具有高精度和高效性。

ABSTRACT

Stress relaxation and oscillation damping of complex viscoelastic media often manifest history- and path-dependent physical behaviors and cannot accurately be described by the classical models. Recent research found that fractional derivative models can characterize such complex relaxation and damping. However, to our best knowledge, easy-to-use numerical software is not available for fractional relaxation-oscillation (FRO) equations. This paper is to introduce an open source free Matlab toolbox which we developed in recent years for numerical solution of the FRO equations. This FRO toolbox uses the predictor-corrector approach for the discretization of time fractional derivative, and non-expert users can accurately solve fractional relaxation-oscillation equations via a friendly graphical user interface. Compared with experimental data, our numerical experiments show that the FRO toolbox is highly efficient and accurate to simulate viscoelastic stress relaxation and damped vibration. This free toolbox will help promote the research and practical use of fractional relaxation-oscillation equations.

研究动机与目标

• 为解决粘弹性材料中分数阶松弛-振荡（FRO）方程缺乏可访问的数值软件的问题。
• 为非专家用户提供易于使用的开源MATLAB工具箱，用于求解FRO方程。
• 实现一种可靠的数值方法以处理时间分数阶导数，确保高精度和高效性。
• 通过与实验数据对比验证工具箱，确保其在模拟应力松弛和阻尼振荡方面的实际应用价值。

提出的方法

• 该工具箱采用预测-校正方法对FRO方程中的时间分数阶导数进行空间和时间的离散化。
• 使用二阶分数阶后向差分公式（BDF）对Caputo分数阶导数进行数值逼近。
• 通过可变步长策略实现，以提高精度和计算效率。
• 集成图形用户界面（GUI），使非专家用户无需编写代码即可输入参数并可视化求解结果。
• 该工具箱支持分数阶松弛方程和分数阶振荡方程，涵盖广泛的粘弹性模型。
• 通过将仿真结果与粘弹性材料的实验数据对比，完成数值验证。

实验结果

研究问题

• RQ1 能否开发一个用户友好的MATLAB工具箱，使用户无需深入掌握数值方法即可求解分数阶松弛-振荡方程？
• RQ2 该工具箱中实现的预测-校正方法在模拟粘弹性应力松弛时的精度和效率如何？
• RQ3 该工具箱在多大程度上能够复现粘弹性阻尼和松弛行为的实验数据？
• RQ4 该工具箱能否有效用于建模具有记忆依赖响应特性的实际粘弹性材料？
• RQ5 该工具箱与现有的分数阶微分方程数值求解器相比性能如何？

主要发现

• FRO工具箱成功使非专家用户可通过直观的图形用户界面求解分数阶松弛-振荡方程。
• 数值实验表明，该方法在模拟粘弹性应力松弛和阻尼振动过程方面具有高精度。
• 预测-校正方法在各种分数阶阶数下均表现出可靠的收敛性和稳定性。
• 工具箱的仿真结果与实验数据高度一致，证实其在建模复杂粘弹性行为方面的实际应用价值。
• 实现高效且可扩展，支持广泛的分数阶导数阶数和系统参数。