function DrawECG
% DrawECG.m
% Author: Alexander Khaustov; alexander dot khaustov at gmail dot com
% Copyright (C) 2008 St.-Petersburg Institute of Cardiological Technics (Incart), www.incart.ru
% This software is released under the terms of the GNU General
% Public License (http://www.gnu.org/copyleft/gpl.html).
%
% Draws the graphs associated with the ECG data (subplot 1), and
% the corresponding landmarks; q, s, Tend and fiducial points
global graph
subplot(graph(1));
global ecg q s tend Align lead Param
SaveAxisLimits;
plot(ecg(:, lead)); hold on;
amp = zeros(1,length(q));
for i=1:length(q)
amp(i) = (max(ecg(q(i):s(i))) - min(ecg(q(i):s(i)))) / 5;
plot([s(i) s(i)], [ecg(s(i), lead) - amp(i), ecg(s(i), lead) + amp(i)], 'k');
plot([q(i) q(i)], [ecg(q(i), lead) - amp(i) / 2, ecg(q(i), lead) + amp(i) / 2], 'k');
end;
if ~isempty(Align)
% for i=1:min([length(Align.fidBase), length(amp)])
% plot([Align.fidBase(i) Align.fidBase(i)], [ecg(Align.fidBase(i), lead) - amp(i), ecg(Align.fidBase(i), lead) + amp(i)], 'g');
% end;
for i=1:min([length(amp), length(Align.fid)])
if (Align.fid(i))
plot([Align.fid(i) Align.fid(i)], [ecg(Align.fid(i), lead) - amp(i), ecg(Align.fid(i), lead) + amp(i)], 'm');
end;
end;
end;
for i=1:min(length(tend), length(amp))
plot([tend(i) tend(i)], [ecg(tend(i), lead) - amp(i), ecg(tend(i), lead) + amp(i)], 'r');
end;
if ~isempty(Param)
if (strcmp(Param.Metric, 'MMA'))
global CurResIndex q s TWARes
qrs = TWARes.res(CurResIndex).StartQRSInd;
x = (q(qrs) + s(qrs)) / 2;
else
global StartQRSInd QRSRelInd s
x = s(StartQRSInd + QRSRelInd);
end;
global hecg;
if (ishandle(hecg))
delete(hecg);
end;
hecg = plot([x x], get(gca, 'YLim'), 'k');
end;
ylabel('amplitude (mCV)')
hold off;
RestoreAxisLimits;
return;