参考文献
薛定宇《分数阶微积分学与分数阶控制》
吴强、黄建华《分数阶微积分》
数值实现
Matlab 2019a 主要基于薛定宇开发的FOTF工具箱
蜜酒厅通讯社 固体地球物理学部
封面及文中照片感谢 @CycleUser 友情提供
前情回顾
形式主义的居士:分数阶微积分和分数阶微分方程数值实验(1)——误差函数与补误差函数zhuanlan.zhihu.com![c2778f055ba33d81dd442d4e409b5b31.png](https://i-blog.csdnimg.cn/blog_migrate/0d27708ae6cfeabcb0ce833519db209a.jpeg)
Gamma函数
定义:
性质:
-
-
-
-
-
-
- 【常见函数】
- 【常见函数】
- 【常见函数】
Matlab实现
1.试绘制
% 薛定宇《分数阶微积分学与分数阶控制》 例2.3
图像
![04c84e31594e07509faab0bb4143cce6.png](https://i-blog.csdnimg.cn/blog_migrate/61d9a2007b739c805dcd1e84cd3e1fc0.png)
2.试用Matlab证明
% 薛定宇《分数阶微积分学与分数阶控制》 p12 例2.4
窗口
Gamma_1_2 =
pi^(1/2)
Gamma_3_2 =
pi^(1/2)/2
Gamma_5_2 =
(3*pi^(1/2))/4
Gamma_7_2 =
(15*pi^(1/2))/8
3.试用Matlab证明
% 薛定宇《分数阶微积分学与分数阶控制》 p12 例2.5
窗口
I1 =
pi/sin(pi*z)
I21 =
-pi/(z*sin(pi*z))
I22 =
gamma(z + 1)*gamma(1 - z)
I3 =
pi/cos(pi*z)
I4 =
1
4.一些积分问题可以由Gamma函数表示,试求出下面不定积分的解析解
% 薛定宇《分数阶微积分学与分数阶控制》 p13 例2.6
窗口
I1 =
piecewise(0 < real(m) & 0 < real(n), beta(m, n)/2, real(m) <= 0 | real(n) <= 0, int(t^(2*m - 1)*(- t^2 + 1)^(n - 1), t, 0, 1))
I2 =
piecewise(x < 1, (2^x*pi^(1/2)*gamma(x/2))/(2*gamma(1/2 - x/2)), x == 1 | 2 <= x, int(t^(x - 1)*cos(t), t, 0, Inf), x < 1 | x in Dom::Interval(1, 2), cos((pi*x)/2)*gamma(x))
5.试求出复参数Gamma函数值
gamma_complex.m
% 薛定宇《分数阶微积分学与分数阶控制》 p13 例2.7
example_2_7.m
% 薛定宇《分数阶微积分学与分数阶控制》 p13 例2.7
窗口
I =
0.112294242282773 + 0.323612885476767i 0.112294242282773 - 0.323612885476767i
不完全Gamma函数
定义:
Matlab实现
6.试绘制出不同
% 薛定宇《分数阶微积分学与分数阶控制》 p13 例2.7
图像
![c3960370cd409513c889dbf97454ed40.png](https://i-blog.csdnimg.cn/blog_migrate/3ca7bbb52f52f33b3d011ad5062f215b.png)
![87fe1d0fe55f72b5e8e95d6c948002c5.png](https://i-blog.csdnimg.cn/blog_migrate/8da01a77373fe5f7fdb00d74ce3ab026.jpeg)