View difference between Paste ID: zPLGBTpn and C3LqLFAV
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
#