function [y] = mulaw_c(s,mu)% mu-law compressor% -----------------% [y] = mulaw_c(s,mu)%  y = compressed output signal%  s = zero-mean nomalized signal between -1 and 1% mu = parameter mu%if mu == 0    y = s;else   y = (log(1+mu*abs(s))/(log(1+mu))).*sign(s);end