close all;
clear all;
a=4;
f0=0.02;
L=50;
n=1:(L-1);
x=a.*cos(2.*pi.*f0.*n);
bit=5;
range=2^5;
plot(x);
interval=(a*2)/(range-1);
partition = [-a:interval:a];
interval2 = (a*2)/(range)
codebook=[-a:interval2:a];
[index quantized]= quantiz(x,partition,codebook);
plot(n,x,'x',n,quantized,'.');