www.gusucode.com > wavelet 源码程序 matlab案例代码 > wavelet/WaveletTimeFreqRidgeExtractionDefaultExample.m

    %% Extract Time-Frequency Ridge from Chirp Signal
% Obtain the wavelet synchrosqueezed transform of a quadratic chirp and 
% extract the maximum time-frequency ridge, in |fridge|, and the associated
% row indices, in |iridge|. 
%%
% Load the chirp signal and obtain its synchrosqueezed transform.
load quadchirp;
[sst,f] = wsst(quadchirp);
%%
% Extract the maximum time-frequency ridge.
[fridge,iridge] = wsstridge(sst);
%%
% Plot the synchrosqueezed transform.
pcolor(tquad,f,abs(sst))
shading interp
title('Synchrosqueezed Transform')
%%
% Overlay the plot of the maximum energy frequency ridge.
hold on
plot(tquad,fridge)
title('Synchrosqueezed Transform with Overlaid Ridge')