Finding random normally distributed number above zero

1 次查看（过去 30 天）
Tishan Anantharajah
评论： Paul 2023-4-15
I am currently using randn to produce a normally distributed number however would like to only take the positive side for the normally distributed curve. How could this be done? Would I have to use a loop?

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

回答（3 个）

Torsten
Torsten 2023-4-15
编辑：Torsten 2023-4-15
Ran in:
Maybe taking abs(randn()) will give you something useful. It's the folded normal distribution (which in this case equals the truncated normal distribution):
https://en.wikipedia.org/wiki/Folded_normal_distribution
https://en.wikipedia.org/wiki/Truncated_normal_distribution
n = 100000;
y = abs(randn(n,1));
hold on
histogram(y,'Normalization','pdf')
fun = @(x)sqrt(2/pi)*exp(-x.^2/2);
x = 0:0.01:4;
plot(x,fun(x))
hold off
  2 个评论
Paul
Paul 2023-4-15
Ran in:
Can also use the HalfNormalDistribution if it's needed to do more than just generate random numbers.
n = 100000;
y = abs(randn(n,1));
hold on
histogram(y,'Normalization','pdf')
x = 0:0.01:4;
plot(x,pdf('HalfNormal',x,0,1))
hold off

请先登录，再进行评论。

Star Strider
Star Strider 2023-4-15
The truncate function is another option.
Image Analyst
Image Analyst 2023-4-15
You could just throw away any less than zero:
r = stDev * randn(100000, 1) + desiredMeanValue;
r(r<=0) = [];
% Or equivalently
r = r(r > 0);
To learn other fundamental concepts, invest 2 hours of your time here:
MATLAB Academy - Free 2 hour training

类别

在 Help Center 和 File Exchange 中查找有关 Triangular Distribution 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by