teratail
回答率
85.48%
質問するログイン新規登録
トップに関する質問PythonでFFTの際にローパスフィルタをかけたいです。

編集履歴

質問編集履歴

3

コード修正

2021/01/01 14:25

投稿

退会済みユーザー
test CHANGED
File without changes
test CHANGED
@@ -22,7 +22,7 @@
22
22
 
23
23
  from scipy import integrate
24
24
 
25
-
25
+ # -*- coding: utf-8 -*-
26
26
 
27
27
  import numpy as np
28
28
 
@@ -82,6 +82,8 @@
82
82
 
83
83
  SpecHeader = 'f[Hz],E'
84
84
 
85
+ # データのパラメータ
86
+
85
87
  SN = 524288 # サンプル数
86
88
 
87
89
  fs = 10000 # カットオフ周波数 S/ sec
@@ -130,7 +132,11 @@
130
132
 
131
133
 
132
134
 
135
+ # calculate energy spectrum
136
+
133
137
  freq1, P1 = signal.welch(f, fs, nperseg=N / 8)
138
+
139
+ # freq2, P2 = signal.welch(y, fs, nperseg=n/64)
134
140
 
135
141
  SpectrumResult = np.array([freq1, P1]).T
136
142
 

2

インデントの修正

2021/01/01 14:25

投稿

退会済みユーザー
test CHANGED
File without changes
test CHANGED
@@ -5,6 +5,8 @@
5
5
  以下にそのコードを示しました(同時にPDFも求めるコードとなっています)。コードを書いていただけると嬉しい限りです。
6
6
 
7
7
 
8
+
9
+ ```python
8
10
 
9
11
  import numpy as np
10
12
 
@@ -153,3 +155,5 @@
153
155
  plt.ylabel("Power")
154
156
 
155
157
  plt.title("Energy spectrum")
158
+
159
+ ```

1

誤字

2021/01/01 13:43

投稿

退会済みユーザー
test CHANGED
File without changes
test CHANGED
@@ -20,7 +20,7 @@
20
20
 
21
21
  from scipy import integrate
22
22
 
23
- # -*- coding: utf-8 -*-
23
+
24
24
 
25
25
  import numpy as np
26
26
 
@@ -80,8 +80,6 @@
80
80
 
81
81
  SpecHeader = 'f[Hz],E'
82
82
 
83
- # データのパラメータ
84
-
85
83
  SN = 524288 # サンプル数
86
84
 
87
85
  fs = 10000 # カットオフ周波数 S/ sec
@@ -130,11 +128,7 @@
130
128
 
131
129
 
132
130
 
133
- # calculate energy spectrum
134
-
135
131
  freq1, P1 = signal.welch(f, fs, nperseg=N / 8)
136
-
137
- # freq2, P2 = signal.welch(y, fs, nperseg=n/64)
138
132
 
139
133
  SpectrumResult = np.array([freq1, P1]).T
140
134
 
@@ -159,7 +153,3 @@
159
153
  plt.ylabel("Power")
160
154
 
161
155
  plt.title("Energy spectrum")
162
-
163
-
164
-
165
- #plt.show()