質問編集履歴
1
ご指摘にあったとおりにコードを貼り付けました
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -8,27 +8,40 @@
|
|
|
8
8
|
要求された 262144x262144 (512.0GB) 配列は、最大配列サイズの設定を超えています。この制限より大きい配列を作成すると、処理に時間がかかり、MATLAB が反応しなくなることがあります。詳細については、配列サイズの制限または設定パネルを参照してください。
|
|
9
9
|
|
|
10
10
|
### 該当のソースコード
|
|
11
|
-
|
|
11
|
+
|
|
12
12
|
```
|
|
13
|
+
Fs = 5000;
|
|
13
|
-
|
|
14
|
+
``` > % サンプリング周波数は5000Hz
|
|
14
15
|
|
|
15
16
|
```
|
|
17
|
+
N = 2^18;
|
|
16
18
|
|
|
17
|
-
|
|
19
|
+
``` > % データの個数
|
|
18
20
|
|
|
21
|
+
```
|
|
22
|
+
i = (0 : N-1);
|
|
19
23
|
|
|
24
|
+
``` > % データの個数に対応するiを準備
|
|
25
|
+
|
|
26
|
+
|
|
27
|
+
```
|
|
20
28
|
C = zeros(N,1);
|
|
21
29
|
|
|
22
30
|
for k = 0 : 1: N-1
|
|
23
|
-
|
|
24
|
-
|
|
31
|
+
ex = exp(complex(0,-(2*pi/N)*k*i))';
|
|
32
|
+
```> %exp(-j*(2*pi/N)*k*i)に対応
|
|
33
|
+
|
|
25
34
|
|
|
35
|
+
```
|
|
26
|
-
|
|
36
|
+
C(k+1,1) = 1/N * sum(cur1.*ex);
|
|
37
|
+
end
|
|
38
|
+
```
|
|
39
|
+
>% 離散フーリエ変換を実行, cur1は測定で得られた電流値
|
|
27
40
|
|
|
28
41
|
|
|
29
|
-
end
|
|
30
42
|
|
|
31
43
|
|
|
44
|
+
|
|
32
45
|
### 試したこと
|
|
33
46
|
|
|
34
47
|
直接法を用いた計算(k*i行列)をしたところメモリ不足で計算できなかったため
|