clear;
clc;

MET1600_ERR = zeros(1, 15);
MET1450_ERR = zeros(1, 15);
MET2250_ERR = zeros(1, 15);
METZ250_ERR = zeros(1, 15); % Met. 2*
MET2150_ERR = zeros(1, 15);

for i=1:15
    pack = load("ecg/A" + sprintf("%05d", i) + ".mat");
    Y = pack.val;

    [~, MET1600_ERR(i)] = spike_count(Y, 600);
    [~, MET1450_ERR(i)] = spike_count(Y, 450);
    
    Y100 = Y(1:3:end);
    MET2250_ERR(i) = diff_count(Y100, 250);
    
    Y50 = Y100(1:2:end);
    METZ250_ERR(i) = diff_count(Y50, 250);
    
    MET2150_ERR(i) = diff_count(Y100, 150);
end