This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English version of the page.

Note: This page has been translated by MathWorks. Click here to see
To view all translated materials including this page, select Country from the country navigator on the bottom of this page.

waveletsupport

CWT filter bank time supports

Syntax

spsi = waveletsupport(fb)
spsi = waveletsupport(fb,thresh)

Description

example

spsi = waveletsupport(fb) returns the wavelet time supports defined as the time interval in which all of the wavelet's energy occurs. The default tolerance is 99.99% of the energy. The time supports are returned in the MATLAB® table spsi. The wavelets are normalized to have unit energy.

spsi = waveletsupport(fb,thresh) specifies the threshold for the integrated energy. The time support of the wavelet is defined as the first instant the integrated energy exceeds thresh and the last instant the integrated energy is less than 1-thresh. If unspecified, thresh defaults to 1e–4.

Examples

collapse all

Create a continuous wavelet transform filter bank. Set the sampling frequency to 1000 Hz and the frequency limits to range from 100 Hz to 200 Hz. Obtain the time-supports of the wavelets in the filter bank.

fb = cwtfilterbank('SamplingFrequency',1000,'FrequencyLimits',[100 200]);
spsi = waveletsupport(fb)
spsi=11×5 table
      CF      IsAnalytic    TimeSupport    Begin      End 
    ______    __________    ___________    ______    _____

       200    "Analytic"       0.032       -0.016    0.016
    186.61    "Analytic"       0.034       -0.017    0.017
    174.11    "Analytic"       0.038       -0.019    0.019
    162.45    "Analytic"        0.04        -0.02     0.02
    151.57    "Analytic"       0.042       -0.021    0.021
    141.42    "Analytic"       0.046       -0.023    0.023
    131.95    "Analytic"       0.048       -0.024    0.024
    123.11    "Analytic"       0.052       -0.026    0.026
    114.87    "Analytic"       0.056       -0.028    0.028
    107.18    "Analytic"        0.06        -0.03     0.03
       100    "Analytic"       0.064       -0.032    0.032

Obtain the time domain wavelets from the filter bank and plot their magnitudes. Use the table to set the x-limits of the plot equal to the wavelet with largest time support.

[psi,t] = wavelets(fb);
plot(t,abs(psi))
xlim([spsi.Begin(end) spsi.End(end)])
title('Time Domain Wavelets')

Input Arguments

collapse all

Continuous wavelet transform (CWT) filter bank, specified as a cwtfilterbank object.

Time-support threshold for the wavelet, specified as a positive real number between 0 and 0.05. The time support of the wavelet is defined as the first instant the integrated energy of the wavelet exceeds thresh and the last instant the integrated energy is less than 1-thresh.

Data Types: double

Output Arguments

collapse all

Wavelet time supports, returned as a Ns-by-5 MATLAB table, where Ns is the number of wavelet bandpass frequencies (equal to the number of scales). The table has the following five variables:

Wavelet center frequency, returned as a positive real number.

Data Types: double

Wavelet designation, returned as a string. Wavelets that do not decay to 5% of their peak value at the Nyquist frequency are not considered analytic. The time support information for those wavelets are returned as NaNs.

Data Types: string

Wavelet time support, returned in samples, seconds, or MATLAB durations. The units of TimeSupport depend on whether you specify SamplingFrequency or SamplingPeriod. If you specify a SamplingFrequency, the units are in seconds. If you specify a SamplingPeriod, the units are the same as the SamplingPeriod. If no SamplingFrequency or SamplingPeriod is specified, the units are in samples.

Data Types: double

Beginning of the wavelet time support, returned as an integer. Begin is defined as the first instant the wavelet integrated energy exceeds the default threshold, 1e–4. Begin has the same units as TimeSupport.

Data Types: double

End of the wavelet time support, returned as an integer. End is defined as the last instant the wavelet integrated energy is less than 1–1e–4. End has the same units as TimeSupport.

Data Types: double

Data Types: table

Introduced in R2018a