31 lines
619 B
Mathematica
31 lines
619 B
Mathematica
|
#!/usr/bin/octave
|
||
|
|
# load data
|
||
|
|
if (length(who("calibdata")) == 0)
|
||
|
|
calibdata = load("calibrate-50a.log");
|
||
|
|
endif
|
||
|
|
|
||
|
|
# split it into each keithley value
|
||
|
|
if (length(who("a")) == 0)
|
||
|
|
for i = 1 : (length(calibdata) / 4000)
|
||
|
|
block=1:4000;
|
||
|
|
a(i, block, 1) = calibdata((i-1) * 4000 + block, 2);
|
||
|
|
a(i, block, 2) = calibdata((i-1) * 4000 + block, 3);
|
||
|
|
a(i, block, 3) = calibdata((i-1) * 4000 + block, 1);
|
||
|
|
endfor
|
||
|
|
endif
|
||
|
|
|
||
|
|
# plot them all
|
||
|
|
#clf;
|
||
|
|
#figure(1);
|
||
|
|
#hold on
|
||
|
|
|
||
|
|
out = [];
|
||
|
|
for i = 1:2:21
|
||
|
|
range = 1:4000;
|
||
|
|
x = a(i, range, 1);
|
||
|
|
y = a(i, range, 2);
|
||
|
|
out = [ out x' y' ];
|
||
|
|
endfor
|
||
|
|
|
||
|
|
save("-ascii","calib.dat","out");
|