as.pulse | R Documentation |
Allows to convert instantaneous frequency determination results
into a single 'pulse' object. This is the format generated by inst.pulse
(and gzc if output = 2
)
as.pulse( dt, f, a = NULL, m = NULL, idt = NULL, mode = NULL, repl = 1, order = NA ) is.pulse(pulse)
dt |
a vector of length n for the depth or time reference |
f |
a data.frame or matrix of n rows of the instantaneous frequencies |
a |
a data.frame or matrix of n rows of the instantaneous amplitudes |
m |
a data.frame or matrix of n rows of the components from which the frequencies and amplitudes were computed from |
idt |
data.frame or matrix of n rows of identity tuning: new dt coordinates to remove the frequency modulation |
mode |
the mode sequence index to give to each replicated IMFs |
repl |
a vector for the number of replicates or a matrix, indicating in which replicate set each point is |
order |
the order of the imf, typically from higher frequency to lower frequency |
pulse |
a pulse object to check |
a list made of $dt (depth/time), $f (instantaneous frequency), $a (instantaneous amplitude) if a is provided, $repl (the replicate id of each point) and $mode (the mode id of each point).
set.seed(42) n <- 600 dt <- seq_len(n) p1 <- 30 p2 <- 240 s30 <- (1 + 0.6 * sin(dt*2*pi/p2)) * sin(dt*2*pi/p1) s240 <- 2 * sin(dt*2*pi/p2) xy <- s30 + s240 dec <- as.emd(xy = xy, dt = dt, imf = matrix(c(s30, s240), ncol = 2)) plot_emd(dec, pdf = FALSE, style = 1) pulse <- inst.pulse(dec, last = TRUE, breaks = 200, bins = 40, cut = 10) is.pulse(pulse) simp.pulse <- as.pulse(pulse$dt, pulse$f) str(simp.pulse)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.