Mathematica中文论坛-非官方

标题: 豆粑粑 matlab scatter high dimensions tensor visualization, multi colorbar [打印本页]

---

作者: meatball1982    时间: 2017-8-29 15:31
标题: 豆粑粑 matlab scatter high dimensions tensor visualization, multi colorbar


matlab scatter high dimensions tensor visualization

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. %% outline
   
6. % plot tensor with value in one figure.
   
7. 
   
8. %% main
   
9. 
   
10. % point number --------------------------------
    
11. n=400;         
    
12. % generate x,y,z  -----------------------------
    
13. z = linspace(0,4*pi,n)+rand(1,n)*0.2;
    
14. x = 2*cos(z) + [1:n]/50.*rand(1,n);
    
15. y = 2*sin(z) + [1:n]/80.*rand(1,n);
    
16. 
    
17. % point properties v1,v2,...,v8  --------------
    
18. v1 = (10*abs(z)+10).*rand(size(z));
    
19. v2 = z + 4  ;
    
20. v3 = sin(z/2*0.2*pi)+2;
    
21. v4 = log(z+2);
    
22. v5 = cos((pi*0.1*y).^0.5) +2  ;
    
23. v6 = y+2;
    
24. v7 = sin(pi*log(z+1))+2;
    
25. v8 = cos(pi*log(z+1))+2;
    
26. 
    
27. % plot parameter -------------------------------
    
28. xmin= -10;xmax = 15;
    
29. ymin= -10;ymax = 15;
    
30. zmin=  0;zmax = 15;
    
31. 
    
32. % view point ----------------------------------
    
33. vx = -50; vy = 20;
    
34. 
    
35. % fontsize ------------------------------------
    
36. mm_fz = 12;
    
37. 
    
38. % set fig size --------------------------------
    
39. fig_pos = [0.09 0.15 0.53 0.75];
    
40. axis_pos= [xmin xmax ymin ymax zmin zmax];
    
41. 
    
42. % plot 3D points ------------------------------
    
43. h=figure(1);
    
44. set(h, 'Position', [100, 100, 800, 500]);
    
45. 
    
46. % plot zeros plane ------------------------------
    
47. ax1 = axes;
    
48. set(ax1,'position',fig_pos );
    
49. h_1 = scatter3(ax1,x,y,zeros(size(x)),3*v3,v4,'o','filled');
    
50. box on
    
51. hidden off
    
52. axis(axis_pos );
    
53. view([vx,vy])
    
54. colormap(ax1,summer);
    
55. xlabel('IamX')
    
56. ylabel('IamY')
    
57. zlabel('IamZ')
    
58. title('IamTitle')
    
59. 
    
60. 
    
61. % plot x plane ------------------------------
    
62. ax2 = axes;
    
63. set(ax2,'position',fig_pos);
    
64. h_2 = scatter3(ax2,xmax*ones(size(x)),y,z,7*v5,v6,'s','filled');
    
65. box off
    
66. axis off
    
67. hidden off
    
68. axis(axis_pos );
    
69. view([vx,vy])
    
70. colormap(ax2,'pink');
    
71. 
    
72. 
    
73. % plot y plane ------------------------------
    
74. ax3 = axes;
    
75. set(ax3,'position',fig_pos);
    
76. h_3 = scatter3(ax3,x,ymax*ones(size(y)),z,8*v7,v8,'^','filled');
    
77. box off
    
78. axis off
    
79. hidden off
    
80. axis(axis_pos );
    
81. view([vx,vy])
    
82. colormap(ax3,'cool');
    
83. 
    
84. % plot xyz  ------------------------------
    
85. ax4 = axes;
    
86. set(ax4,'position',fig_pos);
    
87. h_4 = scatter3(ax4,x,y,z,1.1*v1,v2,'o','filled');
    
88. % box off
    
89. axis off
    
90. hidden off
    
91. axis(axis_pos );
    
92. view([vx,vy])
    
93. colormap(ax4,'jet');
    
94. 
    
95. % set colorbar --------------------------------
    
96. cb1 = colorbar(ax1,'Position',[.65  .15 .05 .6]);
    
97. cb2 = colorbar(ax2,'Position',[.74  .15 .05 .6]);
    
98. cb3 = colorbar(ax3,'Position',[.83   .15 .05 .6]);
    
99. cb3 = colorbar(ax4,'Position',[.92   .15 .05 .6]);
    
100. 
     
101. % set fontsize --------------------------------
     
102. set(ax1,'fontsize',mm_fz)
     
103. set(ax2,'fontsize',mm_fz)
     
104. set(ax3,'fontsize',mm_fz)
     
105. set(ax4,'fontsize',mm_fz)
     
106. 
     
107. 
     
108. h=gcf;
     
109. fi_na=['./fig_8_val_3D'];
     
110. % fun_work_li_035_myfig_out(h,fi_na,3);
     
111. 
     
112. %% plot point properties ---------------------------
     
113. % clf
     
114. % sbp_width=0.8;
     
115. % sbp_heig=0.75;
     
116. % n_row = 4;
     
117. % n_col = 3;
     
118. % [out_pos]=fun_mm_subplot_pos(n_row,n_col,sbp_width,sbp_heig);
     
119. % h=figure(2)
     
120. % set(h, 'Position', [100, 100, 800, 600]);
     
121. %
     
122. % ax=axes('position',out_pos(1,:));
     
123. % plot(x,'.')
     
124. % title('x');
     
125. % ax=axes('position',out_pos(2,:));
     
126. % plot(y,'.')
     
127. % title('y');
     
128. % ax=axes('position',out_pos(3,:));
     
129. % plot(z,'.')
     
130. % title('z');
     
131. %
     
132. % ax=axes('position',out_pos(4,:));
     
133. % plot(v1,'.')
     
134. % title('v1');
     
135. % ax=axes('position',out_pos(5,:));
     
136. % plot(v2,'.')
     
137. % title('v2');
     
138. % ax=axes('position',out_pos(6,:));
     
139. % plot(v3,'.')
     
140. % title('v3');
     
141. %
     
142. % ax=axes('position',out_pos(7,:));
     
143. % plot(v4,'.')
     
144. % title('v4');
     
145. % ax=axes('position',out_pos(8,:));
     
146. % plot(v5,'.')
     
147. % title('v5');
     
148. % ax=axes('position',out_pos(9,:));
     
149. % plot(v6,'.')
     
150. % title('v6');
     
151. %
     
152. % ax=axes('position',out_pos(10,:));
     
153. % plot(v7,'.')
     
154. % title('v7');
     
155. % ax=axes('position',out_pos(11,:));
     
156. % plot(v8,'.')
     
157. % title('v8');
     
158. %
     
159. %
     
160. % h=gcf;
     
161. % fi_na=['./fig_8_val_split'];
     
162. % fun_work_li_035_myfig_out(h,fi_na,3);
     
163. %

复制代码

---

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

Powered by Discuz! X3.2