爱的方程

520 又是一波秀恩爱的季节~

今天看到一个砰砰跳的红心，看了一下代码，原来爱的公式是这样的：

被称为是爱的方程：

17*x^2-16*abs(x)*y+17*y^2-225=0

这是普遍认为比较好看的一个公式，一般写作：

x^2 - abs(x)*y + y^2 = 1

那么，三维的公式：

f=(x^2+9/4*y^2+z^2-1)^3-x^2*z^3-9/80*y^2*z^3;

代码里这么写的。

想把砰砰跳的代码放到这里。。。放不放呢？

反正不是我写的。。。。。。逃

。。。。。。还是放这儿吧~

syms x y z
f=(x^2+9/4*y^2+z^2-1)^3-x^2*z^3-9/80*y^2*z^3;
f=matlabFunction(f);

[x,y,z] = meshgrid(-1.5:.02:1.5,-1:.02:1,-1.5:.02:1.5);    % 画图范围
v = f(x,y,z);
h = patch(isosurface(x,y,z,v,0));
isonormals(x,y,z,v,h)
set(h,'FaceColor','r','EdgeColor','none');
xlabel('x');ylabel('y');zlabel('z');
alpha(1)
grid on; view([1,1,1]); axis equal; camlight; lighting gouraud
%axis off
pos1=get(gca,'position')
pos2=pos1;
pos2(2)=pos2(2)+0.01;
pos2(1)=pos2(1)-0.03;
pos2(3)=pos2(3)+0.08;
pos2(4)=pos2(4)+0.08;

for ii=1:10
    pause(1)
    set(gca,'position',pos1)
    pause(0.1)
    set(gca,'position',pos2)
    pause(0.1)
    set(gca,'position',pos1)
    pause(0.1)
    set(gca,'position',pos2)
    pause(0.1)
    set(gca,'position',pos1)
    pause(0.1)
    set(gca,'position',pos2)
    pause(0.1)
    set(gca,'position',pos1)
    pause(0.1)
    set(gca,'position',pos2)
    pause(0.1)
    set(gca,'position',pos1)
end