Mathematica中文论坛-非官方

标题: 豆粑粑 matlab 断yaxis ，修改 breakYAxes. [打印本页]

---

作者: meatball1982    时间: 2017-5-31 10:46
标题: 豆粑粑 matlab 断yaxis ，修改 breakYAxes.
别人的问题：
http://www.ilovematlab.cn/forum. ... p;extra=#pid3224636

1. clear all
   
2. clc
   
3. 
   
4. %% outline
   
5. % break y axis
   
6. 
   
7. %% main
   
8. x=linspace(0,10*pi,100);
   
9. % y1 = sqrt(3)/3*x +8;
   
10. y=exp(x/7)+sin(x)-30;
    
11. 
    
12. 
    
13. y_break_start = 10;
    
14. y_break_end = 25;
    
15. break_type_w = 'Wave';
    
16. % h=BreakPlot(x,y,y_break_start,y_break_end,break_type);
    
17. break_type_R = 'RPatch';
    
18. 
    
19. % figure;
    
20. subplot(1,2,1)
    
21. plot(x,y,'.')
    
22. % subplot(2,2,2)
    
23. % % BreakPlot(rand(1,21),[1:10,40:50],10,40,'Line');
    
24. % h=BreakPlot(x,y,y_break_start,y_break_end,break_type_R);
    
25. subplot(1,2,2)
    
26. h=BreakPlot(x,y,y_break_start,y_break_end,break_type_w);
    
27. % hold on
    
28. % plot(x,y1,'r-')
    
29. 
    
30. 
    
31. h=gcf;
    
32. fi_na=['./fig_split_y_axis'];
    
33. fun_work_li_035_myfig_out(h,fi_na,3);
    
34. 
    
35. %% logs
    
36. % mod : 31-May-2017 09:07:50
    
37. %

复制代码

1. function h=BreakPlot(x,y,y_break_start,y_break_end,break_type)
   
2. % BreakPlot(x,y,y_break_start,y_break_end,break_type)
   
3. % Produces a plot who's y-axis skips to avoid unecessary blank space
   
4. %
   
5. % INPUT
   
6. % x
   
7. % y
   
8. % y_break_start
   
9. % y_break_end
   
10. % break_type
    
11. %    if break_type='RPatch' the plot will look torn
    
12. %       in the broken space
    
13. %    if break_type='Patch' the plot will have a more
    
14. %       regular, zig-zag tear
    
15. %    if break_plot='Line' the plot will merely have
    
16. %       some hash marks on the y-axis to denote the
    
17. %       break
    
18. %
    
19. % USAGE:
    
20. % figure;
    
21. % BreakPlot(rand(1,21),[1:10,40:50],10,40,'Line');
    
22. % figure;
    
23. % BreakPlot(rand(1,21),[1:10,40:50],10,40,'Patch');
    
24. % figure;
    
25. % BreakPlot(rand(1,21),[1:10,40:50],10,40,'RPatch');
    
26. % figure;
    
27. % x=rand(1,21);y=[1:10,40:50];
    
28. % subplot(2,1,1);plot(x(y>=40),y(y>=40),'.');
    
29. % set(gca,'XTickLabel',[]);
    
30. % subplot(2,1,2);plot(x(y<=20),y(y<=20),'.');
    
31. %
    
32. % IT'S NOT FANCY, BUT IT WORKS.
    
33. 
    
34. % Michael Robbins
    
35. % robbins@bloomberg.net
    
36. % michael.robbins@bloomberg.net
    
37. 
    
38. % data
    
39. if nargin<5 break_type='RPatch'; end;
    
40. if nargin<4 y_break_end=40; end;
    
41. if nargin<3 y_break_start=10; end;
    
42. if nargin<2 y=[1:10,40:50]; end;
    
43. if nargin<1 x=rand(1,21); end;
    
44. 
    
45. 
    
46. y_break_mid=(y_break_end-y_break_start)./2+y_break_start;
    
47. [x_mid_val,ind] = min(abs(y-y_break_mid));
    
48. x_break_mid = x(ind);
    
49. b_shift= y_break_start-sqrt(3)/3*x_break_mid;
    
50. b_shift
    
51. x_break_mid
    
52. y_break_mid
    
53. % erase useless data
    
54. x(y>y_break_start & y <y_break_end)=[];
    
55. y(y>y_break_start & y <y_break_end)=[];
    
56. 
    
57. % leave room for the y_break_end
    
58. % % [junk,i]=min(y>=y_break_end);
    
59. % % if i>y_break_end
    
60. % %     x=[x(1:i-1) NaN x(i:end)];
    
61. % %     y=[y(1:i-1) y_break_mid y(i:end)];
    
62. % % end;
    
63. 
    
64. % remap
    
65. y2=y;
    
66. y2(y2>=y_break_end)=y2(y2>=y_break_end)-(y_break_end-y_break_start);
    
67. 
    
68. % plot
    
69. h=plot(x,y2,'.');
    
70. % h=plot(x,y,'.');
    
71. 
    
72. % make break
    
73. xlim=get(gca,'xlim');
    
74. ytick=get(gca,'YTick');
    
75. [junk,i]=min(ytick<=y_break_start);
    
76. y=(ytick(i)-ytick(i-1))./2+ytick(i-1);
    
77. dy=(ytick(2)-ytick(1))./7;
    
78. xtick=get(gca,'XTick');
    
79. x=xtick(1);
    
80. dx=(xtick(2)-xtick(1))./2;
    
81. n_int=100;
    
82. switch break_type
    
83.     case 'Wave',
    
84.          dx=(xlim(2)-xlim(1))./10;
    
85. %         xx = xlim(1)+dx.*[0:101];
    
86.         xx = linspace(xlim(1),xlim(end),n_int);
    
87. %         yy = sin(xx)+sqrt(3)/3*xx+ytick(i-1)+(ytick(i)-ytick(i-1))./2;
    
88. %        yy =sqrt(3)/3*xx+b_shift;
    
89.          yy = dx*sin(xx/6*pi)+ sqrt(3)/3*xx+b_shift;
    
90.         patch([xx(:);flipud(xx(:))], ...
    
91.         [yy(:)+2*dy;flipud(yy(:)-2.*dy)], ...
    
92.             [1 1 1])
    
93.         
    
94.     case 'Patch',
    
95.         % this can be vectorized
    
96.         dx=(xlim(2)-xlim(1))./10;
    
97.         yy=repmat([y-2.*dy y-dy],1,6);
    
98.         xx=xlim(1)+dx.*[0:11];
    
99.         patch([xx(:);flipud(xx(:))], ...
    
100.             [yy(:);flipud(yy(:)-2.*dy)], ...
     
101.             [1 1 1])
     
102.     case 'RPatch',
     
103.                 % this can be vectorized
     
104.         dx=(xlim(2)-xlim(1))./100;
     
105.         yy=y+rand(101,1).*2.*dy;
     
106.         xx=xlim(1)+dx.*(0:100);
     
107.                 patch([xx(:);flipud(xx(:))], ...
     
108.             [yy(:);flipud(yy(:)-2.*dy)], ...
     
109.             [.8 .8 .8])
     
110.     case 'Line',
     
111.                 line([x-dx x   ],[y-2.*dy y-dy   ]);
     
112.                 line([x    x+dx],[y+dy    y+2.*dy]);
     
113.                 line([x-dx x   ],[y-3.*dy y-2.*dy]);
     
114.                 line([x    x+dx],[y+2.*dy y+3.*dy]);
     
115. end;
     
116. set(gca,'xlim',xlim);
     
117. 
     
118. % map back
     
119. ytick(ytick>y_break_start)=ytick(ytick>y_break_start)+y_break_mid;
     
120. for i=1:length(ytick)
     
121.    yticklabel{i}=sprintf('%4.3d',ytick(i));
     
122. end;
     
123. set(gca,'yticklabel',yticklabel);

复制代码

---

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

Powered by Discuz! X3.2