function yq = quantize(y,min,max,N)

% function yq = quantize(y,min,max,N) performs quantization of vector y to
% N bits over the range "min" to "max."  Quantization is done by rounding.

range = max-min;    % Find range of signal

scale = 2^N/range;  % Scale factor for quantization

yq = round(y*scale);    % Scale up, then round.

yq = yq/scale;  % Scale back down to original.

return;