Generation of random numbers

3 次查看（过去 30 天）
SURBHIT
SURBHIT 2014-5-27
编辑： Mahdi 2014-5-27
How can I generate uniformly distributed random numbers in the range 10^(-4) to 10^(+4)? I tried using rand function but it does not give me the required range. This is what I tried:
dd = (10^(4) - 10^(-4)).*rand(1000,1) + 10^(-4); plot(dd);

请先登录，再回答此问题。

采纳的回答

Mahdi
Mahdi 2014-5-27
编辑：Mahdi 2014-5-27
I think what you're trying to do is done using randn, so that
dd=(10^(4) - 10^(-4)).*randn(1000,1) + 10^(-4)
hist(dd,1000) # To see it in a histogram with 1000 bins
Note the range of x-axis
  3 个评论
显示 1更早的评论
Mahdi
Mahdi 2014-5-27
Then what you did is correct. It looks quite normal to me.
John D'Errico
John D'Errico 2014-5-27
NO. randn is NOT the solution here, since it generates numbers that can fall ANYWHERE in theory.

请先登录，再进行评论。

更多回答（3 个）

Shashank Prasanna
Shashank Prasanna 2014-5-27
That sounds about right. Take a look at the first example in the documentation:
http://www.mathworks.com/help/matlab/ref/rand.html
  1 个评论
SURBHIT
SURBHIT 2014-5-27
but the output is only in the range 10^(-1) to 10^(+4) how can i get values below 10^(-1)? i am not able to generate the smaller values

请先登录，再进行评论。

Roger Wohlwend
Roger Wohlwend 2014-5-27
You use the correct method. When your range is from 10^(-4) to 10^4 it is quite normal that values below 10^(-1) are very, rare.
John D'Errico
John D'Errico 2014-5-27
I think the problem is that while you SAY you want numbers that are uniformly distributed over that interval, you don't really WANT them uniformly distributed.
The probability is actually quite small that for an interval like that, that you will see numbers less than say 1. In fact, suppose we did generate a uniform randm number in the interval [1e-4,1e4]?
What is the probability that any given such deviate will be in the sub-interval [1e-4,1]? This is easy to compute.
(1 - 1e-4)/(1e4 - 1e-4)
ans =
9.999e-05
Essentially, you need to have many such samples before you ever see a SINGLE such deviate in the sub-interval.
So while you SAY you want uniform, I think you really want numbers that are uniformly distributed in a log domain. That is, the logs of your deviates will be uniformly distributed.
% these numbers will be uniform over the interval [-4,4]
E = 8*rand(1,10000) - 4;
% The numbers in R will have the distribution you desire
R = 10.^E;

类别

在 Help Center 和 File Exchange 中查找有关 Uniform Distribution (Continuous) 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by