matlab繪制動態(tài)三維心形代碼(蛋疼的情人節(jié)奉獻)

Matlab繪制三維動態(tài)心形It’sOKtosendapictoyourgirlfriendonValentine'sDay情人節(jié)蛋疼玩意效果圖：原始代碼：%僅供參考，自助修改，原則上自己動手，要是非常強烈的要幫忙%可以聯(lián)系我的QQ865802870，但愿我還在上面.Sourcecode:%構(gòu)造體積方程和坐標軸，畫出圖形;linspace(a,b,c)均勻生成介于a到b的c個值,c的默認為100

。Meshgrid生成矩陣網(wǎng)格。[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',[111],...'XLim',[30120],'YLim',[3565],'ZLim',[3075]);view([-3930]);axisoff%制作出動態(tài)的隱形效果;hiddenon%畫出網(wǎng)格，制作網(wǎng)格動態(tài)效果;%快渲染心得背面:p=patch(isosurface(F,-0.001));set(p,'FaceColor','w','EdgeColor','w');%構(gòu)造Y-Z平面,,描完函數(shù)在該平面的點:foriX=[3538414548515457616467]plane=reshape(F(:,iX,:),101,101);cData=contourc(plane,[00]);xData=iX.*ones(1,cData(2,1));plot3(hAxes,xData,cData(2,2:end),cData(1,2:end),'r');pause(.1),drawnowend%構(gòu)造X-Z平面,描完函數(shù)在該平面的點:foriY=[41444751555861]plane=reshape(F(iY,:,:),101,101);cData=contourc(plane,[00]);yData=iY.*ones(1,cData(2,1));plot3(hAxes,cData(2,2:end),yData,cData(1,2:end),'r');pause(.1),drawnowend%構(gòu)造X-Y平面,描完函數(shù)在該平面的點:foriZ=[363840424446485052545658606264666971]plane=F(:,:,iZ);cData=contourc(plane,[00]);startIndex=1;ifsize(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');endzData=iZ.*ones(1,cData(2,startIndex));plot3(hAxes,cData(1,(startIndex+1):end),...cData(2,(startIndex+1):end),zData,'r');pause(.1),drawnowend%給三維心著色isosurface(x,y,z,val,0);%View視角colormap色圖0-1axisequal;view(-10,10);colormap([10.20.2])復(fù)制代碼%花心;clear;clc;closeallc=5;t=linspace(-c,c);[x,y]=meshgrid(t);z=17*x.^2-16*abs(x).*y+17*y.^2-225;pcolor(x,y,z);shadinginterppause(2);spinmap(10)復(fù)制代碼%構(gòu)造體積方程和坐標軸，畫出圖形;[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',[111],...'XLim',[30120],'YLim',[3565],'ZLim',[3075]);view([-3930]);axisoff%制作出動態(tài)的隱形效果;hiddenon%畫出網(wǎng)格，制作網(wǎng)格動態(tài)效果;%快渲染心得背面:p=patch(isosurface(F,-0.001));set(p,'FaceColor','w','EdgeColor','w');%構(gòu)造Y-Z平面,,描完函數(shù)在該平面的點:foriX=[3538414548515457616467]plane=reshape(F(:,iX,:),101,101);cData=contourc(plane,[00]);xData=iX.*ones(1,cData(2,1));plot3(hAxes,xData,cData(2,2:end),cData(1,2:end),'k');pause(.1),drawnowend%構(gòu)造X-Z平面,描完函數(shù)在該平面的點:foriY=[41444751555861]plane=reshape(F(iY,:,:),101,101);cData=contourc(plane,[00]);yData=iY.*ones(1,cData(2,1));plot3(hAxes,cData(2,2:end),yData,cData(1,2:end),'k');pause(.1),drawnowend%構(gòu)造X-Y平面,描完函數(shù)在該平面的點:foriZ=[363840424446485052545658606264666971]plane=F(:,:,iZ);cData=contourc(plane,[00]);startIndex=1;ifsize(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(