问题描述
我是 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.
谁能帮我解决这个问题。我将不得不对其他发行版进行这样的实现。所以理解这一点会对我有很大帮助。提前致谢
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)