Matlab：从定制的概率密度函数生成随机数

Question

我有一个1977年至1983年1月份3个月降水量的数据集（见附件）。但是，我想根据这些数据生成1984 - 1990年期间的降水数据。因此，我想知道是否有可能创建降水量（1977-1983）的定制概率密度函数，并由此产生所需时期（1984-1990）的随机数（降水数据）。

这可能在Matlab中，有人可以通过这样做帮助我吗？

提前致谢！

Answer 1

直方图将为您提供PDF估算 - 只需将箱数除以样本总数。从那里你可以通过整合来估计CDF。最后，您可以选择介于0和1之间的均匀分布的随机数，并估计将产生该数字的CDF的参数。也就是说，如果y是您选择的随机数，那么您希望找到x，使得CDF（x）= y。 x的值将是具有所需PDF的随机数。

Answer 2

如果您有“统计和机器学习工具箱”，则可以使用'Kernel Distribution'方法评估数据的PDF：

Percip_pd = fitdist(Percip,'Kernel');

然后用它从同一分布中生成N随机数：

y = random(Percip_pd,N,1);

Answer 3

引用@ AnonSubmitter85：

“通过积分来估计CDF。最后，你可以选择一个介于0和1之间的均匀分布的随机数，并估计产生该数字的CDF的参数。也就是说，如果y是你选择的随机数，那么你想要找到x使得CDF（x）= y .x的值将是一个带有所需PDF的随机数。“

%random sampling
 N=10; %number of resamples

 pdf = normrnd(0, 1, 1,100); %your pdf
 s = cumsum(pdf); %its cumulative distribution

 r = rand(N,1); %random numbers between 0 and 1
 for ii=1:N
   inds = find(s>r(ii));
   indeces(ii)=inds(1); %find first value greater than the random number
 end
 resamples = pdf(indeces) %the resamples