matlab绘制‘心’形-动画版的

functiondemoDisplayclearclccloseall;[x,y,z]=heart;bottomZ=min(z(:));radius=abs(bottomZ);axisvis3daxisoffdaspect([1.6,1,1.875]);holdon;markerColor=[0.5,0.5,0.5];lineColor=[0.7,0.7,0.7];figWidth=800figHeight=600;screenSize=get(0,'ScreenSize');X0=(screenSize(3)-figWidth)/2;Y0=(screenSize(4)-figHeight)/2;set(gcf,'Position',[X0,Y0,figWidth,figHeight]);camtarget([0,0,0]);set(gcf,'Renderer','zbuffer');set(gcf,'DoubleBuffer','on');%RGB(255105180)isforpinkcolorset(gcf,'color',[255255255]/255);startAZ=30;satrtEL=20;view(startAZ,satrtEL);%drawpointsforiLoop=1:size(x,1)plot3(x(iLoop,:),y(iLoop,:),z(iLoop,:),'Marker','.','MarkerEdgeColor',markerColor,'LineStyle','none');pause(0.2);end%drawwireframesforiLoop=1:size(x,1)plot3(x(iLoop,:),y(iLoop,:),z(iLoop,:),'color',lineColor,'LineStyle','-');pause(0.1);endforiLoop=1:size(x,2)plot3(x(:,iLoop),y(:,iLoop),z(:,iLoop),'color',lineColor,'LineStyle','-');pause(0.1);end%rotatewireframefori=1:36camorbit(10,0,'data');pause(0.1);end%cleanthecurentaxesdelete(get(gca,'Children'));camlightleft;newLineColor=[0.7,0.4,0.4];%re-drawthewholeheartwithsurffunctionh=surf(x,y,z,'EdgeColor',newLineColor,'FaceColor','r');fori=1:36camorbit(10,0,'data');pause(0.1);endmediumColor=[0.8,0.3,0.3];set(h,'EdgeColor',mediumColor,'FaceLighting','gouraud');fori=1:36camorbit(10,0,'data');drawnow;endset(h,'EdgeColor','none');drawnow;%definetheZ-pathforheartre-shapedeltaZ=linspace(-.3,0,5);deltaZ=[fliplr(deltaZ),deltaZ];deltaZ=repmat(deltaZ,1,3);zlim('manual');foriLoop=1:length(deltaZ)%drawthebouncingheartcurDeltaZ=deltaZ(mod(iLoop,length(deltaZ))+1);ratio=(radius+curDeltaZ)*radius;newZ=z*ratio+curDeltaZ;ifexist('h','var')&&ishandle(h),delete(h);end;h=surf(x,y,newZ,'EdgeColor','none','FaceColor','r','FaceLighting','gouraud');camorbit(-1,0,'data');drawnow;end%drawthebouncingheartdeltaZ=generateSimulatedBouncingPath;zlim([-1,3]);foriLoop=1:length(deltaZ)ifexist('h','var')&&ishandle(h),delete(h);end;curDeltaZ=deltaZ(iLoop);ifcurDeltaZ0%whendeltaZisunderthezerobar,it'shittingthegroundratio=(radius+curDeltaZ)*radius;newZ=z*ratio+curDeltaZ;else%whendeltaZisbiggerthanzero,it'safreeobjectnewZ=z+curDeltaZ;endh=surf(x,y,newZ,'EdgeColor','none','FaceColor','r','FaceLighting','gouraud');%zoominandoutformorevivideffectifiLoop=length(deltaZ)/2camzoom(0.99);elsecamzoom(1.005);endcamorbit(390/length(deltaZ),0,'data');drawnow;endtext(0.3,2.4,'X','FontSize',25,'FontName','TimesNewRoman')%%%%%%%text(0.3,1.8,'X','FontSize',30,'FontName','TimesNewRoman')%%%%%%%endfunctionz=generateSimulatedBouncingPath(iterationNum,roundTimes)%definethezpathofthecenterofballforrountTimesofbouncingalongz-axisifnargin2roundTimes=3;endifnargin1iterationNum=40;endradius=1;deltaH=1.5;%thehighestpointcomapredwithradiusdeltaL=-0.3;%thelowestpointcomapredwithradiusrangeY=[-2,1];%gravityacceleratorconstant%Note:thisissmallertothereal9.8value,butthismakesthebouncingslowerg=5;%definetheacceleratorforthez0%kissimilartospringcoefficient,butit'smorecomplicatedthanthat.%kjustneedstosatisfyoneconditionhere,i.e.itwillstopatlocationdefinedby%deltaL*radiusk=g*deltaH/(-deltaL);%startpointstartZ=deltaH*radius;%speedwhenz=0V0=sqrt(2*startZ*g);%calcualtethetotaltimeforrountTimesofbouncing%timeforz0plusZtime=sqrt(2*deltaH*radius/g);%timeforz0minusZtime=sqrt(-2*deltaL*radius/k);%halfperiodqPeriod=plusZtime+minusZtime;%definetotalTimeqPeriodTime=linspace(0,qPeriod,iterationNum);%zpathwhenz0zPlus=startZ-0.5*g*qPeriodTime(qPeriodTimeplusZtime).^2;%zpathwhenz0zMinus=-V0*(qPeriodTime(qPeriodTime=plusZtime)-plusZtime)+0.5*k*(qPeriodTime(qPeriodTime=plusZtime)-plusZtime).^2;z=[zPlus,zMinus,fliplr(zMinus),fliplr(zPlus)];%thelastzPlusisforthestopactionz=[repmat(z,1,roundTimes),zPlus];y=linspace(rangeY(1),rangeY(2),numel(z));endfunction[cordX,cordY,cordZ]=heart(sizeTheta,sizeFai)ifnargin==0sizeTheta=30;sizeFai=40;elseifnargin2error('heart:NotEnoughArugments','Pleasegiveatleasttwoinputarugments');end;theta=linspace(0,pi,sizeTheta)';nudge=0.0001;%usedtoavoidtheoperlappingfai=linspace(0+nudge,pi/2-nudge,round(sizeFai/4));a=9/4;b=9/80;A=1+(a-1)*(sin(theta).^2)*(sin(fai).^2);B=(sin(theta).^2.*cos(theta).^3)*(1+(b-1)*(sin(fai).^2));rou=zeros(size(A));foriLoop=1:numel(A);curA=A(iLoop);curB=B(iLoop);%thisisthepolarcoordinatesversionofthesexticequationfoundon%=[curA^3,-curB,-3*curA^2,0,3*curA,0,-1];solutions=roots(polyFactors);realRou=real(solutions(abs(imag(solutions))1e-9));rou(iLoop)=realRou(realRou0);end%x,y,zforthequaterofwholeheartx=rou.*(sin(theta)*cos(fai));y=rou.*(sin(theta)*sin(fai));z=rou.*(cos(theta)*ones(size(fai)));%x,y,zforthewholeheartcordX=[x,-fliplr(x),-x,fliplr(x)];cordY=[y,fliplr(y),-y,-fliplr(y)];cordZ=[z,fliplr(z),z,fliplr(z)];end