Mathematica中文论坛-非官方

标题: 豆粑粑　matlab tensor scatter3 multi sup figure [打印本页]

---

作者: meatball1982    时间: 2017-8-24 20:35
标题: 豆粑粑　matlab tensor scatter3 multi sup figure
Got a tensor, H(ff,m,k,n), the value is x+yi.
The plot should be (x,y,ff), with maker size (k)  and color (n)
There are 6 values of m.
So 6 subplot

1. clear all
   
2. clc
   
3. 
   
4. % load Data.mat
   
5. %
   
6. % [num_ff,num_mm,num_k,num_n]=size(H);
   
7. %
   
8. % % % ff_line = [1:num_ff];
   
9. % % % mm_line = [1:num_mm];
   
10. % % % k_line  = [1:num_k ];
    
11. % % % n_line  = [1:num_n ];
    
12. % % %
    
13. % % % x{1} = ff_line;
    
14. % % % x{2} = mm_line;
    
15. % % % x{3} = k_line;
    
16. % % % x{4} = n_line;
    
17. % % % n=4;
    
18. % % %
    
19. % % % [c{1:n}] = ndgrid(x{:});
    
20. % % % c = cell2mat(cellfun(@(a)a(:),c,'un',0));
    
21. % % %
    
22. % % % [H_line]=H(:);
    
23. % % X=real(H_line);
    
24. % % Y=imag(H_line);
    
25. %
    
26. %
    
27. %
    
28. % colormap(jet)
    
29. % all=[];
    
30. % for i=1:num_ff
    
31. %     for j=1:num_mm
    
32. %         for k= 1:num_k
    
33. %             for l=1:num_n
    
34. %                 
    
35. %                 all=[all;real(H(i,j,k,l)),imag(H(i,j,k,l)),i,j,k,l];
    
36. %             end
    
37. %         end
    
38. %     end
    
39. % end
    
40. %
    
41. % save Mat_6d.mat
    
42. 
    
43. load Mat_6d.mat
    
44. 
    
45. 
    
46. % ind=1:7:120000;
    
47. 
    
48. col_dat={'hot';'cool';'summer';'winter';'gray';'pink'};
    
49. view_mat=[-130 40
    
50.           -130 40
    
51.           -130 40
    
52.           -130 40
    
53.           -130 40
    
54.           -130 40];
    
55. sbp_width=0.8;
    
56. sbp_heig=0.8;
    
57. n_row = 1;
    
58. n_col = 1;
    
59. 
    
60. [out_pos]=fun_mm_subplot_pos(n_row,n_col,sbp_width,sbp_heig);
    
61. 
    
62. 
    
63. 
    
64. 
    
65. cho_int=1;
    
66. 
    
67. for i=1:6
    
68.     ind= all(:,4)==i;
    
69. %     tm = all;
    
70.     clf
    
71.     h=figure(1);
    
72.     set(h, 'Position', [100, 100, 600, 400]);
    
73.     ax=axes('position',out_pos(1,:));
    
74.     scatter3(all(ind(1:cho_int:end),1),all(ind(1:cho_int:end),2),all(ind(1:cho_int:end),3),...
    
75.         8+0.3*all(ind(1:cho_int:end),5),sqrt(all(ind(1:cho_int:end),6)),'o','filled')
    
76.     colormap(gca,flipud(col_dat{i}))
    
77.     axis tight
    
78.     xlabel('real')
    
79.     ylabel('imag')
    
80.     zlabel('freq')
    
81.    
    
82.     view(-130, 40)
    
83.     colorbar
    
84.     title(['mm',mat2str(i)])
    
85.     h=gcf;
    
86.     fi_na=['./fig_dif_mm',mat2str(i)];
    
87.      fun_work_li_035_myfig_out(h,fi_na,3);
    
88. end
    

复制代码

Data.mat.zip

2017-8-24 20:34 上传

点击文件名下载附件

1.75 MB, 下载次数: 0

data

---

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

Powered by Discuz! X3.2