Compute wsst, and plot its magnitude:
load mtlb
dt = 1/Fs;
t = 0:dt:numel(mtlb)*dt-dt;
[sst,f] = wsst(mtlb,Fs);
%disp([size(f),size(f(f>2000 & f<=2500)),size(sst)])
pcolor(t,f,abs(sst))
shading interp; colorbar
xlabel("Seconds")
ylabel("Frequency (Hz)")
title("Synchrosqueezed Transform")
Sum up the wavelet amplitudes in band from f=2000 to f=2500:
S1=sum(abs(sst((f>2000 & f<=2500),:)),'all');
fprintf('Sum of wsst amplitudes, from f=2000 to f=2500: %.1f.\n',S1)
S1=58 .1 the sum of about 40,000 array elements. Therefore the mean value of the elements is about 58/4e4~=0.0015. This seems like it could be consistent with the amplitudes, as indicated by the colorbar. At least it is not wildly inconsistent.
