zoom/pc/data/20091104/makehist.m
jim 9c5f811ae2 use unnormalized histogram output
git-svn-id: https://bucket.mit.edu/svn/nilm/zoom@8125 ddd99763-3ecb-0310-9145-efcb8ce7c51f
2009-11-25 16:27:24 +00:00

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)