%% Dirac frequency generator

function [freq,phase] = Dirac(fmin, fmax,fstep)

freq = linspace(fmin,fmax,(fmax-fmin)/fstep);
phase= linespace(0,0,(fmax-fmin)/fstep);