Draw surface of a football (American)

6 次查看(过去 30 天)
Hi everyone!
I got the following assignment:
Draw the surface of an american football which can be described with the following equation:
within the interval
The surface should look like this:
I have the following code:
r = linspace(-1,1,21);
a = linspace(0,2*pi,63);
[R,A] = meshgrid(r,a);
X = R.*cos(A);
Y = R.*sin(A);
Z = acos(R);
surf(X,Y,Z);
axis equal
But every time I run the code, I get the following surface:
What is it that I am doing wrong?
Thank you for any help that anyone could provide, in advance!

采纳的回答

DGM
DGM 2021-11-5
编辑:DGM 2021-11-5
Two things to notice: note the way your points are spaced unevenly along z compared to the reference image. Also note the restricted domain of acos(). This problem would be easier approached by using z as the parameter instead of r.
z = linspace(-pi/2,pi/2,21);
th = linspace(0,2*pi,63).';
X = cos(z).*cos(th);
Y = cos(z).*sin(th);
Z = repmat(z,numel(th),1);
surf(X,Y,Z);
axis equal
colormap(jet)
If you really wanted to do it your way, you could ...
clf; % reset web-plot
r = linspace(0,1,21);
a = linspace(0,2*pi,63);
[R,A] = meshgrid(r,a);
X = R.*cos(A);
Y = R.*sin(A);
Z = acos(R);
surf(X,Y,Z); hold on
surf(X,Y,-Z)
axis equal
  1 个评论
Szabolcs Simon-Guth
Thank you very much for your guidance! I'm new to MATLAB and I didn't think of the fact that the restricted domain, because I concentrated too much on the coding aspect, instead of the maths. Thank you again for all the help! :)

请先登录,再进行评论。

更多回答(1 个)

Yongjian Feng
Yongjian Feng 2021-11-5
编辑:Yongjian Feng 2021-11-5
Hint: your assignment states z should be -Pi/2 < z < Pi/2.
But your r range is (-1, 1), which corresponds to z range 0<z< Pi
In other words, you plot the wrong range of Z. The first plot has Z from about -1.5 to 1.5.
Your plot has Z between 0 < Z < 3.
  1 个评论
Szabolcs Simon-Guth
Alright! I understand. So, I'm going to include these parameters as well. Thank you very much for your help! I did not think of that.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Annotations 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by