Gamma函数相关matlab代码

1、Gamma函数：

Gamma函数matlab代码：

x=:0.5:5
syms t
y=int(t.^(x-)*exp(-t),,inf)
y=double(y)
plot(x,y,'r-o','linewidth',)

图像如下：

2、lgΓ(x)函数

matlab代码：

x=:0.1:5
syms t
y=int(t.^(x-)*exp(-t),,inf)
y=log10(double(y))
plot(x,y,'r-o','linewidth',)

图形：

Gamma分布：

matlab代码：

a=,b=0.5
syms c
d=int(c^(a-)*exp(-c),,inf)
t=:0.1:
g=b.^a*(t.^(a-)).*exp(-b*t)/d
g=double(g)
plot(t,g,'r-.','linewidth',)

图形：

matlab代码：

t=:0.1:
syms c

a=,b=0.5
d=int(c^(a-)*exp(-c),,inf)
g=b.^a*(t.^(a-)).*exp(-b*t)/d
g=double(g)
plot(t,g,'r-.','linewidth',)
hold on;
a=,b=0.5
d=int(c^(a-)*exp(-c),,inf)
g=b.^a*(t.^(a-)).*exp(-b*t)/d
g=double(g)
plot(t,g,'g-.','linewidth',)
hold on
a=,b=0.5
d=int(c^(a-)*exp(-c),,inf)
g=b.^a*(t.^(a-)).*exp(-b*t)/d
g=double(g)
plot(t,g,'b-.','linewidth',)
hold on
a=,b=
d=int(c^(a-)*exp(-c),,inf)
g=b.^a*(t.^(a-)).*exp(-b*t)/d
g=double(g)
plot(t,g,'c-.','linewidth',)
hold on
a=,b=
d=int(c^(a-)*exp(-c),,inf)
g=b.^a*(t.^(a-)).*exp(-b*t)/d
g=double(g)
plot(t,g,'m-.','linewidth',)
hold on;
legend('a=1,b=0.5','a=2,b=0.5','a=3,b=0.5','a=5,b=1','a=9,b=2');

图形：