
%Program na vykreslovanie CSE signalov
clc; clear all; 
% close all;
fvz=500;

typ_trendu=2


[filename,PATH] = uigetfile('*.mat','Choose signal');
file=[PATH filename];
load(file);
time=linspace(0,length(x)/fvz,length(x));

if size(x,1)==3
%     figure
%     subplot(3,1,1)
%     plot(time,x(1,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('X')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,x(2,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('Y')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,x(3,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('Z')
%     xlabel('time [s]')
%     ylabel('[uV]')
%        
else
%     figure
%     subplot(3,1,1)
%     plot(time,x(1,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('I')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,x(2,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('II')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,x(3,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('III')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,x(4,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('aVR')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,x(5,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('aVL')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,x(6,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('aVF')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,x(7,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('V1')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,x(8,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('V2')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,x(9,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('V3')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,x(10,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('V4')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,x(11,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('V5')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,x(12,:))
%     line(0:1/fvz:length(x)/fvz,0,'Color','r')
%     title('V6')
%     xlabel('time [s]')
%     ylabel('[uV]')
end

for i=1:1:size(x,1)
    xZ(i,:)=x(i,:)-mean(x(i,:));
end
%% posun signálu na nulu
if size(xZ,1)==3
%     figure
%     subplot(3,1,1)
%     plot(time,xZ(1,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('X')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xZ(2,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('Y')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xZ(3,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('Z')
%     xlabel('time [s]')
%     ylabel('[uV]')
%        
else
%     figure
%     subplot(3,1,1)
%     plot(time,xZ(1,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('I')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xZ(2,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('II')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xZ(3,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('III')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,xZ(4,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('aVR')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xZ(5,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('aVL')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xZ(6,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('aVF')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,xZ(7,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('V1')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xZ(8,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('V2')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xZ(9,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('V3')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,xZ(10,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('V4')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xZ(11,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('V5')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xZ(12,:))
%     line(0:1/fvz:length(xZ)/fvz,0,'Color','r')
%     title('V6')
%     xlabel('time [s]')
%     ylabel('[uV]')
 end

for i=1:1:size(x,1)
    xT(i,:) = trend_ECG(xZ(i,:),fvz,typ_trendu);
end
% signmál s trendem
if size(xT,1)==3
%     figure
%     subplot(3,1,1)
%     plot(time,xT(1,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('X')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xT(2,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('Y')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xT(3,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('Z')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%        
else
%     figure
%     subplot(3,1,1)
%     plot(time,xT(1,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('I')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xT(2,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('II')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xT(3,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('III')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,xT(4,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('aVR')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xT(5,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('aVL')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xT(6,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('aVF')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,xT(7,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('V1')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xT(8,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('V2')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xT(9,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('V3')
%     xlabel('time [s]')
%     ylabel('[uV]')
%     
%     figure
%     subplot(3,1,1)
%     plot(time,xT(10,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('V4')
%     ylabel('[uV]')
%     subplot(3,1,2)
%     plot(time,xT(11,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('V5')
%     ylabel('[uV]')
%     subplot(3,1,3)
%     plot(time,xT(12,:))
%     line(0:1/fvz:length(xT)/fvz,0,'Color','r')
%     title('V6')
%     xlabel('time [s]')
%     ylabel('[uV]')
end
