大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全系列IDE稳定放心使用
昨天发表了一篇用python教你画心形图表白的文章：
想要表白的看这里，教你用python画不同类型的心形图虏获芳心，值得收藏！！
里面详细介绍了各种心形图的画法以及最终的表白神器，值得点赞收藏！！
同样matlab也可以实现相同的功能并且还可以做得更好，今天就用教你用matlab画动态心形曲线图，不信请看下面：
matlab画动态心形曲线图
利用数学上的格式f(x)=x^2^/^3+e/3*(π-x^2)^1^/^2*sin(a*π*x)可以为你展现非常完美的图像。
% matlab动态心形曲线（保存GIF格式）
clc;clear all;
a=10;
x=-2:0.01:2;
figure(1);
set(gcf,'position',[0,0,800,600],'color','w');
for i=1:100
str_title = strcat('\color{red}a=',num2str(a));
y = abs(x).^(2/3)+(0.9*sqrt((3.3-x.^2))).*sin(a*pi*x);
a=i/10;
figure(i);
set(gcf,'position',[0,0,800,600],'color','w');
p1 = plot(x,y,'r','LineWidth',3);
hold on;
title('\color{red}f(x)=x^2^/^3+e/3*(π-x^2)^1^/^2*sin(a*π*x)','fontsize',15);
text(-0.2,2.3,str_title,'FontName','Times New Roman','FontSize',20);
xlim([-2 2]);
ylim([-1.5 2.5]);
frame = getframe(gcf);
im = frame2im(frame);
[I,map] = rgb2ind(im,256);
if i==1
imwrite(I,map,'心形图.gif','gif','Loopcount',Inf,'DelayTime',0.05);
else
imwrite(I,map,'心形图.gif','gif','WriteMode','append','DelayTime',0.05);
end
close(figure(i));
end
成果图：
这是其中的一张（a=4.6）
怎么样是不是很震撼，喜欢的话感觉点赞收藏起来吧
matlab画动态心形曲线图（基础版）
clc;clear all;
t=0:0.1:2*pi;
x=16*sin(t).^3;
y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);
figure;
h = plot(x(1),y(1),'r');
axis([-20,20,-20,20]);
axis off
[A,map] = rgb2ind(frame2im(getframe),256);
imwrite(A,map,'love.gif','LoopCount',inf,'DelayTime',0.1);
for i = 2:length(t)
h.XData(i) = x(i);
h.YData(i) = y(i);
drawnow
title('\color{red}I LOVE YOU','fontsize',15);
%在本地保存的图片也能显示出来，保存地址默认
[A,map] = rgb2ind(frame2im(getframe),256);
imwrite(A,map,'love.gif','WriteMode','append','DelayTime',0.1);
end
matlab画3D心形图
f=@(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;
f=@(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(-1.5,1.5));
val=f(x,y,z);
isosurface(x,y,z,val,0);
axis equal;view(3);colormap([1 0.2 0.2])
备注
1. 觉得喜欢的小伙伴们来个三连击（点赞+关注+收藏）；
2. 让更多的人可以看到这篇文章；
3. 最重要的是希望可以帮助那些羞涩的程序猿们表白成功。
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/181019.html原文链接：https://javaforall.net
