Reassigned smoothed pseudo Wigner-Ville distribution.
[TFR,RTFR,HAT] = tfrrspwv(X) [TFR,RTFR,HAT] = tfrrspwv(X,T) [TFR,RTFR,HAT] = tfrrspwv(X,T,N) [TFR,RTFR,HAT] = tfrrspwv(X,T,N,G) [TFR,RTFR,HAT] = tfrrspwv(X,T,N,G,H) [TFR,RTFR,HAT] = tfrrspwv(X,T,N,G,H,TRACE) [TFR,RTFR,HAT] = tfrrspwv(...,'plot')
A Nx elements vector or a Nx by 2 array signal.
the time instant(s) with elements in [1 Nx] (default : 1:length(X))
number of frequency bins (default : length(X)). For faster computation N should be a power of 2.
a real vector of odd length: the time smoothing window (default : Hamming(N/10)).
a real vector of odd length: the frequency smoothing window, (default : Hamming(N/4)).
A boolean (or a real scalar) if true (or nonzero),the progression of the algorithm is shown (default : %f).
if one input parameter is 'plot', tfrqview is called and the time-frequency representation will be plotted.
a real N by Nt array: the time-frequency representation.
A real N by Nt array: the reassigned time-frequency representation.
a N by Nt complex matrix: the reassignment vectors.
tfrrspwv computes the smoothed pseudo Wigner-Ville distribution and its reassigned version.
N=128; sig=fmlin(N,0.05,0.15)+fmlin(N,0.3,0.4); t=1:2:N; g=window("kr",15,3*%pi); h=window("kr",63,3*%pi); [tfr,rtfr,hat]=tfrrspwv(sig,t,64,g,h); clf;gcf().color_map= jetcolormap(128); subplot(121) grayplot(t,linspace(0,0.25,N/2),tfr(1:N/2,:)') title("TFR") subplot(122) grayplot(t,linspace(0,0.25,N/2),rtfr(1:N/2,:)') title("RTFR") | ![]() | ![]() |