Generate random numbers from a non-parametric PDF

Question

TheBeginner 2014-11-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/164586-generate-random-numbers-from-a-non-parametric-pdf

回答： Torsten 2014-11-28

Hello,

I know similar questions have already been asked, but I can't seem to solve my problem.

I want to generate random numbers from a probality density function having two blobs, for instance the sum of gaussians. My PDF would therefore look like this :

p(x) = N(x;mu1,sigma1) + N(x;mu2,sigma2);

I have tried the following procedures :

1. Generate random numbers using uniform distribution. Take the PDFs of a vector and invert it using "icdf"

2.Generate random numbers using uniform distribution. Write the CDF (analytic in this case) and invert it.

My code looks like this (for the first idea) :

%Uniform random numbers
  u = rand(1,1000);
%PDF
  x = linspace(0,xmax,N)';
  values = gauss(x,mu1,sigma1) + gauss(x,mu2,sigma2);
  pd = fitdist(values,'kernel');
%Inversion
  y = icdf(pd,u);