圆内均匀分布 matlab,利用matlab中rand函数生成圆和球中均匀分布的数据点

一. rand 函数

rand函数是生成(0,1)之间均匀分布的数据点。

rand(m,n): 生成m行n列的在(0,1)之间的数据点，这些数据服从均匀分布。

例如生成1x2的服从(0,1)之间均匀分布的数据点：

>> rand(1,2)

ans =

0.8147 0.9058

二. 生成半径为r的圆内均匀分布的数据点

%生成圆中均匀分布的随机数据点

angle=rand(1,1000)*2*pi; %(0,2*pi)之间均匀分布数据点

r=sqrt(rand(1,1000));%(0,1)之间r^2均匀分布数据点

x=r.*cos(angle);

y=r.*sin(angle);

figure

plot(x,y,'r.');

axis square

三. 生成半径为r的球内均匀分布数据点

%生成球中均匀分布的随机数据点

angle1=rand(1,1000)*2*pi;

angle2=acos(rand(1,1000)*2-1);;

r=power(rand(1,1000),1/3);

x=r.*cos(angle1).*sin(angle2);

y=r.*sin(angle1).*sin(angle2);

z=r.*cos(angle2);

figure

plot3(x,y,z,'r.');

axis square;

四. 注意事项

1. 在圆中生成均匀分布的随机数据点的时候，半径需要开平方，而不能直接rand生成随机半径值。这是因为均匀分布在这种情形下是针对圆的面积而言，在极坐标下，面积元是,通过积分

不变。面积元是均匀分布的，从而

服从均匀分布，并不是

不服从均匀分布，故开平方。而角度

仍服从

均匀分布，不用开平方。

2. 在球中生成均匀分布随机数据点，半径需要开3次方，angle2用到反余弦。这是由于均匀分布此时指的是体积，即在球中任取一点的概率都相等。在球坐标下，球的体积元是

,通过重积分得体积，

变成

服从均匀分布，故开三次方，

变成

，所以

服从均匀分布，生成rand(-1,1)再反余弦即可的

。