Question

matlab如何画出下面的图像 (x^2+9/4*y^2+z^2-1)^3-x^2*z^3-9/80*y^2*z^3=0

Answer 1

<p>f=@(x,y,z)(x.^2+9/4*y.^2+z.^2-1).^3-x.^2.*z.^3-9/80*y.^2.*z.^3</p>
<p>
implicitmesh(f,[-5
5],150)</p>
<p>-----------------------------</p>
<p>function
h=implicitmesh(f,xlimit,ylimit,zlimit,gd)</p>
<p>%implicitmesh(f,span,gd):画隐函数曲面f(x,y,z)=0的网格图,</p>
<p>%
各坐标范围均限定在span=[lb,ub],</p>
<p>%
网格数为gd，默认为25</p>
<p>if
nargin==2</p>
<p>
ylimit=xlimit;zlimit=xlimit;gd=25;</p>
<p>elseif
nargin==3</p>
<p>
gd=ylimit;ylimit=xlimit;zlimit=xlimit;</p>
<p>elseif
nargin==4</p>
<p>
gd=25;</p>
<p>elseif
nargin==5</p>
<p>else</p>
<p>
error('Error
in
input
arguments')</p>
<p>end</p>
<p>x=linspace(xlimit
(1),xlimit
(2),gd);</p>
<p>y=linspace(ylimit(1),ylimit(2),gd);</p>
<p>z=linspace(zlimit(1),zlimit(2),gd);</p>
<p>[x,y,z]=meshgrid(x,y,z);val=f(x,y,z);</p>
<p>[f,v]=isosurface(x,y,z,val,0);</p>
<p>if
isempty(f)</p>
<p>
warning('There
is
no
graph
in
the
range.');</p>
<p>
p=[];</p>
<p>else</p>
<p>
newplot;</p>
<p>
p=patch('Faces',f,'Vertices',v,'CData',v(:,3),'facecolor','w','EdgeColor','flat');</p>
<p>
isonormals(x,y,z,val,p);view
(3);grid
on</p>
<p>end</p>
<p>if
nargout==0</p>
<p>else</p>
<p>
h=p;</p>
<p>end</p>
<p></p>