% generate unit step function
function f = f_delta(N)

f = zeros(1,N);
f(N/2:N) = 1;
