How to do this figure in MatLab?
Partial code:
M = 256; sigma = 6.0;h = sin(+2*pi*sigma*log((1-([0:M-1]/M))-0.2));h(1:37) = 0; h(193:M) = 0;H= [fft(h) fft(h) fft(h)];k = ([0:floor(M/2)-1 -floor(M/2):-1]);S(1,:)=sum(ifft(H(1:M)))/M*ones(1,M);for n=1:floor(M/2)-1; W=(1./M)*exp(-n^2*k.^2/(2*M^2))... .*exp(-2*pi*i*(n*k/M-sign(n) ... *sigma*log(sigma+abs(n)*k/M))); W(M/2+1:M-ceil(sigma*M/abs(n))+1)=0; W = W/sum(W); S(n+1,:) = ifft(H(n+M+1:n+M+M).* fft(W));end figure(10), contour(abs(S),'ShowText','off'); colorbar
Source code and picture, article DOI: 10.1016/j.sigpro.2004.03.015
Title: "Time-local Fourier analysis with a scalable, phase-modulated analyzing function: the S-transform with a complex window"
Best Answer