使用 Matlab 模拟 Gamma 分布式 RV

如何解决使用 Matlab 模拟 Gamma 分布式 RV

我是 Matlab 的初学者，我需要对此进行解释。我正在尝试生成伽马分布的随机变量，并将它们与 gamcdf 的输出进行比较。

 n = 10000;
% inline function: gamma distributed random variables %generated by the sum of exponentials
randGamma = inline('-sum(log(rand(1,n)))/a','n','a');
for a = [0.2,1.5]
x = randGamma(n,a);
t = linspace(0,5/a,500);
figure
plot(sort(x),(1:n)/n,'rx'),hold on
plot(t,gamcdf(t,a),'k')
axis([0,max(t),1.1])
title(['Gammadistribution a = ',num2str(a)])
end

在我得到的数字中，x 没有出现。我其实有疑问

plot(sort(x),hold on.

谁能帮我解决这个问题。我将不得不对其他发行版进行这样的实现。所以理解这一点会对我有很大帮助。提前致谢