用MATLAB情人节表白

用MATLAB情人节表白% volume dataX,Y,Z = meshgrid(linspace(-3,3,101);F = -X.2.*Z.3 - (9/80).*Y.2.*Z.3 + (X.2 + (9/4).*Y.2 + Z.2 - 1).3;% initialize figurehFig = figure(Menubar,none, Color,w);pos = get(hFig, Position);set(hFig, Position, pos(1)-0.15*pos(3) pos(2) 1.3*pos(3) pos(4);% initialize axeshAxes = axes(Parent,hFig, DataAspectRatio,1 1 1, . XLim,30 120, YLim,35 65, ZLim,30 75);view(-39,30);axis off% Fill the inside of the mesh with an isosurface to% block rendering of the back side of the heartpatch(isosurface(F,-1e-3), FaceColor,w, EdgeColor,none)hidden on % hidden surface removal% contours in the y-z planefor iX = 35 38 41 45 48 51 54 57 61 64 67 plane = reshape(F(:,iX,:), 101 101); cData = contourc(plane, 0 0); xData = iX.*ones(1,cData(2,1); line(xData, cData(2,2:end), cData(1,2:end), . Color,r, Parent,hAxes) pause(.1)end% contours in the x-z planefor iY = 41 44 47 51 55 58 61 plane = reshape(F(iY,:,:), 101 101); cData = contourc(plane, 0 0); yData = iY.*ones(1,cData(2,1); line(cData(2,2:end), yData, cData(1,2:end), . Color,r, Parent,hAxes) pause(.1)end% contours in the x-y planefor iZ = 36 38 40 42 44 46 48 50 52 54 56 58 60 62 64 66 69 71 plane = F(:,:,iZ); cData = contourc(plane, 0 0); startIndex = 1; if size(cData,2) (cData(2,1)+1) startIndex = cData(2,1)+2; zData = iZ.*ones(1,cData(2,1); line(cData(1,2:(startIndex-1), cData(2,2:(startIndex-1), zData, . Color,r, Parent,hAxes) end zData = iZ.*ones(1,cData(2,startIndex); line(cData(1,(startIndex+1):end), cData(2,(startIndex+1):end), zData, . Color,r, Parent,hAxes) pause(.1)end% textprops = FontWeight,bold, FontAngle,italic, FontSize,100;pause(.2)text(7,50,70, I, props:)pause(.5)text(80,50,43, YOU, props:)pause(.2)%