Pseudo Wigner-Ville time-frequency distribution.
[TFR,T,F]=tfrpwv(X) [TFR,T,F]=tfrpwv(X,T) [TFR,T,F]=tfrpwv(X,T,N) [TFR,T,F]=tfrpwv(X,T,N,H) [TFR,T,F]=tfrpwv(X,T,N,H,TRACE) [TFR,T,F]=tfrpwv(...,'plot')
A Nx elements vector (auto-PWV) or a Nx by 2 array signal (cross-PWV).
a real Nt vector with elements in [1 Nx] : time instant(s) (default: 1:NX).
a positive integer: the number of frequency bins (default:Nx). For faster computation N should be a power of 2.
real vector with odd length: the frequency smoothing window,(default: window("hm",N/4)).
It will be normalized such as the middle point equals 1 to preserve signal energy.
if nonzero,the progression of the algorithm is shown (default : 0).
if the last input parameter value is 'plot', tfrqview is called and the time-frequency representation will be plotted.
A real N by Nt array: the time-frequency representation.
A N vector of normalized frequencies.
tfrpwv computes the Pseudo Wigner-Ville distribution of a discrete-time signal X, or the cross Pseudo Wigner-Ville representation between two signals.
Interactive use
Non interactive use
N=128; sig=fmlin(N,0.1,0.4); [TFR,T,F]=tfrpwv(sig,1:N,N,window("hm",N/4+1)); clf;gcf().color_map= jetcolormap(128); Sgrayplot(T,F,TFR'); | ![]() | ![]() |