Finding random normally distributed number above zero
1 次查看(过去 30 天)
显示 更早的评论
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?
0 个评论
回答(3 个)
Torsten
2023-4-15
编辑:Torsten
2023-4-15
Maybe taking abs(randn()) will give you something useful. It's the folded normal distribution (which in this case equals the 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
2023-4-15
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
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:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1357408/image.png)
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Random Number Generation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!