function [y, b, a] = fbcomb(x, g, d)
max_length = max(d);

a = [1, zeros(1, (max_length-1)), g];
b = [zeros(1, max_length), 1];

imp_resp = impz(b, a);

%y = conv(x, imp_resp);

for i=1:2
    y(:,i) = conv(x(:,i),imp_resp);
end

