Mathematica中文论坛-非官方

标题: 豆粑粑，极坐标 非凸 3D曲面. 从不同角度看，显示全曲面 [打印本页]

---

作者: meatball1982    时间: 2017-3-28 13:49
标题: 豆粑粑，极坐标 非凸 3D曲面. 从不同角度看，显示全曲面


别人的问题
http://www.ilovematlab.cn/forum. ... 1&extra=#pid3168135

最近觉得一个3D图从不同的角度看，效果相当的好。
帖子中只是要画出环的部分，将原曲面通过透明的方式显示。似乎效果更加的好。

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. r=linspace(1,1.5,50);
   
6. the=linspace(0,2*pi,200);
   
7. the=[the,0];
   
8. 
   
9. [R,The]=meshgrid(r,the);
   
10. x=R.*cos(The);
    
11. y=R.*sin(The);
    
12. 
    
13. X=x+2;
    
14. Y=y+2;
    
15. Z=X.*Y;
    
16. 
    
17. x_g=linspace(0.5,3.5,100);
    
18. y_g = linspace(0.5,3.5,100);
    
19. [X_g,Y_g]=meshgrid(x_g,y_g);
    
20. Z_g =X_g.*Y_g;
    
21. 
    
22. sbp_width=0.75;
    
23. sbp_heig=0.75;
    
24. n_row = 2;
    
25. n_col = 2;
    
26. [out_pos]=fun_mm_subplot_pos(n_row,n_col,sbp_width,sbp_heig);
    
27. h=figure(1)
    
28. set(h, 'Position', [100, 100, 400, 400]);
    
29. ax=axes('position',out_pos(1,:));
    
30. hold on
    
31. 
    
32. h=surf(X_g,Y_g,Z_g,'edgecolor','none')
    
33. alpha(h,0.2)
    
34. surf(X,Y,Z,Z,'edgecolor','none')
    
35. 
    
36. colormap(jet)
    
37. % colormap(h2,cool)
    
38. view(-20,20)
    
39. axis tight
    
40. grid on
    
41. 
    
42. ax=axes('position',out_pos(2,:));
    
43. hold on
    
44. h=surf(X_g,Y_g,Z_g,'edgecolor','none')
    
45. alpha(h,0.2)
    
46. surf(X,Y,Z,Z,'edgecolor','none')
    
47. colormap(jet)
    
48. axis equal tight
    
49. view(0,90)
    
50. 
    
51. ax=axes('position',out_pos(3,:));
    
52. hold on
    
53. h=surf(X_g,Y_g,Z_g,'edgecolor','none')
    
54. alpha(h,0.2)
    
55. hsurf(X,Y,Z,Z,'edgecolor','none')
    
56. colormap(jet)
    
57. axis tight
    
58. view(0,0)
    
59. grid on
    
60. 
    
61. ax=axes('position',out_pos(4,:));
    
62. hold on
    
63. h=surf(X_g,Y_g,Z_g,'edgecolor','none')
    
64. alpha(h,0.2)
    
65. surf(X,Y,Z,Z,'edgecolor','none')
    
66. colormap(jet)
    
67. axis tight
    
68. view(-90,0)
    
69. grid on
    
70. 
    
71. 
    
72. 
    

复制代码

其中的

1. fun_mm_subplot_pos

复制代码

可以用subplot函数来代替，不过我觉得我的更加好些。

1. function [out_pos]=fun_mm_subplot_pos(n_row,n_col,sbp_width,sbp_heig);
   
2. 
   
3. % subplot pos
   
4. % input : -----------------
   
5. % n_row     : number of row
   
6. % n_col     : number of column
   
7. % sbp_width : subplot width  0-1
   
8. % sbp_heig  : subplot height 0-1
   
9. % output : ----------------
   
10. % out_pos
    
11. %
    
12. %
    
13. % sbp_width=0.8;
    
14. % sbp_heig=0.85;
    
15. % n_row = 4;
    
16. % n_col = 6;
    
17. % [out_pos]=fun_mm_subplot_pos(n_row,n_col,sbp_width,sbp_heig);
    
18. % h=figure(1)
    
19. % set(h, 'Position', [100, 100, 800, 600]);
    
20. % [x,y,z]=peaks(30);
    
21. % caxis_mm=[-5 30];
    
22. % for i=1:n_row*n_col
    
23. %     ax=axes('position',out_pos(i,:));
    
24. %     surf(x,y,z+i,'edgecolor','none')
    
25. %     caxis(caxis_mm)
    
26. %     view(0,90)
    
27. %     hold on
    
28. %     text([2],[2],[30],mat2str(i),'fontsize',20)
    
29. % end
    
30. 
    
31. for i=1:n_row*n_col
    
32. 
    
33.     k1=ceil(i/n_col);
    
34.     k2=i-(k1-1)*n_col;
    
35.    
    
36.     out_pos(i,:)=[(1-sbp_width)/n_col/2 + (k2-1)/n_col,...
    
37.                   (1-sbp_heig )/n_row/2 + (n_row-k1)/n_row, ...
    
38.                   sbp_width/n_col,...
    
39.                   sbp_heig/n_row];
    
40. end
    

复制代码

---

欢迎光临 Mathematica中文论坛-非官方 (http://ilovemathematica.com/)

Powered by Discuz! X3.2