S-transform

S变换源码nhaf=fix(N/2); odvn=1; if nhaf*2==N; odvn=0; end f=[0:nhaf -nhaf+1-odvn:-1]/N; Hft=fft(h); %Compute all frequency domain Gaussians as one matrix invfk=[1./f(2:nhaf+1)]'; W=2*pi*repmat(f,nhaf,1).*repmat(invfk,1,N);