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 | # |