Number generator for 0.00006 to 30 range
Show older comments
I am attempting to run monte carlo simulations and wanted a random number between the end members of 0.000069 and 30 to populate this equation through each pass of a for loop. After looking at my model's output, I noticed that the Kd is heavily skewed toward larger numbers. I then came across this thread which led me to read about creating one's own random number stream. Is there a number generator out there that gives equal weight to each order of magnitude...ie I'd like the number generator to be equally likely to draw .0005, .05, 5, etc. Originally I didn't realize that unifrnd was not appropriate for 0 to 1 range. Thank you in advance for any guidance.
Kd=unifrnd(0.000069,30);
Answers (2)
weikang zhao
on 14 Sep 2022
1 vote
It looks like you don't want these random numbers to follow a uniform distribution, you want their logarithms to follow a uniform distribution.
try this way:
a=log10(0.000069);
b=log10(30);
Kd=10^(unifrnd(a,b));
% Kd=unifrnd(0.000069,30);
logminmax = log10([0.000069,30]);
for i=1:20
kdlog = unifrnd(logminmax(1), logminmax(2));
kd = 10^kdlog
end
Categories
Find more on Parallel Computing in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!