SHOW:
|
|
- or go back to the newest paste.
1 | from scipy import signal | |
2 | ||
3 | x = [0.49, 0.58, 0.49, 0.29, 0.1, -0.19, -0.35, -0.39, -0.29, -0.15, 0.2, | |
4 | 0.49, 0.68, 0.83, 0.78, 0.59, 0.29, -0.1, -0.43, -0.79, -0.92, -0.93, | |
5 | -0.73, -0.44, -0.05, 0.24, 0.59, 0.63, 0.64, 0.39, 0.05, -0.25, -0.58, | |
6 | -0.78, -0.79, -0.68, -0.39, -0.05, 0.29, 0.59, 0.68, 0.74, 0.58, 0.34, | |
7 | 0.0, -0.24, -0.49, -0.58, -0.54, -0.34, -0.05, 0.29, 0.49, 0.63, 0.64, | |
8 | 0.44, 0.14, -0.29, -0.58, -0.93, -1.03, -0.97, -0.69, -0.39, 0.05, 0.49, | |
9 | 0.73, 0.88, 0.88, 0.68, 0.34, -0.04, -0.4, -0.58, -0.73, -0.64, -0.44, | |
10 | -0.14, 0.24, 0.44, 0.68, 0.64, 0.58, 0.3, 0.05, -0.3, -0.53, -0.59, | |
11 | -0.59, -0.34, -0.14, 0.19, 0.39, 0.49, 0.49, 0.34, 0.05, -0.29, -0.59, | |
12 | -0.78, -0.88, -0.78, -0.49, -0.15, 0.3, 0.63, 0.88, 1.03, 0.97, 0.73, | |
13 | 0.44, 0.05, -0.29, -0.49, -0.58, -0.49, -0.3, 0.0, 0.3, 0.53, 0.64, 0.58, | |
14 | 0.4, 0.14, -0.24, -0.59, -0.83, -0.97, -0.98, -0.88, -0.54, -0.29, -0.05, | |
15 | 0.1, 0.1, -0.1, -0.34, -0.74, -1.12, -1.41, -1.47, -1.51, -1.22, -0.88, | |
16 | -0.44, 0.0, 0.34, 0.49, 0.58, 0.39, 0.2, -0.15, -0.44, -0.58, -0.64, | |
17 | -0.53, -0.25, 0.05, 0.44, 0.73, 0.88, 0.88, 0.68, 0.49, 0.1, -0.2, -0.39, | |
18 | -0.48, -0.44, -0.2, 0.15, 0.49, 0.83, 1.02, 0.98, 0.97, 0.74, 0.39, 0.05, | |
19 | -0.25, -0.39, -0.44, -0.29, -0.05, 0.24, 0.49, 0.64, 0.68, 0.54, 0.24, | |
20 | -0.24, -0.64, -0.97, -1.22, -1.18, -1.12, -0.83, -0.39, 0.05, 0.39, 0.63, | |
21 | 0.74, 0.73, 0.49, 0.34, 0.1, -0.05, -0.15, 0.0, 0.2, 0.44, 0.78, 0.97, | |
22 | 1.13, 1.12, 0.98, 0.73, 0.44, 0.15, -0.05, -0.05, 0.05, 0.39, 0.73, 1.12, | |
23 | 1.47, 1.61, 1.61, 1.42, 1.12, 0.63, 0.15, -0.24, -0.49, -0.59, -0.49, | |
24 | -0.24, 0.1, 0.39, 0.58, 0.69, 0.58, 0.25, -0.1, -0.54, -0.92, -1.27, | |
25 | -1.37, -1.37, -1.12, -0.88, -0.54, -0.24, -0.1, 0.0, -0.15, -0.39, -0.68, | |
26 | -0.98, -1.17, -1.22, -1.17, -0.88, -0.49, -0.05, 0.3, 0.58, 0.73, 0.69, | |
27 | 0.44, 0.09, -0.29, -0.68, -0.83, -0.98, -0.88, -0.63, -0.3, 0.05, 0.35, | |
28 | 0.48, 0.54, 0.44, 0.24, -0.04, -0.25, -0.44, -0.49, -0.39, -0.14, 0.14, | |
29 | 0.49, 0.83, 0.93, 0.98, 0.83, 0.53, 0.15, -0.2, -0.53, -0.59, -0.58, | |
30 | -0.4, -0.04, 0.24, 0.59, 0.73, 0.78, 0.63, 0.25, -0.2, -0.68, -1.07, | |
31 | -1.42, -1.47, -1.46, -1.22, -0.93, -0.63, -0.35, -0.19, -0.15, -0.29, | |
32 | -0.49, -0.73, -0.88, -1.03, -0.92, -0.74, -0.43, 0.0, 0.43, 0.69, 0.93, | |
33 | 0.92, 0.78, 0.49, 0.05, -0.24, -0.59, -0.73, -0.68, -0.54, -0.25, 0.1, | |
34 | 0.44, 0.64, 0.78, 0.63, 0.44, 0.1, -0.24, -0.59, -0.68, -0.74, -0.48, | |
35 | -0.25, 0.15, 0.58, 0.83, 0.98, 0.98, 0.83, 0.58, 0.2, -0.15, -0.44, | |
36 | -0.58, -0.54, -0.39, -0.1, 0.15, 0.44, 0.53, 0.54, 0.29, 0.05, -0.34, | |
37 | -0.73, -0.98, -1.02, -1.03, -0.73, -0.39, 0.0, 0.39, 0.68, 0.78, 0.78, | |
38 | 0.49, 0.2, -0.2, -0.49, -0.73, -0.73, -0.63, -0.35, 0.05, 0.34, 0.64, | |
39 | 0.83, 0.73, 0.59, 0.34, 0.0, -0.25, -0.43, -0.54, -0.39, -0.15, 0.2, | |
40 | 0.53, 0.83, 0.98, 0.98, 0.88, 0.58, 0.29, -0.04, -0.35, -0.39, -0.34, | |
41 | -0.14, 0.19, 0.49, 0.78, 1.03, 1.12, 0.98, 0.68, 0.39, -0.05, -0.29, | |
42 | -0.49, -0.49, -0.39, -0.05, 0.25, 0.58, 0.83, 0.88, 0.83, 0.59, 0.24, | |
43 | -0.1, -0.49, -0.73, -0.83, -0.83, -0.58, -0.35, 0.0, 0.3, 0.39, 0.44, | |
44 | 0.29, 0.0, -0.34, -0.68, -0.93, -1.03, -0.97, -0.69, -0.39, 0.0, 0.34, | |
45 | 0.59, 0.68, 0.64, 0.44, 0.14, -0.24, -0.44, -0.64, -0.58, -0.34, -0.1, | |
46 | 0.34, 0.64, 0.97, 1.13, 0.97, 0.78, 0.39, -0.04, -0.4, -0.68, -0.83, | |
47 | -0.73, -0.54, -0.15, 0.15, 0.49, 0.63, 0.64, 0.54, 0.24, -0.15, -0.48, | |
48 | -0.79, -0.92, -0.98, -0.78, -0.54, -0.24, -0.05, 0.19, 0.25, 0.14, -0.09, | |
49 | -0.44, -0.83, -1.03, -1.22, -1.22, -1.02, -0.69, -0.29, 0.05, 0.39, 0.49, | |
50 | 0.53, 0.35, 0.09, -0.19, -0.44, -0.69, -0.68, -0.63, -0.35, -0.04, 0.29, | |
51 | 0.58, 0.79, 0.78, 0.63, 0.29, -0.04, -0.35, -0.63, -0.73, -0.69, -0.39, | |
52 | -0.09, 0.24, 0.63, 0.83, 0.93, 0.83, 0.59, 0.29, -0.1, -0.34, -0.54, | |
53 | -0.58, -0.44, -0.2, 0.15, 0.39, 0.64, 0.68, 0.63, 0.44, 0.1, -0.24, | |
54 | -0.49, -0.73, -0.69, -0.53, -0.25, 0.2, 0.63, 0.98, 1.22, 1.22, 1.07, | |
55 | 0.69, 0.39, -0.05, -0.34, -0.49, -0.44, -0.34, -0.05, 0.24, 0.54, 0.83, | |
56 | 0.83, 0.73, 0.49, 0.1, -0.2, -0.53, -0.64, -0.68, -0.59, -0.24, 0.05, | |
57 | 0.34, 0.63, 0.74, 0.73, 0.54, 0.29, 0.05, -0.3, -0.44, -0.53, -0.49, | |
58 | -0.29, 0.0, 0.19, 0.49, 0.54, 0.49, 0.34, 0.0, -0.34, -0.64, -0.88, | |
59 | -0.93, -0.78, -0.58, -0.2, 0.2, 0.53, 0.74, 0.83, 0.78, 0.53, 0.25, | |
60 | -0.15, -0.44, -0.58, -0.64, -0.49, -0.24, 0.05, 0.39, 0.63, 0.69, 0.63, | |
61 | 0.39, 0.15, -0.24, -0.49, -0.74, -0.78, -0.73, -0.44, -0.19, 0.04, 0.2, | |
62 | 0.39, 0.34, 0.15, -0.15, -0.49, -0.73, -1.02, -1.08, -0.97, -0.69, -0.29, | |
63 | 0.1, 0.53, 0.78, 0.98, 0.88, 0.73, 0.39, 0.05, -0.29, -0.44, -0.44, | |
64 | -0.39, -0.1, 0.15, 0.44, 0.58, 0.64, 0.58, 0.4, 0.09, -0.19, -0.49, | |
65 | -0.54, -0.63, -0.49, -0.2, 0.25, 0.58, 0.83, 1.08, 1.02, 0.79, 0.43, 0.0, | |
66 | -0.39, -0.73, -0.93, -0.92, -0.64, -0.39, 0.0, 0.39, 0.54, 0.68, 0.59, | |
67 | 0.34, 0.05, -0.34, -0.69, -0.83, -0.78, -0.73, -0.44, -0.15, 0.2, 0.49, | |
68 | 0.58, 0.54, 0.34, 0.1, -0.2, -0.39, -0.58, -0.59, -0.49, -0.19, 0.1, | |
69 | 0.39, 0.58, 0.69, 0.58, 0.39, 0.1, -0.24, -0.59, -0.68, -0.83, -0.69, | |
70 | -0.39, -0.1, 0.3, 0.58, 0.74, 0.78, 0.53, 0.25, -0.1, -0.49, -0.83, | |
71 | -1.07, -1.03, -0.88, -0.58, -0.2, 0.1, 0.44, 0.54, 0.58, 0.34, 0.1, | |
72 | -0.19, -0.49, -0.69, -0.78, -0.63, -0.44, -0.15, 0.2, 0.49, 0.53, 0.54, | |
73 | 0.39, 0.05, -0.29, -0.59, -0.78, -0.83, -0.63, -0.4, 0.05, 0.44, 0.88, | |
74 | 1.12, 1.23, 1.17, 0.97, 0.64, 0.29, 0.0, -0.14, -0.15, 0.0, 0.29, 0.59, | |
75 | 0.93, 1.22] | |
76 | ||
77 | f, Pxx = signal.periodogram(x, 5, nfft=512) | |
78 | print ['%.4f' % xx for xx in Pxx[0:10]] | |
79 | # the output will be: | |
80 | - | # ['0.0000', '0.0105', '0.5567', '0.8369', '0.3156', '0.4476', '0.0380', '0.4475', '3.1456', '1.0430'] |
80 | + | # ['0.0000', '0.0105', '0.5567', '0.8369', '0.3156', '0.4476', '0.0380', '0.4475', '3.1456', '1.0430'] |
81 | # |