Generate random polar angles

Joakim

2014 5 28

1 个回答

回答已采纳

6 次查看（30 天）

0 个投票

I want to generate random polar angles given a uniform distribution on the unit sphere. In other words, I want to generate random numbers from the probability density function p(theta) = sin(theta)/2 in the interval [0,pi]. How do I do this?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Follow Question

采纳的回答

Roger Stafford 2014-5-29

编辑：Roger Stafford 2014-5-29

在 MATLAB Online 中打开

0 个投票

The requested distribution of theta is generated below. Consequently the following code will places yellow dots on the unit sphere surface in a statistically uniform manner.

n = 8192;
theta = acos(1-2*rand(n,1)); % <-- The requested distribution
phi = 2*pi*rand(n,1);
x = sin(theta).*cos(phi);
y = sin(theta).*sin(phi);
z = cos(theta);
plot(x,y,z,'y.')