How to generate random points inside sphere using MATLAB

Question

I tried to generate random points inside sphere by using following commands with center at origin by using following code

no_of_spots = 3000
radius=20
rvals = (2)*rand(no_of_spots,1)-(1);
elevation = asin(rvals);
azimuth = 2*pi*rand(no_of_spots,1);
radii = rand(no_of_spots,1)*radius;
[point_x,point_y,point_z] = sph2cart(azimuth,elevation,radii);

I got results as

From figure it is random points are concentrated near origin, ie looking like Gaussian distribution. I need random points distributed inside sphere uniformly or need to shift concentration of points from center to another point. How can I do that/ Can anyone help/ Thanks in advance, Manu

Possible duplicate of Sampling uniformly distributed random points inside a spherical volume — buzjwa
– buzjwa, Commented Oct 19, 2017 at 7:45

frslm · Accepted Answer · 2017-10-19 02:37:54Z

5

If you tweak your radii line from:

radii = rand(no_of_spots,1)*radius;

To:

radii = (rand(no_of_spots,1).^(1/3))*radius;

You should get a more uniform-looking distribution.

This is what Knuth described in The Art of Computer Programming. Vol. 2 and is referenced here.

edited Oct 19, 2017 at 2:37

answered Oct 19, 2017 at 2:23

frslm

2,9773 gold badges16 silver badges26 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

How to generate random points inside sphere using MATLAB

1 Answer 1

Comments

Your Answer

Linked

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Linked

Related