% generate delta function with specified timebase
function f = f_delta(N)

f = zeros(1,N);
f(1) = 1;
