function [z]= resam(x, r1,r2)
%x-wektor
%r-współczynnik

%%%%NADPRÓBKOWANIE%%%%
dl = length(x); %sprawdzamy długość wektora wejściowego x
wsp1 = r1-1;
w = 0;
w(1) = x(1);
for i=2:dl
    gz = [zeros(1,wsp1),x(i)];
w = [w,gz];
end

%%%%FILTR%%%%
R=50; %dB tłumienie  w paśmie zaporowym
fs=8000; % częstotliwość próbkowania
fp=1250; % pasmo przejściowe
fr=1750; %        -||-
f0=fp+(fr-fp)/2; % częstotliwość odcięcia
N=ceil(3.3/((fr-fp)/fs)); %rząd filtru
n=0:N;

hx = 2*f0/fs*sinc(2*f0/fs*(n-N/2));
% h = hx.*hamming(N+1)';
w=0.54-0.46*cos((2*pi*n)/N); % okno Hamminga
h = hx.*w;

xfiltrowane = conv(h,w,'same');

%%%%NADPRÓBKOWANIE%%%%
z=xfiltrowane(1);
dlf=length(xfiltrowane);
for i=r2+1:r2:dlf
    y=xfiltrowane(i);
    z=[z,y];

end

end