function [X1,Y1] = sinc_inter(Y, q)
% SINC_INTER - interpolacja sygnału funkcją sinc
% Wejście:
%   Y - wektor próbek
%   q - skala (np. dla skali 2 otrzymamy 2x więcej próbek
% Wyjście:
%   X1 - czas dla 'nowych próbek',
%   Y1 - interpolacja.

p = width(Y);
X = 1:p;

X1 = voffrept(X, q - 1);
Y1 = zeros(size(X1));
for i=1:1:width(X)
    Y1 = Y1 + Y(i) * sinc(X1 - i); % ???
end
X1 = (X1 - 1) * q + 1;
end