(精品word)matlab动态三维心形.doc

(精品word)matlab动态三维心形
(精品word)matlab动态三维心形
(精品word)matlab动态三维心形

％构造体积方程和坐标轴，画出图形;
［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

％制作出动态的隐形效果;

hiddenon

%画出网格，制作网格动态效果;

%快渲染心得背面:

p=patch（isosurface(F，—0.001)）;
set(p,’FaceColor’,'w’,’EdgeColor',’w'）；

％构造Y—Z平面,,描完函数在该平面的点:
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）,drawnow
end

%构造X-Z平面,描完函数在该平面的点：
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)，drawnow
end

％构造X-Y平面,描完函数在该平面的点:
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’);
end
zData=iZ。＊ones(1,cData（2，startIndex));
plot3（hAxes，cData(1,(startIndex+1）:end），...
cData(2,(startIndex+1）:end)，zData,'r’）；
pause(.1),drawnow
end

%给三维心着色
set(p，’FaceColor'，’r’，'EdgeColor’,'w’）；
pause(。2）；
set(p,’FaceColor’,'w’，’EdgeColor’,’r’）；

%函数已经画完,接下来为文字部分；
pause(.2)

%设置字体大小,粗细,位置等,以下是打出I（心型图）Wendy;

text（7，50,70,'I’，'fontWeight’,'bold'，'FontAngle',’italic’,'FontName','TrebuchetMS'，’fontsize’，60，'Color’，’r')
pause（.5)
text（80,50，43,'YOU','fontWeight',’bold’，'FontAngle’，'italic',’FontName'，'TrebuchetMS’，’fontsize',60,’Color’，’r')
pause(.2）

line([2080］,［5050］,［52。552。5]，’color’，'r'）
line(［5050],[2080］,［52。552。5],'color','r')
line([5050],［5050