Smoothed Pseudo Affine Wigner time-frequency distributions.
[TFR,T,F] = tfrspaw(X) [TFR,T,F] = tfrspaw(X, T) [TFR,T,F] = tfrspaw(X, T, K) [TFR,T,F] = tfrspaw(X, T, K, NH0) [TFR,T,F] = tfrspaw(X, T, K, NH0, NG0) [TFR,T,F] = tfrspaw(X, T, K, NH0, NG0, FMIN,FMAX) [TFR,T,F] = tfrspaw(X, T, K, NH0, NG0, FMIN,FMAX, N) [TFR,T,F] = tfrspaw(X, T, K, NH0, NG0, FMIN,FMAX, N, TRACE) [TFR,T,F] = tfrspaw(...,'plot')
A Nx elements vector (signal) or a Nx by 2 array signal (cross-Smoothed Pseudo Affine Wigner distribution.).
a real Nt vector with elements in [1 Nx] : time instant(s) on which the TFR is evaluated (default: 1:NX).
label of the K-Bertrand distribution. The
distribution with parameterization function
lambdak(u,K) =
(K(exp(-u)-1)/(exp(-Ku)-1))^(1/(K-1))
is computed
(default is 0).
K=-1: Smoothed pseudo (active) Unterberger distribution
K=0: Smoothed pseudo Bertrand distribution
K=1/2: Smoothed pseudo D-Flandrin distribution
K=2: Affine smoothed pseudo Wigner-Ville distribution.
half length of the analyzing wavelet at coarsest scale. A Morlet wavelet is used. NH0 controles the frequency smoothing of the smoothed pseudo Affine Wigner distribution. (default is sqrt(Nx)).
half length of the time smoothing window. If NG0 is set to zero the time smoothing window is a rectangular one and the length is automatically determined.
a positive scalar in ]0 0.5], the normalized lower frequency bound in (Hz) of the analyzed signal. When unspecified, you have to enter it at the command line from the plot of the spectrum.
a positive scalar in ]0 0.5], the normalized upper frequency bound (in Hz) of the analyzed signal. When unspecified, you have to enter it at the command line from the plot of the spectrum.
positive integer: number of analyzed voices. When unspecified, you have to enter it at the command line from the plot of the spectrum.
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.
N by Nt real array: the time-frequency matrix containing the coefficients of the decomposition (abscissa correspond to uniformly sampled time, and ordinates correspond to a geometrically sampled frequency). First row of TFR corresponds to the lowest frequency.
vector of normalized frequencies (geometrically sampled from FMIN to FMAX).
tfrspaw generates the auto- or cross- Smoothed Pseudo Affine Wigner distributions.
Interactive use
Non interactive use
N = 128; sig = altes(N,0.1,0.45); T = 1:N; K = 0; [tfr,t,f] = tfrspaw(sig,1:N,0,2*sqrt(N),0,0.1,0.35,32); clf;gcf().color_map = jetcolormap(128); grayplot(t,f,tfr'); | ![]() | ![]() |