<< tfrrsp Reassigned Time-Frequency Processing tfrrstan >>

stftb >> stftb > Reassigned Time-Frequency Processing > tfrrspwv

tfrrspwv

Reassigned smoothed pseudo Wigner-Ville distribution

Calling Sequence

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

Parameters

X :

A Nx elements vector or a Nx by 2 array signal.

T :

the time instant(s) with elements in [1 Nx] (default : 1:length(X))

N :

number of frequency bins (default : length(X)). For faster computation N should be a power of 2.

G :

a real vector of odd length: the time smoothing window (default : Hamming(N/10)).

H :

a real vector of odd length: the frequency smoothing window, (default : Hamming(N/4)).

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 :

a real N by Nt array: the time-frequency representation.

RTFR :

A real N by Nt array: the reassigned time-frequency representation.

HAT :

a N by Nt complex matrix: the reassignment vectors.

Description

tfrrspwv computes the smoothed pseudo Wigner-Ville distribution and its reassigned version.

Examples

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

Authors


Report an issue
<< tfrrsp Reassigned Time-Frequency Processing tfrrstan >>