matlab gamma函数_分数阶微积分和分数阶微分方程数值实验（2）——Gamma函数

参考文献

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

吴强、黄建华《分数阶微积分》

数值实现

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

蜜酒厅通讯社 固体地球物理学部

封面及文中照片感谢 @CycleUser 友情提供

前情回顾

形式主义的居士：分数阶微积分和分数阶微分方程数值实验（1）——误差函数与补误差函数​zhuanlan.zhihu.com

---

Gamma函数

定义：

性质：

• 【常见函数】
• 【常见函数】
• 【常见函数】

Matlab实现

1.试绘制

区间上的Gamma函数曲线

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

图像

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.试绘制出不同

值下的不完全Gamma函数曲线

% 薛定宇《分数阶微积分学与分数阶控制》 p13 例2.7

图像

---

拦——@CycleUser