質問編集履歴
1
15sの窓でパワースペクトル算出の部分を追加しました。
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -38,6 +38,20 @@
|
|
|
38
38
|
f = np.arange(0,fs/2,fs/n)
|
|
39
39
|
return power , f
|
|
40
40
|
|
|
41
|
+
|
|
42
|
+
|
|
43
|
+
######################################################追加部分です
|
|
44
|
+
#15sの窓でパワースペクトル算出
|
|
45
|
+
time_interval = 15
|
|
46
|
+
loop_len = list(range(round(len(img_mean)/(fs*time_interval))))
|
|
47
|
+
fft_result = []
|
|
48
|
+
for n,i in enumerate(loop_len):
|
|
49
|
+
if n ==0:
|
|
50
|
+
img = img_mean[:time_interval*fs]
|
|
51
|
+
fft_result.append(signal_fft(img,fs))
|
|
52
|
+
else :
|
|
53
|
+
img = img_mean[time_interval*fs*n : time_interval*fs*(n+1)]
|
|
54
|
+
fft_result.append(signal_fft(img,fs))
|
|
41
55
|
# パワースペクトルグラフ化
|
|
42
56
|
#print(fft_result)
|
|
43
57
|
time = np.arange(0, len(fft_result[0][0]), 1)
|