最后活跃于 1 month ago

Paras aulas de física experimental.

修订 fc51c7fc73474a7c40ebf95d182a4317f6055684

generate.py 原始文件
1#!/usr/bin/env python
2# -*- coding: utf-8 -*-
3from __future__ import unicode_literals
4
5import sys
6import numpy as np
7import matplotlib.mlab as mlab
8import matplotlib.pyplot as plt
9
10width = 20 # cm
11height = 15 # cm
12
13if len(sys.argv) < 4:
14 print "Usage: generate.py data.txt title output"
15 exit(1)
16
17output = sys.argv[3]
18title = sys.argv[2].decode(encoding='UTF-8',errors='strict')
19
20data = np.fromfile(sys.argv[1], sep='\n')
21
22num_bins = 28
23avg = np.average(data)
24std = np.std(data)
25
26fig, ax = plt.subplots()
27n, bins, patches = ax.hist(data, num_bins, color='burlywood', histtype='stepfilled')
28
29print "Média: %s Desvio Padrão: %s" %(avg, std)
30
31ax.axvline(avg, color='red', label="Média = %.2f s" %avg)
32ax.axvline(avg-std, color='lightblue', label="Média - Desvio Padrão = %.2f s" %(avg-std))
33ax.axvline(avg+std, color='darkblue', label="Média + Desvio Padrão = %.2f s" % (avg+std))
34
35ax.set_xlabel('Tempo (s)')
36ax.set_ylabel('Contagem (n)')
37ax.set_title(title)
38ax.set_ylim( None, n.max() * 1.2)
39
40legend = ax.legend(loc='upper left', shadow=True, prop={'size':10})
41fig.tight_layout()
42
43fig.set_size_inches(width / 2.54, height / 2.54)
44plt.savefig('%s.png' %output, dpi=100)