分享一个Matlab老祖宗520代码，有需要可以一用。
一、静态心形
代码如下：
%静态爱心clear allclcf=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;[x,y,z]=meshgrid(linspace(-3,3));val=f(x,y,z);[p,v]=isosurface(x,y,z,val,0);patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat');view(3);grid on;axis equal;
二、动态心形
代码如下：
%动态爱心clear; clc; close all; % NOTICE: Your MATLAB version should be at least R2019a !
% --------------------------- functions -------------------------------f = @(x, y, z)(x.^2 + 2.25*y.^2 + z.^2 - 1).^3 - ... x.^2.* z.^3 - 0.1125*y.^2.*z.^3;g = @(x, y, z)(sqrt(x.^2+y.^2)-2.5).^2 + z.^2 - 0.4^2;
% -------------------------- generate data ----------------------------t = linspace(-5, 5);
[x1, y1, z1] = meshgrid(t);[x2, y2, z2] = meshgrid(t);
val1 = f(x1, y1, z1);val2 = g(x2, y2, z2);
[p1, v1] = isosurface(x1, y1, z1, val1, 0);[p2, v2] = isosurface(x2, y2, z2, val2, 0);
% --------------------------- basic plot ------------------------------figure()subplot(1, 1, 1)
h = patch('faces',p1,'vertices',v1,'facevertexcdata',jet(size(v1,1)),... 'facecolor','w','edgecolor','flat'); hold on;patch('faces',p2,'vertices',v2,'facevertexcdata',jet(size(v2,1)),... 'facecolor','w','edgecolor','flat');grid on; axis equal; axis([-3,3,-3,3,-1.5,1.5]); view(3)title(["程式解说",... "给个关注呗！"],'Interpreter','latex','position',[3.3,4])warning('off');% 大声说出来！% 大声说出来！T = suptitle("$I\ Love\ U\ !$"); % 大声说出来！% 大声说出来！set(T,'Interpreter','latex','FontSize',24)
% -------------------------- generate gif -----------------------------pic_num = 1;for i = 1:20 v1 = 0.98 * v1; set(h, 'vertices', v1); drawnow; F = getframe(gcf); I = frame2im(F); [I,map]=rgb2ind(I,256); if pic_num == 1 imwrite(I,map,'BeatingHeart.gif','gif','Loopcount',inf,'DelayTime',0.05); else imwrite(I,map,'BeatingHeart.gif','gif','WriteMode','append','DelayTime',0.05); end pic_num = pic_num + 1;endfor i = 1:20 v1 = v1 / 0.98; set(h, 'vertices', v1); drawnow; F = getframe(gcf); I = frame2im(F); [I,map] = rgb2ind(I,256); imwrite(I,map,'BeatingHeart.gif','gif','WriteMode','append','DelayTime',0.05); pic_num = pic_num + 1;end
% ----------------------------- dynamic -------------------------------while true for i = 1:20 v1 = 0.98 * v1; set(h, 'vertices', v1); drawnow; end for i = 1:20 v1 = v1 / 0.98; set(h, 'vertices', v1); drawnow; endend
（获取更多知识，前往wx 公z号 程式解说）
原文来自https://mp.weixin.qq.com/s/o3uTht1I-QsHY76TjNQMsQ
