Frercuencia cardiaca en matlab

REPORTE DE FRECUENCIA CARDAICA
1. [s,fs]=wavread('Normal Split S1.wav')
En esta línea se abre el audio a analizar, o formalmente es la lectura de la bioseñal
2. fs=2756.5;
Aquí se mesuriza lafrecuencia en Hz, para obtener este dato, solo se analiza el archivo de audio para determinar la frecuencia. En sí, el dato que nos arrojaba Matlab es de 11025, pero este valor corresponde a dos ciclosde sístole y diástole. Lo dividi entre 4 para obtener el valor de una sola pulsación, ya sea sístole o diástole de un solo ciclo.
3. max_freq=2;
Este es el rango máximo en Hz deseado para graficar.4. N=length(s);
En esta línea, obtenemos la duración de la señal obtenida en segundos.
5. f=(1:N)*fs/N;
En esta línea, se construye el vector de frecuencia.
6. t=(1:N)/fs
Aquí se construye elvector de tiempo, con respecto a la frecuencia en Hz.
7. plot(t,s)
Aqui se grafica el tiempo de la señal obtenida.
8. X=fft(s);
Aqui hacemos uso de la transformada de Fourier para trabajar con loscomponentes de la frecuencia; de ahí obtenemos el espectro de la señal “signal spectrum” que son los componentes de la frecuencia de una señal.
9. m_plot=round(2/(fs/N));
Con esta linea, buscamos mpara dos Hz; es decir, la frecuencia especifica relacionada a un componente dado, se da por f=m/T, donde m representa el numero harmónico.
10. subplot(2,1,1);
11.plot(f(1:m_plot-1),abs(X(2:m_plot)));
Aquí se grafica la magnitud del espectro de nuestra señal.
12. phase=unwrap(angle(X));
En esta línea se calcula la fase del espectro
13. phase=phase*360/(2*pi);
Aquí descomponemos la fase delespectro
14. subplot(2,1,2);
15. plot(f(1:m_plot-1),phase(2:m_plot));
En esta línea se grafica la fase del espectro.
16. [peak m_peak]=max(abs(X(2:m_plot)));
Aquí se busca m (el numero harmónico)como magnitud máxima del pico.
17. max_freq=f(m_peak)
Se calcula y arroja la frecuencia.
18. max_time=1/max_freq
Calcula y arroja el tiempo máximo.
19. latidos_min=60/max_time
Se calcula y...