豆粑粑 matlab zoom 局部放大

meatball1982

这两天只画图了。
一个图，局部放大。
之前作过，但这次的要求是要画两条线，标注两个关键的点。在matlab 中，通过一个axes, 一个新的axes, 然后在第一个axes上　axis off实现。
麻烦是麻烦了些，但效果还不错。对数据的展示，有不错的效果。


两个不同的版本，分别对应的着框的小角和线的端点。
程序在附件，主要麻烦的地方在于把小图中点对应到大图中的作标中。

meatball1982

更新一下，通过设置第二个图为透明，直接显示第一个图中的线。

效果同附件。

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. %% outline
   
6. % zoom in with line cross different axis
   
7. 
   
8. %% main
   
9. x=linspace(0,3*pi,1000);
   
10. rnd = rand(size(x));
    
11. 
    
12. y=x+exp(sin(2*x))+exp(x./3)+exp(x).*rnd/100;
    
13. seedm(3)
    
14. 
    
15. zoom =[3.8 4.3 8 14]
    
16. pos1 = [0.1 0.15 0.85 0.7];
    
17. pos2 = [0.2 0.4 0.3 0.25];
    
18. 
    
19. 
    
20. % find line1 ends
    
21. 
    
22. tm = x-zoom(1);
    
23. [va,inbeg] = min(abs(tm));
    
24. 
    
25. tm = x-zoom(2);
    
26. [va,inend] = min(abs(tm));
    
27. 
    
28. 
    
29. 
    
30. h=figure(1)
    
31. set(h, 'Position', [1000, 1000, 550, 500]);
    
32. ax11=axes('position',pos1);
    
33. hold on
    
34. plot(x,y,'.-')
    
35. axis([0 10 -20 200])
    
36. set(gca,'fontsize',14,'linewidth',2)
    
37. 
    
38. % % % line : from corner to corner -----------------------------------------
    
39. % lin1_x = [zoom(1) (0.2-0.1)/0.85 * (10-0)];
    
40. % lin1_y = [zoom(3) (0.4-0.15)/0.7 * (200-(-20))-20];
    
41. % lin2_x = [zoom(2) (0.2+0.3-0.1)/0.85*(10-0)];
    
42. % lin2_y = [zoom(3) (0.4-0.15)/0.7 * (200-(-20))-20];
    
43. 
    
44. % % line : from line to line -----------------------------------------
    
45. lin1_x = [zoom(1) (0.2-0.1)/0.85 * (10-0)];
    
46. lin1_y = [y(inbeg) ((y(inbeg)-8)*0.25/(14-8)+(0.4-0.15))/0.7 * (200-(-20))-20 ];
    
47. lin2_x = [zoom(2) (0.2+0.3-0.1)/0.85*(10-0)];
    
48. lin2_y = [y(inend) ((y(inend)-8)*0.25/(14-8)+(0.4-0.15))/0.7 * (200-(-20))-20];
    
49. plot(lin1_x,lin1_y,'r-','linewidth',2)
    
50. plot(lin2_x,lin2_y,'g-','linewidth',2)
    
51. axis([0 10 -20 200])
    
52. % axis off
    
53. box on
    
54. 
    
55. rectangle('Position',[zoom(1) zoom(3) zoom(2)-zoom(1) zoom(4)-zoom(3)],'linewidth',2)
    
56. 
    
57. ax11=axes('position',pos2);
    
58. plot(x,y,'.-','linewidth',2,'markersize',12)
    
59. axis(zoom)
    
60. set(gca,'linewidth',2)
    
61. 
    
62. set(gcf, 'Color','none')
    
63. set(gca,'color','none')
    
64. 
    
65. 
    
66. % ax11=axes('position',pos1);
    
67. % hold on
    
68. % plot(lin1_x,lin1_y,'r-','linewidth',2)
    
69. % plot(lin2_x,lin2_y,'g-','linewidth',2)
    
70. % axis([0 10 -20 200])
    
71. % axis off

复制代码