Signal-Processing/testcode_by_conv.m at main · VibAuth/Signal-Processing · GitHub

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
94
95
96
97
98
99
100
101
102
103
close all;
clearvars;
% Vibration parameter setting
rate = 1500;
vibLength = 1.44;
vibSectionLength = vibLength + 0.3;
coarseInterval = rate * vibSectionLength;

% File parameter setting
path = './../Vib-Data/0405data_widegap/';
signalfile = './chirp.csv';
signal = csvread(signalfile);

%% Dabin data
filename_list = {'db1.csv','db2.csv','db3.csv','db4.csv','db5.csv', ...
    'db7.csv','db8.csv','db9.csv','db10.csv','db11.csv'};

cnt = 1;
for i = 1:length(filename_list)
    filename = filename_list{i};
    [db_x(cnt:cnt+2,:), db_y(cnt:cnt+2,:),db_z(cnt:cnt+2,:), db_xf(cnt:cnt+2,:),db_yf(cnt:cnt+2,:), db_zf(cnt:cnt+2,:)] ...
    = func_signalcut_by_conv(path, filename, signal, coarseInterval);
    cnt = cnt+3;
end

%% Heesu data
filename_list = {'hs1.csv','hs2.csv','hs3.csv','hs5.csv','hs6.csv', ...
    'hs7.csv','hs8.csv','hs9.csv','hs11.csv','hs12.csv'};

cnt = 1;
for i = 1:length(filename_list)
    filename = filename_list{i};
    [hs_x(cnt:cnt+2,:), hs_y(cnt:cnt+2,:),hs_z(cnt:cnt+2,:), hs_xf(cnt:cnt+2,:),hs_yf(cnt:cnt+2,:), hs_zf(cnt:cnt+2,:)] ...
    = func_signalcut_by_conv(path, filename, signal, coarseInterval);
    cnt = cnt+3;
end

%% Jinseon Data
filename_list = {'js1.csv','js2.csv','js3.csv','js4.csv','js5.csv', ...
    'js6.csv','js7.csv','js8.csv','js9.csv','js10.csv'};

cnt = 1;
for i = 1:length(filename_list)
    filename = filename_list{i};
    [js_x(cnt:cnt+2,:), js_y(cnt:cnt+2,:),js_z(cnt:cnt+2,:), js_xf(cnt:cnt+2,:),js_yf(cnt:cnt+2,:), js_zf(cnt:cnt+2,:)] ...
    = func_signalcut_by_conv(path, filename, signal, coarseInterval);
    cnt = cnt+3;
end

% figure('Name','Jinseon Data','NumberTitle','off');
% for cnt = 1:45
%     subplot(5,9,cnt)
%     plot(js_z(cnt,:))
% end

%% Correlation

% gather all data
all_zf(1:30,:) = db_zf;
all_zf(31:60,:) = hs_zf;
all_zf(61:90,:) = js_zf;
figure()
imagesc(corr(all_zf'))
daspect([1 1 1])
xticks(0:30:90);
yticks(0:30:90);
title('zaxis fft corr (140~190)')

all_xf(1:30,:) = db_xf;
all_xf(31:60,:) = hs_xf;
all_xf(61:90,:) = js_xf;
figure()
imagesc(corr(all_xf'))
daspect([1 1 1])
xticks(0:30:90);
yticks(0:30:90);
title('xaxis fft corr (140~190)')

all_yf(1:30,:) = db_yf;
all_yf(31:60,:) = hs_yf;
all_yf(61:90,:) = js_yf;
figure()
imagesc(corr(all_yf'))
daspect([1 1 1])
xticks(0:30:90);
yticks(0:30:90);
title('yaxis fft corr (140~190)')

avg_all_axis = (all_xf + all_yf + all_zf)/3;
figure()
imagesc(corr(avg_all_axis'))
daspect([1 1 1])
xticks(0:30:90);
yticks(0:30:90);
title('fft value avg corr (140~190)')

avg_all_axis_corr = (corr(all_zf') + corr(all_xf') + corr(all_yf'))/3;
figure()
imagesc(avg_all_axis_corr)
daspect([1 1 1])
xticks(0:30:90);
yticks(0:30:90);
title('fft corr value avg corr (140~190)')