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

参考文献

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

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

数值实现

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

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

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

前情回顾

形式主义的居士:分数阶微积分和分数阶微分方程数值实验(1)——误差函数与补误差函数​zhuanlan.zhihu.com
c2778f055ba33d81dd442d4e409b5b31.png

Gamma函数

定义:

性质:

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

Matlab实现

1.试绘制

区间上的Gamma函数曲线
% 薛定宇《分数阶微积分学与分数阶控制》 例2.3

图像

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

值下的不完全Gamma函数曲线
% 薛定宇《分数阶微积分学与分数阶控制》 p13 例2.7

图像

c3960370cd409513c889dbf97454ed40.png

87fe1d0fe55f72b5e8e95d6c948002c5.png
拦——@CycleUser
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值