9c5f811ae2
git-svn-id: https://bucket.mit.edu/svn/nilm/zoom@8125 ddd99763-3ecb-0310-9145-efcb8ce7c51f
30 lines
619 B
Matlab
30 lines
619 B
Matlab
# Make a histogram of the error (in milliamps)
|
|
|
|
a=load("log-locked");
|
|
len=size(a,1)
|
|
keith=a(:,3);
|
|
calib=a(1,4)
|
|
dac=a(:,5);
|
|
adc=a(:,6);
|
|
lookup=load("../../../firmware/lookup.inc");
|
|
meas=-(lookup(dac + 1)/64 - calib * adc);
|
|
keith_amps = keith * 50; # 50 turns
|
|
p = polyfit(meas, keith_amps, 1)
|
|
meas_amps = meas * p(1) + p(2); # calibration fit
|
|
err_amps = keith_amps - meas_amps;
|
|
|
|
[n, x] = hist(err_amps*1000, [-15:0.5:15]);
|
|
|
|
for i = 0:30
|
|
cumulative(i+1)=sum(n((31-i):(31+i))/len);
|
|
endfor
|
|
figure(1);
|
|
plot(x(31:61),cumulative);
|
|
|
|
|
|
savedata = [ x' n' ];
|
|
save("-ascii", "hist.dat", "savedata");
|
|
|
|
figure(2);
|
|
bar(x,n)
|