你是否准备好表白了呢?可在好多人的眼里我们理工男是无聊、木讷、不善交往的代名词,尤其是在跟女孩子相处过程中总被埋怨不解风情,可是我们理科生也有爱啊,只不过我们的表达方式很有特点而已。本文就为理工男支招表爱神器。
%构造体积方程和坐标轴,画出图形;[X,Y,Z] = meshgrid(linspace(-3,3,101));%3D心型图方程如下;F = -X.^2.*Z.^3-(9/80).*Y.^2.*Z.^3+(X.^2+(9/4).*Y.^2+Z.^2-1).^3;hFigure = figure;sz = get(hFigure, 'Position');set(hFigure, 'Position', [sz(1)-0.15*sz(3) sz(2) 1.3*sz(3) sz(4)]);set(hFigure,'color','w', 'menu','none')hAxes = axes('Parent',hFigure,'NextPlot','add', 'DataAspectRatio',... [1 1 1], 'XLim',[30 120],'YLim',[35 65],'ZLim',[30 75]);view([-39 30]);axis off% 制作出动态的隐形效果;hidden on% 画出网格,制作网格动态效果;% 快渲染心得背面:p = patch(isosurface(F,-0.001));set(p,'FaceColor','w','EdgeColor','w');% 构造Y-Z平面,,描完函数在该平面的点: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)); plot3(hAxes,xData,cData(2,2:end),cData(1,2:end),'r'); pause(.1), drawnowend% 构造X-Z平面,描完函数在该平面的点: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)); plot3(hAxes,cData(2,2:end),yData,cData(1,2:end),'r'); pause(.1), drawnowend% 构造X-Y平面,描完函数在该平面的点: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)); plot3(hAxes,cData(1,2:(startIndex-1)),... cData(2,2:(startIndex-1)),zData,'r'); end zData = iZ.*ones(1,cData(2,startIndex)); plot3(hAxes,cData(1,(startIndex+1):end),... cData(2,(startIndex+1):end),zData,'r'); pause(.1), drawnowend%函数已经画完,接下来为文字部分;pause(.2)%设置字体大小,粗细,位置等,以下是打出I (心型图) YOU;text(7,50,70,'I', 'fontWeight','bold','FontAngle','italic','fontsize',60)pause(.5)text(80,50,43,'You', 'fontWeight','bold','FontAngle','italic','fontsize',60)pause(.2)line([20 80],[50 50],[52.5 52.5], 'color','w')line([50 50],[20 80],[52.5 52.5], 'color','w')line([50 50],[50 50],[30 80], 'color','w')%制作者签名;text(30,60,30,'Made By 天天Matlab 2020/01/04', 'fontsize',10)text(35,45,30,'', 'fontsize',8);
f=@(x,y,z)(x.^2+(9./4).*y.^2+z.^2-1).^3-x.^2.*z.^3-(9./80).*y.^2.*z.^3;[x,y,z]=meshgrid(linspace(-1.5,1.5));val=f(x,y,z);isosurface(x,y,z,val,0);axis equal;view(30,0);colormap([1 0.20 0.2])axis off
scrsz = get(0,'ScreenSize'); h=figure('Position',[scrsz(3)*0.1 scrsz(4)*0.1 scrsz(3)*0.8 scrsz(4)*0.8]); x=-20:0.1:-18; y=0*x; plot(x,y,'y'); axis([-20,24,-15,10]) axis off; title('\fontsize{18}\rm To My Love') hold on; for i=0:1:400 a=-18+i*0.1; b=-18+i*0.1+1; x=a-0.01:0.01:b-0.01; y=3*(x-16).*(x-18)./abs((x-16).*(x-18))-sqrt((4-(abs(x-16)-abs(x-18)).^2)./16)-sqrt((36-(abs(x-14)-abs(x-20)).^2)./16)+3/4*(x+14).*(x+16)./abs((x+14).*(x+16))+abs(x)-abs(x/2+2)-abs(x/2-2)+(abs(x/2+4)+abs(x/2-4)-abs(x/2+2)-abs(x/2-2)+4.5-4.5*(x+14).*(x+16)./abs((x+14).*(x+16))-3*(x-16).*(x-18)./abs((x-16).*(x-18))+sqrt((4-(abs(x-16)-abs(x-18)).^2)./16)+3*(x-14).*(x-20)./abs((x-14).*(x-20))-sqrt((36-(abs(x-14)-abs(x-20)).^2)./16)).*sin(3*pi*x)-15/4; t='r'; plot(x,y,t); pause(0.05); end x=22:0.1:24; y=0*x; plot(x,y,'y'); hold off;