回答編集履歴
4
spanはNの間違いのため修正
test
CHANGED
|
@@ -116,7 +116,7 @@
|
|
|
116
116
|
|
|
117
117
|
```
|
|
118
118
|
|
|
119
|
-
ここで注意する必要があるのはK[0]は単純な複素数ではなく、フーリエ変換済みの結果であり、要素
|
|
119
|
+
ここで注意する必要があるのはK[0]は単純な複素数ではなく、フーリエ変換済みの結果であり、要素をN個を持つ配列であることです(numpyのndarray)。numpyの様々な関数は引数に配列を渡しても演算可能で、結果を配列で返してくれます。すなわち
|
|
120
120
|
|
|
121
121
|
```python
|
|
122
122
|
|
|
@@ -134,7 +134,7 @@
|
|
|
134
134
|
|
|
135
135
|
...
|
|
136
136
|
|
|
137
|
-
amp[0][
|
|
137
|
+
amp[0][N-1] = np.sqrt(K[0][N-1].real ** 2 + K[0][N-1].imag ** 2)
|
|
138
138
|
|
|
139
139
|
|
|
140
140
|
|
3
より正確にリストと配列を区別
test
CHANGED
|
@@ -52,7 +52,7 @@
|
|
|
52
52
|
|
|
53
53
|
|
|
54
54
|
|
|
55
|
-
関数fourier()で行っていることは、|<---N個--->|を先頭から順々にフーリエ変換し、結果を
|
|
55
|
+
関数fourier()で行っていることは、|<---N個--->|を先頭から順々にフーリエ変換し、結果をリストに格納しています。FFTのアルゴリズム上、Nは2のn乗(オリジナルのサイトはN=2の8乗=256)である必要があります。
|
|
56
56
|
|
|
57
57
|
|
|
58
58
|
|
|
@@ -86,7 +86,7 @@
|
|
|
86
86
|
|
|
87
87
|
```
|
|
88
88
|
|
|
89
|
-
先に修正した関数fourier()に、このmonoを渡してフーリエ変換を行
|
|
89
|
+
先に修正した関数fourier()に、このmonoを渡してフーリエ変換を行うと、span個のフーリエ変換済みのリストが返ってきます。
|
|
90
90
|
|
|
91
91
|
```python
|
|
92
92
|
|
|
@@ -94,7 +94,7 @@
|
|
|
94
94
|
|
|
95
95
|
```
|
|
96
96
|
|
|
97
|
-
print(len(K))とやると結果を格納した変数`K`がspan個の要素を持つ
|
|
97
|
+
print(len(K))とやると結果を格納した変数`K`がspan個の要素を持つことがわかります。そうすると振幅を計算する部分
|
|
98
98
|
|
|
99
99
|
```python
|
|
100
100
|
|
2
リストを配列に変更
test
CHANGED
|
@@ -52,7 +52,7 @@
|
|
|
52
52
|
|
|
53
53
|
|
|
54
54
|
|
|
55
|
-
関数fourier()で行っていることは、|<---N個--->|を先頭から順々にフーリエ変換し、結果を
|
|
55
|
+
関数fourier()で行っていることは、|<---N個--->|を先頭から順々にフーリエ変換し、結果を配列に格納しています。FFTのアルゴリズム上、Nは2のn乗(オリジナルのサイトはN=2の8乗=256)である必要があります。
|
|
56
56
|
|
|
57
57
|
|
|
58
58
|
|
|
@@ -94,7 +94,7 @@
|
|
|
94
94
|
|
|
95
95
|
```
|
|
96
96
|
|
|
97
|
-
print(len(K))とやると結果を格納した変数`K`がspan個の
|
|
97
|
+
print(len(K))とやると結果を格納した変数`K`がspan個の要素を持つ配列であることがわかります。そうすると振幅を計算する部分
|
|
98
98
|
|
|
99
99
|
```python
|
|
100
100
|
|
|
@@ -116,7 +116,7 @@
|
|
|
116
116
|
|
|
117
117
|
```
|
|
118
118
|
|
|
119
|
-
ここで注意する必要があるのはK[0]は単純な複素数ではなく、フーリエ変換済みの結果であり、
|
|
119
|
+
ここで注意する必要があるのはK[0]は単純な複素数ではなく、フーリエ変換済みの結果であり、要素数span個を持つ配列であることです(numpyのndarray)。numpyの様々な関数は引数に配列を渡しても演算可能で、結果を配列で返してくれます。すなわち
|
|
120
120
|
|
|
121
121
|
```python
|
|
122
122
|
|
1
typoを編集
test
CHANGED
|
@@ -106,13 +106,13 @@
|
|
|
106
106
|
|
|
107
107
|
```python
|
|
108
108
|
|
|
109
|
-
amp[0] = np.sqrt(K[0].real ** 2 + K[0].
|
|
109
|
+
amp[0] = np.sqrt(K[0].real ** 2 + K[0].imag ** 2)
|
|
110
|
-
|
|
110
|
+
|
|
111
|
-
amp[1] = np.sqrt(K[1].real ** 2 + K[1].
|
|
111
|
+
amp[1] = np.sqrt(K[1].real ** 2 + K[1].imag ** 2)
|
|
112
112
|
|
|
113
113
|
...
|
|
114
114
|
|
|
115
|
-
amp[span-1] = np.sqrt(K[span-1].real ** 2 + K[span-1].
|
|
115
|
+
amp[span-1] = np.sqrt(K[span-1].real ** 2 + K[span-1].imag ** 2)
|
|
116
116
|
|
|
117
117
|
```
|
|
118
118
|
|
|
@@ -120,7 +120,7 @@
|
|
|
120
120
|
|
|
121
121
|
```python
|
|
122
122
|
|
|
123
|
-
amp[0] = np.sqrt(K[0].real ** 2 + K[0].
|
|
123
|
+
amp[0] = np.sqrt(K[0].real ** 2 + K[0].imag ** 2)
|
|
124
124
|
|
|
125
125
|
```
|
|
126
126
|
|
|
@@ -128,13 +128,13 @@
|
|
|
128
128
|
|
|
129
129
|
```python
|
|
130
130
|
|
|
131
|
-
amp[0][0] = np.sqrt(K[0][0].real ** 2 + K[0][0].
|
|
131
|
+
amp[0][0] = np.sqrt(K[0][0].real ** 2 + K[0][0].imag ** 2)
|
|
132
|
-
|
|
132
|
+
|
|
133
|
-
amp[0][1] = np.sqrt(K[0][1].real ** 2 + K[0][1].
|
|
133
|
+
amp[0][1] = np.sqrt(K[0][1].real ** 2 + K[0][1].imag ** 2)
|
|
134
134
|
|
|
135
135
|
...
|
|
136
136
|
|
|
137
|
-
amp[0][span-1] = np.sqrt(K[0][span-1].real ** 2 + K[0][span-1].
|
|
137
|
+
amp[0][span-1] = np.sqrt(K[0][span-1].real ** 2 + K[0][span-1].imag ** 2)
|
|
138
138
|
|
|
139
139
|
|
|
140
140
|
|