简单的思考（一）：MATLAB实现心形线

今天刷B站的时候看见了：
在这里插入图片描述
于是想着自己能不能也做出来

clc;clear;
% 定义x的范围
x = -2:0.01:2;% 初始化图形
figure;
set(gcf,'position',[0,0,800,600],'color','w');
h1 = plot(x,abs(x).^(2/3) + (0.9*sqrt((3.3-x.^2))).*sin(0.1*pi*x),'r','LineWidth',3);
hold on
xlim([-2 2]);
ylim([-1.5 2.5]);
title('\color{red}f(x)=x^2^/^3+e/3*(π-x^2)^1^/^2*sin(a*π*x)','fontsize',20)
str_title= strcat ('\color{red}a=',num2str(0.1));
text(-0.2,2.3,str_title,'FontName','Times New Roman', 'FontSize',20);filename = '心形线.gif';  % GIF文件名% 创建动画
for a = 0.2:0.1:20% 更新数据y1 = abs(x).^(2/3) + (0.9*sqrt((3.3-x.^2))).*sin(a*pi*x);% 删除旧的文本对象delete(findobj('Type','text'));% 添加新的文本str_title= strcat ('\color{red}a=',num2str(a));text(-0.2,2.3,str_title,'FontName','Times New Roman', 'FontSize',20);% 更新图形set(h1, 'YData', real(y1));% 更新图形窗口drawnow;% 捕获当前帧frame = getframe(gcf);im = frame2im(frame);[imind,cm] = rgb2ind(im,256);% 写入GIF文件if a == 0.2imwrite(imind,cm,filename,'gif', 'Loopcount',inf, 'DelayTime',0.01);elseimwrite(imind,cm,filename,'gif','WriteMode','append', 'DelayTime',0.01);end% 暂停pause(0.01);
end% 关闭hold状态
hold off;