心用matlab画出来的,用matlab画3D心(surf,mesh)

% 1.核心程序%

%designed by CONST from BUAA%

% 可根据电脑性能更改参数

th=linspace(-pi,pi,100);

ph=linspace(0,2*pi,100);

r=linspace(0,1.5,1000);

[t,p]=meshgrid(th,ph);

f=zeros(1,1000);

sf=zeros(100,100);

for i=1:100

for j=1:100

x=sin(ph(i))*cos(th(j));

y=cos(ph(i))*cos(th(j));

z=sin(th(j));

f=((x*r).^2 + (9/4)*(y*r).^2 + (z*r).^2 - 1).^3 - (x*r).^2.*(z*r).^3 - (9/80)*(y*r).^2.*(z*r).^3; [minf,ind]=min(abs(f));

sf(i,j)=r(ind);

end

end

[x,y,z]=sph2cart(p,t,sf);

% 2.mesh红心%

mesh(x,y,z)

axis equal

Colormap([1,0,0])

% 3.surf红心%

surf(x,y,z)