File Exchange

image thumbnail

Simple Rejection Sampling

version 1.0 (1.68 KB) by

Return values sampled from a user defined distribution. Samples are not guaranteed IID.

5 Downloads

Updated

View License

SAMPLEDIST Sample from an arbitrary distribution
    sampleDist(f,M,N,b) retruns an array of size X of random values sampled from the distribution defined by the probability density function referred to by handle f, over the range b = [min, max]. M is the threshold value for the proposal distribution, such that f(x) < M for all x in b.
 
    sampleDist(...,true) also generates a histogram of the results with an overlay of the true pdf.
 
    Examples:
    %Sample from a step function over [0,1]:
    X = sampleDist(@(x)1.3*(x>=0&x<0.7)+0.3*(x>=0.7&x<=1),...
                   1.3,1e6,[0,1],true);
    %Sample from a normal distribution over [-5,5]:
    X = sampleDist(@(x) 1/sqrt(2*pi) *exp(-x.^2/2),...
                   1/sqrt(2*pi),1e6,[-5,5],true);

Comments and Ratings (0)

MATLAB Release
MATLAB 7.9 (R2009b)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video