samer@32: function S=phasespec(X)
samer@32: % PHASESPEC(X): get phase spectra of columns of x
samer@32: %
samer@32: % phasespec :: [[N,T]->real] -> [[M,T]->complex]] :- M is N/2+1.
samer@32: 
samer@32: m=size(X,1)/2;
samer@32: S=fft(X);
samer@32: S=S(1:m+1,:)./abs(S(1:m+1,:)); % get complex phases of components.
samer@32: % first and last rows should be ones
samer@32: