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;
xlimit=[-1.5 1.5];
ylimit=[-.8 .8];
zlimit=[-1.5 1.5];
x=linspace(xlimit(1),xlimit(2),50);
y=linspace(ylimit(1),ylimit(2),50);
z=linspace(zlimit(1),zlimit(2),50);
[x,y,z]=meshgrid(x,y,z);val=f(x,y,z);
[f,v]=isosurface(x,y,z,val,0);
newplot;
p=patch('Faces',f,'Vertices',v,'CData',v(:,3),'facecolor','flat','EdgeColor','k');
isonormals(x,y,z,val,p);view(3);grid on
set(p,'AmbientStrength',.5);
title('Add you title here~~','FontSize',16,'Color','m');
hold on;
f=@(x,y,z)(sqrt(x.^2+y.^2)-2).^2+z.^2-.09;
xlimit=[-2.3,2.3];
ylimit=[-2.3,2.3];
zlimit=[-.3,.3];
x=linspace(xlimit(1),xlimit(2),25);
y=linspace(ylimit(1),ylimit(2),25);
z=linspace(zlimit(1),zlimit(2),25);
[x,y,z]=meshgrid(x,y,z);val=f(x,y,z);
[f,v]=isosurface(x,y,z,val,0);
newplot;
p=patch('Faces',f,'Vertices',v,'CData',v(:,3),'facecolor','flat','EdgeColor','k');
isonormals(x,y,z,val,p);view(3);grid on
colormap hsv;set(p,'AmbientStrength',.8,'FaceAlpha',.5);
axis off;axis equal;shading interp;camlight;lighting gouraud;
view([2 14]);

复制到Matlab命令行>>后面运行即可