06efd7b5c2
git-svn-id: https://bucket.mit.edu/svn/nilm/zoom@8679 ddd99763-3ecb-0310-9145-efcb8ce7c51f
34 lines
891 B
Matlab
34 lines
891 B
Matlab
a = load("dactest.log");
|
|
|
|
for i = 0:0x3ff;
|
|
amin(i+1) = min(a(find(a(:,2) == i),3));
|
|
amax(i+1) = max(a(find(a(:,2) == i),3));
|
|
amean(i+1) = mean(a(find(a(:,2) == i),3));
|
|
astd(i+1) = std(a(find(a(:,2) == i),3));
|
|
end
|
|
|
|
# First line is the maximum deviation that we saw
|
|
# (ie. 14 bits means "1 part in 2^14, over the full 10V scale")
|
|
# Second line is 3 standard deviations, so 99.9% for a normal distribution.
|
|
|
|
t = 1:1024;
|
|
|
|
figure(1);
|
|
plot(t, log2(10./(amax(t)-amin(t))), ';Max;', t, log2(10./(3*astd(t))), ';99.9%;');
|
|
xlabel("DAC command")
|
|
ylabel("Deviation from mean (bits)");
|
|
|
|
figure(2);
|
|
p = polyfit(t, amean(t), 1)
|
|
plot(t, log2(10 ./ ((t * p(1) + p(2)) - amean(t))));
|
|
xlabel("DAC command")
|
|
ylabel("Average error from linear (bits)");
|
|
|
|
if input('Print? [y/N] ', 's') == 'y'
|
|
disp('Printing');
|
|
figure(1);
|
|
print('-landscape');
|
|
figure(2);
|
|
print('-landscape');
|
|
end
|