t=(0:1/1000:1);
f=4;
A=1;
x=A*sin(2*pi*f*t);
B=1; %ilosc bitow kwantowania
MN=2^(B-1);
xq=((ceil(0.999999999*x/max(x)*MN)/MN)-(0.5^B))*max(x);
%sygnał jest przybliżany do środka przedziału
blad=x-xq;
dokladnosc_przyblizenia=xq/x
plot(t,x,t,xq,t,blad)