首页分享利用MATLAB绘制各种表白心形

利用MATLAB绘制各种表白心形

来源：花匠小妙招时间：2024-12-01 08:28

利用MATLAB绘制各种表白心形♥

跳动三维心形

表白情人节必备
面向工科屌丝男

固定心形的各种画法见下面连接
链接: https://blog.csdn.net/weixin_44044161/article/details/115342172.

在这里插入图片描述

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(["$(x^2+frac{9}{4}y^2+z^2-1)^3-x^2z^3-frac{9}{80}y^2z^3=0$",... "$(sqrt{x^2+y^2}-R)^2 +z^2 = r^2$"],'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; end for 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; end end % ------------------------------- end --------------------------------- 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879