caputo分数阶导数程序_基于Caputo导数的分数阶非线性振动系统响应计算

基于

Caputo

导数的分数阶非线性振动系统响应计算

李亚杰

1,2

，

吴志强

*1,2

，

章国齐

1,2

【摘

要】

摘

要

:

研究了含分数阶

Caputo

导数的非线性振动系统响应的数值计

算方法。首先，由

Caputo

分数阶导数算子的叠加关系，得到含分数阶导数项

非线性振动系统状态方程的标准形式。其次，基于

Caputo

导数与

Riemann-

Liouville

导数和

Grunwald-Letnikov

导数间的关系，推导计算了

Caputo

导

数的一般数值迭代格式。本文方法不要求状态方程中各分数阶导数阶数相等，

弱化了已有算法中对分数阶导数阶数的限制，并可推广到多自由度的情形。随

后，选择若干有解析解的算例验证了本文方法的正确性。最后，以多吸引子共

存的分数阶

Duffing

振子系统为例，比较

Caputo

和

GL

两种算法所得结果，

说明了用

GL

算法求解存在的问题。

【期刊名称】

计算力学学报

【年

(

卷

),

期】

2018(035)004

【总页数】

7

【关键词】

关键词

:

分数阶导数；算法；数值迭代格式；导数阶数；系统初值

修改稿收到日期

: 2017-07-08.

基金项目

:

国家自然科学基金

(11672349,11372211)

；

国家

973

计划

(2014CB046805)

资助项目

.

1

引

言

分数阶微积分是对整数阶微积分的一种推广，将微积分运算的阶次从传统的整

数阶扩展到了分数阶和复数阶，至今已有三百多年的发展历史，已成为研究反

常扩散、多孔介质力学、非牛顿流体力学、粘弹性力学以及软物质物理等学科