matlab绘制Gamma分布的概率密度曲线！

首先简单给出gamma函数的 x>=0,且当x<0时，。

其中,其中afa>0.

称X服从形状参数α，反尺度参数β的g分 布，记为X ∼ Gamma(α, β)

function s=gramma
x=0:0.1:20;
s1=grammma(x,1,0.5);

s2=grammma(x,2,0.5);

s3=grammma(x,3,0.5);

s4=grammma(x,5,1.0);

s5=grammma(x,9,2.0);

plot(x,s1,x,s2,x,s3,x,s4,x,s5)
legend('afa=1,beta=0.5','afa=2,beta=0.5','afa=3,beta=0.5','afa=5,beta=1.0','afa=9,beta=2.0')
function s=grammma(x,afa,beta)

s=beta^afa.*x.^(afa-1).*exp(-beta*x)/fun(afa);
function s=fun(afa)
syms x
y=exp(-x)*x^(afa-1);
s=int(y,x,0,inf);

注意matlab中本身有一个名字叫gamma的内置函数，所以避开此名字。

运行结果如下图，可以观察到只有afa=1，beta=0.5的趋势不同，其余只有集中点的差别，得以验证gamma分布。