wpspectrum - Wavelet packet spectrum

Syntax

[SPEC,TIMES,FREQ] = wpspectrum(WPT,Fs) [...] = wpspectrum(WPT,Fs,'plot') [...,TNFO] = wpspectrum(...)

Description

[SPEC,TIMES,FREQ] = wpspectrum(WPT,Fs) returns a matrix of wavelet packet spectrum estimates, SPEC, for the binary wavelet packet tree object, WPT. Fs is the sampling frequency in Hertz. SPEC is a 2^J-by-N matrix where J is the level of the wavelet packet transform and N is the length of the time series. TIMES is a 1-by-N vector of times and FREQ is a 1-by-2^J vector of frequencies.

[...] = wpspectrum(WPT,Fs,'plot') displays the wavelet packet spectrum.

[...,TNFO] = wpspectrum(...) returns the terminal nodes of the wavelet packet tree in frequency order.

Input Arguments

WPT

WPT is a binary wavelet packet tree of class wptree.

Fs

Sampling frequency in Hertz as a scalar of class double.

Default: 1

plot

The string 'plot' displays the wavelet packet spectrum. Enter 'plot' after Fs to produce a plot of the wavelet packet spectrum.

Output Arguments

`SPEC`	Wavelet packet spectrum. `SPEC` is a 2^J-by-N matrix where J is the level of the wavelet packet transform and N is the length of node 0 in the wavelet packet tree object. The frequency spacing between the rows of `SPEC` is Fs/2^J+1.
`TIMES`	Time vector. `TIMES` is a vector of times in seconds equal in length to node 0 of the wavelet packet tree object. The time spacing between elements is 1/Fs.
`FREQ`	Frequency vector. `FREQ` is a vector of frequencies of length 2^J where J is the level of the wavelet packet tree object. The frequency spacing in `FREQ` is Fs/2^J+1.
`TNFO`	Terminal nodes. `TNFO` is a vector of the terminal nodes of the wavelet packet tree object in frequency order.

Definitions

Wavelet Packet Spectrum

The wavelet packet spectrum contains the absolute values of the coefficients from the frequency-ordered terminal nodes of the input binary wavelet packet tree. The terminal nodes provide the finest level of frequency resolution in the wavelet packet transform. If J denotes the level of the wavelet packet transform and Fs is the sampling frequency, the terminal nodes approximate bandpass filters of the form:

At the terminal level of the wavelet packet tree, the transform divides the interval from 0 to the Nyquist frequency into bands of approximate width

Examples

Wavelet packet spectrum for signal consisting of two sinusoids with disjoint support:

fs = 500;
t = 0:1/fs:4;
y = sin(32*pi*t).*(t<2) + sin(128*pi*t).*(t>=2);
subplot(2,1,1);
plot(t,y); 
axis tight
title('Analyzed Signal');

% Wavelet packet spectrum
level = 6;
wpt = wpdec(y,level,'sym6');
subplot(2,1,2);
[S,T,F] = wpspectrum(wpt,fs,'plot');

Wavelet packet spectrum of chirp:

fs = 1000;
t = 0:1/fs:2; 
% create chirp signal        
y = sin(256*pi*t.^2);

% Plot the analyzed signal
subplot(2,1,1);
plot(t,y); 
axis tight
title('Analyzed Signal');

% Wavelet packet spectrum
level = 6;
wpt = wpdec(y,level,'sym8');
subplot(2,1,2);
[S,T,F] = wpspectrum(wpt,fs,'plot');

Algorithms

wpspectrum computes the wavelet packet spectrum as follows:

Extract the wavelet packet coefficients corresponding to the terminal nodes. Take the absolute value of the coefficients.
Order the wavelet packet coefficients by frequency ordering.
Determine the time extent on the original time axis corresponding to each wavelet packet coefficient. Repeat each wavelet packet coefficient to fill in the time gaps between neighboring wavelet packet coefficients and create a vector equal in length to node 0 of the wavelet packet tree object.