-
Notifications
You must be signed in to change notification settings - Fork 43
/
stopTorqueControlBalancing.m
93 lines (66 loc) · 2.93 KB
/
stopTorqueControlBalancing.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% RUN THIS SCRIPT TO REMOVE LOCAL PATHS ADDED WHEN RUNNING THE
% CONTROLLER.
%
% In the Simulink model, this script is run every time the user presses
% the 'terminate' button.
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
rmpath('./src/')
% Create a folder for collecting data
if Config.SAVE_WORKSPACE
if (~exist(['experiments',date],'dir'))
mkdir(['experiments',date]);
end
matFileList = dir(['./experiments',date,'/*.mat']);
c = clock;
save(['./experiments',date,'/exp_',num2str(c(4)),'-',num2str(c(5)),'.mat'])
end
% Compare the Yarp time with the Simulink time to catch if the required
% integration time step is respected during the simulation
if Config.CHECK_INTEGRATION_TIME && exist('yarp_time','var')
sim_time = yarp_time.time;
% normalize the yarp time over the first value (at t_sim = 0)
yarp_time0 = yarp_time.signals.values - yarp_time.signals.values(1);
% fast check of yarp time vs. Simulink time
if Config.PLOT_INTEGRATION_TIME
figure
hold on
plot(yarp_time0,0:length(yarp_time0)-1,'o')
plot(sim_time,0:length(sim_time)-1,'o-')
xlabel('Time [s]')
ylabel('Iterations')
grid on
legend('Yarp Time','Simulink Time')
title('Yarp time vs. Simulink time')
end
% number of times the real time step was bigger than twice the
% desired time step value
numOfTimeStepViolations = sum(diff(yarp_time0) > 2*Config.tStep);
if numOfTimeStepViolations > 1 && numOfTimeStepViolations <= 50
warning(['The time step tolerance of ', num2str(Config.tStep), '[s] has been violated at least once.'])
elseif numOfTimeStepViolations > 50
warning(['The time step tolerance of ', num2str(Config.tStep), '[s] has been violated more than 50 times.'])
end
end
% If a joint hits the limits or an encoder spike is detected, print a
% warning message displaying the name of the joint
if Config.EMERGENCY_STOP_WITH_ENCODER_SPIKES && exist('res_check_spikes','var')
if ~isempty(res_check_spikes.signals.values)
checkEachJoint = res_check_spikes.signals.values(end,:);
for k = 1:length(checkEachJoint)
if checkEachJoint(k)
warning(['A spike occurred on joint ', WBTConfigRobot.ControlledJoints{k}])
end
end
end
end
if Config.EMERGENCY_STOP_WITH_JOINTS_LIMITS && exist('res_check_range','var')
if ~isempty(res_check_range.signals.values)
checkEachJoint = res_check_range.signals.values(end,:);
for k = 1:length(checkEachJoint)
if checkEachJoint(k)
warning(['Joint ', WBTConfigRobot.ControlledJoints{k}, ' exited the range'])
end
end
end
end