Difference between revisions of "Ltf matlab source"
From DIDEAS Wiki
m |
m |
||
Line 31: | Line 31: | ||
ylabel('Ak motor rotations') | ylabel('Ak motor rotations') | ||
</pre> | </pre> | ||
+ | =test2= |
Revision as of 15:32, 17 December 2009
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')