极坐标绘图

wusemm

请教要如何才能绘制出如下的效果？
我用polar效果只能会出部分效果，角度有点差异

[ 本帖最后由 ChaChing 于 2010-6-17 23:28 编辑 ]

ChaChing

印象中官网中好像有, LZ先搜下!

wusemm

请问什么是官网，网址？

ChaChing

参考下
http://www.mathworks.com/matlabc ... -function-halfpolar

ChaChing

好像直接把4F程式列出较为完整!

1. function h = halfPolar(phi,gain,linestyle,xtickval)
   
2. % HALFPOLAR function performs the polar plot in radian angle range [0 pi]
   
3. % using half polar coordinates
   
4. %
   
5. % HALFPOLAR(phi,gain) makes a plot with phi in radians angle range [0 pi]
   
6. % and gain in half polar coordinates. Phi and gain could be vector or matrix
   
7. % with the same dimensions. When Phi is a vector(1XN) and vector(1XN) and
   
8. % matrix(MXN)gain values are allowable. When Phi is a matrix(MXN) and then
   
9. % the gain should only be matrix(MXN) with the same row numbers.
   
10. %
    
11. % HALFPOLAR(phi,gain, linestyle) uses specified linestyle defined in cell
    
12. % variable linestyle to plot phi and gain, like: linestyle = {'ko-','b--'}
    
13. %
    
14. % HALFPOLAR(phi,gain, linestyle, xtickval) plots phi and gain with the
    
15. % prescribed tick value xtickval. like: xtickval = [15 25 45 75 105]
    
16. %
    
17. % IMPORTANT NOTE:
    
18. % HALFPOLAR is designed especially for the radiated sould pressure level
    
19. % demonstration. Accordingly the gain value is designed to be dB value.
    
20. % the default xtickvalue is also designed to match the dB value range.
    
21. % So you are permitted to modify the function in order to correctly
    
22. % demonstrate data for your specified applications.
    
23. %
    
24. % HALFPOLAR is coded straightforward without variables checking and fault
    
25. % testing.The author don't take the responsibility for the program bug.
    
26. %
    
27. % Example:
    
28. % phi = linspace(0,pi,20);
    
29. % gain = 80 * rand(3,length(phi)) + 40;
    
30. % h = halfPolar(phi,gain,{'k-','g-.','r--','k-'},[30 40 65 80 95 120])
    
31. %
    
32. % See also polar;
    
33. %
    
34. % Don't heisitate to let me know the problems about the function:
    
35. % lhd06@mails.thu.edu.cn
    
36. % Departmen of Electrical Engineering,
    
37. % Tsinghua University, Beijing, China.
    
38. %
    
39. % Reference function: polar.m provided by the matlab graphics toolbox
    
40. %
    
41. % If the program is modified and redistributed, the description above
    
42. % should be included.
    
43. min_gain = min(min(gain)); max_gain = max(max(gain));
    
44. space_val = 4; spoke_val = 6; spoketickval = linspace(0,pi,spoke_val+1);
    
45. if nargin < 2, error('Not enough input arguments.'); end
    
46.   
    
47. % if nargin < 3, rangedb = [min_gain max_gain]; end
    
48. if nargin < 3  linestyle = 'auto';xtickval = linspace(round(min_gain) - 10,round(max_gain) + 10,space_val);  end
    
49. if nargin < 4, xtickval = linspace(round(min_gain) - 10,round(max_gain) + 10,space_val); end
    
50. % normalization the input value
    
51. gain_range = xtickval(end) - xtickval(1);   %obtain the gain range
    
52. gain_scale = abs((gain - xtickval(1)))/gain_range; %[0,1]
    
53. cax = newplot; tc = get(cax,'xcolor'); ls = get(cax,'gridlinestyle');
    
54.    
    
55. th = linspace(0,pi,200); xunit = cos(th); yunit = sin(th);
    
56. % plot background if necessary
    
57. if ~isstr(get(cax,'color'))
    
58.        patch('xdata',xunit,'ydata',yunit, 'edgecolor',tc,'facecolor',get(gca,'color'));
    
59. end
    
60. hold on; [m1,n1] =size(phi); [m2,n2] =size(gain);
    
61. if(m1==1)
    
62.     for i=1:m2
    
63.         x = gain_scale(i,:) .* cos(phi); y = gain_scale(i,:) .* sin(phi);
    
64.         if(strcmp(linestyle,'auto')), h(i) = plot(x,y,'linewidth',2);
    
65.         else h(i) = plot(x,y,linestyle{i},'linewidth',2);
    
66.         end
    
67.     end
    
68. else
    
69.     if(m1 ~= m2), error('Matrix should be have the same dimension!');
    
70.     else      
    
71.         for i=1:m2
    
72.          x = gain_scale(i,:) .* cos(phi(i,:)); y = gain_scale(i,:) .* sin(phi(i,:));
    
73.             if(strcmp(linestyle,'auto')), h(i) = plot(x,y,'linewidth',2);  %,'LineWidth', 2
    
74.             else h(i) = plot(x,y,linestyle{i},'linewidth',2);
    
75.             end
    
76.         end
    
77.     end
    
78. end
    
79. set(gca,'dataaspectratio',[1 1 1]);axis off;
    
80. % legend_val = legend(h(1:end),'Controlled','Uncontrolled');
    
81. % set(legend_val,'location','SouthOutside');
    
82. hold on;
    
83. %define the circle
    
84. contour_val = zeros(1,length(xtickval));
    
85. contour_val = abs((xtickval - xtickval(1))/gain_range);
    
86. for k=2:length(contour_val) %gain circles
    
87.     plot(xunit*contour_val(k), yunit*contour_val(k),ls,'color','black');
    
88.     text(contour_val(k),-0.05, sprintf('%.3g',xtickval(k)),'horiz', 'center', 'vert', 'middle');
    
89.     text(- contour_val(k),-0.05, sprintf('%.3g',xtickval(k)),'horiz', 'center', 'vert', 'middle');
    
90. end
    
91. text(contour_val(1),-0.05,sprintf('%.3g',xtickval(1)),'horiz', 'center', 'vert', 'middle');
    
92. text(-0.45,-0.15,'Sound Pressure Level(SPL) /dB');
    
93. %plot the spokes and  % annotate spokes in degrees
    
94. cst = cos(spoketickval); snt = sin(spoketickval);
    
95. for k = 1:length(spoketickval) - 1
    
96.     plot(cst(k) * contour_val,snt(k) * contour_val,ls,'color',tc,'linewidth',1, 'handlevisibility','off');
    
97.     text(1.07*cst(k),1.07*snt(k), sprintf('%.3g^{o}',spoketickval(k)/pi*180, 'horiz', 'center', 'vert', 'middle');
    
98. end
    
99. text(1.08*cst(end),1.07*snt(end), sprintf('%.3g^{o}',spoketickval(end)/pi*180),'horiz', 'center', 'vert', 'middle');
    
100. hold off;
     

复制代码