我在一个程序中验证过,不是一条线,正常,如:
[x,y,z]=meshgrid(linspace(-1.3,1.3));
val=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (1/9)*y.^2.*z.^3;%isosurface(x,y,z,val,0)
p = patch(isosurface(x,y,z,val,0));
isonormals(x,y,z,val,p)
set(p,'FaceColor','red','EdgeColor','none');
axis equal
view(-10,24)%colormap(hot)
axis off
camlight
rotate3d on
colormap(jet)
colorbar('location','EastOutside')
%colorbar('East')
ylabel(colorbar,'Emission intensity (a.u.)','fontsize',18)
set(gca,'fontsize',14,'linewidth',2)