-
Notifications
You must be signed in to change notification settings - Fork 0
/
iec60063.m
executable file
·64 lines (63 loc) · 2.44 KB
/
iec60063.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
%
% iec60063.m
%
% Created by Léa Strobino.
% Copyright 2013 hepia. All rights reserved.
%
function R = iec60063(R,e)
if nargin < 2
e = 12;
end
switch round(e)
case 3
e = [220 470 1000];
case 6
e = [150 220 330 470 680 1000];
case 12
e = [120 150 180 220 270 330 390 470 560 680 820 1000];
case 24
e = [110 120 130 150 160 180 200 220 240 270 300 330 ...
360 390 430 470 510 560 620 680 750 820 910 1000];
case 48
e = [105 110 115 121 127 133 140 147 154 162 169 178 ...
187 196 205 215 226 237 249 261 274 287 301 316 ...
332 348 365 383 402 422 442 464 487 511 536 562 ...
590 619 649 681 715 750 787 825 866 909 953 1000];
case 96
e = [102 105 107 110 113 115 118 121 124 127 130 133 ...
137 140 143 147 150 154 158 162 165 169 174 178 ...
182 187 191 196 200 205 210 215 221 226 232 237 ...
243 249 255 261 267 274 280 287 294 301 309 316 ...
324 332 340 348 357 365 374 383 392 402 412 422 ...
432 442 453 464 475 487 499 511 523 536 549 562 ...
576 590 604 619 634 649 665 681 698 715 732 750 ...
768 787 806 825 845 866 887 909 931 953 976 1000];
case 192
e = [101 102 104 105 106 107 109 110 111 113 114 115 ...
117 118 120 121 123 124 126 127 129 130 132 133 ...
135 137 138 140 142 143 145 147 149 150 152 154 ...
156 158 160 162 164 165 167 169 172 174 176 178 ...
180 182 184 187 189 191 193 196 198 200 203 205 ...
208 210 213 215 218 221 223 226 229 232 234 237 ...
240 243 246 249 252 255 258 261 264 267 271 274 ...
277 280 284 287 291 294 298 301 305 309 312 316 ...
320 324 328 332 336 340 344 348 352 357 361 365 ...
370 374 379 383 388 392 397 402 407 412 417 422 ...
427 432 437 442 448 453 459 464 470 475 481 487 ...
493 499 505 511 517 523 530 536 542 549 556 562 ...
569 576 583 590 597 604 612 619 626 634 642 649 ...
657 665 673 681 690 698 706 715 723 732 741 750 ...
759 768 777 787 796 806 816 825 835 845 856 866 ...
876 887 898 909 920 931 942 953 965 976 988 1000];
otherwise
error('Second argument value is not supported: %.0f',e)
end
R = real(R);
for i = 1:length(R(:))
if isinf(R(i)) || R(i) < 0
R(i) = NaN;
end
m = ceil(log10(R(i))-3);
[~,l] = min(abs(log10(R(i))-log10(e)-m));
R(i) = e(l)*10^m;
end