<< tfrscalo Bilinear Time-Frequency Processing in the Affine Class tfrunter >>

stftb >> stftb > Bilinear Time-Frequency Processing in the Affine Class > tfrspaw

tfrspaw

Smoothed Pseudo Affine Wigner time-frequency distributions.

Calling Sequence

[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')

Parameters

X :

A Nx elements vector (signal) or a Nx by 2 array signal (cross-Smoothed Pseudo Affine Wigner distribution.).

T:

a real Nt vector with elements in [1 Nx] : time instant(s) on which the TFR is evaluated (default: 1:NX).

K :

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.

NH0 :

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)).

NG0 :

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.

FMIN:

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.

FMAX :

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.

N :

positive integer: number of analyzed voices. When unspecified, you have to enter it at the command line from the plot of the spectrum.

TRACE :

A boolean (or a real scalar) if true (or nonzero),the progression of the algorithm is shown (default : %f).

'plot':

if one input parameter is 'plot', tfrqview is called and the time-frequency representation will be plotted.

TFR :

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.

F :

vector of normalized frequencies (geometrically sampled from FMIN to FMAX).

Description

tfrspaw generates the auto- or cross- Smoothed Pseudo Affine Wigner distributions.

Examples

Interactive use

N = 128;
sig = altes(N,0.1,0.45); 
tfrspaw(sig,"plot");

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');

Authors


Report an issue
<< tfrscalo Bilinear Time-Frequency Processing in the Affine Class tfrunter >>