質問編集履歴
1
15sの窓でパワースペクトル算出の部分を追加しました。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -78,6 +78,34 @@
|
|
|
78
78
|
|
|
79
79
|
|
|
80
80
|
|
|
81
|
+
|
|
82
|
+
|
|
83
|
+
|
|
84
|
+
|
|
85
|
+
######################################################追加部分です
|
|
86
|
+
|
|
87
|
+
#15sの窓でパワースペクトル算出
|
|
88
|
+
|
|
89
|
+
time_interval = 15
|
|
90
|
+
|
|
91
|
+
loop_len = list(range(round(len(img_mean)/(fs*time_interval))))
|
|
92
|
+
|
|
93
|
+
fft_result = []
|
|
94
|
+
|
|
95
|
+
for n,i in enumerate(loop_len):
|
|
96
|
+
|
|
97
|
+
if n ==0:
|
|
98
|
+
|
|
99
|
+
img = img_mean[:time_interval*fs]
|
|
100
|
+
|
|
101
|
+
fft_result.append(signal_fft(img,fs))
|
|
102
|
+
|
|
103
|
+
else :
|
|
104
|
+
|
|
105
|
+
img = img_mean[time_interval*fs*n : time_interval*fs*(n+1)]
|
|
106
|
+
|
|
107
|
+
fft_result.append(signal_fft(img,fs))
|
|
108
|
+
|
|
81
109
|
# パワースペクトルグラフ化
|
|
82
110
|
|
|
83
111
|
#print(fft_result)
|