clear;
clc;

s1 = load("ecg/A00001.mat");
s2 = load("ecg/A00004.mat");

Y = s2.val;
X = 1:width(Y);
DIV = 300; % Pierwotny pomiar był wykonany w cz. 300 Hz

ERR_MAX = zeros(1, 15);
ERR_MEAN = ERR_MAX;
ERR_MEDIAN = ERR_MAX;

prob_freq = 4; % częstotliwość 300 / prob_freq

X1 = 1:prob_freq:width(X);
Y1 = Y(X1);
close all;
hold on;
plot(X, Y);
plot(X1, Y1, 'or', 'LineWidth', 1.5);
[~, Yr] = sinc_inter(Y1, prob_freq); % Xr powinien być równy X.
plot(X,Yr, '--b', 'LineWidth', 1.5);
[~, Ye] = e_inter(Y1, prob_freq);
plot(X,Ye, ':m', 'LineWidth', 1.5);
grid on;
xlim([1000 2000]);
plot([10 10], [-1.5 1.5], '--k', 'LineWidth', 2);
plot([70 70], [-1.5 1.5], '--k', 'LineWidth', 2);
xlabel("t");
ylabel("A");
legend('Sygnał wejściowy', 'Punkty pomiaru', 'sinc-inter', 'e-inter');