MATLAB之心形图绘制

一、静态心形图绘制

（1）效果展示

（2）静态心形原始代码

 clc;
 clear all;
 const=;
 % 均布三位坐标
 x=-:0.05:;
 y=-:0.05:;
 z=-:0.05:;
 [x,y,z]=meshgrid(x,y,z); % 绘制三位坐标点
 % 心形函数
 f=(x.^ + (/)*y.^ + z.^ - ).^ - x.^.*z.^ - (/)*y.^.*z.^-const;
 p=patch(isosurface(x,y,z,f,));  % 连接各点，组成封闭的图形
 set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
 daspect([  ])
 view()
 camlight; lighting phong

二、动态心形图绘制

（1）效果展示

(2) 动态心形源代码

 clc;
 clear;
 filename='heart';
 % 三位坐标均布
 [x,y,z]=meshgrid(linspace(-,));
 % 心形函数
 p=(x.^+(/)*y.^+z.^-).^-x.^.*z.^-(/)*y.^.*z.^;
 [faces,verts,colors] = isosurface(x,y,z,p,,x);
 % 循环绘制心形形成的过程
 for i=:
 figure(i)
 pp=patch('Faces',faces(:+i*,:),'Vertices',verts);
 set(pp,'FaceColor','red','EdgeColor','none');
 view(-,)
 axis off
 axis equal
 axis tight
 camlight
 lighting gouraud
 pause(0.5)
 f(i) = getframe(i);
 imind = frame2im(f(i));
 [imind,cm] = rgb2ind(imind,);
 if i == 

         imwrite(imind,cm,filename,'gif', 'Loopcount',inf,'DelayTime',0.5);
 else
         imwrite(imind,cm,filename,'gif','WriteMode','append','DelayTime',0.5);
 end
 close(i)
 end