为你心动——Matlab代码

效果图

n = 100;
t = linspace(0, 2*pi, n);
x = 16*sin(t).^3;
y = 13*cos(t) - 5*cos(2*t) - 2*cos(3*t) - cos(4*t);
handle = scatter(x, y, '.m');
ax = gca;
ax.Color = [0, 0, 0];
ax.XColor = 'none';
ax.YColor = 'none';
hold on;
handle0 = scatter(x, y, '.m');
hold off;
n = 5000;
t = linspace(0, 2*pi, n);
r = 10 - abs(normrnd(10, 1, [1, n]) - 10);
x = 16*sin(t).^3;
y = 13*cos(t) - 5*cos(2*t) - 2*cos(3*t) - cos(4*t);
x = r.*x;
y = r.*y;
ax.XLim = [-200, 200];
ax.YLim = [-200, 150];
handle.XData = x;
handle.YData = y;
handle0.XData = 1.2*x;
handle0.YData = 1.2*y;
n = 100;
t = linspace(0, 2*pi, n);
x = 16*sin(t).^3;
y = 13*cos(t) - 5*cos(2*t) - 2*cos(3*t) - cos(4*t);
handle = scatter(x, y, '.m');
ax = gca;
ax.Color = [0, 0, 0];
ax.XColor = 'none';
ax.YColor = 'none';
hold on;
handle0 = scatter(x, y, '.m');
hold off;
n = 5000;
t = linspace(0, 2*pi, n);
r = 10 - abs(normrnd(10, 1, [1, n]) - 10);
x = 16*sin(t).^3;
y = 13*cos(t) - 5*cos(2*t) - 2*cos(3*t) - cos(4*t);
x = r.*x;
y = r.*y;
ax.XLim = [-200, 200];
ax.YLim = [-200, 150];
handle.XData = x;
handle.YData = y;
handle0.XData = 1.2*x;
handle0.YData = 1.2*y;
len = 2000;
tt = linspace(0, 2*pi, len);
for i = 1:len
    r0 = 10 - abs(normrnd(10, 2, [1, len]) - 10);
    x0 = r0 * 16.*sin(tt).^3;
    y0 = r0 .* (13*cos(tt) - 5*cos(2*tt) - 2*cos(3*tt) - cos(4*tt));
    handle0.XData = 1.2*x0;
    handle0.YData = 1.2*y0;
    d = sqrt(x.^2 + y.^2);
    l = (d - min(d)) / (max(d) - min(d));
    c = 1 + (2 - 1.8*l)/8.*sin(80*tt(i));
    handle.XData = c.*x;
    handle.YData = c.*y;    
    drawnow;
    pause(0.01);
end