質問編集履歴
3
コード修正
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
インデントの修正
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
誤字
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -20,7 +20,7 @@
|
|
|
20
20
|
|
|
21
21
|
from scipy import integrate
|
|
22
22
|
|
|
23
|
-
|
|
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()
|