Mathematica中文论坛-非官方

标题: 给别人画了一个不错的图 [打印本页]

---

作者: meatball1982    时间: 2017-5-22 16:43
标题: 给别人画了一个不错的图
是统计轨迹，并画一定范围之内数据的高斯分布。
真TMD优雅，我靠。

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. dat=load('../data/combine.dat');
   
6. n=length(dat(:,1));
   
7. x=(1:n)/10;
   
8. dat1=dat(500:end,1);
   
9. dat2=dat(500:end,2);
   
10. 
    
11. bi = linspace(-13,-6,70);
    
12. [hi1,bi1]=hist(dat1,bi);
    
13. [hi2,bi2]=hist(dat2,bi);
    
14. 
    
15. 
    
16. fx=@(b,x)b(1)/sqrt(2*pi*b(2))*exp(-(x-b(3)).^2/(2*b(2)^2));
    
17. x_bin=bi;
    
18. x1=linspace(min(x_bin),max(x_bin),200);
    
19. 
    
20. y1 = hi1;
    
21. b0=[35,2,-8];
    
22. b=b0;
    
23. for l=1:10
    
24.     % Solve nonlinear curve-fitting problems in least-squares sense.
    
25.     % just a guess.
    
26.     b=lsqcurvefit(fx,b,x_bin,y1);
    
27.     % Nonlinear regression,
    
28.     b=nlinfit(x_bin,y1,fx,b);
    
29.    
    
30. end
    
31. b1=b;
    
32. y_hi1 = fx(b1,x1);
    
33. 
    
34. 
    
35. y2 = hi2;
    
36. b0=[30,2,-11];
    
37. b=b0;
    
38. for l=1:10
    
39.     % Solve nonlinear curve-fitting problems in least-squares sense.
    
40.     % just a guess.
    
41.     b=lsqcurvefit(fx,b,x_bin,y2);
    
42.     % Nonlinear regression,
    
43.     b=nlinfit(x_bin,y2,fx,b);
    
44.    
    
45. end
    
46. b2=b;
    
47. y_hi2 = fx(b2,x1);
    
48. 
    
49. 
    
50.             
    
51. % hold on
    
52. % plot(bi1,hi1,'r.')
    
53. % plot(x1,y_hi1,'r-')
    
54. %
    
55. % plot(bi,hi2,'b.')
    
56. % plot(x1,y_hi2,'b-')
    
57. 
    
58. mm_font=14;
    
59. 
    
60. h=figure(1);
    
61. set(h, 'Position', [100, 100, 800, 400]);
    
62. % subplot(1,2,1)
    
63. 
    
64. pos1=[0.1 0.1 0.6 0.8];
    
65. subplot('position',pos1)
    
66. hold on
    
67. x_pat1=[ 50 100 100 50];
    
68. y_pat1=[-13.1 -13.1 -5.9 -5.9];
    
69. h_pat1=patch(x_pat1,y_pat1,'gray')
    
70. alpha(h_pat1,0.1);
    
71. plot(x,dat(:,1),'r-')
    
72. 
    
73. plot(x,dat(:,2),'b-')
    
74. 
    
75. axis tight
    
76. box on
    
77. axis ([0 100 -13.1 -5.9])
    
78. set(gca,'fontsize',mm_font);
    
79. set(gca,'xtick',[0:10:100]);
    
80. 
    
81. 
    
82. % subplot(1,2,2)
    
83. pos2=[0.7 0.1 0.2 0.8];
    
84. subplot('position',pos2)
    
85. hold on
    
86. 
    
87. x_pat2=[ 0 50 50 0];
    
88. y_pat2=[-13.1 -13.1 -5.9 -5.9];
    
89. h_pat2=patch(x_pat2,y_pat2,'gray')
    
90. alpha(h_pat2,0.1);
    
91. 
    
92. plot(y_hi1,x1,'r-','linewidth',2)
    
93. plot(y_hi2,x1,'b-','linewidth',2)
    
94. 
    
95. hi = [hi1;hi2]';
    
96. h_bar=barh(bi,hi);
    
97. h1=h_bar(1);
    
98. h2=h_bar(2);
    
99. 
    
100. set(h1,'facecolor','r','edgecolor','none','barwidth',2.5)
     
101. set(h2,'facecolor','b','edgecolor','none','barwidth',2.5)
     
102. axis([0 50 -13.1 -5.9])
     
103. set(gca,'ytick',[]);
     
104. 
     
105. 
     
106. text(20,-6.5,['\mu=',mat2str(floor(b1(3)*1000)/1000)],'fontsize',mm_font)
     
107. text(20,-7.0,['\sigma=',mat2str(floor(b1(2)*1000)/1000)],'fontsize',mm_font)
     
108. text(20,-12,['\mu=',mat2str(floor(b2(3)*1000)/1000)],'fontsize',mm_font)
     
109. text(20,-12.5,['\sigma=',mat2str(floor(b2(2)*1000)/1000)],'fontsize',mm_font)
     
110. 
     
111. 
     
112. set(gca,'fontsize',mm_font);
     
113. box on
     
114. 
     
115. h=gcf;
     
116. fi_na=['../file_img/fig_traj_hist'];
     
117. fun_work_li_035_myfig_out(h,fi_na,3);
     
118. 
     
119. 
     
120. 
     
121. % plot(dat(:,1),dat(:,2),'.')

复制代码

fig_traj_hist.png (432.17 KB, 下载次数: 1600)

下载附件

2017-5-22 16:43 上传

---

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

Powered by Discuz! X3.2