function y = cfilt(x, lag);
% y=cfilt(x,lag) - cancellation comb filter - lag is in samples

lag = fix(lag);
n = size(x,2);
y = x(:,lag+1:n) -  x(:,1:n-lag);