mulaw_e(yq,B,mu)

function [sq] = mulaw_e(yq,B,mu)

% mu-law expander

% ---------------

% [sq] = mulaw_e(yq,B,mu)

% sq = expanded quantized signal

% yq = quantized input in integers: 0 <=yq <= 2^B-1

%  B = bits/sample

% mu = compander parameter

%

swing = (2^B-1)/2;

yq = ((yq-swing)/swing);

if mu == 0

  sq = yq;

else

  sq = (((1+mu).^(abs(yq))-1)/mu).*sign(yq);

end