How to plot a surface curve with inequalities?

Question

Aditya Zade 2023-5-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1975784-how-to-plot-a-surface-curve-with-inequalities

评论： Aditya Zade 2023-5-30

采纳的回答： Sulaymon Eshkabilov

M = 0.85 : 1/200 : 1 ;

Psi = 0 : 1/200 : 30 ;

dp = 0.6 : 1/200 : 1 ;

[ X, Y, Z ] = meshgrid( dp, Psi, M ) ;

ineq = ( X > 1 - ( 4 * Y / 360 ) ) & ( X > (2/pi)*asin(Z/2) ) & ( X < (2/pi)*asin(Z) ) ;

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

请先登录，再进行评论。

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

Answer 1

Sulaymon Eshkabilov 2023-5-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1975784-how-to-plot-a-surface-curve-with-inequalities#answer_1247184

在 MATLAB Online 中打开

Note that 1st of all, in order to compare X vs. Y vs. Z, their size "MUST" match. Therefore, to make their size compatible, e.g. 200.

Here is one of the possible solutions of this exercise:

N = 200;

M = linspace(0.85 ,1, N) ;

Psi = linspace(0,30, N) ;

dp = linspace(0.6,1, N) ;

[ X, Y] = meshgrid(dp, Psi) ;

[~, ZZ] = meshgrid(M, M);

IND = ( X > 1 - ( 4 * Y / 360 )) & ( X > (2/pi)*asin(ZZ/2) ) & ( X < (2/pi)*asin(ZZ)) ;

Xs = X.*(IND);

Ys = Y.*(IND);

Zs = ZZ.*(IND);

meshc(Xs, Ys, Zs)

xlabel('dp')

ylabel('psi')

zlabel('M')

Check your specified inequality conditions whether they are correct or not.

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

Aditya Zade 2023-5-30

Hi Sulaymon,

Thanks for the answer. But it seems the shape of the plot is wrong. The inequalities are correct in your script. For your reference, the plot drawn by John is correct. The plot should only be valid between M = 0.87 to 1.

请先登录，再进行评论。

Answer 2

John D'Errico 2023-5-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1975784-how-to-plot-a-surface-curve-with-inequalities#answer_1247174

编辑：John D'Errico 2023-5-30

在 MATLAB Online 中打开

And exactly what is the problem? I think what you do not understand is, what you have created is NOT a surface.

M = 0.85 : 1/200 : 1 ;

Psi = 0 : 1/200 : 30 ;

dp = 0.6 : 1/200 : 1 ;

[ X, Y, Z ] = meshgrid( dp, Psi, M ) ;

ineq = ( X > 1 - ( 4 * Y / 360 ) ) & ( X > (2/pi)*asin(Z/2) ) & ( X < (2/pi)*asin(Z) );

X(~ineq) = NaN;

Y(~ineq) = NaN;

Z(~ineq) = NaN;

plot3(X(:),Y(:),Z(:),'.')

The result is a set of points that live in a triangular region in X and Y, with one curved surface, where Z is a function of X and Y.

But there is no "surface" that you have created. You might think of what you have done is to create a volume, the set of points delineated by those constraints.

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

Aditya Zade 2023-5-30

编辑：Aditya Zade 2023-5-30

Hi John,

The problem here is its hard to understand the shape of the 3d set of points. Is there a way to have boundaries or change of color of the set of points as the value of M (Z axis) increases? Something similar to what we have in a surface plot.

请先登录，再进行评论。