wdf.zip_WDF_WDF matlab_Wigner_wigner function_zip

function y = wdf(x, t, f) T = length(t); F = length(f); dt = t(2)-t(1); df = f(2)-f(1); N = round(1/(df*dt*2)); n0 = t(1)/dt; n1 = t(T)/dt; m0 = f(1)/df; m1 = f(F)/df; m = f/df; M = round(mod(m,N)+1); y = zeros(F,T); for n = n0:n1 Q = round(min(n-n0,n1-n)); c1= zeros(1,N); for q = 0:(2*Q) c1(q+1)=x((n+q-Q)-n0+1)*conj(x((n-q+Q)-n0+1)); end C = fft(c1,N); y(:,n-n0+1)=2*dt*exp(j*2*pi*m*Q/N).*C(M); end end