HHSpecPlot: Display Hilbert periodogram

In hht: The Hilbert-Huang Transform: Tools and Methods

Display Hilbert periodogram

Description

This function displays the Hilbert periodogram, with options to plot individual IMFs and also the Fourier periodogram for comparison.

Usage

HHSpecPlot(hspec, freq.span = NULL, scaling = "none", imf.list = NULL, 
    show.total = TRUE, show.fourier = FALSE, scale.fourier = FALSE, 
    show.imfs = FALSE, legend = TRUE, ...)

Arguments

hspec	Data structure returned by HHSpectrum
freq.span	Frequency range to plot, NULL plots all of them
scaling	Amplitude scaling, can be "log" (log 10), "sqrt" (square root), defaults to "none".
imf.list	Which IMFs to plot, requires show.imfs = TRUE.
show.total	Show the ensemble Hilbert spectrogram
show.fourier	Show the Fourier periodogram
scale.fourier	Scale Fourier and Hilbert spectra to each other for easier comparison
show.imfs	Plot individual IMF spectra
legend	Determines whether or not a legend is shown
...	This function supports some optional parameters as well: xlab - X axis label ylab - Y axis label legend.location - where to put the legend total.col - color of ensemble Hilbert periodogram total.lwd - lwd of ensemble Hilbert periodogram total.lty - lty of ensemble Hilbert periodogram imf.cols - colors of IMF periodogram imf.lwd - lwds of IMF periodogram imf.lty - ltys of IMF periodogram fourier.col - color of Fourier periodogram fourier.lwd - lwd of Fourier periodogram fourier.lty - lty of Fourier periodogram main - figure title

Details

This function plots the Hilbert periodogram of a signal, with options to show periodograms of individual IMFs. You can also plot a simple Fourier periodogram for comparison.

Author(s)

Daniel Bowman [e-mail]

See Also

HHSpectrum, HHGramImage

Examples

#Here we see how the EMD produces a dyadic filter bank for uniform random noise
#The frequency distributions of all but the first IMF display a Chi-Square distribution
#See Huang, N. E. & Wu, Z. 
#A review on Hilbert-Huang Transform: Method and its applications to geophysical studies.
#Reviews of Geophysics, 2008, 46, RG2006

#The EMD of this signal may take a couple of minutes to run

set.seed(628)
sig  <-  runif(10000)
tt  <-  seq_len(length(sig)) * 0.01

## Not run: emd.result  <-  Sig2IMF(sig, tt)

dfreq  <-  0.1
## Not run: hspec  <-  HHSpectrum(emd.result, dfreq)

## Not run: HHSpecPlot(hspec, show.imfs = TRUE, 
imf.list = 1:10, show.total = TRUE, scaling = "sqrt", 
imf.lwd = rep(2, 10), total.lty = 3)
## End(Not run)

hht documentation built on March 31, 2023, 10:08 p.m.