function y=compexp(x,comp,a,Fs)
% Compressor/expander
% comp	- compression: 0>comp>-1, expansion: 0<comp<1
% a		    - filter parameter <1
h=filter([(1-a)^2],[1.0000 -2*a a^2],abs(x)); 
h=h/max(h);
h=h.^comp;

y=x.*h;

y=y*max(abs(x))/max(abs(y));

