load wyniki;
N=lenght(x);
n=(0:N-1)';
%modf=abs(f)/(N/2); %bez skl
%modf=abs(f)/N; %ze skl
plot(n, 20*log10(modf), '.');
axis([0 N-1 -100 30]);
%f=(x prazek/N)*fs