function [segmentos tfins pulsacaomedia]=trocos(sinal,batidas,inicio,freqamostragem,NUMbatidas)
#trocos.m returns segment of 'NUMbatidas' beats, end of T-wave on each beat and average heart rithm
#Copyright (C) 2008 Rui Rodrigues <rapr@fct.unl.pt>
#This software is released under the terms of the GNU General Public License (http://www.gnu.org/copyleft/gpl.html)
%sinal é o ecg a considerar mas só a coluna de um dos 'leads'
%batidas é o ficheiro por exemplo da forma:
% 'ihr -r edb/e0105 -a atr -i -v>batidas_edb_e0105.txt'
%inicio é o numero da batida onde se vão iniciar a 128 batidas a considerar
%determinar pulsacao media
soma=sum(batidas(inicio:inicio+NUMbatidas-1,2));
pulsacaomedia=soma/NUMbatidas;
%dados do artigo ??
pulsacaoref=120;
aux=floor(225/(1000/freqamostragem)*(pulsacaoref/pulsacaomedia));
%este valor é demasiado grande (será que estou a fazer algo errado?)
numamostraspsegmento=2*ceil(0.6*aux/2);
%vou substituir por
%numamostraspsegmento=2*ceil(0.4*aux/2)
%delay=floor(55/(1000/freqamostragem)*(pulsacaoref/pulsacaomedia));
%a partir daqui esta diferente do programa inicial(artigo)
%vai ser usado o outro artigo para localizar o fim da onda T
segmentos=zeros(numamostraspsegmento,NUMbatidas);
tfins=zeros(1,NUMbatidas);
for i=1:NUMbatidas
fim=tfim(sinal(batidas(inicio+i-1,1):batidas(inicio+i,1)),freqamostragem,batidas(inicio+i-1,2));
tfins(1,i)=batidas(inicio+i-1,1)+fim;
comeco=batidas(inicio+i-1,1)+fim-numamostraspsegmento+1;
segmentos(:,i)=sinal(comeco:comeco+numamostraspsegmento-1);
endfor
endfunction
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%funcao para localizar, numa batida, o fim da onda T
%parte daquilo que esta feito em 'tlimites.m'
function Tfim=tfim(sinal_batida,freqamostragem,batida)
w=ceil(32*freqamostragem/250);
l=length(sinal_batida);
%inicio da procura do tfim
if (batida>68)
inicio=ceil(0.15*l+freqamostragem/250*37);
else
inicio=ceil(freqamostragem/250*70);
endif
%fim da procura do tfim
if (batida>68)
fim=ceil(0.7*l-freqamostragem/250*9);
else
fim=ceil(freqamostragem/250*(0.2*l+101));
endif
Tfim=procuratfins(sinal_batida,inicio,fim,w);
endfunction
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%copia da versao de 'tlimites'
function res=procuratfins(sinal_batida,inicio,fim,w)
%threshold
L=6;
A=zeros(fim-inicio+1,1);
for j=inicio:fim
A(j-inicio+1)=sum(sinal_batida(j-w+1:j))-w*sinal_batida(j);
endfor
[m1 res1]=max(A);
[m2 res2]=min(A);
if ((abs(m1/m2)>1/L)&&(abs(m1/m2)<L))
res=max([res1,res2])+inicio-1;
else
[a res]=max([abs(m1),abs(m2)]);
res=res+inicio-1;
endif
endfunction
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%