% -------------------------------------------------------------------------------------------------
% TSmerge.m:
%
% Copyright (C) 2004 Federico Cantini, Clinical Physiology Institute, CNR, Pisa, Italy
%
% 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 2 of the License, or (at your option) any later version.
%
% This program is distributed "as is" and "as available" in the hope that it will be useful,
% but WITHOUT ANY WARRANTY of any kind; 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, write to the Free Software Foundation, Inc.,
% 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
%
% For any comment or bug report, please send e-mail to: federico.cantini@ifc.cnr.it
% -------------------------------------------------------------------------------------------------
function out = TSmerge (main_dir, CasiDir, record_list)
if nargin < 1
record_learn ={'t01','t02','t03','t04','t05','t06','t07','t08','t09','t10','s01','s02','s03','s04','s05','s06','s07','s08','s09','s10'};
record_list = record_learn;
main_dir = '.\';
end
for ii = 1:length(record_list)
clear beatGroup;
file = fullfile(main_dir, 'DatiQc',CasiDir,[record_list{ii},'_QRSclass2.mat']);
load (file);
max_jj = 1;
max_num = length(beatGroup.vq{1});
for jj = 2:length (beatGroup.vq)
if (length(beatGroup.vq{1}) > max_num)
max_num = length(beatGroup.vq{1});
max_jj = jj;
end
end
QRSclass{ii} = beatGroup.avgC1Q(:,max_jj);
end
for ii = 1:length(QRSclass)/2
for jj = 1:length(QRSclass)
if (ii ~= jj)
d(jj) = dist_p(QRSclass{ii},QRSclass{jj});
else
d(jj) = 99999999;
end
end
[pat(ii).value,pat(ii).ind] = min(d);
end
for ii = 1:length(QRSclass)/2
out {ii,1} = record_list{ii};
out {ii,2} = record_list{pat(ii).ind};
fprintf('%s <--> %s\n',record_list{ii},record_list{pat(ii).ind});
end