1. function [WD] = EVD_tfrwv(x0)
 2. X=fft(x);
 3. X=[X(1:N/2+1);zeros(N,1);X(N/2+2:N+1)];
 4. x=2*ifft(X);
 5. x=[zeros(N,1);x;zeros(N,1)];
 6. X=zeros(2*N+1);
 7. for k=1:2*N+1  $ $X(:,k)=x(k+(0:2*N)); $ $end
 8. ww=X.*conj(flipud(X));
 9. WW=fft(ww([N+1:2*N+1,1:N],:));
 10. WW=real(WW([N+2:2*N+1,1:N+1],:));