[zz] I <3 Math: 用MATLAB画出3D心型

不素偶写的code :P

原帖地址：http://stackoverflow.com/questio ... aped-mesh-in-matlab



MATLAB code：

1. % Initialize the volume data, figure, and axes:
   
2. [X,Y,Z] = meshgrid(linspace(-3,3,101));
   
3. F = -X.^2.*Z.^3-(9/80).*Y.^2.*Z.^3+(X.^2+(9/4).*Y.^2+Z.^2-1).^3;
   
4. hFigure = figure;
   
5. sz = get(hFigure, 'Position');
   
6. set(hFigure, 'Position', [sz(1)-0.15*sz(3) sz(2) 1.3*sz(3) sz(4)]);
   
7. set(hFigure,'color','w', 'menu','none')
   
8. 
   
9. hAxes = axes('Parent',hFigure,'NextPlot','add',...
   
10.     'DataAspectRatio',[1 1 1],...
    
11.     'XLim',[30 120],'YLim',[35 65],'ZLim',[30 75]);
    
12. view([-39 30]);
    
13. axis off
    
14. 
    
15. % hidden surface removal
    
16. hidden on
    
17. 
    
18. % Fill the inside of the mesh with an isosurface to
    
19. % block rendering of the back side of the heart
    
20. p = patch(isosurface(F,-0.001));
    
21. set(p,'FaceColor','w','EdgeColor','none');
    
22. 
    
23. % Create and plot contours in the y-z plane:
    
24. for iX = [35 38 41 45 48 51 54 57 61 64 67]
    
25.     plane = reshape(F(:,iX,:),101,101);
    
26.     cData = contourc(plane,[0 0]);
    
27.     xData = iX.*ones(1,cData(2,1));
    
28.     plot3(hAxes,xData,cData(2,2:end),cData(1,2:end),'k');
    
29.     pause(.1), drawnow
    
30. end
    
31. 
    
32. % Create and plot contours in the x-z plane:
    
33. for iY = [41 44 47 51 55 58 61]
    
34.     plane = reshape(F(iY,:,:),101,101);
    
35.     cData = contourc(plane,[0 0]);
    
36.     yData = iY.*ones(1,cData(2,1));
    
37.     plot3(hAxes,cData(2,2:end),yData,cData(1,2:end),'k');
    
38.     pause(.1), drawnow
    
39. end
    
40. 
    
41. % Create and plot contours in the x-y plane:
    
42. for iZ = [36 38 40 42 44 46 48 50 52 54 56 58 60 62 64 66 69 71]
    
43.     plane = F(:,:,iZ);
    
44.     cData = contourc(plane,[0 0]);
    
45.     startIndex = 1;
    
46.     if size(cData,2) > (cData(2,1)+1)
    
47.         startIndex = cData(2,1)+2;
    
48.         zData = iZ.*ones(1,cData(2,1));
    
49.         plot3(hAxes,cData(1,2:(startIndex-1)),...
    
50.                 cData(2,2:(startIndex-1)),zData,'k');
    
51.     end
    
52.     zData = iZ.*ones(1,cData(2,startIndex));
    
53.     plot3(hAxes,cData(1,(startIndex+1):end),...
    
54.         cData(2,(startIndex+1):end),zData,'k');
    
55.     pause(.1), drawnow
    
56. end
    
57. pause(.2)
    
58. 
    
59. text(7,50,70,'I', 'fontWeight','bold','FontAngle','italic','fontsize',100)
    
60. pause(.5)
    
61. text(80,50,43,'Math', 'fontWeight','bold','FontAngle','italic','fontsize',100)
    
62. pause(.2)
    
63. 
    
64. line([20 80],[50 50],[52.5 52.5], 'color','k')
    
65. line([50 50],[20 80],[52.5 52.5], 'color','k')
    
66. line([50 50],[50 50],[30 80], 'color','k')
    
67. text(20,50,50,'x')
    
68. text(48,20,50,'y')
    
69. text(45,50,80,'z')
    
70. 
    
71. text(30,60,30,'(x^2+9/4y^2+z^2-1)^3 - x^2z^3-9/80y^2z^3=0', 'fontsize',8)
    
72. text(35,45,30,'-3<x,y,z<3', 'fontsize',8)

复制代码

####

for Statistics geeks:

端木mm发的水贴都 总是这么有含金量

以前看过mathematica版本的，只要一句话……比MATLAB方便好多啊……
3D:
http://goo.gl/jDxGi
2D:
http://goo.gl/UXePi

回复 3# zach

赞～ 2d版本也不错～

我想知道这些公式是如何推出来的那?

我想知道这些公式是如何推出来的那?
小K 发表于 2011-3-24 23:57

貌似是神马笛卡尔曲线

哈! 我有这个心的t-shirt! 穿到实验室被大家感慨, chinese geek啊.......哈哈, 有人很羡慕也想买, 我说是在中国的网站买的, 于是对方很郁闷, 哈哈哈
是这个, 是matrix67设计的, 这里, 我前些日子才知道这个数学大牛居然是北大中文系的!

哈! 我有这个心的t-shirt! 穿到实验室被大家感慨, chinese geek啊.......哈哈, 有人很羡慕也想买, 我说是在 ...
RehsuR 发表于 2011-3-25 00:18

太简单了有没有%……
自己拿MATLAB画出来，然后去：
http://www.customink.com/
把图片upload上去，还可以写你想写的whatever内容

太简单了有没有%……
自己拿MATLAB画出来，然后去：

把图片upload上去，还可以写你想写的whatever内容 ...
zach 发表于 2011-3-24 23:33

那哥们儿天天用mathematica做modeling, 肯定是可以做, 但是还没geek到那个份儿上吧...懒得追细节了, 只是觉得特cool, 呵呵~

哈! 我有这个心的t-shirt! 穿到实验室被大家感慨, chinese geek啊.......哈哈, 有人很羡慕也想买, 我说是在 ...
RehsuR 发表于 2011-3-25 13:18

    国内买了带过去？
   读Phd前就这么Geek啊....
     仰慕....