-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathvis.m
79 lines (70 loc) · 1.92 KB
/
vis.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
%% EMG filtered visualization
close all
subjects = [4];
file='T1';
eval(['T = ',file,';']);
for i_sub = subjects
for i_trial=1:length(T(i_sub).EMG.Right)
figure
for i=1:16
subplot(4,4,i)
plot(T(i_sub).EMG.Right(i_trial).M(i,:))
end
end
end
%% KIN visualization
close all
subjects = [4];
file='T1';
eval(['T = ',file,';']);
for i_sub = subjects
for i_trial=1:length(T(i_sub).KIN.Right)
figure
for i=1:9
subplot(3,3,i)
plot(T(i_sub).KIN.Right(i_trial).M(i,:))
end
end
end
%% COM visulaization
subjects = [3];
file='T1';
eval(['T = ',file,';']);
figure
for i_sub = subjects
numTrials = length(T(i_sub).KIN.Right);
for i_trial=1:numTrials
for i=1:3
sig = T(i_sub).ForcePlate(i_trial).COM(:,i);
subplot(numTrials,4,4*(i_trial-1)+i)
plot(sig)
rmsVal = rms(sig-mean(sig));
stdVal = std(sig);
title(['rms: ',num2str(rmsVal),' /std: ',num2str(stdVal)])
end
subplot(numTrials,4,4*(i_trial-1)+4)
plot(T(i_sub).ForcePlate(i_trial).COM(:,1), T(i_sub).ForcePlate(i_trial).COM(:,2))
axis equal
end
end
%% COP visulaization
subjects = [3];
file='T1';
eval(['T = ',file,';']);
figure
for i_sub = subjects
numTrials = length(T(i_sub).KIN.Right);
for i_trial=1:numTrials
for i=1:2
sig = T(i_sub).ForcePlate(i_trial).COP(:,i);
subplot(numTrials,3,3*(i_trial-1)+i)
plot(sig)
rmsVal = rms(sig-mean(sig));
stdVal = std(sig);
title(['rms: ',num2str(rmsVal),' /std: ',num2str(stdVal)])
end
subplot(numTrials,3,3*(i_trial-1)+3)
plot(T(i_sub).ForcePlate(i_trial).COP(:,1), T(i_sub).ForcePlate(i_trial).COP(:,2))
axis equal
end
end