% Segnali e Sistemi a.a. 2021-22
% Es1 Lab 6, ECG
clear all
clc
close all

%%%%%%%%%%%%%%%% CARICA I DATI 

load data_ecg.mat %carica i dati ecg ed ecg_dist
fc = 125; %Hz
wc = 2*pi*fc; %rad/s 
Tc = 1/fc; %s
%sottrazione della componente continua
s = ecg-mean(ecg);
N = length(s);
t = [0:N-1]*Tc; %asse dei tempi

figure(1)
plot(t,s,'b-')
title('ECG')
xlabel('t(s)')
ylabel('microV')
grid on
axis([0 t(end) ylim])
% zoom
figure(2)
plot(t,s,'b-')
title('ECG')
xlabel('t(s)')
ylabel('microV')
grid on
axis([0 5 ylim])


%%%%%%%%%%%%%%%% TRASFORMATA DI FOURIER 

% calcolo della trasformata di Fourier

% plot del modulo della trasformata

% zoom in scala lineare

% identificazione del primo massimo: funzione find