function hh = dtmfdesign(fb, L, fs)

for k = 1 : length(fb)
    hh(k, :) = 2*cos(2*pi*fb(k)*(0:L-1)/fs)/L;
end

hh = hh';


