latex 分数_【重点】分数阶微积分和分数阶微分方程数值实验（9）——R-L分数阶微积分...

参考文献

薛定宇《分数阶微积分学与分数阶控制》

数值实现

Matlab 2019a 主要基于薛定宇开发的FOTF工具箱

地球物理局 地震波动力学实验室 分数阶组

# 本文优先满足个人需求
# 欢迎批评指正，禁止转载
# 侵删

---

Riemann-Liouville分数阶微积分

1. 高阶整数阶积分公式

先考虑整数阶积分，显然，给定函数

的一阶积分可以表示为：

从该结果再求一次积分，则可以得出原函数的二阶积分：

类似地，可以推导出

阶积分公式为：

2.Riemann-Liouville分数阶微积分定义

定义1：函数

的

阶Riemann-Liouville积分定义为：

式中，

表示Riemann-Liouville定义。

两端的下标为积分的上下限。

定义2：函数

的Riemann-Liouville分数阶微分的定义为：

式中，

，且

。

3.常用函数的Riemann-Liouville微积分公式

（1）幂函数

，且

，则

满足：

（2）指数函数

的Riemann-Liouville分数阶导数为：

（3）一些常用函数的Riemann-Liouville分数阶导数公式为：

式中，

为Heaviside函数，

为冲击函数（Dirac函数）。

因为变量

在积分边界与被积函数中同时出现，所以一般

的分数阶导数的解析解很难求出，只能使用数值方法来求解Riemann-Liouville分数阶导数。

4. 初始时刻平移的性质

如果已知

，则Riemann-Liouville分数阶导数可以表示为：

式中，

。

---

数值实现

1.数值算法

（1）由已知的函数获得函数值向量

；

（2）令

，计算出二项式系数

；

（3）由下式计算向量

：

（4）在计算步长

下计算向量

的

阶差分。

• Riemann_Liouville_differential.m

% 薛定宇《分数阶微积分学与分数阶控制》 p47 

2.例子：已知阶跃函数0.75阶导数的解析解为

，试比较Grunwald-Letnikov与Riemann-Liouville分数阶导数数值解的精度。

• example3_9.m

% 薛定宇《分数阶微积分学与分数阶控制》 p47 

• 图像

直接计算的Riemann-Liouville分数阶导数有较大误差。在实际应用中不建议采用该函数计算Riemann-Liouville分数阶微积分。

---