質問編集履歴
1
ご指摘にあったとおりにコードを貼り付けました
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -18,23 +18,39 @@
|
|
|
18
18
|
|
|
19
19
|
### 該当のソースコード
|
|
20
20
|
|
|
21
|
-
|
|
21
|
+
|
|
22
22
|
|
|
23
23
|
```
|
|
24
24
|
|
|
25
|
+
Fs = 5000;
|
|
26
|
+
|
|
25
|
-
|
|
27
|
+
``` > % サンプリング周波数は5000Hz
|
|
26
28
|
|
|
27
29
|
|
|
28
30
|
|
|
29
31
|
```
|
|
30
32
|
|
|
33
|
+
N = 2^18;
|
|
31
34
|
|
|
32
35
|
|
|
36
|
+
|
|
37
|
+
``` > % データの個数
|
|
38
|
+
|
|
39
|
+
|
|
40
|
+
|
|
41
|
+
```
|
|
42
|
+
|
|
43
|
+
i = (0 : N-1);
|
|
44
|
+
|
|
45
|
+
|
|
46
|
+
|
|
33
|
-
|
|
47
|
+
``` > % データの個数に対応するiを準備
|
|
34
48
|
|
|
35
49
|
|
|
36
50
|
|
|
37
51
|
|
|
52
|
+
|
|
53
|
+
```
|
|
38
54
|
|
|
39
55
|
C = zeros(N,1);
|
|
40
56
|
|
|
@@ -42,19 +58,29 @@
|
|
|
42
58
|
|
|
43
59
|
for k = 0 : 1: N-1
|
|
44
60
|
|
|
61
|
+
ex = exp(complex(0,-(2*pi/N)*k*i))';
|
|
45
62
|
|
|
63
|
+
```> %exp(-j*(2*pi/N)*k*i)に対応
|
|
46
64
|
|
|
47
|
-
|
|
65
|
+
|
|
48
66
|
|
|
49
67
|
|
|
50
68
|
|
|
69
|
+
```
|
|
70
|
+
|
|
51
|
-
|
|
71
|
+
C(k+1,1) = 1/N * sum(cur1.*ex);
|
|
72
|
+
|
|
73
|
+
end
|
|
74
|
+
|
|
75
|
+
```
|
|
76
|
+
|
|
77
|
+
>% 離散フーリエ変換を実行, cur1は測定で得られた電流値
|
|
52
78
|
|
|
53
79
|
|
|
54
80
|
|
|
55
81
|
|
|
56
82
|
|
|
57
|
-
|
|
83
|
+
|
|
58
84
|
|
|
59
85
|
|
|
60
86
|
|