用 MATLAB情人节表白% volume data [X,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 figure hFig = 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 axes hAxes = 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 heart patch(isosurface(F,-1e-3), 'FaceColor','w', 'EdgeColor','none') hidden on % hidden surface removal % contours in the y-z plane for 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 plane for 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 plane for 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 % text props = {'FontWeight','bold', 'FontAngle','italic', 'FontSize',100}; pause(.2) text(7,50,70, 'I', props{:}) pause(.5) text(80,50,43, 'YOU', props{:}) pause(.2) % equation props = {'FontSize',10, 'Interpreter','latex'}; text(30,45,30, 'I LOVE YOU', props{:}); drawnow
