MATLAB画爱心

一个简单的画图程序

f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;
[x,y,z]=meshgrid(linspace(-3,3));
val=f(x,y,z);
[p,v]=isosurface(x,y,z,val,0);
Hnd1=patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat');
f1=@(x,y,z) ((x/2).^2+ (9./4).*(y/2).^2 + (z/2).^2 - 1).^3 - (x/2).^2.*(z/2).^3 - (9./80).*(y/2).^2.*(z/2).^3;
val1=f1(x,y,z);
[p1,v1]=isosurface(x,y,z,val1,0);
Hnd2=patch('faces',p1,'vertices',v1,'facevertexcdata',jet(size(v1,1)),'facecolor','w','edgecolor','flat');
view(3);
grid on;
axis equal;
%set函数对属性进行修改
set(Hnd2,'FaceColor','none','EdgeColor',[0.5,0.5,0.5]);