-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathview_signals.m
More file actions
86 lines (75 loc) · 2.58 KB
/
view_signals.m
File metadata and controls
86 lines (75 loc) · 2.58 KB
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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% ECGAUG -- view_signals
% Version 1.0
%
% Views template signals provided with ECGAug software
%
% Copyright 2020 Jonathan W. Waks and Hans F. Stabenau
% Beth Israel Deaconess Medical Center, Boston MA
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program. If not, see <https://www.gnu.org/licenses/>.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
% Folder where signals are present
input_folder = 'C:\ECGAug\templates';
% Load list of files in load directory
file_list_struct = dir(fullfile(input_folder, '*.mat'));
num_files = length(file_list_struct);
file_list = cell(num_files, 1);
for i = 1:num_files
file_list{i} = fullfile(file_list_struct(i).folder, file_list_struct(i).name);
end
orig_directory = file_list_struct(i).folder;
file_list = sort(file_list');
for i=1:length(file_list)
load(file_list{i})
sig{i} = vm_sig.vm;
Q(i) = vm_sig.Q;
S(i) = vm_sig.S;
Tend(i) = vm_sig.Tend;
end
% Number of figures with 5x5 subplots
pgs = ceil(length(file_list)/25)
% Plot
for j=1:pgs
if j<pgs
figure(j)
for i=1+(25*(j-1)):25+(25*(j-1))
subplot(5,5,i-(25*(j-1)))
plot(sig{i})
hold on
scatter(Q(i),sig{i}(Q(i)))
scatter(S(i),sig{i}(S(i)))
scatter(Tend(i),sig{i}(Tend(i)))
line([0 length(sig{i})],[0 0])
[~,fname,~] = fileparts(file_list{i});
title(fname)
end
else
remainder = length(file_list)-(25*(pgs-1))
figure(j)
for i=1+(25*(j-1)):remainder+(25*(j-1))
subplot(5,5,i-(25*(j-1)))
plot(sig{i})
hold on
scatter(Q(i),sig{i}(Q(i)))
scatter(S(i),sig{i}(S(i)))
scatter(Tend(i),sig{i}(Tend(i)))
line([0 length(sig{i})],[0 0])
[~,fname,~] = fileparts(file_list{i});
title(fname)
end
end
end