代码如下,时间紧促,下次详细说明
syms x y z;
f=(x^2+9/4*y^2+z^2-1)^3-x^2*z^3-9/80*y^2*z^3;
f=matlabFunction(f);
[x,y,z] = meshgrid(-1.5:.02:1.5,-1:.02:1,-1.5:.02:1.5); % 画图范围,设置[-1.5,1.5]区间
v = f(x,y,z);
h = patch(isosurface(x,y,z,v,0)); %创建补片对象的句柄
isonormals(x,y,z,v,h) %计算法向量
set(h,'FaceColor','r','EdgeColor','None');
alpha(1) %设置透明度
grid on; view([1,1,1]); axis equal; camlight; lighting gouraud %设置光照
axis off;
pos1=get(gca,'position')
pos2=pos1;
pos2(2)=pos2(2)+0.01;
pos2(1)=pos2(1)-0.03;
pos2(3)=pos2(3)+0.08;
pos2(4)=pos2(4)+0.08;
for ii=1:10
pause(1)
set(gca,'position',pos1)
pause(0.1)
set(gca,'position',pos2)
pause(0.1)
set(gca,'position',pos1)
pause(0.1)
set(gca,'position',pos2)
pause(0.1)
set(gca,'position',pos1)
pause(0.1)
set(gca,'position',pos2)
pause(0.1)
set(gca,'position',pos1)
pause(0.1)
set(gca,'position',pos2)
pause(0.1)
set(gca,'position',pos1)
end