How to generate uniformly distributed points inside the volume of frustrum with base radius R and tip radius r and with a height of h.

Question

anuradha verma 2022-11-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1849073-how-to-generate-uniformly-distributed-points-inside-the-volume-of-frustrum-with-base-radius-r-and-ti

评论： anuradha verma 2022-11-12

I want to generate points inside the volume of a frustrum which is at a certain height h1 from (0,0,0) and the axis of frustrum is parallel to z-axis.

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

请先登录，再进行评论。

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

Answer 1

Bruno Luong 2022-11-12

编辑：Bruno Luong 2022-11-12

在 MATLAB Online 中打开

Here we go

r = 1;

R = 2; % must be > r

h = 3;

if R <= r

error('Non valid parameter')

end

Zmin = r*h/(R-r); % Position where the frustrum starts

Zmax = R*h/(R-r); % Position where the frustrum ends

N = 1e4; % Number of point

zr = (Zmin/Zmax).^3;

z = (zr + (1-zr)*rand(1,N)).^(1/3);

rho = R*sqrt(rand(1,N)).*z;

theta = (2*pi)*rand(1,N);

x = rho.*cos(theta);

y = rho.*sin(theta);

z = Zmax*(1-z);

scatter3(x,y,z,'.');

axis equal

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

anuradha verma 2022-11-12

Thank you so much sir for your help. Iam highly oblighed.

请先登录，再进行评论。

Answer 2

Image Analyst 2022-11-12

在 MATLAB Online 中打开

OK, but do you have a question? I'm sure you used

n = 10000; % Whatever.  It's the number of points
x = 2*R*rand(n, 1)
y = 2*R*rand(n, 1);
z = h1 * rand(n, 1);

to generate uniformly distributed points in the rectangular volume. And then you probably threw out points outside the frustrum (truncated cone) volume by looking at each point's radius and the radius of the frustrum at that z level. But what is your question? Is this homework, or is there a real world use case for this?

How do I get help on homework questions on MATLAB Answers? - MATLAB Answers - MATLAB Central