-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathexeSG.m
34 lines (32 loc) · 2.21 KB
/
exeSG.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
% Execute Savitzky-Golay Filter on All Data
% polynominal degree
p = 3;
% window size based on half eye blink lenght (200ms) (still dunno why, that value give better result T.T)
% actual value is 102 data point, but SG window should be odd number, so +1
l = 103;
% Do Savitzky-Golay filter on Pre Condition
sgPre_Adrian = sgolayfilt(Pre_Adrian, p, l); Pre_Adrian = Pre_Adrian - sgPre_Adrian;
sgPre_Aji = sgolayfilt(Pre_Aji, p, l); Pre_Aji = Pre_Aji - sgPre_Aji;
sgPre_Aris = sgolayfilt(Pre_Aris, p, l); Pre_Aris = Pre_Aris - sgPre_Aris;
% sgPre_Bella = sgolayfilt(Pre_Bella, p, l); Pre_Bella = Pre_Bella - sgPre_Bella;
% sgPre_Gafur = sgolayfilt(Pre_Gafur, p, l); Pre_Gafur = Pre_Gafur - sgPre_Gafur;
sgPre_Ima = sgolayfilt(Pre_Ima, p, l); Pre_Ima = Pre_Ima - sgPre_Ima;
sgPre_Jodi = sgolayfilt(Pre_Jodi, p, l); Pre_Jodi = Pre_Jodi - sgPre_Jodi;
% sgPre_Mas_Fadil = sgolayfilt(Pre_Mas_Fadil, p, l); Pre_Mas_Fadil = Pre_Mas_Fadil - sgPre_Mas_Fadil;
% sgPre_Rizky = sgolayfilt(Pre_Rizky, p, l); Pre_Rizky = Pre_Rizky - sgPre_Rizky;
sgPre_Sely = sgolayfilt(Pre_Sely, p, l); Pre_Sely = Pre_Sely - sgPre_Sely;
% sgPre_Supri = sgolayfilt(Pre_Supri, p, l); Pre_Supri = Pre_Supri - sgPre_Supri;
sgPre_Wega = sgolayfilt(Pre_Wega, p, l); Pre_Wega = Pre_Wega - sgPre_Wega;
% Do Savitzky-Golay filter on Pos Condition
sgPos_Adrian = sgolayfilt(Pos_Adrian, p, l); Pos_Adrian = Pos_Adrian - sgPos_Adrian;
sgPos_Aji = sgolayfilt(Pos_Aji, p, l); Pos_Aji = Pos_Aji - sgPos_Aji;
sgPos_Aris = sgolayfilt(Pos_Aris, p, l); Pos_Aris = Pos_Aris - sgPos_Aris;
% sgPos_Bella = sgolayfilt(Pos_Bella, p, l); Pos_Bella = Pos_Bella - sgPos_Bella;
% sgPos_Gafur = sgolayfilt(Pos_Gafur, p, l); Pos_Gafur = Pos_Gafur - sgPos_Gafur;
sgPos_Ima = sgolayfilt(Pos_Ima, p, l); Pos_Ima = Pos_Ima - sgPos_Ima;
sgPos_Jodi = sgolayfilt(Pos_Jodi, p, l); Pos_Jodi = Pos_Jodi - sgPos_Jodi;
% sgPos_Mas_Fadil = sgolayfilt(Pos_Mas_Fadil, p, l); Pos_Mas_Fadil = Pos_Mas_Fadil - sgPos_Mas_Fadil;
% sgPos_Rizky = sgolayfilt(Pos_Rizky, p, l); Pos_Rizky = Pos_Rizky - sgPos_Rizky;
sgPos_Sely = sgolayfilt(Pos_Sely, p, l); Pos_Sely = Pos_Sely - sgPos_Sely;
% sgPos_Supri = sgolayfilt(Pos_Supri, p, l); Pos_Supri = Pos_Supri - sgPos_Supri;
sgPos_Wega = sgolayfilt(Pos_Wega, p, l); Pos_Wega = Pos_Wega - sgPos_Wega;