基于matlab的动态心形图案

心形图案的函数式为

$$y=\sqrt[3]{x^2}+0.9\times \sqrt{3.3-x^2}\times sin(a\times \pi \times x)$$

其中x为自变量，a决定图形的形状，故在编写程序时，可以先定义x的范围，然后使用for循环来控制a的数值，配合上延时指令，实现心形不断变化的效果。源程序如下所示：

x=-2:0.01:2;
for a=1:50
y=(x.^2).^(1/3) + (0.9*sqrt((3.3-x.^2))).*sin(a*pi*x);
plot(x,y)
pause(0.1)
end

for a=50:-1:1
y=(x.^2).^(1/3) + (0.9*sqrt((3.3-x.^2))).*sin(a*pi*x);
plot(x,y)
pause(0.1)
end

程序运行后的效果为变化的心形图案。