回答編集履歴
1
fixed code
test
CHANGED
|
@@ -24,7 +24,7 @@
|
|
|
24
24
|
w, h = signal.freqz(b, a, fs = fs)
|
|
25
25
|
return signal.filtfilt(b, a, x)
|
|
26
26
|
|
|
27
|
-
def fft(x
|
|
27
|
+
def fft(x):
|
|
28
28
|
spec = fftpack.fft(x)
|
|
29
29
|
amp = abs(spec)
|
|
30
30
|
amp = amp / (len(x) / 2)
|
|
@@ -36,7 +36,7 @@
|
|
|
36
36
|
axs[0, 0].set_title('original wave')
|
|
37
37
|
axs[0, 0].set_ylabel('amp')
|
|
38
38
|
|
|
39
|
-
amp, freq = fft(data
|
|
39
|
+
amp, freq = fft(data)
|
|
40
40
|
axs[0, 1].plot(freq[:fs // 20], amp[:fs // 20])
|
|
41
41
|
axs[0, 1].set_title('original spectrum')
|
|
42
42
|
axs[0, 1].set_ylabel('amp')
|
|
@@ -47,7 +47,7 @@
|
|
|
47
47
|
axs[1, 0].set_xlabel('time[s]')
|
|
48
48
|
axs[1, 0].set_ylabel('amp')
|
|
49
49
|
|
|
50
|
-
amp, freq = fft(hp
|
|
50
|
+
amp, freq = fft(hp)
|
|
51
51
|
axs[1, 1].plot(freq[:fs // 20], amp[:fs // 20])
|
|
52
52
|
axs[1, 1].set_title('high passed spectrum')
|
|
53
53
|
axs[1, 1].set_xlabel('frequency[Hz]')
|