cz_prob = 8000;
czestotliwosc = 3800;

ts = 0:1/20000:(2*pi*(1/czestotliwosc));
s = sin(czestotliwosc*ts);
t = 0:1/cz_prob:(2*pi*(1/czestotliwosc));
x = sin(czestotliwosc*t);

%kwantyzacja
kw = quantizer([8 1], 'float');
X1 = quantize(kw, x);
X2 = compand(X1,87.6,max(X1),'a/compressor');


%Sygnał zkwantowany i po kodowaniu
kwant_kod = num2bin(kw,X2);

%Sygnał zkwantowany bez kodowania
kwant_bez = num2bin(kw,X1);
%bez kompresji
zamiana_bez=bin2num(kw,kwant_bez);


%zamiana z binarnych na analog
zamiana_z = bin2num(kw,kwant_kod);
%kompresja/ekspander
Y = compand(zamiana_z,87.6,max(zamiana_z),'a/expander');

%sygnał po kodeku
y = quantize(kw, Y);

%sygnał bez kodeku
z= quantize(kw, zamiana_bez);


figure(1)
stem(t,z)
hold on
plot(ts,s,'black')
title('Wykres sygnału bez użycia kodeku')
figure(2)
stem(t,y)
hold on
plot(ts,s,'green')
title('Wykres sygnału z użyciem kodeku')