Meshgrid multi-valued function data

Question

Ahmed Zankoor 2018-6-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/405304-meshgrid-multi-valued-function-data

编辑： Ahmed Zankoor 2018-6-12

If I have data z = f(x,y) which have at some (or all) (x,y) coordinates more than a value for z. How can I represent the corresponding z value for the X and Y mesh grid, where [X Y] = meshgrid(x,y). I also wonder how this is handled in sphere function. For example: >> [x y z] = sphere(50) , surf(x,y,z) Since x,y,z meshgrids are not in the normal format of mesh grids?

2 个评论
显示无隐藏无

Rik 2018-6-12

You can best think of meshgrid as returning coordinate grids. That means that you can use it as an input to a function that returns an non-scalar, but you'll have to use cellfun or a similar function.

sphere does not use meshgrid internally. The reason the vector outputs of sphere are not in a meshgrid format seems simple: you'dd have way too many points that are useless. You can have a peek by typing edit sphere, although I'dd urge you never to edit internal functions.

Ahmed Zankoor 2018-6-12

编辑：Ahmed Zankoor 2018-6-12

在 MATLAB Online 中打开

Thank you Rik, I got your point for sphere. What I can not understand is the meshgrid part, how can we express (x1,y1,z1) and (x1,x2,z2) using the same mishgrid X,Y matrices. Take this for example:

{radius = [3];
x = linspace(-radius,radius,100);
y = linspace(-radius,radius,100);
[X Y] = meshgrid(x,y);
Z = - sqrt(radius^2 - X.^2 - Y.^2); 
Z(imag(Z)~= 0) = 0 ;
surf(X,Y,Z)
Z(imag(Z)~= 0) = 0 ;
surf(X,Y,Z)
end}

What if I want to surf both halves of the sphere ??

请先登录，再进行评论。

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

Answer 1

Rik 2018-6-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/405304-meshgrid-multi-valued-function-data#answer_324405

在 MATLAB Online 中打开

You should use hold on to plot multiple parts. The code below is a fix to your code.

radius=3;
x = linspace(-radius,radius,100);
y = linspace(-radius,radius,100);
[X,Y] = meshgrid(x,y);
Z = sqrt(radius^2 - X.^2 - Y.^2);
Z(imag(Z)~= 0) = 0 ;
%set coordinates outside of the sphere to the edge
L=(X.^2+Y.^2)>(radius^2);
phi=atan2(Y,X);
X(L)=radius*cos(phi(L));
Y(L)=radius*sin(phi(L));
%put the two halves together
Z = [Z,-Z];
X = [X,X];
Y = [Y,Y];
figure(1),clf(1)
surf(X,Y,Z)

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

Ahmed Zankoor 2018-6-12

编辑：Ahmed Zankoor 2018-6-12

Thank you so much Rik. So this is my point, meshgrid 2D matrices must be single-valued (single value of z corresponding to each (x,y)), otherwise, we should use multiple surfaces (different Z matrices for the same X and Y).

请先登录，再进行评论。