Mathematica中文论坛-非官方

标题: 豆粑粑　matlab radarplot [打印本页]

---

作者: meatball1982    时间: 2018-2-19 15:29
标题: 豆粑粑　matlab radarplot
看着效果还不错。
其它的语言似乎都有直接用的函数，matlab，落后了。
修改一下别人的函数。挺好用。

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. 
   
6. %% outline
   
7. % test radarplot
   
8. 
   
9. %% main
   
10. % P=rand(13,2);
    
11. % P(:,1)=P(:,1)*180/pi;
    
12. % P(:,2)=20*P(:,2);
    
13. % radarPlot(P)
    
14. 
    
15. nPoints = 4;
    
16. nDimensions = 6;
    
17. P = rand(nDimensions, nPoints);
    
18. P
    
19. pointNames = arrayfun( @(i)sprintf('p_{%d}', i),...
    
20. 1:nPoints, 'UniformOutput', false);
    
21. 
    
22. radarPlot(P, 'o-','LineWidth', 4, 'MarkerFaceColor', [0,0,0])
    
23. % radarPlot(P, 'o-','LineWidth', 4)
    
24. 
    
25. legend(pointNames{:}, 'Location', 'Best');
    
26. title('Radar Plot Demo');

复制代码

1. % RADARPLOT spiderweb or radar plot
   
2. % radarPlot(P) Make a spiderweb or radar plot using the columns of P as datapoints.
   
3. %  P is the dataset. The plot will contain M dimensions(or spiderweb stems)
   
4. %  and N datapoints (which is also the number of columns in P). Returns the
   
5. %  axes handle
   
6. %
   
7. % radarPlot(P, ..., lineProperties) specifies additional line properties to be
   
8. % applied to the datapoint lines in the radar plot
   
9. %
   
10. % h = radarPlot(...) returns the handles to the line objects.
    
11. function varargout = radarPlot( P, varargin )
    
12. %
    
13. % demo
    
14. %%%%%%%%
    
15. % nPoints = 4;
    
16. % nDimensions = 6;
    
17. % P = rand(nDimensions, nPoinhts);
    
18. % P
    
19. % pointNames = arrayfun( @(i)sprintf('p_{%d}', i),...
    
20. % 1:nPoints, 'UniformOutput', false);
    
21. %
    
22. % radarPlot(P, 'o-','LineWidth', 3, 'MarkerFaceColor', [0,0,0])
    
23. % legend(pointNames{:}, 'Location', 'Best');
    
24. % title('Radar Plot Demo');
    
25. %%%%%%%
    
26. 
    
27. %%% Get the number of dimensions and points
    
28. [M, N] = size(P);
    
29. 
    
30. %%% Plot the axes
    
31. % Radial offset per axis
    
32. th = (2*pi/M)*(ones(2,1)*(M:-1:1));
    
33. % Axis start and end
    
34. r = [0;1]*ones(1,M);
    
35. % Conversion to cartesian coordinates to plot using regular plot.
    
36. [x,y] = pol2cart(th, r);
    
37. hLine = line(x, y,...
    
38.     'LineWidth', 1.5,...
    
39.     'Color', [1, 1, 1]*0.7  );
    
40. 
    
41. for i = 1:numel(hLine)
    
42.     set(get(get(hLine(i),'Annotation'),'LegendInformation'),...
    
43.         'IconDisplayStyle','off'); % Exclude line from legend
    
44. end
    
45. 
    
46. toggle = ~ishold;
    
47. 
    
48. if toggle
    
49.     hold on
    
50. end
    
51. 
    
52. %%% Plot axes isocurves
    
53. n_cir=4;
    
54. % Radial offset per axis
    
55. th = (2*pi/M)*(ones(n_cir,1)*(M:-1:1));
    
56. % Axis start and end
    
57. r = (linspace(0.1, 0.9, n_cir)')*ones(1,M);
    
58. % Conversion to cartesian coordinates to plot using regular plot.
    
59. [x,y] = pol2cart(th, r);
    
60. hLine = line([x, x(:,1)]', [y, y(:,1)]',...
    
61.     'LineWidth', 1,...
    
62.     'Color', [1, 1, 1]*0.5  );
    
63. for i = 1:numel(hLine)
    
64.     set(get(get(hLine(i),'Annotation'),'LegendInformation'),...
    
65.         'IconDisplayStyle','off'); % Exclude line from legend
    
66. end
    
67. 
    
68. 
    
69. %%% Insert axis labels
    
70. 
    
71. % Compute minimum and maximum per axis
    
72. minV = min(P,[],2);
    
73. maxV = max(P,[],2);
    
74. for j = 1:M
    
75.     % Generate the axis label
    
76.     msg = sprintf('x_{%d} = %5.2f ... %5.2f',...
    
77.         j, minV(j), maxV(j));
    
78.     [mx, my] = pol2cart( th(1, j), 1.1);
    
79.     text(mx, my, msg);
    
80. end
    
81. axis([-1,1,-1,1]*1.5)
    
82. 
    
83. % Hold on to plot data points
    
84. hold on
    
85. 
    
86. % Radius
    
87. R = 0.8*((P - (minV*ones(1,N)))./((maxV-minV)*ones(1,N))) + 0.1;
    
88. R = [R; R(1,:)];
    
89. Th = (2*pi/M) * ((M:-1:0)'*ones(1,N));
    
90. 
    
91. % polar(Th, R)
    
92. [X, Y] = pol2cart(Th, R);
    
93. 
    
94. h = plot(X, Y, varargin{:});
    
95. axis([-1,1,-1,1])
    
96. axis square
    
97. axis off
    
98. 
    
99. if toggle
    
100.     hold off
     
101. end
     
102. 
     
103. if nargout > 0
     
104.     varargout{1} = h;
     
105. end

复制代码

---

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

Powered by Discuz! X3.2