Ltf matlab source

From DIDEAS Wiki
Revision as of 15:32, 17 December 2009 by Ceb (talk | contribs)
Jump to: navigation, search

PF Users Navigation:

Edit

potential function

fname = 'static_test_experiment_log-20091216_1.txt';
s = importdata(fname);

% compute the JR3 DC offset from the first 2 log entries which are +-0.2 amps
dc = (s.data(1,:) + s.data(2,:) )/2;

ndx = 1:length(s.data);   % use full data set
bd_pos = s.data(ndx,1);
current = s.data(ndx,2);
qiq = s.data(ndx,3);
ank_pos = s.data(ndx,4);
lw = s.data(ndx,5);
jr3z = s.data(ndx,8) - dc(8);  % apply DC correction, grab the Fz component
jr3x = s.data(ndx,6) - dc(6);  
jr3y = s.data(ndx,7) - dc(7);  


ndx = find( bd_pos == -11000)
subplot(2,1,1);
plot(jr3z(ndx), ank_pos(ndx)/2048,'.'); grid on
title('Ak motor rotations, vs JR3 force at Bw=-11000, Imot=-9 : 9 amps')
ylabel('Ak motor rotations')
xlabel('JR3z force N')

subplot(2,1,2);
plot(lw(ndx), ank_pos(ndx)/2048,'.'); grid on
title('Ak motor rotations, vs LW force at Bw=-11000, Imot=-9 : 9 amps')
xlabel('LW force N (over 12 hrs)')
ylabel('Ak motor rotations')

test2